南阳理工学院OJ - 0003 - 多边形重心问题

题目要求我们给出按顺序给出的N点组成的多边形的面积和重心横纵坐标的和，考查计算几何。
我们有任意N边形的有向面积计算公式：

$$S'=\sum_{i=1}^{n-1}\frac{(x_i-x_0)(y_{i+1}-y_0)-(y_i-y_0)(x_{i+1}-x_0)}{2}$$
记按顺序排列的连续两点与选取某一点的三角形的有向面积：
$$S_{i,i+1}=\frac{(x_i-x_0)(y_{i+1}-y_0)-(y_i-y_0)(x_{i+1}-x_0)}{2},\text{i=1,2,...,n-1}$$
重心的X坐标：
$$X_{i,i+1}=\frac{x_0+x_i+x_{i+1}}{3}$$
重心的Y坐标：
$$Y_{i,i+1}=\frac{y_0+y_i+y_{i+1}}{3}$$
则任意N边形的重心计算公式：
$$X=\sum_{i=1}^{n-1}\left[\frac{S_{i,i+1}}{S'}X_{i,i+1}\right],Y=\sum_{i=1}^{n-1}\left[\frac{S_{i,i+1}}{S'}Y_{i,i+1}\right]$$
根据题意，我们变形方便计算：
$$S'=\sum_{i=1}^{n-1}S_{i,i+1}$$
$$sgm=\sum_{i=1}^{n-1}S_{i,i+1}(x_0+x_i+x_{i+1}+y_0+y_i+y_{i+1})$$
则：
$$S=|S'|，X+Y=\frac{sgm}{3S'}$$
在题目定义中，如果S=0，则重心公式的分母为0，并且题目规定这种情况为(0,0)，这是需要注意的。

#include <cstdio>
#include <cmath>
using namespace std;
const double eps=1e-7;
int main(){
    int n,m,j;
    double s,x0,y0,x1,y1,x2,y2,sgm,tems;
    scanf("%d",&n);
    for(int i=0;i<n;i++){
        scanf("%d",&m);
        if(m<3){
            for(j=0;j<m;j++)scanf("%lf%lf",&x1,&y1);
            printf("0.000 0.000\n");
        }
        else{s=0;sgm=0;
            scanf("%lf%lf%lf%lf",&x0,&y0,&x1,&y1);
            for(j=2;j<m;j++){scanf("%lf%lf",&x2,&y2);
                s+=(tems=(x1-x0)*(y2-y0)-(x2-x0)*(y1-y0));
                sgm+=tems*(x0+x1+x2+y0+y1+y2);
                x1=x2;y1=y2;
            }
            if(fabs(s)<eps)printf("0.000 0.000\n");
            else{
                printf("%.3lf %.3lf\n",fabs(s)/2,sgm/(3*s));
            }
        }
    }
}