[SCOI2007] 蜥蜴

link

之前写过，今天考到最后挂了，原因竟然是数组开小了，就很玄妙。

就把每个点拆成入点和出点，中间连边边权是桩子上的最大数量。然后正常到达的点之间连无线边，建立超级源点给每个蜥蜴的位置连1边，跑网络流即可。数据范围挺小的。

好久没写网络流了……

#include<bits/stdc++.h>
//#define feyn
const int N=40010;
const int M=25;
const int maxn=1e7;
using namespace std;
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;
}
inline int get(){
	char w=getchar();
	while(w<'0'||w>'9')w=getchar();
	return w-'0';
}
inline bool gget(){
	char w=getchar();
	while(w!='.'&&w!='L')w=getchar();
	return w=='L';
}
inline int min(int s1,int s2){
	return s1<s2?s1:s2;
}

int m,n,dis,cnt,r[M][M],c[M][M],a[M][M];
int ss,tt,num;

struct edge{
	int t,v,next;
}e[N<<4];
int head[N],esum=1;
inline void add(int fr,int to,int val){
	e[++esum]=(edge){to,val,head[fr]};head[fr]=esum;
}
inline void adde(int fr,int to,int val){
	
	add(fr,to,val);
	add(to,fr,0);
}

int q[N],ll,rr;
int t[N],nt,d[N];
inline bool check(){
	nt++;t[ss]=nt;d[ss]=1;
	q[ll=rr]=ss;
	while(ll<=rr){
		int wh=q[ll++];
		for(int i=head[wh],th;i;i=e[i].next){
			if(e[i].v==0)continue;th=e[i].t;
			if(t[th]^nt){t[th]=nt;q[++rr]=th;d[th]=d[wh]+1;}
		}
	}
	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)continue;
		if(d[wh]+1!=d[th=e[i].t])continue;
		int now=dinic(th,min(val,e[i].v));
		used+=now;val-=now;
		e[i].v-=now;e[i^1].v+=now;
		if(val==0)break;
	}
	if(val)d[wh]=-1;
	return used;
}
void solve(){
	int ans=0;
	while(check()){
		ans+=dinic(ss,maxn);
	}
	printf("%d\n",num-ans);
	return;
}

signed main(){
	
	#ifdef feyn
	freopen("in.txt","r",stdin);
	#endif
	
	read(m);read(n);read(dis);
	for(int i=1;i<=m;i++){
		for(int j=1;j<=n;j++){
			a[i][j]=get();
			r[i][j]=++cnt;c[i][j]=++cnt;
			if(a[i][j])adde(r[i][j],c[i][j],a[i][j]);
		}
	}
	ss=++cnt;tt=++cnt;
	for(int i=1;i<=m;i++){
		for(int j=1;j<=n;j++){
			if(a[i][j]==0)continue;
			for(int s1=-3;s1<=3;s1++){
				for(int s2=-3;s2<=3;s2++){
					if(s1==0&&s2==0)continue;
					if(s1*s1+s2*s2>dis*dis)continue;
					int x=i+s1,y=j+s2;
					if(x<1||y<1||x>m||y>n)adde(c[i][j],tt,maxn);
					else adde(c[i][j],r[x][y],maxn);
				}
			}
		}
	}
	for(int i=1;i<=m;i++){
		for(int j=1;j<=n;j++){
			if(gget())adde(ss,r[i][j],1),num++;
		}
	}
	solve();
	
	return 0;
}