ZOJ 2587 Unique Attack

判断最小割是否唯一，先求一次最大流，然后在残留网络中分别从源汇开始dfs一次，找出最小割[S,T]，如果SUT不包含所有点，那么最小割不唯一

 

假设点i不被SUT包含，那么残留网络中s不能到达i,i不能到达t,即进入i的边和从i出去的边都满流，假设某条进入i的边x满流，这些流量从若干条边y流出i，那么，如果选x为割边，或者选所有对应的y为割边，不会影响最大流，即最小割容量不变，最小割也就不唯一

 

如果最小割不唯一，假设有两个最小割[S,T],[S2,T2]，则存在点i,所有进入i的边和从i出去的边都满流，否则最小割唯一，dfs得到的S包含于S2，i属于S2-S，又由于i不属于T,且i不属于T2,SUT就不包含i,也就不包含所有点

 

 

代码：

#include<iostream> #include<algorithm> #include<memory.h> #include<stdio.h> using namespace std; const int MAX=805; const int inf=1<<30; struct node { int v,c,next; }g[MAX*MAX]; int dis[MAX],num[MAX],pre[MAX],adj[MAX],cur[MAX],vis1[MAX],vis2[MAX]; int e,n,m,s,t,vn; void add(int u,int v,int c) { g[e].v=v; g[e].c=c; g[e].next=adj[u]; adj[u]=e++; g[e].v=u; g[e].c=0; g[e].next=adj[v]; adj[v]=e++; } void dfs1(int u) { vis1[u]=1; for(int i=adj[u];i!=-1;i=g[i].next) { if(g[i].c&&!vis1[g[i].v]) dfs1(g[i].v); } } void dfs2(int u) { vis2[u]=1; for(int i=adj[u];i!=-1;i=g[i].next) { if(g[i^1].c&&!vis2[g[i].v]) dfs2(g[i].v); } } int sap() { int i,u,v,flag,aug=inf,flow=0; for(i=1;i<=vn;i++) { cur[i]=adj[i]; num[i]=dis[i]=0; } num[0]=vn; pre[s]=u=s; while(dis[s]<vn) { flag=0; for(i=adj[u];i!=-1;i=g[i].next) { v=g[i].v; if(g[i].c&&dis[u]==dis[v]+1) { flag=1; aug=min(aug,g[i].c); pre[v]=u; cur[u]=i; u=v; if(u==t) { flow+=aug; while(u!=s) { u=pre[u]; g[cur[u]].c-=aug; g[cur[u]^1].c+=aug; } aug=inf; } break; } } if(flag) continue; if(--num[dis[u]]==0) break; for(dis[u]=inf,i=adj[u];i!=-1;i=g[i].next) { v=g[i].v; if(g[i].c&&dis[v]<dis[u]) { dis[u]=dis[v]; cur[u]=i; } } dis[u]++; num[dis[u]]++; u=pre[u]; } return flow; } int main() { int i,j,a,b,k; while(scanf("%d%d%d%d",&n,&m,&s,&t)!=EOF) { if(n==0&&m==0&&s==0&&t==0) break; e=0; vn=n; memset(adj,-1,sizeof(adj)); while(m--) { scanf("%d%d%d",&a,&b,&k); add(a,b,k); add(b,a,k); } sap(); memset(vis1,0,sizeof(vis1)); memset(vis2,0,sizeof(vis2)); dfs1(s); dfs2(t); int flag=1; for(i=1;i<=vn;i++) { if((vis1[i]&&vis2[i])||(!vis1[i]&&!vis2[i])) flag=0; } if(flag) puts("UNIQUE"); else puts("AMBIGUOUS"); } return 0; }