BZOJ1458 士兵占领-最大流

传送门

题目大意：自行参考

题解：首先一看方格图最值问题就是网络流（尽管我是知道了网络流才做的QwQ）

肯定是对每一行建一排点，每一列建一排点，建成一个二分图。

显然行和列连边容量为1.

首先注意到要求最小值，发现最小割没办法处理（因为从s对每一行的连边割去是没有意义的，因为可能只“割去”一部分）

然后注意到，如果s和每一行连边的容量表示“最多/最少能放多少个士兵的时候”，这就是个有上下界的“最小网络流”。

但是很明显我这么蒟蒻肯定不会做。

所以挖掘这个题的性质，首先是个二分图，但是并没有什么卵用。

最重要的一点是，按照刚才那个方法建图，是一定能跑出“满流”这种状态的！而且是全部的边（不包括反向弧）都是满流的！（这就对应了所有非障碍的地方都放士兵的情况）

所以显然我们想要从s“减少”最多的流，使得满足下界；

“减少”想到了什么？退流！没错，从s减少多少流，等价于从t出发多少流到达s！

所以可以这么做：先跑一边最大流（或者直接建图的时候设置一下容量和流量即可，其实不跑也可以，只是说出来大家体会一下精神），然后对于每条边，用上界减下界得到这条边最多能退多少流，并把这条边“反过来”（就是说原来是(u,v)现在是(v,u)），从t出发跑一边最大流即可！

但是这么做明显还不够清真，因为太麻烦了。

其实我们还可以找到另一种方法，得到相同的模型，其实上文说了，第一遍最大流没有必要跑，“意会”一下就好了。

所以可以直接依旧正着建边，对每条边，以行i为例，从s向i建边(s,i,m-l[i]-k1[i]),k1[i]就是i行有多少障碍。

其实就是刚才的“上界减去下界”，这样用“意会出来的最大流”，其实大家已经遇见了，答案就是nm-k，再减去跑出来的最大流就是答案。

这个是什么意思呢？

其实是一种很常用的思想，既然求的是最小值，那么就是一个数（相对定值）减去另一个数的最大值。

在这个例子中，相对定值的就是所有格子都放士兵的情况，减去的，就是把那些“多余”的士兵踢走，踢得越多，剩下的就越少，而这个是可以用网络最大流做的。

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<queue>
#include<climits>
#define MAXN 300
#define MAXM 50000
#define INF INT_MAX
#define build_edge(u,v,c) add_edge(u,v,c),add_edge(v,u,0)
using namespace std;
struct edges{
	int to,pre,resf;
}e[MAXM];
int h[MAXN],a[MAXN],b[MAXN],cur[MAXN],lev[MAXN];
bool vis[MAXN],has[MAXN][MAXN];queue<int> q;int etop;
int add_edge(int u,int v,int c)
{
	etop++;
	e[etop].to=v;
	e[etop].resf=c;
	e[etop].pre=h[u];
	h[u]=etop;
	return 0;
}
bool bfs(int s,int t)
{
	memset(vis,false,sizeof(vis));
	memset(lev,0,sizeof(lev));
	while(!q.empty()) q.pop();
	q.push(s);lev[s]=0;vis[s]=true;
	while(!q.empty())
	{
		int x=q.front();q.pop();
		for(int i=h[x];i;i=e[i].pre)
			if(!vis[e[i].to]&&e[i].resf)
			{
				vis[e[i].to]=true;
				lev[e[i].to]=lev[x]+1;
				q.push(e[i].to);
			}
	}
	return vis[t];
}
int dfs(int s,int t,int a)
{
	if(s==t||!a) return a;
	int flow=0,f;
	for(int &i=cur[s];i;i=e[i].pre)
		if(lev[e[i].to]==lev[s]+1&&(f=dfs(e[i].to,t,min(e[i].resf,a)))>0)
		{
			a-=f;flow+=f;e[i].resf-=f;
			e[((i-1)^1)+1].resf+=f;
			if(!a) break;
		}
	return flow;
}
int main()
{
	int n,m,k;scanf("%d%d%d",&n,&m,&k);
	for(int i=1;i<=n;i++) scanf("%d",&a[i]);
	for(int i=1;i<=m;i++) scanf("%d",&b[i]);
	for(int i=1;i<=n;i++) a[i]=m-a[i];
	for(int i=1;i<=m;i++) b[i]=n-b[i];
	for(int i=1;i<=k;i++)
	{
		int x,y;scanf("%d%d",&x,&y);
		a[x]--;b[y]--;has[x][y]=true;
	}
	int s=0,t=n+m+1,totans=n*m-k;
	for(int i=1;i<=n;i++)
		if(a[i]<0)
		{
			printf("JIONG\n");
			return 0;
		}
		else build_edge(s,i,a[i]);
	for(int i=1;i<=m;i++)
		if(b[i]<0)
		{
			printf("JIONG\n");
			return 0;
		}
		else build_edge(i+n,t,b[i]);
	for(int u=1;u<=n;u++)
		for(int v=1;v<=m;v++)
			if(!has[u][v]) build_edge(u,v+n,1);
//	for(int i=1;i<=n;i++) cout<<a[i]<<" ";cout<<endl;
//	for(int i=1;i<=m;i++) cout<<b[i]<<" ";cout<<endl;
	while(bfs(s,t))
	{
		for(int i=s;i<=t;i++)
			cur[i]=h[i];
		totans-=dfs(s,t,INF);
	}
	printf("%d\n",totans);
	return 0;
}