bzoj1137: [POI2009]Wsp 岛屿【半平面交】

Description

Byteotia岛屿是一个凸多边形。城市全都在海岸上。按顺时针编号1到n。任意两个城市之间都有一条笔直的道路相连。道路相交处可以自由穿行。有一些道路被游击队控制了，不能走，但是可以经过这条道路与未被控制的道路的交点。问从城市1到n的最短距离。

Input

第一行正整数n m表示城市数和被控制的岛屿数（3≤n≤10^5 1≤m≤10^6）接下来n行每行两个整数x y表示每个城市的坐标。（|x|,|y|≤10^6）接下来m行描述一条不能走的道路（起点和终点）。数据保证有解。

Output

输出一个实数，最短距离，误差10^-5以内均算正确。

Sample Input

6 9

-12 -10

-11 6

-4 12

6 14

16 6

18 -2

3 4

1 5

2 6

2 3

4 5

3 5

1 3

3 6

1 6

Sample Output

42.000000000

解题思路：

注意可以在道路相交处换向，且所有点是顺时针给出的，所以1和n是相邻的。
容易发现，最短路径就是所有可行道路右侧（小连大）的半平面交，而注意一个点有用的道路就只有它连向编号最大的一条。这样边数就是$O(n)$的了。
为了方便，连上n->1,最后减去即可。

#include<bits/stdc++.h>
using namespace std;

int getint()
{
    int i=0,f=1;char c;
    for(c=getchar();(c!='-')&&(c<'0'||c>'9');c=getchar());
    if(c=='-')f=-1,c=getchar();
    for(;c>='0'&&c<='9';c=getchar())i=(i<<3)+(i<<1)+c-'0';
    return i*f;
}

const int N=200005,M=2000005;
struct point
{
    double x,y;
    point(){}
    point(double _x,double _y):x(_x),y(_y){}
    inline friend point operator - (const point &a,const point &b)
    {return point(a.x-b.x,a.y-b.y);}
    inline friend point operator + (const point &a,const point &b)
    {return point(a.x+b.x,a.y+b.y);}
    inline friend point operator * (const point &a,const double &b)
    {return point(a.x*b,a.y*b);}
    inline friend double operator * (const point &a,const point &b)
    {return a.x*b.y-a.y*b.x;}
    inline double dis(){return sqrt(x*x+y*y);}
}p[N];
struct line
{
    point st,v;double ang;
    line(){}
    line(point _st,point _v):st(_st),v(_v){ang=atan2(v.y,v.x);}
    inline friend bool operator < (const line &a,const line &b)
    {return a.ang<b.ang;}
}l[N];
int n,m,cnt;double ans;
pair<int,int>e[M];

inline bool onleft(const point &a,const line &b)
{
    return b.v*(a-b.st)>0;
}

inline point get_inter(const line &a,const line &b)
{
    double s1=(a.st-b.st)*b.v;
    double s2=b.v*(a.st+a.v-b.st);
    return a.st+a.v*(s1/(s1+s2));
}

int main()
{
    //freopen("lx.in","r",stdin);
    n=getint(),m=getint();
    for(int i=1;i<=n;i++)p[i].x=getint(),p[i].y=getint();
    for(int i=1;i<=m;i++)
    {
        e[i].first=getint(),e[i].second=getint();
        if(e[i].first>e[i].second)swap(e[i].first,e[i].second);
    }
    sort(e+1,e+m+1);
    for(int i=1,j=1;i<n;i++)
    {
        while(i==e[j+1].first)j++;
        if(i==e[j].first)
        {
            int to=n,k=j;
            while(e[k].second==to&&e[k].first==i)k--,to--;
            if(to>i)l[++cnt]=line(p[to],p[i]-p[to]);
        }
        else l[++cnt]=line(p[n],p[i]-p[n]);
    }m=cnt;
    l[++m]=line(p[1],p[n]-p[1]),ans-=(p[n]-p[1]).dis();
    sort(l+1,l+m+1);
    int head,tail=1;
    for(int i=2;i<=m;i++)
    {
        if(l[i].ang!=l[tail].ang)l[++tail]=l[i];
        else if(onleft(l[i].st,l[tail]))l[tail]=l[i];
    }
    m=tail,head=tail=1;
    for(int i=2;i<=m;i++)
    {
        while(head<tail&&!onleft(p[tail],l[i]))tail--;
        while(head<tail&&!onleft(p[head+1],l[i]))head++;
        l[++tail]=l[i];
        if(head<tail)p[tail]=get_inter(l[tail],l[tail-1]);
    }
    while(head<tail&&!onleft(p[tail],l[head]))tail--;
    p[head]=get_inter(l[head],l[tail]);
    ans+=(p[tail]-p[head]).dis();
    for(int i=head;i<tail;i++)ans+=(p[i+1]-p[i]).dis();
    printf("%0.8lf",ans);
    return 0;
}