【bzoj 1066&poj 2711】 蜥蜴 【SCOI2007】

Description

　　在一个r行c列的网格地图中有一些高度不同的石柱，一些石柱上站着一些蜥蜴，你的任务是让尽量多的蜥蜴逃
到边界外。 每行每列中相邻石柱的距离为1，蜥蜴的跳跃距离是d，即蜥蜴可以跳到平面距离不超过d的任何一个石
柱上。石柱都不稳定，每次当蜥蜴跳跃时，所离开的石柱高度减1（如果仍然落在地图内部，则到达的石柱高度不
变），如果该石柱原来高度为1，则蜥蜴离开后消失。以后其他蜥蜴不能落脚。任何时刻不能有两只蜥蜴在同一个
石柱上。

Input

　　输入第一行为三个整数r，c，d，即地图的规模与最大跳跃距离。以下r行为石竹的初始状态，0表示没有石柱
，1~3表示石柱的初始高度。以下r行为蜥蜴位置，“L”表示蜥蜴，“.”表示没有蜥蜴。

Output

　　输出仅一行，包含一个整数，即无法逃离的蜥蜴总数的最小值。

Sample Input

5 8 2
00000000
02000000
00321100
02000000
00000000
........
........
..LLLL..
........
........

Sample Output

1

HINT

 

100%的数据满足：1<=r, c<=20, 1<=d<=4

这道题是网络流，将每个石柱拆成两个点，中间连一条容量为该石柱高度的边，设置超级源和汇，若石柱上有蜥蜴，则将与超级源连一条容量为1的边，若从该石柱可以跳出，则将与超级汇连一条容量为INF的边，对于相互可达的两个石柱，在和间连一条边，然后跑一遍dinic即可，下面是程序：

#include<stdio.h>
#include<string.h>
#include<iostream>
using namespace std;
const int N=805,M=500005,INF=0x7fffffff;
struct queue{
	int l,r,a[N];
	void clear(){
		l=0,r=1;
	}
	bool empty(){
		return l+1==r;
	}
	void push(int x){
		a[r]=x;
		r=(r+1)%N;
	}
	void pop(){
		l=(l+1)%N;
	}
	int front(){
		return a[(l+1)%N];
	}
}q;
struct edge{
	int v,w,next;
}e[M];
int n,m,c,head[N],k,d[N],map[25][25],mark[25][25];
void add(int u,int v,int w){
	e[k]=(edge){v,w,head[u]};
	head[u]=k++;
	e[k]=(edge){u,0,head[v]};
	head[v]=k++;
}
bool BFS(){
	int i;
	q.clear();
	memset(d,0,sizeof(d));
	d[0]=1;
	q.push(0);
	while(!q.empty()){
		int u=q.front();
		q.pop();
		for(i=head[u];i!=-1;i=e[i].next){
			if(!d[e[i].v]&&e[i].w){
				d[e[i].v]=d[u]+1;
				q.push(e[i].v);
			}
		}
	}
	return d[N-1];
}
int dfs(int u,int f){
	if(f==0||u==N-1){
		return f;
	}
	int i,t,use=0;
	for(i=head[u];i!=-1;i=e[i].next){
		if(d[e[i].v]==d[u]+1&&(t=dfs(e[i].v,min(f-use,e[i].w)))){
			e[i].w-=t;
			e[i^1].w+=t;
			use+=t;
			if(use==f){
				return f;
			}
		}
	}
	if(!use){
		d[u]=-1;
	}
	return use;
}
int dinic(){
	int ans=0;
	while(BFS()){
		ans+=dfs(0,INF);
	}
	return ans;
}
char ch(){
	char c=getchar();
	while((c<'0'||c>'4')&&c!='L'&&c!='.'){
		c=getchar();
	}
	return c;
}
int min(const int a,const int b){
	return a<b?a:b;
}
int abs(const int n){
	return n>0?n:-n;
}
bool check(int &x1,int &y1,int &x2,int &y2){
	return map[x1][y1]&&map[x2][y2]&&abs(x1-x2)+abs(y1-y2)<=c;
}
int main(){
	int t=0,i,j,k,l,cnt=0;
	scanf("%d%d%d",&n,&m,&c);
	memset(head,-1,sizeof(head));
	for(i=1;i<=n;i++){
		for(j=1;j<=m;j++){
			mark[i][j]=++t;
		}
	}
	for(i=1;i<=n;i++){
		for(j=1;j<=m;j++){
			map[i][j]=ch()-'0';
		}
	}
	for(i=1;i<=n;i++){
		for(j=1;j<=m;j++){
			if(ch()=='L'){
				add(0,mark[i][j],1);
				cnt++;
			}
			if(map[i][j]){
				add(mark[i][j],mark[i][j]+t,map[i][j]);
				if(min(min(i-1,j-1),min(n-i,m-j))<c){
					add(mark[i][j]+t,N-1,INF);
				}
			}
			for(k=1;k<=n;k++){
				for(l=1;l<=m;l++){
					if((i!=k||j!=l)&&check(i,j,k,l)){
						add(mark[i][j]+t,mark[k][l],INF);
					}
				}
			}
		}
	}
	printf("%d\n",cnt-dinic());
	return 0;
}