skywalker — Wed, 01 Aug 2007 14:13:00 GMT

#include
#include
#include
using namespace std;

struct point
{
double x;
double y;
};

//求多边�Ş的重心算�?br>//说明:
//求多边�Ş重心�q�不是简单的把求三角形的重心公式推广��p��?br>//我的��法是在�q�面上取一�?一般取原点, �q�样可以减少很多计算, 而且使思�\更清晰^_^)
//�q�样��得��C��N个三角�ŞOP[i]P[i+1](其中点的��序要�ؓ逆时针的),
//分别求出�q�N个三角�Ş的重心Ci和面�U�Ai(注意此处面积是又向面�U? ��是用叉乘求面积时保留其正负�?
//在求出A = A1+A2+...+AN(同样保留正负��L(f��ng)��代数相加)
//最�l�重心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;
}

skywalker 2007-08-01 22:13 发表评论

skywalker — Wed, 01 Aug 2007 12:12:00 GMT

求两个或N个数的最大公�U�数(gcd)和最��公倍数(lcm)的较优算�?/span>

#include
#include
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;
}

//两个数的最大公�U�数--�Ƨ几里得��法
int gcd(int a, int b)
{
if (a < b)
swap(a, b);

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

//n个数的最大公�U�数��法
//说明:
//把n个数保存��Z��个数�l?br>//参数为数�l�的指针和数�l�的大小(需要计��的数的个数)
//然后先求出gcd(a[0],a[1]), 然后��所求的gcd与数�l�的下一个元素作为gcd的参数��l�求gcd
//�q�样��׃�生一个递归的求ngcd的算�?br>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个数的最��公倍数��法
//��法�q�程和n个数的最大公�U�数求法�c�M��
//求出头两个的最��公倍数,再将�ƺ和大三个数求最��公倍数直到数组末尾
//�q�样产生一个递归的求nlcm的算�?br>int nlcm(int *a, int n)
{
if (n == 1)
return *a;
else
return lcm(a[n-1], nlcm(a, n-1));
}

skywalker 2007-08-01 20:12 发表评论

日韩视频在线观看免费,亚洲欧洲日本国产,国产欧美一区二区三区另类精品