题目大意
给定一个凸多边形(没有三点共线),每个端点逆时针编号为0——n-1,现在在多边形内随机一个点,求这个点与0、1号点围成三角形面积小于等于与其它相邻点围成的三角形的概率。
n≤100000 坐标值在int范围
分析
设这个点为(x,y),两个端点为(x1,y1),(x2,y2),然后用叉积来表示三角形面积。
那么对于除0,1之外的相邻点,都可以列出一个不等式来。整理可以得到ax+by+c≤0的形式。
那么就可以上半平面交了!并且这道题的特殊性使得可以不用判无解。。
注意:(x,y)必须在凸多边形内,所以要以点1、0又作为一个半平面加进来。而且可能有b等于0的情况。
时间复杂度O(nlogn)
人生第一道半平面交,所以代码有点丑
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
const int N=100010;
typedef long long LL;
typedef double db;
const db eps=1e-10,inf=2e9;
int n,h,t,m;
db s;
struct P
{
db x,y;
P (db _x=0,db _y=0)
{
x=_x; y=_y;
}
}p1[N],p2[N];
P operator + (P a,P b)
{
return P(a.x+b.x,a.y+b.y);
}
P operator - (P a,P b)
{
return P(a.x-b.x,a.y-b.y);
}
P operator * (P a,db b)
{
return P(a.x*b,a.y*b);
}
db operator * (P a,P b)
{
return a.x*b.x+a.y*b.y;
}
db operator ^ (P a,P b)
{
return a.x*b.y-a.y*b.x;
}
db arg(P a)
{
return atan2(a.y,a.x);
}
bool equ(db x,db y)
{
return fabs(x-y)<=eps;
}
struct L
{
P p,v;
L () {}
L (P _p,P _v)
{
p=_p; v=_v;
}
}pl[N],D[N];
bool cmp(L a,L b)
{
return arg(a.v)<arg(b.v);
}
bool in(P a,L b)
{
return ((a-b.p)^b.v)>=0;
}
P ict(L a,L b)
{
return b.p+b.v*((a.v^(a.p-b.p))/(a.v^b.v));
}
char c;
int read()
{
int x=0,sig=1;
for (c=getchar();c<'0' || c>'9';c=getchar()) if (c=='-') sig=-1;
for (;c>='0' && c<='9';c=getchar()) x=x*10+c-48;
return x*sig;
}
db calc()
{
pl[m++]=L(P(-inf,-inf),P(0,1));
pl[m++]=L(P(-inf,inf),P(1,0));
pl[m++]=L(P(inf,inf),P(0,-1));
pl[m++]=L(P(inf,-inf),P(-1,0));
sort(pl,pl+m,cmp);
n=1;
for (int i=1;i<m;i++)
{
if (!equ(arg(pl[n-1].v),arg(pl[i].v))) pl[n++]=pl[i];
else if (in(pl[i].p,pl[n-1])) pl[n-1]=pl[i];
}
h=1; t=2;
D[1]=pl[0]; D[2]=pl[1];
for (int i=2;i<n;i++)
{
for (;h<t && !in(ict(D[t-1],D[t]),pl[i]);t--);
for (;h<t && !in(ict(D[h],D[h+1]),pl[i]);h++);
D[++t]=pl[i];
}
for (;h<t && !in(ict(D[t],D[t-1]),D[h]);t--);
m=0;
p2[0]=ict(D[t],D[h]);
for (int i=h;i<t;i++) p2[++m]=ict(D[i],D[i+1]);
p2[++m]=ict(D[t],D[h]);
db ret=0.0;
for (int i=0;i<m;i++) ret+=(p2[i]^p2[i+1]);
return ret;
}
int main()
{
n=read();
for (int i=0;i<n;i++)
{
int x=read(),y=read();
p1[i]=P(x,y);
}
p1[n]=p1[0];
for (int i=0;i<n;i++) s+=(p1[i]^p1[i+1]);
pl[m++]=L(p1[1],p1[0]-p1[1]);
for (int i=1;i<n;i++)
{
db a=p1[0].y-p1[1].y-p1[i].y+p1[i+1].y,b=p1[1].x-p1[0].x-p1[i+1].x+p1[i].x,c=(p1[0]^p1[1])-(p1[i]^p1[i+1]);
pl[m].p=P(b?0:-c/a,b?-c/b:0); pl[m++].v=P(b,-a);
}
printf("%.4lf\n",calc()/s);
return 0;
}