题目大意:按顺时针顺序给出一个N边形,求N边形的核的面积
其实是做POI2009的一道题时被告知要用半平面交于是才找这道题练练手的..
不由得又想起2016年北大夏令营第二场被这道题支配的恐惧.....
当时我们学校好像只有alone_wolf一个人调出来了..
不由得又想起2016年北大夏令营第二场被这道题支配的恐惧.....
当时我们学校好像只有alone_wolf一个人调出来了..
#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();
}