【BZOJ1458】士兵占领 最小流

【BZOJ1458】士兵占领

Description

有一个M * N的棋盘，有的格子是障碍。现在你要选择一些格子来放置一些士兵，一个格子里最多可以放置一个士兵，障碍格里不能放置士兵。我们称这些士兵占领了整个棋盘当满足第i行至少放置了Li个士兵, 第j列至少放置了Cj个士兵。现在你的任务是要求使用最少个数的士兵来占领整个棋盘。

Input

第一行两个数M, N, K分别表示棋盘的行数，列数以及障碍的个数。 第二行有M个数表示Li。 第三行有N个数表示Ci。 接下来有K行，每行两个数X, Y表示(X, Y)这个格子是障碍。

Output

输出一个数表示最少需要使用的士兵个数。如果无论放置多少个士兵都没有办法占领整个棋盘，输出”JIONG!” (不含引号)

Sample Input

4 4 4
1 1 1 1
0 1 0 3
1 4
2 2
3 3
4 3

Sample Output

4
数据范围

M, N <= 100, 0 <= K <= M * N

题解：再说一遍最小流怎么写吧

设源点s，汇点t，新建源点S，汇点T

从s向行连下界li，上界∞的边，从列向t连下界ci，上界∞的边，从行向列连上界1的边，从t到s连上界∞的边，跑从S到T的可行流，记录t到s这条边的反向边流量x1

在删掉S和T，t到s，跑从t到s的最大流x2，x1-x2即为答案

 

#include <cstdio>
#include <cstring>
#include <iostream>
#include <queue>
using namespace std;
int m,n,K,S,T,t1,t2,cnt,ans;
int map[110][110],d[210],to[30000],next[30000],val[30000],head[210];
queue<int> q;
void add(int a,int b,int c)
{
	to[cnt]=b;
	val[cnt]=c;
	next[cnt]=head[a];
	head[a]=cnt++;
}
int bfs()
{
	int i,u;
	while(!q.empty())	q.pop();
	memset(d,0,sizeof(d));
	d[S]=1,q.push(S);
	while(!q.empty())
	{
		u=q.front(),q.pop();
		for(i=head[u];i!=-1;i=next[i])
		{
			if(!d[to[i]]&&val[i])
			{
				d[to[i]]=d[u]+1;
				if(to[i]==T)	return 1;
				q.push(to[i]);
			}
		}
	}
	return 0;
}
int dfs(int x,int mf)
{
	if(x==T)	return mf;
	int k,i,temp=mf;
	for(i=head[x];i!=-1;i=next[i])
	{
		if(d[to[i]]==d[x]+1&&val[i])
		{
			k=dfs(to[i],min(temp,val[i]));
			if(!k)	d[to[i]]=0;
			val[i]-=k,val[i^1]+=k,temp-=k;
			if(!temp)	break;
		}
	}
	return mf-temp;
}
int main()
{
	memset(head,-1,sizeof(head));
	scanf("%d%d%d",&m,&n,&K);
	S=m+n+2,T=m+n+3;
	int i,j,a,b;
	for(i=1;i<=m;i++)
	{
		scanf("%d",&a);
		add(0,i,1<<30),add(i,0,0);
		add(S,i,a),add(i,S,0);
		t1+=a;
	}
	add(0,T,t1),add(T,0,0);
	for(i=1;i<=n;i++)
	{
		scanf("%d",&a);
		add(m+i,m+n+1,1<<30),add(m+n+1,m+i,0);
		add(m+i,T,a),add(T,m+i,0);
		t2+=a;
	}
	for(i=1;i<=K;i++)
	{
		scanf("%d%d",&a,&b);
		map[a][b]=1;
	}
	for(i=1;i<=m;i++)
		for(j=1;j<=n;j++)
			if(!map[i][j])
				add(i,m+j,1),add(m+j,i,0);
	add(S,n+m+1,t2),add(n+m+1,S,0);
	while(bfs())	ans+=dfs(S,1<<30);
	if(ans<max(t1,t2))
	{
		printf("JIONG");
		return 0;
	}
	for(i=head[S];i!=-1;i=next[i])	val[i]=val[i^1]=0;
	for(i=head[T];i!=-1;i=next[i])	val[i]=val[i^1]=0;
	val[cnt-1]=val[cnt-2]=0;
	S=m+n+1,T=0;
	while(bfs())	ans-=dfs(S,1<<30);
	printf("%d",ans);
	return 0;
}