#include <iostream>
#include <cmath>
#include <iomanip>
using namespace std;

struct point
{
 double x;
 double y;
};

//求多邊形的重心算法
//說明:
//求多邊形重心并不是簡單的把求三角形的重心公式推廣就行了
//我的算法是在平面上取一點(一般取原點, 這樣可以減少很多計算, 而且使思路更清晰^_^)
//這樣就得到了N個三角形OP[i]P[i+1](其中點的順序要為逆時針的),
//分別求出這N個三角形的重心Ci和面積Ai(注意此處面積是又向面積, 就是用叉乘求面積時保留其正負號)
//在求出A = A1+A2+...+AN(同樣保留正負號的代數相加)
//最終重心C = sigma(Ai+Ci)/A;
point gravity(point *p, int n)
{
 double area = 0;
 point center;
 center.x = 0;
 center.y = 0;

 for (int i = 0; i < n-1; i++)
 {
  area += (p[i].x*p[i+1].y - p[i+1].x*p[i].y)/2;
  center.x += (p[i].x*p[i+1].y - p[i+1].x*p[i].y) * (p[i].x + p[i+1].x);
  center.y += (p[i].x*p[i+1].y - p[i+1].x*p[i].y) * (p[i].y + p[i+1].y);
 }

 area += (p[n-1].x*p[0].y - p[0].x*p[n-1].y)/2;
 center.x += (p[n-1].x*p[0].y - p[0].x*p[n-1].y) * (p[n-1].x + p[0].x);
 center.y += (p[n-1].x*p[0].y - p[0].x*p[n-1].y) * (p[n-1].y + p[0].y);

 center.x /= 6*area;
 center.y /= 6*area;

 return center;
}

                          求兩個或N個數的最大公約數(gcd)和最小公倍數(lcm)的較優算法

#include <iostream>
#include <cmath>
using namespace std;

int gcd(int a, int b);
int ngcd(int *a, int n);
int lcm(int a, int b);
int nlcm(int *a, int n);

int main()
{
 //int a,b;
 //cin >> a >> b;
 //cout << lcm(a, b) << endl;
 int *a = new int[3];
 a[0] = 3;
 a[1] = 4;
 a[2] = 5;
 cout << nlcm(a, 3) << endl;

 return 0;
}

//兩個數的最大公約數--歐幾里得算法
int gcd(int a, int b)
{
 if (a < b)
  swap(a, b);

 if (b == 0)
  return a;
 else
  return gcd(b, a%b);
}

//n個數的最大公約數算法
//說明:
//把n個數保存為一個數組
//參數為數組的指針和數組的大小(需要計算的數的個數)
//然后先求出gcd(a[0],a[1]), 然后將所求的gcd與數組的下一個元素作為gcd的參數繼續求gcd
//這樣就產生一個遞歸的求ngcd的算法
int ngcd(int *a, int n)
{
 if (n == 1)
  return *a;

 return gcd(a[n-1], ngcd(a, n-1));
}

//兩個數的最小公倍數(lcm)算法
//lcm(a, b) = a*b/gcd(a, b)
int lcm(int a, int b)
{
 return a*b/gcd(a, b);
}

//n個數的最小公倍數算法
//算法過程和n個數的最大公約數求法類似
//求出頭兩個的最小公倍數,再將欺和大三個數求最小公倍數直到數組末尾
//這樣產生一個遞歸的求nlcm的算法
int nlcm(int *a, int n)
{
 if (n == 1)
  return *a;
 else
  return lcm(a[n-1], nlcm(a, n-1));
}