2018.10.11 bzoj1458: 士兵占领（最大流）

传送门
最大流简单题。

---

正向做不是很好做。
考虑反过来求出最多可以删去几个点。
这时可以定义行列每匹配一个表示删去了一个点。
然后利用原题的限制定义容量就可以了。
代码：

#include<bits/stdc++.h>
#define N 205
#define M 200005
using namespace std;
inline int read(){
	int ans=0;
	char ch=getchar();
	while(!isdigit(ch))ch=getchar();
	while(isdigit(ch))ans=(ans<<3)+(ans<<1)+(ch^48),ch=getchar();
	return ans;
}
int n,m,k,l[N],h[N],ans,L[N],H[N];
struct edge{int v,next,c;};
bool matrix[N][N];
struct Dinic{
	edge e[M<<1];
	int first[N],cur[N],cnt,s,t,d[N];
	inline void addedge(int u,int v,int c){e[++cnt].v=v,e[cnt].c=c,e[cnt].next=first[u],first[u]=cnt;}
	inline void add(int u,int v,int c){addedge(u,v,c),addedge(v,u,0);}
	inline void init(){memset(first,-1,sizeof(first)),cnt=-1,s=0,t=n+m+1;}
	inline bool bfs(){
		queue<int>q;
		for(int i=1;i<=t;++i)d[i]=-1;
		q.push(s);
		while(!q.empty()){
			int x=q.front();
			q.pop();
			for(int i=first[x];~i;i=e[i].next){
				int v=e[i].v;
				if(~d[v]||!e[i].c)continue;
				d[v]=d[x]+1,q.push(v);
			}
		}
		return ~d[t];
	}
	inline int dfs(int x,int f){
		if(x==t||!f)return f;
		int flow=f;
		for(int&i=cur[x];~i;i=e[i].next){
			int v=e[i].v;
			if(!flow)break;
			if(e[i].c&&d[v]==d[x]+1){
				int tmp=dfs(v,min(flow,e[i].c));
				if(!tmp)d[v]=-1;
				e[i].c-=tmp,e[i^1].c+=tmp,flow-=tmp;
			}
		}
		return f-flow;
	}
	inline int solve(){
		int ret=0;
		while(bfs())memcpy(cur,first,sizeof(first)),ret+=dfs(s,0x3f3f3f3f);
		return ret;
	}
}dinic;
int main(){
	n=read(),m=read(),k=read(),ans=n*m,dinic.init();
	for(int i=1;i<=n;++i)h[i]=read(),H[i]=m;
	for(int i=1;i<=m;++i)l[i]=read(),L[i]=n;
	while(k--){
		int x=read(),y=read();
		matrix[x][y]=1,--ans,--H[x],--L[y];
		if(H[x]<h[x]||L[y]<l[y])return puts("JIONG!"),0;
	}
	for(int i=1;i<=n;++i)dinic.add(dinic.s,i,H[i]-h[i]);
	for(int i=1;i<=m;++i)dinic.add(i+n,dinic.t,L[i]-l[i]);
	for(int i=1;i<=n;++i)for(int j=1;j<=m;++j)if(!matrix[i][j])dinic.add(i,j+n,1);
	cout<<ans-dinic.solve();
	return 0;
}