BZOJ1137 洛谷3493 [POI2009]WSP-Island 岛屿 半平面交

题目链接

首先说明，我在BZOJ上过这题的时候洛谷由于没有SPJ，所以没法过，于是我不确定这份代码在洛谷上能不能过，标题写上洛谷题号只是为了有了SPJ之后其他人搜题解更方便。

题意：有一个n个点的凸多边形，任意两点之间都有一条路，道路长度为两点间距离，可以在道路的交界处走到另一条路，现在有m条路不能走，按照在凸多边形上的顺序给出n给点的坐标，求1号点到n号点的最短距离。n<=1e5,m<=1e6

题解：
好题。这个可以在任意交点拐弯的条件是做出此题的关键。

由于道路长度就是两点间距离，而且点的编号是按照凸多边形上的顺序给出的，所以对于每个点，它要向n走肯定直接会选择编号最大的那个点，如果1直接能到n那么答案肯定是1到n两点的距离，原因是两点间线段最短。由于可以直接在任意的路的交点拐弯，所以我们并不需要选择了一条路就走到底，画几个图会发现，只要到了与当前路相交的某条路，如果相交的那条路能通往编号更大的点，那么在这里拐弯一定会比继续沿着原路往下走更优，原因还是两点之间线段最短。这样我们会发现在原多边形内现在这些有用的路其实形成了一个凸壳，再加上n到1这条边便可以形成一个凸包，而这个凸包可以看作是由一些边的一部分拼在一起组成的，这些边都是与该点所能到达的编号最大的点所连的边，所以虽然总的边数很多，这样一来有用的边就很少了。由于n在每条边的哪一侧是确定的，所以我们拐弯只会向靠近n的那一侧拐弯，那么其实每条直线便是限制了一个行走的范围，不可能越过直线往远离n的方向走，那么就是一个一堆直线的半平面交问题了。最后求出由各种半平面交得到的凸多边形，用这个凸多边形的周长减去n到1的那条人为加进去的边的边长，就是最后最短距离的答案了。

代码：

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

int n,m,lst[100010],ji,h,t,gg;
double ans;
vector <int> v[100010];
struct point
{
    double x,y;
    point(double X=0,double Y=0)
    {
        x=X;
        y=Y;
    }
}a[200010],p[400010];
point operator + (point x,point y)
{
    return point(x.x+y.x,x.y+y.y);
}
point operator - (point x,point y)
{
    return point(x.x-y.x,x.y-y.y);
}
point operator * (point x,double y)
{
    return point(x.x*y,x.y*y);
}
struct line
{
    point x,y;
    double ang;
}l[400010],q[400010];
inline double dis(point x,point y)
{
    return sqrt((x.x-y.x)*(x.x-y.x)+(x.y-y.y)*(x.y-y.y));
} 
inline double chaji(point x,point y)
{
    return x.x*y.y-x.y*y.x;
}
inline int cmp(line x,line y)
{
    if(x.ang!=y.ang)
    return x.ang<y.ang;
    return chaji(y.y-x.x,y.x-x.x)>0;
}
point jiao(line x,line y)
{
    double z=chaji(x.y-x.x,y.x-x.x)/chaji(y.y-y.x,x.y-x.x);
    return y.x+(y.y-y.x)*z;
}
int main()
{
    scanf("%d%d",&n,&m);
    gg=n;
    for(int i=1;i<=n;++i)
    scanf("%lf%lf",&a[i].x,&a[i].y);
    for(int i=1;i<=m;++i)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        v[x].push_back(y);
        v[y].push_back(x);
    }
    for(int i=1;i<=n;++i)
    {
        sort(v[i].begin(),v[i].end());
        ji=v[i].size();
        for(int j=ji-1;j>=0;--j)
        {
            if(v[i][j]!=n-ji+j+1)
            {
                lst[i]=n-ji+j+1;
                break;
            }
        }
        if(!lst[i])
        lst[i]=n-ji;
    }
    if(lst[1]==n)
    {
        printf("%.10lf\n",dis(a[1],a[n]));
        return 0;
    }
    ji=0;
    for(int i=1;i<=n;++i)
    {
        if(lst[i]>i)
        {
            l[++ji].x=a[i];
            l[ji].y=a[lst[i]];
            l[ji].ang=atan2(l[ji].y.y-l[ji].x.y,l[ji].y.x-l[ji].x.x);
        }
    }
    l[++ji].x=a[n];
    l[ji].y=a[1];
    l[ji].ang=atan2(l[ji].y.y-l[ji].x.y,l[ji].y.x-l[ji].x.x);
    sort(l+1,l+ji+1,cmp);
    l[0].ang=-2e9;
    n=0;
    for(int i=1;i<=ji;++i)
    {
        if(l[i].ang!=l[i-1].ang)
        {
            l[++n]=l[i];
        }
    }
    h=1;
    t=1;
    q[1]=l[1];
    for(int i=2;i<=n;++i)
    {
        while(h<t&&chaji(l[i].x-p[t-1],l[i].y-p[t-1])>=0)
        --t;
        while(h<t&&chaji(l[i].x-p[h],l[i].y-p[h])>=0)
        ++h;
        q[++t]=l[i];
        if(h<t)
        p[t-1]=jiao(q[t-1],q[t]);
    }
    while(h<t&&chaji(q[h].x-p[t-1],q[h].y-p[t-1])>=0)
    --t;
    p[t]=jiao(q[t],q[h]);
    p[++t]=p[h];
    for(int i=h;i<t;++i)
    ans+=dis(p[i],p[i+1]);
    ans-=dis(a[1],a[gg]);
    printf("%.10lf\n",ans);
    return 0;
}