POJ2125：Destroying The Graph

题目链接

【题意】

给一张有向图，现在要选择一些点，删掉图中的所有边。具体操作为：选择点i，可
以选择删除从i 出发的所有有向边或者进入i 的所有有向边，分别有个代价ini 和
outi，求最小的代价删掉所有边。并输出删除方案。

【分析】

一条边不是被入点删掉，就是出点。就可以拆点构二分图，然后最小割。
方案？
考虑残量网络，从S做BFS，那么对于左侧被选中的点，就是level值为0，及由S连向它的边满流了，对于右侧被选中的点，就是level值不为0的点，及由S可以流向它，那么一定没有删掉这个点所在的某条边中的另外一个点。

注意一下边的方向。。

【代码】

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <queue>
#include <cmath>
#include <cstring>
#define N 205
#define M 10405
#define INF 1e9+1
using namespace std;
typedef long long ll;

int read()
{
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
    return x*f; 
}

int n,m,cnt=1,S,T,ans,sum;
int b[M],p[N],nextedge[M],w[M],cur[N];
int Level[N];
bool Flag[N];

class Node{
    public:
        int x,y;
}E[N];

void Add(int x,int y,int z)
{
    cnt++;
    b[cnt]=y;
    nextedge[cnt]=p[x];
    p[x]=cnt;
    w[cnt]=z;
}

void Anode(int x,int y,int z){
    Add(x,y,z);Add(y,x,0);
}

bool Bfs()
{
    queue<int>q;
    q.push(S);
    for(int i=0;i<=T;i++) Level[i]=0;
    Level[S]=1;
    while(!q.empty())
    {
        int k=q.front();q.pop();
        for(int i=p[k];i;i=nextedge[i])
        {
            int v=b[i],f=w[i];
            if(!Level[v]&&f)
            {
                Level[v]=Level[k]+1;
                q.push(v);
            }
        }
    }
    return Level[T];
}

int Dfs(int x,int maxf)
{
    if(x==T||!maxf) return maxf;
    int rtn=0;
    for(int i=cur[x];i&&maxf>rtn;i=nextedge[i])
    {
        int v=b[i],f=w[i];
        if(Level[v]==Level[x]+1&&f)
        {
            f=Dfs(v,min(maxf-rtn,f));
            w[i]-=f;w[i^1]+=f;
            if(w[i]>0) cur[x]=i;
            rtn+=f;
        }
    }
    if(!rtn) Level[x]=0;
    return rtn;
}

void Dinic()
{
    while(Bfs())
    {
        for(int i=0;i<=T;i++) cur[i]=p[i];
        ans+=Dfs(0,INF);
    }
}

void Input_Init()
{
    n=read(),m=read();T=n<<1|1;
    for(int i=1;i<=n;i++) Anode(0,i,read());
    for(int i=1;i<=n;i++) Anode(i+n,T,read());
    for(int i=1;i<=m;i++)
    {
        static int x,y;
        x=read(),y=read();
        Anode(y,x+n,INF);
    }
}

void Solve()
{
    Bfs();
    printf("%d\n",ans);
    for(int i=1;i<=n;i++) if(!Level[i]) sum++,E[sum].x=i,E[sum].y=1;
    for(int i=n+1;i<=n<<1;i++) if(Level[i]) sum++,E[sum].x=i-n,E[sum].y=0;
    printf("%d\n",sum);
    for(int i=1;i<=sum;i++)
    {
        printf("%d ",E[i].x);
        if(E[i].y) printf("+\n");
        else printf("-\n");
    }
}

int main()
{
    Input_Init();
    Dinic();
    Solve();
    return 0;
}