hihoCoder1083积分

简述

　　我太蒻了，这题交了好几遍才$\color{green}AC$。
　　首先要知道重心$(x_g,y_g)$的求法。
　　

$$x_g=\frac{\iint_D xdxdy}{\iint_D dxdy}\\y_g=\frac{\iint_D ydxdy}{\iint_D dxdy}$$
　　而 $\iint_D dxdy=S$
　　而积分是可加的，所以把原多边形剖分成三角形来做，原式子可以化为
　　 $$\iint(x+y)dxdy=\iint xdxdy+\iint ydxdy\\=\sum(xg+yg)S$$
　　那么这题最后成了每个三角形的重心横纵坐标之和乘以其面积。
　　凹多边形也是可以剖成三角形的，做法是先把所有的有向面积加起来，最后如果有向面积的和为负，就要取个相反数。

代码

//计算几何
#include <cstdio>
#include <algorithm>
#include <cmath>
#define maxn 100000
using namespace std;
struct point{double x, y;}pt[maxn];
struct vec{point a; double x, y;};
double operator*(vec v1, vec v2){return v1.x*v2.y-v2.x*v1.y;}
vec operator-(point p1, point p2)
{return (vec){(point){0.0,0.0},p2.x-p1.x,p2.y-p1.y};}
double x[maxn], y[maxn];
int main()
{
    int N, i, T;
    vec v1, v2;
    double xg, yg, ans=0, t=0;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&N);
        ans=t=0;
        for(i=1;i<=N;i++)scanf("%lf%lf",x+i,y+i);
        for(i=3;i<=N;i++)
        {
            v1=(point){x[i-1],y[i-1]}-(point){x[1],y[1]};
            v2=(point){x[i],y[i]}-(point){x[1],y[1]};
            xg=(x[1]+x[i-1]+x[i])/3, yg=(y[1]+y[i-1]+y[i])/3;
            t+=v1*v2/2.0;
            ans+=v1*v2/2.0*(xg+yg);
        }
        if(t<0)ans=-ans;
        printf("%.2lf\n",ans);
    }
    return 0;
}