[bzoj4445] [Scoi2015]小凸想跑步

题目大意

给定一个凸多边形（没有三点共线），每个端点逆时针编号为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;
}