[BZOJ1797][Ahoi2009]Mincut 最小割 做题笔记

题目来源：http://www.lydsy.com/JudgeOnline/problem.php?id=1797
这题考察判断网络流的可行割与必须割。
引用一下lydrainbowcat的结论

• 最小割的必须边
  
  • 一定在最小割中的边、扩大容量后能增大最大流的边， ① 满流；② 残余网络中S能到入点、出点能到T。 从S开始DFS、T开始反向DFS，标记到达的点，然后枚举满流边即可。
• 最小割的可行边
  
  • 被某一种最小割的方案包含的边， ① 满流；② 删掉之后在残余网络中找不到u到v的路径。 在残余网络中tarjan求SCC，(u,v)两点在同一SCC中说明残余网络中存在u到v路径。

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <algorithm>
#include <queue>
using namespace std;
const int N=5005,M=120005,inf=0x3fffffff;
int head[N],d[N],ver[M<<1],e[M<<1],next[M<<1];
int dfn[N],low[N],hav[N],bel[N],q[N<<1];
bool inq[N];
int tot=1,n,m,s,t,cnt=0,maxflow=0,top=0,scc=0;

void add (int u,int v,int w) {
    ver[++tot]=v;e[tot]=w;next[tot]=head[u];head[u]=tot;
    ver[++tot]=u;e[tot]=0;next[tot]=head[v];head[v]=tot;
}

bool bfs () {
    queue<int> q;
    memset(d,0,sizeof(d));
    d[s]=1;q.push(s);
    while (!q.empty()) {
        int now=q.front(); q.pop();
        for (int i=head[now];i;i=next[i]) 
            if (!d[ver[i]]&&e[i]) {
                d[ver[i]]=d[now]+1;
                q.push(ver[i]);
                if (ver[i]==t) return 1;
            }
    }
    return 0;
}

int dinic (int x,int f) {
    int rest=f;
    if (x==t) return f;
    for (int i=head[x];i&&rest;i=next[i]) 
        if (e[i]&&d[ver[i]]==d[x]+1) {
            int now=dinic(ver[i],min(rest,e[i]));
            if (!now) d[ver[i]]=0;
            rest-=now;
            e[i]-=now;
            e[i^1]+=now;
        }
    return f-rest;
}

void Tarjan (int x) {
    dfn[x]=low[x]=++cnt;
    q[++top]=x;inq[x]=1;
    int v,now=0;
    for (int i=head[x];i;i=next[i]) 
        if (e[i]) {//注意是在残量网络跑scc
            if (!dfn[v=ver[i]])
                Tarjan(v),low[x]=min(low[x],low[v]);
            else if (inq[v])
                low[x]=min(low[x],dfn[v]);
        }
    if (low[x]==dfn[x]) {
        scc++;
        while (top&&x!=now) {//
            now=q[top--];
            inq[now]=0;
            hav[scc]++;
            bel[now]=scc;
        }
    }
}

int main () {
    int u,v,w,tmp;
    scanf("%d%d%d%d",&n,&m,&s,&t);
    for (int i=1;i<=m;i++) {
        scanf("%d%d%d",&u,&v,&w);
        add(u,v,w);
    }
    while (bfs()) 
        while (tmp=dinic(s,inf)) maxflow+=tmp;
    for (int i=1;i<=n;i++) 
        if (!dfn[i]) Tarjan(i);
    for (int i=2;i<=tot;i+=2) {
        if (bel[ver[i]]!=bel[ver[i^1]] && !e[i]) printf("%d ",1); else printf("%d ",0);
        if (bel[ver[i]]!=bel[ver[i^1]] && !e[i] && bel[ver[i]]==bel[t] && bel[ver[i^1]]==bel[s]) printf("%d\n",1); else printf("%d\n",0);
    }
}