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

多边形求重心：

       得先知道三角形的重心： 设三个顶点为A（x1,y1） B(x2,y2)  C(x3,y3)，则重心为（(x1+x2+x3)/3,(y1+y2+y3)/3）；

      三角形的面积为：向量AB^AC/2，(肯定有正有负)，最后加绝对值就行。但是在本题中要求的是多边形的面积，考虑到有凹变形这种特殊情况，所以在算的过程中不要加绝对值，

在最后输出面积的时候再加绝对值。

      接下来就开始攻多边形了：

   将多边形划分成三角形。

          多边形的重心 =  每个三角形的重心*每个三角形的面积（要带入正负）的和/多边形的面积（也要带入正负）

           重心的坐标不会因为正负号而出现误差，在算的时候会抵消掉。

        多边形的面积 = 每个三角形的面积之和 取绝对值。

PS：

          我在做题的时候发现在求多边形重心的时候，分子分母都有多边形的面积和那直接不用不就好了，这是不可以的，因为本题中会给出不是顶点的点（点在直线上）。

         N变形的重心坐标为：(各个顶点的横坐标之和/N，各个顶点的纵坐标之和/N）。.

#include <stdio.h>
#include <math.h>
struct point
{
    double x;
    double y;
}p[10005];
double mul(struct point p1,struct point p2,struct point p3)
{
   // printf("%lf\n",(p2.x-p1.x)*(p3.y-p1.y)-(p3.x-p1.x)*(p2.y-p1.y));
    return (p2.x-p1.x)*(p3.y-p1.y)-(p3.x-p1.x)*(p2.y-p1.y);
}
int main (void)
{
   // freopen("题目3.txt","r",stdin);
    int t,n,i,j;
    double ssum,xsum,area,ysum;
    while(scanf("%d",&t)==1)
    {
        for(i=1;i<=t;i++)
        {
            xsum=0;
            ssum=0;
            ysum=0;
            scanf("%d",&n);
            for(j=1;j<=n;j++)
            {
                scanf("%lf%lf",&p[j].x,&p[j].y);
            }
            for(j=2;j<n;j++)
            {
                area=mul(p[1],p[j],p[j+1])/2.0;
                ssum+=area;
                         xsum+=((p[1].x+p[j].x+p[j+1].x)*area);
                         ysum+=((p[1].y+p[j].y+p[j+1].y)*area);


            }
            xsum/=ssum;
            ysum/=ssum;
            printf("%.3lf ",fabs(ssum));
            if(ssum==0)
                printf("0.000\n");
            else
                printf("%.3lf\n",(xsum+ysum)/3.0);
        }
    }
    return 0;
}