POJ1279 Art Gallery

題目大意：按順時針順序給出一個N邊形，求N邊形的核的面積

其實是做POI2009的一道題時被告知要用半平面交於是才找這道題練練手的..
不由得又想起2016年北大夏令營第二場被這道題支配的恐懼.....
當時我們學校好像只有alone_wolf一個人調出來了..

那這道題的話...就是一個半平面交模板題咯
半平面交知識點->這裏
學習+寫代碼+調代碼 基本花了我一天的時間QAQ
我還是太弱了

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 10010
using namespace std;
struct node
{
	double x,y;
	node(){}
	node(double xx,double yy){x=xx;y=yy;}
	friend node operator +(node x,node y){
		return node(x.x+y.x,x.y+y.y);
	}
	friend node operator -(node x,node y){
		return node(x.x-y.x,x.y-y.y);
	}
	friend node operator *(node x,double y){
		return node(x.x*y,x.y*y);
	}
};
struct line{node a,b;double ang;};
double xj(const node &x,const node &y){return x.x*y.y-x.y*y.x;}
node a[N];
line L[N];
line q[N];int h,t;node p[N];
bool cmp(line x,line y)
{
	if(x.ang!=y.ang) return x.ang<y.ang;
	return xj(y.b-x.a,y.a-x.a)>0;
}
bool left(node x,line l)
{
	return xj(l.a-x,l.b-x)>0; 
}
node jiao(line l1,line l2)
{
	double b=xj(l1.b-l1.a,l2.a-l1.a)/xj(l2.b-l2.a,l1.b-l1.a);
	return l2.a+(l2.b-l2.a)*b;
}
void print(node x){cout<<x.x<<' '<<x.y<<' ';}
void print(line x){print(x.a);cout<<'*';print(x.b);cout<<endl;}
void doit()
{
	int n;
	scanf("%d",&n);
	int i,j,x,y;
	for(i=1;i<=n;i++)
	scanf("%lf%lf",&a[i].x,&a[i].y);
	a[n+1]=a[1];
	for(i=1;i<=n;i++)
	{
		L[i].a=a[i];L[i].b=a[i+1];
		L[i].ang=atan2(L[i].b.y-L[i].a.y,L[i].b.x-L[i].a.x);
	}
	sort(L+1,L+n+1,cmp);
	int nn=1;
	for(i=2;i<=n;i++)
	if(L[i].ang!=L[i-1].ang) nn++,L[nn]=L[i];
	n=nn;
	h=t=1;q[1]=L[1];
	for(i=2;i<=n;i++)
	{
		while(h<t&&left(p[t-1],L[i])) t--;
		while(h<t&&left(p[h],L[i])) h++;
		t++;q[t]=L[i];
		if(h<t) p[t-1]=jiao(q[t-1],q[t]);
	}
	while(h<t&&left(p[t-1],q[h])) t--;
	p[t]=jiao(q[h],q[t]);
	double ans=0;
	for(i=h+1;i<t;i++)
	ans+=xj(p[i]-p[h],p[i+1]-p[h]);
	printf("%.2lf\n",ans/2);
}
int main()
{
	int T;
	scanf("%d",&T);
	while(T--)
	doit();
}