BIT1035 远古石碑

这题调了N久啊有木有，交到hdu都AC了的代码，在BIT就是WA，我调了三个小时了有木有

问题：

已知多边形的各点的坐标

计算一个多边形的重心的坐标，多边形的质量分布均匀

按逆时针顺序给出各个点的坐标

如图

将多边形作如图切割，红点是各个三角形的重心，连起来就是另一个多边形

将各个三角形看成一个点，即整个三角形的重量都集中在一个点上

然后求各个重心所构成的多边形的重心

例如，各个重心的坐标为(x1,y2)(x2,y2)...(xn,yn)

各个重心的重量为w1,w2,,,,,wn

多边形的重心为X,Y

X=(x1*w1+x2*w2+....+xn*wn)/n

Y=(y1*w1+y2*w2+.....+yn*wn)/n

接下来上代码

#include<iostream>#include<algorithm>#include<cstdio>using namespace std;class point//坐标{public:double x,y;};class Vector//向量{public:double x,y;double cross(Vector a)//叉积  {  return x*a.y-a.x*y;  }  };Vector operator-(point a,point b){Vector res;res.x=a.x-b.x;res.y=a.y-b.y;return res;}point CalcCenter(point a,point b,point c){point res;res.x=(a.x+b.x+c.x);res.y=(a.y+b.y+c.y);return res;}int main(){int n;while (~scanf("%d",&n)){point ans;ans.x=0;ans.y=0;double sumWeight=0;point a,b,c;scanf("%lf %lf",&a.x,&a.y);scanf("%lf %lf",&b.x,&b.y);scanf("%lf %lf",&c.x,&c.y);Vector t1,t2;t1=b-a;t2=c-a;point p=CalcCenter(a,b,c);double area=t1.cross(t2);ans.x+=p.x*area;ans.y+=p.y*area;sumWeight+=area;for (int i = 3; i < n; i++){point temp;scanf("%lf %lf",&temp.x,&temp.y);t1=t2;t2=temp-a;b=c;c=temp;p=CalcCenter(a,b,c);area=t1.cross(t2);ans.x+=area*p.x;ans.y+=area*p.y;sumWeight+=area;}ans.x/=sumWeight;ans.y/=sumWeight;printf("%.2lf %.2lf\n",ans.x/3+0.00000001,ans.y/3+0.00000001);//不加这个0.00000000001就会出现输出-0.0的情况，在hdu可以ac，在bit就wa了}return 0;}