Description
A,B两个国家正在交战,其中A国的物资运输网中有N个中转站,M条单向道路。设其中第i (1≤i≤M)条道路连接了vi,ui两个中转站,那么中转站vi可以通过该道路到达ui中转站,如果切断这条道路,需要代价ci。现在B国想找出一个路径切断方案,使中转站s不能到达中转站t,并且切断路径的代价之和最小。
小可可一眼就看出,这是一个求最小割的问题。但爱思考的小可可并不局限于此。现在他对每条单向道路提出两个问题:
问题一:是否存在一个最小代价路径切断方案,其中该道路被切断?
问题二:是否对任何一个最小代价路径切断方案,都有该道路被切断?
现在请你回答这两个问题。
Input
第一行有4个正整数,依次为N,M,s和t。第2行到第(M+1)行每行3个正整数v,u,c表示v中转站到u中转站之间有单向道路相连,单向道路的起点是v,终点是u,切断它的代价是c(1≤c≤100000)。
注意:两个中转站之间可能有多条道路直接相连。
同一行相邻两数之间可能有一个或多个空格。
Output
对每条单向边,按输入顺序,依次输出一行,包含两个非0即1的整数,分别表示对问题一和问题二的回答(其中输出1表示是,输出0表示否)。
同一行相邻两数之间用一个空格隔开,每行开头和末尾没有多余空格。
Sample Input
6 7 1 6
1 2 3
1 3 2
2 4 4
2 5 1
3 5 5
4 6 2
5 6 3
Sample Output
1 0
1 0
0 0
1 0
0 0
1 0
1 0
Data Constraint
Hint
【样例说明】
设第(i+1)行输入的边为i号边,那么{1,2},{6,7},{2,4,6}是仅有的三个最小代价切割方案。它们的并是{1,2,4,6,7},交是 。
【数据规模和约定】
测试数据规模如下表所示:
数据编号 1 2 3 4 5 6 7 8 9 10
N= 10 20 200 200 1000 1000 1000 2000 3000 4000
M= 50 200 2000 2000 20000 20000 40000 50000 60000 60000
分析:
首先对原图跑最大流,得到残余网络。
很容易得知若一条边u->v存在u->v的增广路,则这条边一定不是割边,反之亦然则可能是割边。
那么我们就把点分成了s集和t集还有一部分不确定。若一条边连接了s集和t集,则这条边一定是必然存在度最小割中的边。
那么我们就得到了这样一个算法:
把残余网络中每个SCC缩成一个点,这样图中剩下的边则必然是满流边,且图中的每一个割都对应了原图中的一个最小割。设点i所在的强连通分量为id[i].那么一条边可能存在于割集中当且仅当id[u]!=id[v]。一条边一定存在于割集中当且仅当id[u]==id[s]且id[v]==id[t]。
(by ymw%%%) dalao表示已经忘记证明了= =
代码:
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int N=4005;
const int M=60005;
int n,m,s,t,cnt,last[N],cur[N],dfn[N],low[N],tim,sta[N],top,dis[N],tot,id[N];
bool ins[N];
struct edge{int to,next,c;}e[M*2];
queue<int> que;
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;
}
void addedge(int u,int v,int c)
{
e[++cnt].to=v;e[cnt].c=c;e[cnt].next=last[u];last[u]=cnt;
e[++cnt].to=u;e[cnt].c=0;e[cnt].next=last[v];last[v]=cnt;
}
bool bfs()
{
for (int i=1;i<=n;i++) dis[i]=0;
while (!que.empty()) que.pop();
dis[s]=1;que.push(s);
while (!que.empty())
{
int u=que.front();que.pop();
for (int i=last[u];i;i=e[i].next)
if (e[i].c&&!dis[e[i].to])
{
dis[e[i].to]=dis[u]+1;
if (e[i].to==t) return 1;
que.push(e[i].to);
}
}
return 0;
}
int dfs(int x,int maxf)
{
if (x==t|!maxf) return maxf;
int ret=0;
for (int &i=cur[x];i;i=e[i].next)
if (e[i].c&&dis[e[i].to]==dis[x]+1)
{
int f=dfs(e[i].to,min(e[i].c,maxf-ret));
e[i].c-=f;
e[i^1].c+=f;
ret+=f;
if (maxf==ret) break;
}
return ret;
}
void dinic()
{
while (bfs())
{
for (int i=1;i<=n;i++) cur[i]=last[i];
dfs(s,1000000000);
}
}
void tarjan(int x)
{
dfn[x]=low[x]=++tim;
ins[x]=1;sta[++top]=x;
for (int i=last[x];i;i=e[i].next)
{
if (!e[i].c) continue;
if (!dfn[e[i].to])
{
tarjan(e[i].to);
low[x]=min(low[x],low[e[i].to]);
}
else if (ins[e[i].to]) low[x]=min(low[x],dfn[e[i].to]);
}
if (low[x]==dfn[x])
{
tot++;int y=sta[top--];id[y]=tot;ins[y]=0;
while (y!=x) y=sta[top--],id[y]=tot,ins[y]=0;
}
}
int main()
{
n=read();m=read();s=read();t=read();cnt=1;
for (int i=1;i<=m;i++)
{
int x=read(),y=read(),z=read();
addedge(x,y,z);
}
dinic();
for (int i=1;i<=n;i++)
if (!dfn[i]) tarjan(i);
for (int i=2;i<=cnt;i+=2)
if (e[i].c) puts("0 0");
else
{
if (id[e[i].to]!=id[e[i^1].to]) printf("1 ");
else printf("0 ");
if (id[e[i].to]==id[t]&&id[e[i^1].to]==id[s]) puts("1");
else puts("0");
}
return 0;
}