【BZOJ2960】跨平面（对偶图+朱刘算法）

Description

Input

第一行两个整数n和m，表示点与线段的数目。
接下来n行，每行两个整数x和y，表示第i个点的坐标，点从1到n编号。
接下来m行，每行四个整数p,q,V1和V2，表示存在一条从第p个点连向第q个点的线段，激活p->q这个方向的费用为V1，另一个方向费用为V2。
保证若两条线段相交，则交点是它们的公共端点。

Output

输出一行一个正整数，表示最小总激活费用。

Sample Input

4 5

0 0

1 0

0 1

1 1

1 2 0 0

1 3 0 3

2 3 1 0

2 4 2 0

4 3 0 0

Sample Output

3

HINT

　　对于100%的数据，n≤3000，区域数不超过1000，点坐标绝对值不超过1W，每条边激活费用不超过100。

题解：真是一道神题，学到了不少新姿势√
首先，不知道对偶图的请自行百度。因为这道题要用到平面图转对偶图。首先，找一条没有经过的单向边，开始深搜，直到回到自己的点为止，这样就围成了一个区域，在回溯时给经过的每条边都标上这个区域的编号。然后求它的对偶图。这就转成了一个有向图，在上面求最小树形图即可。

代码如下：

#include<iostream>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<queue>
#include<math.h>
#define ll long long
#define inf 0x7f7f7f7f
using namespace std;
namespace Graph
{
    #define N 2005
    #define M 50005
    int X[M],Y[M],val[M],c2;  
    int nex[M],to[M],hd[N],cnt,tot;  
    int in[N],fa[N],n;  
    bool vis[N],inc[N];
    void add(int x,int y,int v)
    {
        X[++c2]=x,Y[c2]=y,val[c2]=v;
    }
    void addedge(int x,int y)
    {
        to[++cnt]=y,nex[cnt]=hd[x],hd[x]=cnt;
    }
    int dfs(int x)
    {
        inc[x]=1;
        if(inc[fa[x]]){tot+=in[x];return fa[x];}
        else 
        {
            int tmp=dfs(fa[x]);
            if(tmp && x!=tmp)
            {
                tot+=in[x];
                return tmp;
            }
            else if(x==tmp)
            {
                tot+=in[x];
                return 0;
            }
        }
        inc[x]=0;
        return 0;
    }
    int check()  
    {  
        for(int i=0;i<=n;i++)  
        if(!vis[i]) in[i]=inf,hd[i]=0,inc[i]=1;  
        for(int i=1;i<=c2;i++)  
        {  
            if(!vis[X[i]] && !vis[Y[i]])  
            if(in[Y[i]]>val[i])  
            {  
                in[Y[i]]=val[i];  
                fa[Y[i]]=X[i];  
            }   
        }  
        int ret=0;  
        for(int i=1;i<=n;i++)  
        if(!vis[i]) ret+=in[i];  
        cnt=0;  
        for(int i=1;i<=n;i++)  
        if(!vis[i]) addedge(fa[i],i);  
        queue<int> q;  
        q.push(0);  
        while(!q.empty())  
        {  
            int x=q.front();q.pop();  
            inc[x]=0;  
            for(int i=hd[x];i;i=nex[i]) q.push(to[i]);   
        }  
        for(int i=0;i<=n;i++)  
        if(!vis[i] && inc[i])  
        {  
            memset(inc,0,sizeof(inc[0])*(n+1));  
            dfs(i);  
            return -1;  
        }  
        return ret;  
    }  
    int solve()  
    {  
        int ans,sum=0;  
        for(int i=1;i<=c2;i++) sum+=val[i];  
        for(int i=1;i<=n;i++) add(0,i,sum+1);  
        tot=0;  
        while((ans=check())==-1)  
        {  
            for(int i=1;i<=c2;i++)  
            {  
                if(inc[X[i]] && !vis[Y[i]] && !inc[Y[i]]) X[i]=n+1;  
                else if(!vis[X[i]] && !inc[X[i]] && inc[Y[i]])  
                {  
                    val[i]-=in[Y[i]];  
                    Y[i]=n+1;  
                }  
            }  
            for(int i=0;i<=n;i++)  
            if(inc[i]) vis[i]=1,inc[i]=0;  
            n++;  
        }  
        return ans+tot-sum-1;  
    }  
    #undef N
    #undef M
}
namespace Plane
{
    #define N 3005
    #define M 15000
    struct Point
    {
        int x,y;
        Point(){}
        Point(int _x,int _y):x(_x),y(_y){}
    }a[N];
    struct Line
    {
        int x,y,id;
        bool v;
        double k;
        inline bool operator<(const Line& rhs) const {return k<rhs.k;}
    }l[M];
    int nex[M],to[M],num[M],hd[N],tot;
    bool v[M],vis[N];
    int val[M][2],X[M][2];
    void add(int x,int y,int nm,bool vv)
    {
        to[++tot]=y,num[tot]=nm,v[tot]=vv,nex[tot]=hd[x],hd[x]=tot;
    }
    int dfs(int x,int fa)
    {
        if(vis[x]) return ++Graph::n;
        for(int i=hd[x];i;i=nex[i])
        {
            if(to[i]==fa)
            {
                i=nex[i];
                if(i) return X[num[i]][v[i]]=dfs(to[i],x);
                else return X[num[hd[x]]][v[hd[x]]]=dfs(to[hd[x]],x); 
            }
        }
        return 0;
    }
    void solve()
    {
        int n,m;
        scanf("%d%d",&n,&m);
        for(int i=1;i<=n;i++) scanf("%d%d",&a[i].x,&a[i].y);
        for(int i=1;i<=m;i++)
        {
            scanf("%d%d%d%d",&l[i*2-1].x,&l[i*2-1].y,&val[i][0],&val[i][1]);
            l[i*2-1].id=i,l[i*2-1].v=0;
            l[i*2-1].k=atan2(a[l[i*2-1].y].y-a[l[i*2-1].x].y,a[l[i*2-1].y].x-a[l[i*2-1].x].x);
            l[i*2].x=l[i*2-1].y,l[i*2].y=l[i*2-1].x,l[i*2].id=i,l[i*2].v=1;
            l[i*2].k=atan2(a[l[i*2-1].x].y-a[l[i*2-1].y].y,a[l[i*2-1].x].x-a[l[i*2-1].y].x);
        }
        sort(l+1,l+m*2+1);
        for(int i=1;i<=m*2;i++) add(l[i].x,l[i].y,l[i].id,l[i].v);  
        for(int i=1;i<=n;i++)  
        {  
            vis[i]=1;  
            for(int j=hd[i];j;j=nex[j])  
            if(!X[num[j]][v[j]]) X[num[j]][v[j]]=dfs(to[j],i);  
            vis[i]=0;  
        }  
        for(int i=1;i<=m;i++)  
        {  
            if(val[i][0]) Graph::add(X[i][0],X[i][1],val[i][0]);  
            if(val[i][1]) Graph::add(X[i][1],X[i][0],val[i][1]);  
        }  
    }
    #undef N
    #undef M
}
int main()
{
    Plane::solve();
    printf("%d\n",Graph::solve()); 
    return 0;
}