AC日记——[CQOI2014]危桥 洛谷 P3163

题目描述

Alice和Bob居住在一个由N座岛屿组成的国家，岛屿被编号为0到N-1。某些岛屿之间有桥相连，桥上的道路是双向的，但一次只能供一人通行。其中一些桥由于年久失修成为危桥，最多只能通行两次。Alice希望在岛屿al和a2之间往返an次（从al到a2再从a2到al算一次往返）。同时，Bob希望在岛屿bl和b2之间往返bn次。这个过程中，所有危桥最多通行两次，其余的桥可以无限次通行。请问Alice和Bob能完成他们的愿望吗？

输入输出格式

输入格式：

本题有多组测试数据。每组数据第一行包含7个空格隔开的整数，分别为N、al、a2、an、bl、b2、bn。接下来是一个N行N列的对称矩阵，由大写字母组成。矩阵的i行j列描述编号i一1和j-l的岛屿间的连接情况，若为”O“则表示有危桥相连：为”N“表示有普通的桥相连：为”X“表示没有桥相连。|

输出格式：

对于每组测试数据输出一行，如果他们都能完成愿望输出”Yes“，否则输出”No“。

输入输出样例

输入样例#1：

4 0 1 1 2 3 1
XOXX
OXOX
XOXO
XXOX
4 0 2 1 1 3 2
XNXO
NXOX
XOXO
OXOX

输出样例#1：

Yes
No
数据范围
4

说明

4<=N<50

0<=a1, a2, b1, b2<=N-1

1 <=an. b<=50

思路：

　　网络流；

　　s连a1，b1，t连a2,b2跑一遍；

　　s连a1,b2，t连a2,b1跑一遍：

　　如果过了就yes负责no；

来，上代码：

#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>

#define maxn 55
#define INF 0x7ffffff

using namespace std;

struct EdgeType {
int v,f,e;
};
struct EdgeType edge[maxn*maxn*maxn];

int deep[maxn],s,t,cnt;
int n,a1,a2,an,b1,b2,bn,head[maxn];

inline void edge_add(int u,int v,int f)
{
edge[++cnt].v=v,edge[cnt].f=f,edge[cnt].e=head[u],head[u]=cnt;
edge[++cnt].v=u,edge[cnt].f=0,edge[cnt].e=head[v],head[v]=cnt;
}

bool bfs()
{
for(int i=0;i<maxn;i++) deep[i]=-1;
queue<int>que;que.push(s),deep[s]=0;
while(!que.empty())
{
int now=que.front();que.pop();
for(int i=head[now];i;i=edge[i].e)
{
if(deep[edge[i].v]<0&&edge[i].f>0)
{
deep[edge[i].v]=deep[now]+1;
if(edge[i].v==t) return true;
que.push(edge[i].v);
}
}
}
return false;
}

int flowing(int now,int flow)
{
if(now==t||flow<=0) return flow;
int oldflow=0;
for(int i=head[now];i;i=edge[i].e)
{
if(edge[i].f<=0||deep[edge[i].v]!=deep[now]+1) continue;
int pos=flowing(edge[i].v,min(flow,edge[i].f));
if(pos>0)
{
flow-=pos;
oldflow+=pos;
edge[i].f-=pos;
edge[i^1].f+=pos;
if(flow==0) return oldflow;
}
}
if(oldflow==0) deep[now]=-1;
return oldflow;
}

int main()
{
char ch[maxn][maxn];
while(scanf("%d%d%d%d%d%d%d",&n,&a1,&a2,&an,&b1,&b2,&bn)==7)
{
int ans_a=0,ans_b=0;cnt=1;
for(int i=0;i<maxn;i++) head[i]=0;
for(int i=0;i<n;i++)
{
scanf("%s",ch[i]);
for(int j=0;j<n;j++)
{
if(ch[i][j]=='O') edge_add(i,j,2);
if(ch[i][j]=='N') edge_add(i,j,INF);
}
}
s=maxn-2,t=maxn-1;
edge_add(s,a1,an*2),edge_add(s,b1,bn*2);
edge_add(a2,t,an*2),edge_add(b2,t,bn*2);
while(bfs()) ans_a+=flowing(s,INF);
cnt=1;for(int i=0;i<maxn;i++) head[i]=0;
for(int i=0;i<n;i++)
{
for(int j=0;j<n;j++)
{
if(ch[i][j]=='O') edge_add(i,j,2);
if(ch[i][j]=='N') edge_add(i,j,INF);
}
}
edge_add(s,a1,an*2),edge_add(s,b2,bn*2);
edge_add(a2,t,an*2),edge_add(b1,t,bn*2);
while(bfs()) ans_b+=flowing(s,INF);
if(ans_a>=(an+bn)*2&&ans_b>=(an+bn)*2) printf("Yes\n");
else printf("No\n");
}
return 0;
}