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

#include<stdio.h>
#include<math.h>
const int N = 10005;

typedef struct
{
	double x;
	double y;
}POINT;

POINT point[N];

inline double CrossProduct(int i)
{
	return point[i].x * point[i + 1].y - point[i + 1].x * point[i].y;
}

int main()
{
	int n,t;
	int i;
	double temp;
	double area,x,y;
	scanf("%d", &t);
	while(t--)
	{
		scanf("%d",&n);
		for(i = 0; i < n; ++i)
			scanf("%lf %lf",&point[i].x, &point[i].y);

		point[n].x = point[0].x;//记下第一个点，形成循环
		point[n].y = point[0].y;

		area = 0.0;
		x = 0.0;
		y = 0.0;

		for(i = 0; i < n; ++i)
		{
			temp = CrossProduct(i) / 2.0;
			area += temp;
			x += temp * (point[i].x + point[i + 1].x) / 3.0;
			y += temp * (point[i].y + point[i + 1].y) / 3.0;
		}

		area = fabs(area);//我说了，叉乘有方向，最后记得取绝对值
		if(area < 0.0000001)
			printf("0.000 0.000\n");
		else
			printf("%.3lf %.3lf\n", area, fabs(x + y) / area);
	}
}