大视野1458--士兵占领（网络流）

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

思路：超源向行连边，边权为该行最多删点数。列向超汇连边，边权为该列最多删点数。行列没障碍则连边权为1的边。

这样求一次最大流则得出最多可删点数。用总点数-k-最大流即可。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <vector>
using namespace std;
#define maxn 308
#define maxm 20080
#define inf 0x3f3f3f3f
int first[maxn],dis[maxn],num[maxn];
int r[maxn],c[maxn];
int vv[maxm],ww[maxm],nxt[maxm];
bool Map[maxn][maxn];
int e,NN;
inline int min(int a,int b)
{
	return a>b?b:a;
}
void addedge(int u,int v,int w)
{
	vv[e] = v;	ww[e] = w;	nxt[e] = first[u];	first[u] = e++;
	vv[e] = u;	ww[e] = 0;	nxt[e] = first[v];	first[v] = e++;
}

int dfs(int u,int s,int d,int cost)
{
	if(u == d)	return cost;
	int ans = 0;
	int _min = NN;
	for(int i = first[u];i != -1;i = nxt[i])
	{
		int v = vv[i];
		if(ww[i])
		{
			if(dis[v] + 1 == dis[u])
			{
				int t = dfs(v,s,d,min(ww[i],cost));
				ww[i] -= t;
				ww[i^1] += t;
				ans += t;
				cost -= t;
				if(dis[s] == NN) return ans;
				if(!cost)	break;
			}
			if(_min > dis[v])_min = dis[v];
		}
	}
	if(!ans)
	{
		if(--num[dis[u]] == 0)	dis[s] = NN;
		dis[u] = _min + 1;
		num[dis[u]]++;
	}
	return ans;
}

int isap(int s,int d)
{
	memset(dis,0,sizeof(dis));
	memset(num,0,sizeof(num));
	memset(Map,0,sizeof(Map));
	num[0] = NN;
	int ans = 0;
	while(dis[s] < NN)
		ans += dfs(s,s,d,inf);
	return ans;
}

int main()
{
	//freopen("in.txt","r",stdin);
	int m,n,k;
	while(scanf("%d%d%d",&m,&n,&k)==3)
	{
		bool flag = false;
		e = 0;
		memset(first,-1,sizeof(first));
		memset(Map,0,sizeof(Map));
		NN = n+m+2;
		for(int i = 1;i <= m;i++)
		{
			scanf("%d",&r[i]);
			r[i] = n-r[i];
		}
		for(int i = 1;i <= n;i++)
		{
			scanf("%d",&c[i]);
			c[i] = m-c[i];
		}
		for(int i = 1;i <= k;i++)
		{
			int x,y;
			scanf("%d%d",&x,&y);
			Map[x][y] = 1;
			r[x] -= 1;	c[y] -= 1;
		}
		for(int i = 1;i <= m;i++)
			if(r[i] < 0)	flag = true;
		for(int i = 1;i <= n;i++)
			if(c[i] < 0)	flag = true;
		if(flag)
		{
			puts("JIONG!");
			continue;
		}
		for(int i = 1;i <= m;i++)
			addedge(0,i,r[i]);
		for(int i = 1;i <= n;i++)
			addedge(i+m,n+m+1,c[i]);
		for(int i = 1;i <= m;i++)
		{
			for(int j = 1;j <= n;j++)
			{
				if(!Map[i][j])
					addedge(i,j+m,1);
			}
		}
		int ans = n*m-isap(0,n+m+1)-k;
		printf("%d\n",ans);
	}
	return 0;
}