bzoj1797 [Ahoi2009]最小割（最小割+tarjan缩点）

拖了一个世纪的题。嗯。
求最小割，问每一条边是否可能是割边，是否一定是割边。
我们先跑一遍dinic，然后对残余网络进行scc缩点。结论如下：

1.如果（u，v）可能被割掉，则scc[u]!=scc[v]
2.如果（u，v）一定被割掉，则scc[u]=scc[S],scc[v]=scc[T]

为什么呢？首先割边一定是满流边。
因为不满流的边一定不是割边，一定有更优的。这个想一想还是可以理解的。

其次，我们可以把所有点分成三部分，S可以到达的点V1，可以到达T的点V2，其他的点V3。割开之后，V1一定属于S集合，V2一定属于T集合。V3既可以属于S集合，也可以属于T集合。
对于结论1，我们肯定有scc[S]!=scc[T],如果scc[u]==scc[v]即u,v两点强连通，则u,v一定在一个点集里,则边(u,v)一定不可能是割边。否则把u放在S集合里，v放在T集合里，一定可以是割边。

对于结论2，从S集合指向T集合，满流，无话可说，一定是割边。

另外对于一些博客所说的，缩完点之后的新图的边一定是满流的，我觉得不太准确，因为还可能存在一些流量为0的边。所以输出答案时还是要判一下这条边是否是满流边的。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
#include <stack>
using namespace std;
#define inf 0x3f3f3f3f
#define N 4010
#define M 60010
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,m,S,T,h[N],num=1,lev[N],dfn[N],low[N],bel[N],scc=0,dfnum=0,cur[N];
struct edge{
    int to,next,val;
}data[M<<1];
struct Edge{
    int x,y;
}e[M];
inline void add(int x,int y,int val){
    data[++num].to=y;data[num].next=h[x];h[x]=num;data[num].val=val;
    data[++num].to=x;data[num].next=h[y];h[y]=num;data[num].val=0;
}
inline bool bfs(){
    queue<int>q;memset(lev,0,sizeof(lev));
    q.push(S);lev[S]=1;
    while(!q.empty()){
        int x=q.front();q.pop();
        for(int i=h[x];i;i=data[i].next){
            int y=data[i].to;if(lev[y]||!data[i].val) continue;
            lev[y]=lev[x]+1;if(y==T) return 1;q.push(y);
        }
    }return 0;
}
inline int dinic(int x,int loww){
    if(x==T) return loww;int tmp=loww;
    for(int &i=cur[x];i;i=data[i].next){
        int y=data[i].to;if(lev[y]!=lev[x]+1||!data[i].val) continue;
        int res=dinic(y,min(tmp,data[i].val));
        if(!res) lev[y]=0;else tmp-=res,data[i].val-=res,data[i^1].val+=res;
        if(!tmp) return loww;
    }return loww-tmp;
}
stack<int>qq;bool inq[N];
inline void tarjan(int x){
    dfn[x]=low[x]=++dfnum;qq.push(x);inq[x]=1;
    for(int i=h[x];i;i=data[i].next){
        if(!data[i].val) continue;int y=data[i].to;
        if(!dfn[y]) tarjan(y),low[x]=min(low[x],low[y]);
        else if(inq[y]) low[x]=min(low[x],dfn[y]);
    }if(low[x]==dfn[x]){
        ++scc;while(1){
            int y=qq.top();qq.pop();inq[y]=0;
            bel[y]=scc;if(y==x) break;
        }
    }
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();S=read();T=read();
    for(int i=1;i<=m;++i){
        e[i].x=read();e[i].y=read();add(e[i].x,e[i].y,read());
    }while(bfs()){memcpy(cur,h,sizeof(h));dinic(S,inf);}
    for(int i=1;i<=n;++i) if(!dfn[i]) tarjan(i);S=bel[S];T=bel[T];
    for(int i=1;i<=m;++i){
        int x=bel[e[i].x],y=bel[e[i].y];
        if(data[i<<1].val){puts("0 0");continue;}//可能会有流量为0的边！ 
        putchar(x!=y?'1':'0');putchar(' ');
        putchar(x==S&&y==T?'1':'0');putchar('\n');
    }return 0;
}