Description
Alice和Bob居住在一个由N座岛屿组成的国家,岛屿被编号为0到N-1。某些岛屿之间有桥相连,桥上的道路是双
向的,但一次只能供一人通行。其中一些桥由于年久失修成为危桥,最多只能通行两次。Alice希望在岛屿al和a2之间往返an次(从al到a2再从a2到al算一次往返)。同时,Bob希望在岛屿bl和b2之间往返bn次。这个过程中,所有危桥最多通行两次,其余的桥可以无限次通行。请问Alice和Bob能完成他们的愿望吗?
Input
本题有多组测试数据。
每组数据第一行包含7个空格隔开的整数,分别为N、al、a2、an、bl、b2、bn。
接下来是一个N行N列的对称矩阵,由大写字母组成。矩阵的i行j列描述编号i一1和j-l的岛屿间的连接情况,若为“O”则表示有危桥相连:为“N”表示有普通的桥相连:为“X”表示没有桥相连。
|
Output
对于每组测试数据输出一行,如果他们都能完成愿望输出“Yes”,否则输出“No”。
Sample Input
4 0 1 1 2 3 1
XOXX
OXOX
XOXO
XXOX
4 0 2 1 1 3 2
XNXO
NXOX
XOXO
OXOX
XOXX
OXOX
XOXO
XXOX
4 0 2 1 1 3 2
XNXO
NXOX
XOXO
OXOX
Sample Output
Yes
No
数据范围
4<=N<50
O<=a1, a2, b1, b2<=N-1
1 <=an. b<=50
No
数据范围
4<=N<50
O<=a1, a2, b1, b2<=N-1
1 <=an. b<=50
题解:
源点向a1,b1连容量为an的边.
a2,b2向汇点连容量为bn的边.
是危桥的边容量为1
普通的桥容量为inf.
验证最大流是否为an+bn即可.
但是这样处理不了a1流到b2的情况.
所以把b1,b2交换一下再跑一遍即可.
代码:
#include<iostream>
#include<cstdio>
#include<cstring>
#define N 100
#define M 20000
#define inf 210000000
using namespace std;
int point[N],next[M<<1],n,m,a1,a2,b1,b2,an,bn;
int T,dis[N],pre[N],gap[N],cur[N],map[N][N],cnt(1);
char ch[N];
bool f;
struct use{
int st,en,v;
}e[M<<1];
void add(int x,int y,int v){
//cout<<x<<' '<<y<<' '<<v<<endl;
next[++cnt]=point[x];point[x]=cnt;
e[cnt].st=x;e[cnt].en=y;e[cnt].v=v;
next[++cnt]=point[y];point[y]=cnt;
e[cnt].st=y;e[cnt].en=x;e[cnt].v=0;
}
int isap(){
int i,u(1),mn,ans(0);
memset(dis,0,sizeof(dis));
memset(gap,0,sizeof(gap));
for(i=1;i<=T;i++) cur[i]=point[i];
gap[0]=T;
while (dis[1]<T){
f=false;
for (i=cur[u];i;i=next[i])
if (e[i].v&&dis[e[i].en]+1==dis[u]){cur[u]=i;f=1;break;}
if (f){
pre[u=e[i].en]=i;
if (u==T){
mn=inf;
for (i=T;i!=1;i=e[pre[i]].st) mn=min(mn,e[pre[i]].v);
ans+=mn;
for (i=T;i!=1;i=e[pre[i]].st) e[pre[i]].v-=mn,e[pre[i]^1].v+=mn;
u=1;
}
}
else{
gap[dis[u]]--;if (!gap[dis[u]]) return ans;
for (mn=T,i=point[u];i;i=next[i])
if (e[i].v) mn=min(mn,dis[e[i].en]);
gap[dis[u]=mn+1]++;cur[u]=point[u];
if (u!=1) u=e[pre[u]].st;
}
}
return ans;
}
int main(){
while (scanf("%d",&n)!=EOF){
memset(point,0,sizeof(point));
memset(map,0,sizeof(map));
scanf("%d%d%d%d%d%d%d",&a1,&a2,&an,&b1,&b2,&bn);
T=n+2;cnt=1;
for (int i=1;i<=n;i++){
scanf("%s",ch+1);
for (int j=1;j<=n;j++){
if (ch[j]=='X') map[i][j]=0;
if (ch[j]=='O') map[i][j]=1,add(i+1,j+1,1);
if (ch[j]=='N') map[i][j]=2,add(i+1,j+1,inf);
}
}
add(1,a1+2,an);add(1,b1+2,bn);
add(a2+2,T,an);add(b2+2,T,bn);
if (isap()!=an+bn){puts("No");continue;}
else{
memset(point,0,sizeof(point));cnt=1;
for (int i=1;i<=n;i++)
for (int j=1;j<=n;j++){
if (map[i][j]==1) add(i+1,j+1,1);
if (map[i][j]==2) add(i+1,j+1,inf);
}
add(1,a1+2,an);add(1,b2+2,bn);
add(a2+2,T,an);add(b1+2,T,bn);
if (isap()==an+bn) puts("Yes");
else puts("No");
}
}
}