acm-多边形重心问题

多边形重心问题

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

难度： 5

描述

在某个多边形上，取n个点，这n个点顺序给出，按照给出顺序将相邻的点用直线连接， （第一个和最后一个连接），所有线段不和其他线段相交，但是可以重合，可得到一个多边形或一条线段或一个多边形和一个线段的连接后的图形；  
如果是一条线段,我们定义面积为0，重心坐标为（0,0）.现在求给出的点集组成的图形的面积和重心横纵坐标的和；

输入

第一行有一个整数0
每组数据第一行有一个整数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

代码：

#include
#include
using namespace std;
typedef struct node
{
  double x;
  double y;
}Node;

Node p[10001];
double cal(int i,int j)
{
  double res;
  res=p[i].x*p[j].y-p[j].x*p[i].y;
  return res;
}
int main()
{
  int n;
  scanf("%d",&n);
  while(n--)
  {
    int m;
    scanf("%d",&m);
    for(int i=0;i
    {
      scanf("%lf%lf",&p[i].x,&p[i].y); 
    }
    double s1=0,s2=0;
    double area,gra;
    for(int i=0;i
    {
      s1+=p[(i+1)%m].y*p[i].x;
      s2+=p[(i+1)%m].y*p[(i+2)%m].x;
    }
    area=fabs(s1-s2)/2;
    if(area==0)
      gra=0;
    else
    {
      p[m].x = p[0].x; 
                p[m].y = p[0].y;
      double A = 0; 
                for (int i = 0; i < m; i++) 
                { 
                        A += cal(i, i+1); 
                } 
                A /= 2; 
                   
                double cx = 0;
      double cy = 0; 
                for (int i = 0; i < m; i++) 
                { 
                        cx += (p[i].x + p[i+1].x)*cal(i, i+1); 
                        cy += (p[i].y + p[i+1].y)*cal(i, i+1); 
                } 
                cx /= 6*A; 
                cy /= 6*A;
      gra=cx+cy;   
    } 
    printf("%.3lf %.3lf\n",area,gra); 
  }
  return 0;
}