1. 程式人生 > >多邊形重心及面積

多邊形重心及面積

#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; 
}