多边形重心问题 NYOJ 3 思路+代码

多边形重心问题

时间限制：3000 ms  |  内存限制：65535 KB
难度：5

描述在某个多边形上，取n个点，这n个点顺序给出，按照给出顺序将相邻的点用直线连接， （第一个和最后一个连接），所有线段不和其他线段相交，但是可以重合，可得到一个多边形或一条线段或一个多边形和一个线段的连接后的图形； 

如果是一条线段,我们定义面积为0，重心坐标为（0,0）.现在求给出的点集组成的图形的面积和重心横纵坐标的和；

输入第一行有一个整数0<n<11,表示有n组数据；

每组数据第一行有一个整数m<10000,表示有这个多边形有m个顶点；
输出输出每个多边形的面积、重心横纵坐标的和，小数点后保留三位；
样例输入

3
3
0 1
0 2
0 3
3
1 1
0 0
0 1
4
1 1
0 0
0 0.5
0 1

样例输出

0.000 0.000
0.500 1.000
0.500 1.000

上传者
iphxer

思路：
首先得知道两个定理。
已知三角形各顶点A（Xi，Yi）（ i=1，2，3）

S=((x2-x1)*(y3-y1)-(x3-x1)*(y2-y1));

重心坐标（Xg，Yg）为

Xg=(x1+x2+x3)/3;

Yg=(y1+y2+y3)/3;

有两种方法，我用的第二种。

算法一：一般适合凸多边形

      n边多边形可以分成n-2个三角形，将这些三角形看做质点（质点的位置是三角形的重心x1,x2,..，质量是面积s1,s2,..），

      那么多边形就由这些质点组成，质点坐标以其质量为权的加权算术平均数即是多边形重心坐标x。

x=(x1*s1+x2*s2+...)/(s1+s2+...)
y=(y1*s1+y2*s2+...)/(s1+s2+...)

s=s1+s2+...

算法二：任意多边形

       将算法一改进，n边多边形中每两个点（有顺序）加上原点可构成n个三角形,

       将这些三角形看做质点（质点的位置是三角形的重心x1,x2...，质量是面积(有正负)s1,s2,...），

       那么多边形就由这些质点组成，质点坐标以其质量为权的加权算术平均数即是多边形重心坐标x。

       多边形的面积s是这n个三角形面积(有正负)的代数和的绝对值。

x=(x1*s1+x2*s2+...)/(s1+s2+...)
x=(x1*s1+x2*s2+...)/(s1+s2+...)

s=|s1+s2+...|

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

double cross(double a[2],double b[2]){   //求向量a,b的叉积大小
return a[0]*b[1]-a[1]*b[0];
}

int main(void)
{
int ncase;
cin>>ncase;
while(ncase--)
{
int i,k;
double S,tS,a[10002][2],sx,sy;
cin>>k;
for(i=1;i<=k;i++)
cin>>a[i][0]>>a[i][1];
S=0.;sx=0.;sy=0.;          //S面积,xy横纵坐标和
for(i=1;i<=k;i++)
{
tS=cross(a[i],a[i%k+1])/2.;//最后一个与第一个相连 取余
S+=tS;
sx+=tS*(a[i][0]+a[i%k+1][0])/3;//s*x1
sy+=tS*(a[i][1]+a[i%k+1][1])/3;
}
if(fabs(S)<1e-7)//精度判断
puts("0.000 0.000");
else
printf("%.3lf %.3lf\n",fabs(S),(sx+sy)/S);
}
return 0;
}