[BZOJ3504]CQOI2014危桥|最大流

源向a1和b1连流量为2*an，2*bn的边，a2和b2向汇连流量为2*an，2*bn的边，其他边非危桥连inf，危桥连2的边，跑最大流，但是这时满流不一定有解，因为可能有a1流到b2的流这样，所以把b1，b2反过来重构图再跑一次最大流，若能满流则有解。证明：假设第一次a1流了x的流到b2，第二次的时候由于是无向图，a1到a2流还是2*an-x，那么a1到b1是x，所以有x的流是b1-a1-b2，所以一样能流到b2。。

#include<cstdio>#include<iostream>#define N 60#define inf 0x3ffffffusing namespace std;struct edge{int e,f,next;}ed[N*N*2];int n,i,j,a1,a2,an,b1,b2,bn,ne,nd,t,a[N],d[N],que[N],u[N],mp[N][N];char s[N];void add(int s,int e,int f){ed[++ne].e=e;ed[ne].f=f;ed[ne].next=a[s];a[s]=ne;}bool bfs(int s,int t){int i,j,head=1,tail=1,get;for (i=0;i<=nd;i++) u[i]=d[i]=0;u[s]=1;que[1]=s;while (head<=tail){get=que[head++];for (j=a[get];j;j=ed[j].next)if (ed[j].f&&!u[ed[j].e]){d[ed[j].e]=d[get]+1;u[ed[j].e]=1;que[++tail]=ed[j].e;}}return d[t]!=0;}int extend(int x,int minf,int t){if (x==t) return minf;int f=minf,del,j;for (j=a[x];j;j=ed[j].next)if (ed[j].f&&d[ed[j].e]==d[x]+1){del=extend(ed[j].e,min(minf,ed[j].f),t);ed[j].f-=del;ed[j^1].f+=del;minf-=del;if (!minf) break;}if (f==minf) d[x]=0;return f-minf;}int dinic(int s,int t){int ans=0;while (bfs(s,t)) ans+=extend(s,inf,t);return ans;}void build(){ne=1;for (int i=0;i<=nd;i++) a[i]=0;add(0,a1,2*an);add(a1,0,0);add(0,b1,2*bn);add(b1,0,0);add(nd,a2,0);add(a2,nd,2*an);add(nd,b2,0);add(b2,nd,2*bn);for (int i=1;i<=n;i++)for (int j=i+1;j<=n;j++)if (mp[i][j]) add(i,j,mp[i][j]),add(j,i,mp[i][j]);}int main(){while (scanf("%d%d%d%d%d%d%d",&n,&a1,&a2,&an,&b1,&b2,&bn)!=EOF){nd=n+1;a1++;a2++;b1++;b2++;for (i=1;i<=n;i++){scanf("%s",s+1);for (j=i+1;j<=n;j++) if (s[j]=='O') mp[i][j]=2;else if (s[j]=='N') mp[i][j]=inf;else mp[i][j]=0;}build();t=dinic(0,nd);if (t<2*an+2*bn) printf("No\n");else{swap(b1,b2);build();t=dinic(0,nd);if (t<2*an+2*bn) printf("No\n");else printf("Yes\n");}}}