多邊形重心及面積
阿新 • • 發佈:2019-02-01
#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; }