这里学到的是一个输出最小割的方案!!!
在跑完最大流后,属于最小割的边肯定cap都为零了。
我们再以源点为起点进行dfs搜索,标记。
然后最小割会将其阻拦住,然后我们再遍历边,看有哪些边一端被标记了另一端没有被标记。
这样就好了。
#include<iostream>
#include<algorithm>
#include<set>
using namespace std;
typedef pair<int,int> pii;
const int max_n = 100;
const int max_m = 1000;
const int inf = 2e9;
struct edge {
int to,cap,next;
}E[max_m<<2];
int n,m;
int head[max_n],cnt=1,gap[max_n],last[max_n],d[max_n],que[max_n],ql,qr;
void add(int u,int to,int cap) {
E[++cnt].to=to;E[cnt].cap=cap;
E[cnt].next=head[u];head[u]=cnt;
E[++cnt].to=u;E[cnt].cap=0;
E[cnt].next=head[to];head[to]=cnt;
}int aug(int x,int s,int t,int mi) {
if (x==t) return mi;
int flow=0;
for (int &i=last[x],v=E[i].to;i;i=E[i].next,v=E[i].to) if (d[x]==d[v]+1) {
int tmp=aug(v,s,t,min(mi,E[i].cap));
flow+=tmp,mi-=tmp,E[i].cap-=tmp,E[i^1].cap+=tmp;
if (!mi) return flow;
}
if (!(--gap[d[x]])) d[s]=n+1;
++gap[++d[x]],last[x]=head[x];
return flow;
}
int maxflow(int s,int t) {
fill(gap,gap+n+10,0);fill(d,d+n+10,0);
++gap[d[t]=1];
for (int i=1;i<=n;++i) last[i]=head[i];
que[ql=qr=1]=t;
while (ql<=qr) {
int x=que[ql++];
for (int i=head[x],v=E[i].to;i;i=E[i].next,v=E[i].to) if (!d[v]) ++gap[d[v]=d[x]+1],que[++qr]=v;
}
int ret=aug(s,s,t,inf);
while (d[s]<=n) ret+=aug(s,s,t,inf);
return ret;
}
int vis[100];
void dfs(int u){
vis[u]=1;
for (int i=head[u];i;i=E[i].next){
if (E[i].cap==0)continue;
if (vis[E[i].to])continue;
dfs(E[i].to);
}
}
int main(){
while(scanf("%d %d",&n,&m)){
if (n==m&&n==0)break;
fill(head,head+n+5,0);cnt=1;
for (int i=1,u,v,c;i<=m;++i){
scanf("%d %d %d",&u,&v,&c);
add(u,v,c);add(v,u,c);
}
maxflow(1,2);
set<pii> res;
fill(vis,vis+max_n,0);
dfs(1);
for (int u=1;u<=n;++u){
for (int i=head[u];i;i=E[i].next){
if (i&1)continue;
if (E[i].cap==0&&(vis[E[i].to]^vis[u]))
res.insert(pii(min(u,E[i].to),max(u,E[i].to)));
}
}
for (auto it=res.begin();it!=res.end();++it)
printf("%d %d\n",(*it).first,(*it).second);
printf("\n");
}
}
京公网安备 11010502036488号