bzoj3051[WC2013]平面图（树上倍增+平面图转对偶图+扫描线）

简要题意：二维平面上n个点，点之间有一些连线，连线不在点之外的地方相交，将平面分为若干个区域。给出一些询问点对，问从这个点所在的区域走到另一个点所在的区域的最小代价。

题解：这道题首先可以把平面图转对偶图，这一点比较容易发现。然后对于从左指向右的线段，运用扫描线的思想，扫到左端点加入平衡树，扫到右端点从平衡树中删除。因为两线互不相交，所以相对位置不变。然后建立平面直角坐标系，y轴可以随意左右平移。对于每一条线段，我们使右端点正好在y轴上，然后选择两线段左端点x坐标比较大的作为比较直线，计算这条直线与两线段的交点的高低。然后就是一堆细节。真实的码农题。我写了3个namespace不然受不了

#include<bits/stdc++.h>
using namespace std;
const int N=1e5+,inf=1e9;
const double eps=1e-;
int n,m,q,num,bel[N<<];
struct edge{int v,nxt,w,id;}e[N<<];
struct Edge{int u,v,w;double o;}d[N<<];
struct point{
    int x,y,id;
    void init(int i)
    {
        double xx,yy;scanf("%lf%lf",&xx,&yy);
        x=(int)(xx*+eps),y=(int)(yy*+eps),id=i;
    }
}p[N],Q[N<<];
struct node{int x,y,t;};
bool operator<(node a,node b){return a.x!=b.x?a.x<b.x:a.t<b.t;}
namespace MST{
    int cnt,ecnt,f[N],hd[N],fa[N][],d[N][],dep[N];
    Edge E[N<<];
    edge e2[N<<];
    void adde(int x,int y,int z)
    {
        e2[++ecnt]=(edge){y,hd[x],z},hd[x]=ecnt;
        e2[++ecnt]=(edge){x,hd[y],z},hd[y]=ecnt;
    }
    int find(int x){return x==f[x]?x:f[x]=find(f[x]);}
    bool cmp(Edge a,Edge b){return a.w<b.w;}
    void add(int u,int v,int w){E[++cnt]=(Edge){u,v,w};}
    void dfs(int x)
    {
        for(int i=;i<=;i++)
        {
            fa[x][i]=fa[fa[x][i-]][i-];
            d[x][i]=max(d[x][i-],d[fa[x][i-]][i-]);
        }
        for(int i=hd[x];i;i=e2[i].nxt)
        if(e2[i].v!=fa[x][])
        dep[e2[i].v]=dep[x]+,fa[e2[i].v][]=x,d[e2[i].v][]=e2[i].w,dfs(e2[i].v);
    }
    int lca(int x,int y)
    {
        if(dep[x]<dep[y])swap(x,y);
        int ret=,t=dep[x]-dep[y];
        for(int i=;(<<i)<=t;i++)if(t&(<<i))ret=max(ret,d[x][i]),x=fa[x][i];
        if(x==y)return ret;
        for(int i=;~i;i--)
        if(fa[x][i]!=fa[y][i])ret=max(ret,max(d[x][i],d[y][i])),x=fa[x][i],y=fa[y][i];
        return max(ret,max(d[x][],d[y][]));
    }
    void work()
    {
        sort(E+,E+cnt+,cmp);
        for(int i=;i<=num;i++)f[i]=i;
        for(int i=;i<=cnt;i++)
        {
            int x=find(E[i].u),y=find(E[i].v);
            if(x!=y)adde(E[i].u,E[i].v,E[i].w),f[x]=y;
        }
        dfs();
        for(int i=;i<=q;i++)
        {
            if(bel[i]==||bel[i+q]==){puts("-1");continue;}
            int ans=lca(bel[i],bel[i+q]);
            if(ans==inf)puts("-1");else printf("%d\n",ans);
        }
    }
}
namespace ScanLine{
    struct cmp{
        bool operator()(int i,int j)
        {
            if(d[i].u==d[j].u)return d[i].o>d[j].o;
            int x=max(p[d[i].u].x,p[d[j].u].x);
            double x1=1.0*(p[d[i].v].y-p[d[i].u].y)*(x-p[d[i].v].x)/(p[d[i].v].x-p[d[i].u].x)+p[d[i].v].y;
            double x2=1.0*(p[d[j].v].y-p[d[j].u].y)*(x-p[d[j].v].x)/(p[d[j].v].x-p[d[j].u].x)+p[d[j].v].y;
            return x1>x2;
        }
    };
    set<int,cmp>S;
    node a[N<<];
    int cnt=;
    void work()
    {
        for(int i=;i<=(m<<)+;i++)
        if(p[d[i].u].x<p[d[i].v].x)a[++cnt]=(node){p[d[i].u].x,i,},a[++cnt]=(node){p[d[i].v].x,i,};
        for(int i=;i<=q;i++)a[++cnt]=(node){Q[i].x,i,},a[++cnt]=(node){Q[i+q].x,i+q,};
        sort(a+,a+cnt+);
        for(int i=;i<=cnt;i++)
        if(!a[i].t)S.erase(a[i].y);
        else if(a[i].t==)S.insert(a[i].y);
        else{
            p[n+]=(point){a[i].x,Q[a[i].y].y,};
            p[n+]=(point){a[i].x-,Q[a[i].y].y,};
            d[(m<<)+]=(Edge){n+,n+,,atan2(p[n+].y-p[n+].y,p[n+].x-p[n+].x)};
            S.insert((m<<)+);
            set<int,cmp>::iterator it=S.find((m<<)+);
            if(it!=S.begin())bel[a[i].y]=e[*--it].id;else bel[a[i].y]=;
            S.erase((m<<)+);
        }
    }
}
namespace Graph{
    int cnt=,hd[N],nxt[N<<];
    vector<pair<double,int> >G[N];
    point s=(point){inf,inf,};
    void adde(int x,int y,int z)
    {
        e[++cnt]=(edge){y,hd[x],z,-},hd[x]=cnt;
        e[++cnt]=(edge){x,hd[y],z,-},hd[y]=cnt;
    }
    void dosort(int x)
    {
        for(int i=hd[x];i;i=e[i].nxt)
        G[x].push_back(make_pair(atan2(p[e[i].v].y-p[x].y,p[e[i].v].x-p[x].x),i));
        sort(G[x].begin(),G[x].end());
        int sz=G[x].size();
        for(int i=;i<sz;i++)nxt[G[x][i].second^]=G[x][(i+)%sz].second;
    }
    void find(int x){for(int i=x;e[i].id<;i=nxt[i])e[i].id=num;}
    void work()
    {
        cnt=;
        for(int i=;i<=n;i++)
        {
            p[i].init(i);
            if(s.x>p[i].x||s.x==p[i].x&&s.y>p[i].y)s=p[i];
        }
        for(int i=,x,y,z;i<=m;i++)
        {
            scanf("%d%d%d",&x,&y,&z);
            d[i<<]=(Edge){x,y,z,atan2(p[y].y-p[x].y,p[y].x-p[x].x)};
            d[(i<<)+]=(Edge){y,x,z,atan2(p[x].y-p[y].y,p[x].x-p[y].x)};
            adde(x,y,z);
        }
        scanf("%d",&q);
        for(int i=;i<=q;i++)Q[i].init(i),Q[i+q].init(i+q);
        for(int i=;i<=n;i++)dosort(i);
        num=;
        find(G[s.id][].second);
        for(int i=;i<=cnt;i++)if(e[i].id<)++num,find(i);
        for (int i=;i<=cnt;i+=)
        if(e[i].id!=&&e[i^].id!=)MST::add(e[i].id,e[i^].id,e[i].w);
        else MST::add(e[i].id,e[i^].id,inf);
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    Graph::work();
    ScanLine::work();
    MST::work();
}