多边形重心及面积

#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
#define INF 0.0000001
const int N=10010;
struct point{
	double x;
	double y;
	point():x(0),y(0){
	}
}p[N];
int main(){
	int t,n;
	scanf("%d",&t);
	while(t--){
		double sum=0;	
		point a; 
		scanf("%d",&n);
	
		for(int i=0;i<n;++i)
			scanf("%lf%lf",&p[i].x,&p[i].y);
			for(int i=1;i<=n;++i){//i==n时回到第一个点，下标为0 
				double s=(p[i%n].x*p[i-1].y-p[i%n].y*p[i-1].x)/2.0;//计算三角形面积
				  
				  sum+=s;//面积累加和
				   a.x+=s*(p[i%n].x+p[i-1].x)/3.0;//三角形横坐标平均值*三角形面积的累加和 
				   a.y+=s*(p[i%n].y+p[i-1].y)/3.0;
			}
			if(fabs(sum-0)<INF) puts("0.000 0.000");
			else
			printf("%.3lf %.3lf\n",fabs(sum),(a.x+a.y)/sum);
	}
	return 0; 
}