BZOJ 1066 POJ 2711 [SCOI2007]蜥蜴

1066: [SCOI2007]蜥蜴

Time Limit: 1 Sec   Memory Limit: 162 MB
Submit: 2944   Solved: 1471
[ Submit][ Status][ Discuss]

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

Source

【题目分析】

   

参考了mhr的ppt才做出来orzzzzzzzz。。。感觉网络流里拆点的方法很有趣啊

【代码】

#include <cstdio>
#include <cstring>
#include <string>
#include <iostream>
#include <queue>
#define INF 10000000
#define M(a) memset(a,0,sizeof a)
#define F(i,j,k) for (int i=j;i<=k;++i)
using namespace std;
int list[51][51],num[51][51];
int n,m,d,en=0;
int S=0,T,sum=0,all=0;
int u[100001],v[100001],h[100001],ne[100001],f[100001],map[100001];
inline void add(int a,int b,int r)
{
	u[en]=a;v[en]=b;ne[en]=h[a];f[en]=r;h[a]=en++;
	u[en]=b;v[en]=a;ne[en]=h[b];f[en]=0;h[b]=en++;
}
inline bool tell()
{
	memset(map,-1,sizeof map);
	queue <int> q;
	q.push(S);
	map[S]=0;
	while (!q.empty()){
		int u=q.front(); q.pop();
		for (int i=h[u];i!=-1;i=ne[i]){
			if (map[v[i]]==-1&&f[i]){
				map[v[i]]=map[u]+1;
				q.push(v[i]);
			}
		}
	}
	if (map[T]!=-1) return true;
	else return false;
}
inline int zeng (int k,int now)
{
	if (k==T) return now;
	int r=0;
	for (int i=h[k];i!=-1&&now>r;i=ne[i]){
		if (map[k]+1==map[v[i]]&&f[i]!=0){
			int t=zeng(v[i],min(now-r,f[i]));
			f[i]-=t;f[i^1]+=t;r+=t;
		}
	}
	if (!r) map[k]=-1;
	return r;
}
inline int dinic()
{
	int r=0,t;
	while (tell()) while (t=zeng(S,INF)) r+=t;
	return r; 
}
inline bool ok(int a,int b,int c,int e)
{
	if ((c-a)*(c-a)+(e-b)*(e-b)<=d*d) return true;
	return false;
}
int main()
{
	memset(h,-1,sizeof h);
	memset(ne,-1,sizeof ne);
	memset(u,-1,sizeof u);
	memset(v,-1,sizeof v);
	scanf("%d%d%d",&n,&m,&d);
	F(i,1,n)F(j,1,m)
	{
		char ch;
		cin>>ch;
		list[i][j]=ch-'0';
		if (list[i][j]!=0){
			num[i][j]=++sum;
		}
	}
	T=sum*2+1;
	F(i,1,n)F(j,1,m)
	{
		char ch;
		cin>>ch;
		if (ch=='L') {
			add(S,num[i][j]*2,1);
			all++;
		}
	}
	F(i,1,n)F(j,1,m) if ((i<=d||j<=d||i+d>n||j+d>m)&&(num[i][j])){
		add(num[i][j]*2+1,T,INF);
	}
	F(i,1,n)F(j,1,m) if (num[i][j]) {
		add(num[i][j]*2,num[i][j]*2+1,list[i][j]);
	}
	F(i,1,n)F(j,1,m){
		if (num[i][j]){
			F(x,min(i-d,1),min(i+d,n)) F(y,min(j-d,1),min(j+d,m))
			  if (num[x][y]&&(i!=x||j!=y)) if (ok(i,j,x,y)) {
			  	add(num[i][j]*2+1,num[x][y]*2,INF);
			  }
		}
	}
	cout<<all-dinic()<<endl;
}