计算几何 - Icebergs - Gym 102501F

计算几何 - Icebergs - Gym 102501F

题意：

$给定n个不重叠的多边形，每个多边形给定m_i个点的坐标，$

$计算出n个多边形的面积之和。$

Examples
Input

1
4
0 0
1 0
1 1
0 1

Output

1

Input

2
5
98 35
79 90
21 90
2 36
50 0
3
0 0
20 0
0 20

Output

6100

数据范围：

$1\le N\le 1000，3\le P\le 50，0\le x,y\le 10^6$

---

注意点：

$模板函数是求的是有向面积，绝对值应当在最后返回的时候取，无论是求凹多边形和凸多边形。$

$且需要保证点的顺序是按照逆时针/顺时针给出的，否则无法唯一确定一个多边形。$

$最后输出四舍五入，应该先用double类型求和，最后再强转为longlong，$

$若直接用longlong型进行求和，每加一次都可能会损失一次精度。$

代码：

#include<iostream>
#include<cmath>
#include<cstdio>
#include<vector>
#include<algorithm>

#define ll long long

using namespace std;

const double eps=1e-10;
const double pi=acos(-1.0);

struct Point 
{
    double x,y;
    Point(double x=0,double y=0) : x(x), y(y) {}
};

//点与向量
typedef Point Vector;

Vector operator + (Vector A,Vector B) { return Vector(A.x+B.x,A.y+B.y); }

Vector operator - (Point A,Point B) { return Vector(A.x-B.x,A.y-B.y); }

Vector operator * (Vector A,double p) { return Vector(A.x*p,A.y*p); }

Vector operator / (Vector A,double p) { return Vector(A.x/p,A.y/p); }

double Cross(Vector A, Vector B) { return A.x*B.y - A.y*B.x; }  //若A到B逆时针则为正，否则为负

double PolygonArea(Point p[], int n) //多边形的有向面积
{
    double area=0;
    for(int i=1;i<n-1;i++) area+=Cross(p[i]-p[0],p[i+1]-p[0]);
    return fabs(area/2);
}

Point V[1010][60];

int main()
{
    int m,n;
    scanf("%d",&m);
    double area=0;
    for(int i=0;i<m;i++)
    {
        scanf("%d",&n);
        for(int j=0;j<n;j++) scanf("%lf%lf",&V[i][j].x,&V[i][j].y);
        area+=PolygonArea(V[i],n);
    }
    printf("%lld\n",(ll)area);
    return 0;
}