[BZOJ3504][CQOI2014]危桥

bzoj
luogu

sol

很显然的建边方式：普通的桥建双向流量(inf)的边，危桥建双向流量(2)的边。
然后分别从(a_1)，(b_1)向(a_2)，(b_2)跑最大流，保证流量分别要大于等于(a_n)，(b_n)。
多源汇且源汇一一对应的网络流是一类NPC问题，但在这里只有两源两汇，可以采用两次(Dinic)的方式，第一次(S)连(a_1),(b_1)，(T)连(a_2),(b_2)跑最大流，第二次(S)连(a_1),(b_2)，(T)连(a_2),(b_1)再跑一遍，只要两遍都能保证满流就行了。

code

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
using namespace std;
const int N = 60;
const int inf = 1e9;
struct edge{int to,nxt,w;}a[N*N<<2];
int n,a1,a2,an,b1,b2,bn,head[N],cnt,dep[N],cur[N];
char map[N][N];queue<int>Q;
void link(int u,int v,int w)
{
	a[++cnt]=(edge){v,head[u],w};
	head[u]=cnt;
	a[++cnt]=(edge){u,head[v],w};
	head[v]=cnt;
}
bool bfs(int s,int t)
{
	memset(dep,0,sizeof(dep));
	dep[s]=1;Q.push(s);
	while (!Q.empty())
	{
		int u=Q.front();Q.pop();
		for (int e=head[u];e;e=a[e].nxt)
			if (a[e].w&&!dep[a[e].to])
				dep[a[e].to]=dep[u]+1,Q.push(a[e].to);
	}
	return dep[t];
}
int dfs(int u,int f,int t)
{
	if (u==t) return f;
	for (int &e=cur[u];e;e=a[e].nxt)
		if (a[e].w&&dep[a[e].to]==dep[u]+1)
		{
			int tmp=dfs(a[e].to,min(a[e].w,f),t);
			if (tmp) {a[e].w-=tmp;a[e^1].w+=tmp;return tmp;}
		}
	return 0;
}
int Dinic(int s,int t)
{
	int res=0;
	while (bfs(s,t))
	{
		for (int i=1;i<=t;++i) cur[i]=head[i];
		while (int tmp=dfs(s,inf,t)) res+=tmp;
	}
	return res;
}
int solve(int s1,int s2,int t1,int t2)
{
	int s=n+1,t=n+2;
	memset(head,0,sizeof(head));cnt=1;
	for (int i=1;i<n;++i)
		for (int j=i+1;j<=n;++j)
			if (map[i][j]=='N') link(i,j,inf);
			else if (map[i][j]=='O') link(i,j,2);
	link(s,s1,an<<1);link(s,s2,bn<<1);
	link(t1,t,an<<1);link(t2,t,bn<<1);
	return Dinic(s,t);
}
int main()
{
	while (scanf("%d%d%d%d%d%d%d",&n,&a1,&a2,&an,&b1,&b2,&bn)!=EOF)
	{
		++a1;++a2;++b1;++b2;
		for (int i=1;i<=n;++i) scanf("%s",map[i]+1);
		puts(solve(a1,b1,a2,b2)==2*(an+bn)&&solve(a1,b2,a2,b1)==2*(an+bn)?"Yes":"No");
	}
	return 0;
}