地址:http://acm.bit.edu.cn/mod/programming/view.php?a=521
题意:求多边形重心。
这里的多边形重心是指质量均匀分布在整个多边形上的重心,而不是质量分布在顶点上的,所以不能直接坐标全加起来再除以n。
三角剖分,因为三角形的重心即等于质量分布在三角形三个顶点上时的重心,即等于三个顶点坐标之和除以3.那么就可以求出这些三角形的重心,
那么多边形质量就相当于质量集中在这些三角形的重心上,并且每个重心所占质量比重不同,取决于这些三角形的面积。
然后就可以用质点系的重心求法,把这些重心坐标分别乘以权(三角形面积)相加再除以多边形总面积,即得。
http://blog.csdn.net/zhangwei1120112119/article/details/8563262
#include<iostream>
#include<cmath>
typedef struct POINT
{
double x,y;
}point;
point a[100005],ans,temp;
double operator ^(point a,point b)
{
return a.x*b.y-a.y*b.x;
}
int main()
{
int n,i;
while(~scanf("%d",&n))
{
ans.x=ans.y=0;
double area=0;
for(i=1;i<=n+1;i++)
{
if(i!=n+1) scanf("%lf%lf",&a[i].x,&a[i].y);
if(i==1) a[n+1]=a[1];
else
{
temp.x=(a[i].x+a[i-1].x)/3,temp.y=(a[i].y+a[i-1].y)/3;
temp.x*=a[i]^a[i-1],temp.y*=a[i]^a[i-1];
ans.x+=temp.x,ans.y+=temp.y;
area+=a[i]^a[i-1];
}
}
ans.x/=area,ans.y/=area;
printf("%.2lf %.2lf\n",ans.x+0.00000001,ans.y+0.00000001);
}
return 0;
}