[bzoj1797] [Ahoi2009]Mincut 最小割

　　问最小割中可能成为割边和一定会成为割边的边有哪些。

　　膜了半天各路题解。

　　比较详细的解释：

　　http://blog.csdn.net/horizon_smz/article/details/50889806

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=4023,maxm=60233;
struct zs{
    int too,pre,flow;
}e[maxm<<1];int tot,last[maxn];
short dis[maxn],dl[maxn];
int st[maxn],top,dfn[maxn],low[maxn],tim,bel[maxn],cnt;
bool ins[maxn];
int i,j,k,n,m,s,t;
 
int ra;char rx;
inline int read(){
    rx=getchar(),ra=0;
    while(rx<'0'||rx>'9')rx=getchar();
    while(rx>='0'&&rx<='9')ra*=10,ra+=rx-48,rx=getchar();return ra;
}
inline bool bfs(){
    memset(dis,0,(n+1)<<1);
    int l=0,r=1,i,now;dl[1]=s,dis[s]=1;
    while(l<r&&!dis[t])
        for(i=last[now=dl[++l]];i;i=e[i].pre)if(e[i].flow&&!dis[e[i].too])
            dl[++r]=e[i].too,dis[e[i].too]=dis[now]+1;
    return dis[t];
}
inline int min(int a,int b){return a<b?a:b;}
int dfs(int x,int mx){
    if(x==t)return mx;
    int i,used=0,w;
    for(i=last[x];i;i=e[i].pre)if(e[i].flow&&dis[e[i].too]==dis[x]+1){
        w=dfs(e[i].too,min(e[i].flow,mx-used));if(w){
            e[i].flow-=w,e[i^1].flow+=w,used+=w;
            if(used==mx)return mx;
        }
    }
    dis[x]=0;return used;
}
 
void tarjan(int x){
    int i;
    dfn[x]=low[x]=++tim,ins[x]=1,st[++top]=x;
    for(i=last[x];i;i=e[i].pre)if(e[i].flow){
        if(!dfn[e[i].too])tarjan(e[i].too),low[x]=min(low[x],low[e[i].too]);
        else if(ins[e[i].too])low[x]=min(low[x],dfn[e[i].too]);
    }
    if(low[x]==dfn[x]){
        cnt++;
        while(st[top]!=x)bel[st[top]]=cnt,ins[st[top]]=0,top--;
        bel[x]=cnt,ins[x]=0,top--;
    }
}
 
inline void insert(int a,int b,int c){
    e[++tot].too=b,e[tot].flow=c,e[tot].pre=last[a],last[a]=tot,
    e[++tot].too=a,e[tot].flow=0,e[tot].pre=last[b],last[b]=tot;
}
 
int main(){
    n=read(),m=read(),s=read(),t=read(),tot=1;
    for(i=1;i<=m;i++)j=read(),k=read(),insert(j,k,read());
    while(bfs())dfs(s,1002333333);
    for(i=1;i<=n;i++)if(!dfn[i])tarjan(i);
    for(i=2;i<=tot;i+=2){
        if(e[i].flow){puts("0 0");continue;}
        if(bel[e[i^1].too]!=bel[e[i].too])printf("1 ");else printf("0 ");
        if(bel[e[i^1].too]==bel[s]&&bel[e[i].too]==bel[t])puts("1");else puts("0");
    }
    return 0;
}