【分析】
刚这个题刚了一早上,终于懂了一些。
结论:
在残量网络中跑tarjan后,对每条满流边(u,v):
①(u,v)能够出现在某个最小割集中,当且仅当Belong[u]!=Belong[v];
②(u,v)必定出现在最小割集中,当且仅当Belong[u]==Belong[S ]且Belong[v]==Belong[T]。
为什么呢?
先说②:这代表从S–>u有通路,v–>T有通路,所以若将该边流量增大,则S–>T又可增广。所以该边在此条链中一定为割边,所以必定出现在最小割集中。
对于①:【引用jcvb】
将每个SCC缩成一个点,得到的新图就只含有满流边了。那么新图的任一s-t割都对应原图的某个最小割,从中任取一个把id[u]和id[v]割开的割即可证明。
自己yy一下,也就是残量网络中还存在u–>v的路径,所以u–>v并没有被割掉,自然边(u,v)就不是割边。(我也不知道我想的对不对。。)
PS:因为少个else。。傻傻的查了半天,调出来把自己蠢哭了╮(╯▽╰)╭
【代码】
#include <cstdio>
#include <iostream>
#include <queue>
#include <vector>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <stack>
#define N 4005
#define M 120005
#define INF 1000000000
using namespace std;
typedef long long ll;
typedef pair<ll,ll> pa;
int read()
{
int x=0,f=1;char ch=getchar();
while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
while(isdigit(ch)){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
int n,m,cnt=1,S,T,ind,ans,scc;
int b[M],p[N],nextedge[M],w[M],cur[N];
int Level[N],Dfn[N],Low[N],Belong[N];
bool instack[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);
}
void Input_Init()
{
n=read(),m=read(),S=read(),T=read();
for(int i=1;i<=m;i++)
{
static int x,y,z;
x=read(),y=read(),z=read();
Anode(x,y,z);
}
}
bool Bfs()
{
queue<int>q;
q.push(S);
for(int i=1;i<=n;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=1;i<=n;i++) cur[i]=p[i];
ans+=Dfs(S,INF);
}
}
stack<int>st;
void Tarjan(int x)
{
Dfn[x]=Low[x]=++ind;
st.push(x);instack[x]=1;
for(int i=p[x];i;i=nextedge[i]) if(w[i])
{
int v=b[i];
if(!Dfn[v])
{
Tarjan(v);
Low[x]=min(Low[x],Low[v]);
}
else if(instack[v]) Low[x]=min(Low[x],Dfn[v]);
}
if(Dfn[x]==Low[x])
{
scc++;
int now=-1;
while(now!=x)
{
now=st.top();st.pop();
instack[now]=0;
Belong[now]=scc;
}
}
}
int main()
{
Input_Init();
Dinic();
for(int i=1;i<=n;i++)
if(!Dfn[i]) Tarjan(i);
for(int i=2;i<=cnt;i+=2)
{
if(w[i]) printf("0 0\n");
else
{
if(Belong[b[i]]!=Belong[b[i+1]])
printf("1 ");
else printf("0 ");
if(Belong[S]==Belong[b[i+1]]&&Belong[T]==Belong[b[i]])
printf("1\n");
else printf("0\n");
}
}
return 0;
}