poj-1279 Art Gallery

题意：

给出一个n个点的多边形，求多边形核的面积；

多边形的核即为多边形内部的点集，使其中的点可以看到多边形的任意顶点和边；

n<=1500；

题解：

这应该也算是模板题的一种吧；

求多边形核的方法现在是O(nlogn)的；

将所有的直线按极角排序(math库中atan2()函数)，然后线性的在双端队列里面搞搞；

这样求出了半平面交的边集，也顺带求出了点集；

求面积的话就直接叉积扫一遍就好了；

但是这道题不一样的地方是没有说顶点按什么序给出；

所以可以先顺时针求一边，如果半平面交无解再逆时针求，都无解那多边形不存在核，答案为0；

半平面交无解就是求出的边集小于三；

为了省事我两次直接粘的代码= =，可以写的函数化一些以缩短码长；

代码：

#include<math.h>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define N 110000
using namespace std;
typedef long long ll;
const double EPS=1e-10;
const double INF=1e100;
struct Point
{
	double x,y;
	Point(){}
	Point(const double &_,const double &__):x(_),y(__){}
	friend Point operator + (const Point &a,const Point &b)
	{
		return Point(a.x+b.x,a.y+b.y);
	}
	friend Point operator - (const Point &a,const Point &b)
	{
		return Point(a.x-b.x,a.y-b.y);
	}
	friend double operator * (const Point &a,const Point &b)
	{
		return a.x*b.x+a.y*b.y;
	}
	friend double operator ^ (const Point &a,const Point &b)
	{
		return a.x*b.y-a.y*b.x;
	}
	friend Point operator * (const double a,const Point &b)
	{
		return Point(a*b.x,a*b.y);
	}
}poi[N],p[N];
struct Line
{
	Point p,v;
	double alpha;
	Line(){};
	Line(const Point &_,const Point &__):p(_),v(__-_){alpha=atan2(v.y,v.x);}
	friend bool operator < (const Line &a,const Line &b)
	{
		return a.alpha<b.alpha;
	}
	friend Point getP(const Line &a,const Line &b)
	{
		Point u=a.p-b.p;
		double temp=(b.v^u)/(a.v^b.v);
		return a.p+temp*a.v;
	}
}l[N];
int q[N],st,en;
bool Onleft(const Line &l,const Point &a)
{
	return (l.v^(a-l.p))>EPS;
}
void HPI(int n)
{
	sort(l+1,l+n+1);
	q[st=en=1]=1;
	for(int i=2;i<=n;i++)
	{
		while(st<en&&!Onleft(l[i],p[en]))
			en--;
		while(st<en&&!Onleft(l[i],p[st+1]))
			st++;
		if(fabs(l[q[en]].alpha-l[i].alpha)>EPS)
			q[++en]=i;
		else
			q[en]=Onleft(l[i],l[q[en]].p)?q[en]:i;
		p[en]=getP(l[q[en-1]],l[q[en]]);
	}
	while(st<en&&!Onleft(l[q[st]],p[en]))
		en--;
}	
int main()
{
	int c,T,n,i,j;
	double ans;
	scanf("%d",&T);
	for(c=1;c<=T;c++)
	for(c=1;c<=T;c++)
	{
		scanf("%d",&n);
		for(i=1;i<=n;i++)
			scanf("%lf%lf",&poi[i].x,&poi[i].y);
		for(i=2;i<=n;i++)
		{
			l[i]=Line(poi[i],poi[i-1]);
		}
		l[1]=Line(poi[1],poi[n]);
		HPI(n);
		bool flag=1;
		if(en-st+1<3)
		{
			for(i=2;i<=n;i++)
				l[i]=Line(poi[i-1],poi[i]);
			l[1]=Line(poi[n],poi[1]);
			HPI(n);
			if(en-st+1<3)
				flag=0;
		}
		if(flag==0)
		{
			puts("0.00");
			continue;
		}
		p[st]=getP(l[q[st]],l[q[en]]);
		for(i=st,ans=0;i<en;i++)
			ans+=p[i]^p[i+1];
		ans+=p[en]^p[st];
		if(ans<0)	ans=-ans;
		printf("%.2f\n",ans/2+EPS);
	}
	return 0;
}