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;
}