求多边形的重心-hdu 1115

//先把所有的坐标按一个方向排序，比如逆时针，这里可以用极坐标排序；

//排好序后，以原点为基点，用n边形的n个点中的每相邻的两个点和基点组成一个三角形，这样就构成了n-1个三角形，并求出相应三角形的重心ci和叉积面积ai（这里保留方向性）；（这里为什么用原点作为基点，因为用原点作为基点时，求叉积时简单），设A=A1+A2....Ai，最后的n多边形重心C=∑(ci×Ai)/A；

//上面的公式详见：http://www.cppblog.com/linqiu/archive/2007/08/01/29153.html

 

//代码如下：

 

#include<stdio.h>
#include<math.h>
#include<algorithm>
using namespace std;
struct point
{
	__int64 x,y;
}pos[1001000];
void gravity(int n,double &x,double &y)
{
	int i;
	double A;
	A=x=y=0;
	for(i=0;i<n-1;i++)
	{
		A+=(pos[i].x*pos[i+1].y-pos[i+1].x*pos[i].y);
		x+=(pos[i].x*pos[i+1].y-pos[i+1].x*pos[i].y)*(pos[i].x+pos[i+1].x);
		y+=(pos[i].x*pos[i+1].y-pos[i+1].x*pos[i].y)*(pos[i].y+pos[i+1].y);
	}
	A+=(pos[n-1].x*pos[0].y-pos[0].x*pos[n-1].y);
	x+=(pos[n-1].x*pos[0].y-pos[0].x*pos[n-1].y)*(pos[n-1].x+pos[0].x);
	y+=(pos[n-1].x*pos[0].y-pos[0].x*pos[n-1].y)*(pos[n-1].y+pos[0].y);
	x/=3*A;
	y/=3*A;
}
int main()
{
	int n,i,cas;
	double x,y;
	scanf("%d",&cas);
	while(cas--)
	{
		scanf("%d",&n);	
		for(i=0;i<n;i++)
			scanf("%I64d%I64d",&pos[i].x,&pos[i].y);
		gravity(n,x,y);
		printf("%.2lf %.2lf\n",x,y);
	}
	return 0;
}