【bzoj3504】【cqoi2014】【危桥】【最大流】

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

Sample Output

Yes
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"); 
	 }	
  }	
}