两遍网络流 - 危桥（CQOI2014）

Analysis

建图很好想，建立源点和汇点
源点向S1，S2建一条流量为次数*2的边
T1，T2向汇点同理建边
最后检查是否满流即可

但只做一遍是有问题的（省选题果然不是盖的）
因为在流的时候可能会存在S1最后流到了T2，S2流到了T1
一个显然的反例
如果要求0->3,1->2
跑出来的最大流是合法的，但实际情况不对

一个巧妙的解决办法：
把一条路径反过来再跑一遍，如果还满流，则说明有解。

正确性：
如果本身是合法的（也就是S1->T1,S2->T2）那么将路径反过来再跑一定也合法。所以不会将本来正确的情况判错
如果本身不合法但将路径反过来后合法，
反过来满流也分为两种情况：
1.S1->S2,T2->T1
因为已经满足S1->T2了，则S1也可到T1，正确
2.S1->T1,T2->S2 这种情况显然可行

Code

#include<bits/stdc++.h>
#define in read()
#define re register
using namespace std;
inline int read(){
	char ch;int f=1,res=0;
	while((ch=getchar())<'0'||ch>'9') if(ch=='-') f=-1;
	while(ch>='0'&&ch<='9'){
		res=(res<<1)+(res<<3)+(ch^48);
		ch=getchar();
	}
	return f==1?res:-res;
}
const int inf=1e9;
const int N=55,M=6000;
int n,S,T,sum=0,a1,a2,a3,b1,b2,b3;
int nxt[M],head[N],to[M],cap[M],ecnt=1;
inline void add(int x,int y,int z){
	nxt[++ecnt]=head[x];head[x]=ecnt;to[ecnt]=y;cap[ecnt]=z;
	nxt[++ecnt]=head[y];head[y]=ecnt;to[ecnt]=x;cap[ecnt]=0;
}
int cur[N],lev[N];
bool bfs(){
	queue<int> q;
	for(re int i=S;i<=T;++i) lev[i]=-1,cur[i]=head[i];
	lev[S]=0;q.push(S);
	while(!q.empty()){
		int u=q.front();q.pop();
		for(re int e=head[u];e;e=nxt[e]){
			int v=to[e];
			if(lev[v]!=-1||cap[e]<=0) continue;
			lev[v]=lev[u]+1;
			if(v==T) return 1;
			q.push(v);
		}
	}
	return 0;
}
int dinic(int u,int flow){
	if(u==T) return flow;
	int delta=0,res=0;
	for(re int &e=cur[u];e;e=nxt[e]){
		int v=to[e];
		if(lev[v]>lev[u]&&cap[e]){
			delta=dinic(v,min(flow-res,cap[e]));
			if(delta){
				cap[e]-=delta;cap[e^1]+=delta;
				res+=delta;if(res==flow) return res;
			}
		}
	}
	return res;
}
int f[55][55];
void build(){
	memset(head,0,sizeof(head));
	ecnt=1;
	for(re int i=1;i<=n;++i)
		for(re int j=1;j<=n;++j){
			if(!f[i][j]) continue;
			add(i,j,f[i][j]==1?2:inf);
		}
}
void init(){
	memset(f,0,sizeof(f));
	S=0;T=n+1;sum=0;
}

int main(){
	while((scanf("%d%d%d%d%d%d%d",&n,&a1,&a2,&a3,&b1,&b2,&b3))!=EOF){
		init();
		a1++;a2++;b1++;b2++;
		sum=(a3+b3)*2;
		char st[55];
		for(re int i=1;i<=n;++i){
			scanf("%s",st+1);
			for(re int j=1;j<=n;++j){
				if(st[j]=='O') f[i][j]=1;
				else if(st[j]=='N') f[i][j]=2;
			}
		}
		build();
		add(S,a1,a3*2);add(a2,T,a3*2);
		add(S,b1,b3*2);add(b2,T,b3*2);
		int maxflow=0;
		while(bfs()) maxflow+=dinic(S,inf);
		bool flag=0;
		if(maxflow!=sum) flag=1;
		if(!flag){
			build();
			add(S,a1,a3*2);add(a2,T,a3*2);
			add(S,b2,b3*2);add(b1,T,b3*2);
			int maxflow=0;
			while(bfs()) maxflow+=dinic(S,inf);
			if(maxflow!=sum) flag=1;
		}
		if(flag) puts("No");
		else puts("Yes");
	}
	return 0;
}