【半平面交】洛谷P4196 模板

题目在这里
题目求解的是很多多边形的交的面积。
其实很多多边形就是很多线段，由于是逆时针给出的点，所以也就是很多线段的左半平面。
故直接半平面交即可。
求出凸壳后用三角剖分求出面积即可。
本代码是复制洛谷题解里面的，当作模板使用。

#include <bits/stdc++.h>
using namespace std;
const int maxn = 505;
const double eps = 1e-8;
int n;
int tot,cnt,head,tail;
struct Point
{
    double x,y;
    friend Point operator - (Point a,Point b)
	{
        Point t;
        t.x=a.x-b.x;t.y=a.y-b.y;
        return t;
    }
    friend Point operator + (Point a,Point b)
	{
        Point t;
        t.x=a.x+b.x;t.y=a.y+b.y;
        return t;
    }
    friend double operator * (Point a,Point b)
	{
        return a.x*b.x+a.y*b.y;
    }
    friend double operator ^ (Point a,Point b) // a^b 表示叉积 
	{
        return a.x*b.y-b.x*a.y;
    }
    friend Point operator * (double k,Point b)
	{
        Point t;
        t.x=k*b.x;t.y=k*b.y;
        return t; 
    }
}b[maxn];
int sign(double x)
{
    return fabs(x)<=eps ?  0 : (x>0 ? 1 : -1);
}
struct Line
{
    Point p,v;
    double poa;
    friend bool operator < (Line x,Line y)
	{
        return sign(x.poa-y.poa)==0 ? sign((x.v-x.p) ^ (y.v-x.p)) > 0 : sign(x.poa-y.poa)<0;
        //因为我是向量左侧求交，所以极角相同时靠左的更优，把优的放在后面，方便之后的操作 ，可以画图体会一下 
    }
}a[maxn],q[maxn];
Point inter(Line a,Line b)//求交点 
{
    Point p1=a.p,p2=b.p,v1=a.v,v2=b.v;
    v1=v1-p1;v2=v2-p2;
    Point u=p2-p1;
    Point p=p2+((u^v1)/(v1^v2))*v2;
    return p;
}
bool pd(Line i,Line j,Line k)
{
    Point p=inter(i,j);
    return sign((k.v-k.p) ^ (p-k.p))<0;
}
void Half_Plane()
{
    sort(a+1,a+tot+1);//排序 
    for(int i=1;i<=tot;i++) 
	{
        if(sign(a[i].poa-a[i-1].poa)!=0) cnt++;
        a[cnt]=a[i];//因为排过序，即使极角相同，后面的也比前面的优 
    }
    head=1;tail=0;
    q[++tail]=a[1];q[++tail]=a[2];
    for(int i=3;i<=cnt;i++)
	{
        while(head<tail&&pd(q[tail-1],q[tail],a[i])) tail--;//维护双端队列 
        while(head<tail&&pd(q[head+1],q[head],a[i])) head++;
        q[++tail]=a[i];
    }
    while(head<tail&&pd(q[tail-1],q[tail],q[head])) tail--;
    while(head<tail&&pd(q[head+1],q[head],q[tail])) head++;
    q[tail+1]=q[head];
    tot=0;
    for(int i=head;i<=tail;i++) b[++tot]=inter(q[i],q[i+1]);
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
	{
        int k;
        scanf("%d",&k);
        for(int j=1;j<=k;j++) scanf("%lf%lf",&b[j].x,&b[j].y);
        b[k+1]=b[1];
        for(int j=1;j<=k;j++) a[++tot].p=b[j],a[tot].v=b[j+1];
    }
    for(int i=1;i<=tot;i++) a[i].poa=atan2(a[i].v.y-a[i].p.y,a[i].v.x-a[i].p.x);
    Half_Plane();
    b[tot+1]=b[1];
    double ans = 0;
    if(tot>2) 
	{
		for(int i=1;i<=tot;i++) ans+=(b[i]^b[i+1]); //b[i]为凸壳上的点 
	}
    ans=fabs(ans)/2.0;
    printf("%.3lf\n",ans);
    return 0;
}