士兵占领

link

很巧妙的题。可以用上下界网络流来写但我不想。

正着推不是很好处理，因为正着那就是限定边的流量下界，还不如用上下界网络流呢。考虑反着来，先把每个不是障碍点的地方摆上士兵，可以计算出每一行和列的士兵人数，这样就可以找出每一行和列最多可以删除的士兵数量，这不就限定了流量的上界了吗。考虑每个流量都会牵涉到两个节点，而这两个节点又分属两个交集为空的集合中，考虑用二分图的模型来跑。由于删掉的士兵越多实际使用到的士兵就越少，跑最大流即可。

#include<bits/stdc++.h>
//#define feyn
using namespace std;
const int N=1010;
const int M=5e6;
const int maxn=1e9;
inline void read(int &wh){
	wh=0;int f=1;char w=getchar();
	while(w<'0'||w>'9'){if(w=='-')f=-1;w=getchar();}
	while(w>='0'&&w<='9'){wh=wh*10+w-'0';w=getchar();}
	wh*=f;return;
}

int m,n,num,aa[N],ab[N],ba[N],bb[N];

int a[N],b[N],ss,tt,cnt;

struct edge{
	int t,v,next;
}e[M];
int head[N],esum=1;
inline void adde(int fr,int to,int val){
	e[++esum]=(edge){to,val,head[fr]};head[fr]=esum;
}
inline void add(int fr,int to,int val){
	adde(fr,to,val);adde(to,fr,0);
}
int q[N],ll,rr,d[N],t[N],nt;
inline bool check(){
	q[ll=rr=1]=ss;t[ss]=++nt;d[ss]=1;
	while(ll<=rr){
		int wh=q[ll++];
		for(int i=head[wh],th;i;i=e[i].next){
			if(e[i].v==0||t[th=e[i].t]==nt)continue;
			t[th]=nt,d[th]=d[wh]+1,q[++rr]=th;
		}
	}
	return t[tt]==nt;
}
inline int dinic(int wh,int val){
	if(wh==tt)return val;
	int used=0;
	for(int i=head[wh],th;i;i=e[i].next){
		if(e[i].v==0||d[th=e[i].t]!=d[wh]+1)continue;
		int now=dinic(th,min(val,e[i].v));
		if(now)used+=now,val-=now,e[i].v-=now,e[i^1].v+=now;
		if(val==0)break;
	}
	if(val)d[wh]=-1;return used;
}

bool del[N][N];

signed main(){
	
	#ifdef feyn
	freopen("in.txt","r",stdin);
	#endif
	
	read(m);read(n);read(num);
	ss=++cnt,tt=++cnt;int s1,s2;
	for(int i=1;i<=m;i++)read(aa[i]),a[i]=++cnt;
	for(int i=1;i<=n;i++)read(ab[i]),b[i]=++cnt;
	for(int i=1;i<=num;i++){
		read(s1);read(s2);del[s1][s2]=true;
		ba[s1]++;bb[s2]++;
	}
	bool ok=false;
	for(int i=1;i<=m;i++)if(n-aa[i]-ba[i]<0)ok=true;
	for(int i=1;i<=n;i++)if(m-ab[i]-bb[i]<0)ok=true;
	if(ok){printf("JIONG!\n");return 0;}
	for(int i=1;i<=m;i++)add(ss,a[i],n-aa[i]-ba[i]);
	for(int i=1;i<=n;i++)add(b[i],tt,m-ab[i]-bb[i]);
	for(int i=1;i<=m;i++){
		for(int j=1;j<=n;j++){
			if(del[i][j])continue;
			add(a[i],b[j],1);
		}
	}
	int ans=0;
	while(check())ans+=dinic(ss,maxn);
	printf("%d\n",m*n-num-ans);
	
	return 0;
}