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..
........
........
00000000
02000000
00321100
02000000
00000000
........
........
..LLLL..
........
........
Sample Output
1
HINT
100%的数据满足:1<=r, c<=20, 1<=d<=3
很裸的拆点网络流。。每个点拆成u和v
v连可以到达的u,容量INF。如果当前点可以跳出边界则v连汇点,容量INF
有蜥蜴的地方从源点连u,容量为1。
T T拆点的序号弄错了结果找了半个多小时错。。
#include<cstdio>
#include<string>
#include<cstring>
using namespace std;
int head[100001];
struct map
{
int f;
int s,t;
int next;
}a[400001];
int b[5001];
int edge;
int p;
int q[400001],d[400001];
inline void add(int s,int t,int f)
{
a[edge].next=head[s];
head[s]=edge;
a[edge].s=s;
a[edge].t=t;
a[edge].f=f;
}
inline bool bfs()
{
int l=0,r=0;
memset(q,0,sizeof(q));
r++;
q[r]=0;
memset(d,-1,sizeof(d));
d[0]=0;
int i,k;
while(l<r)
{
l++;
int k=q[l];
for(i=head[k];i!=0;i=a[i].next)
{
if(a[i].f>0&&d[a[i].t]==-1)
{
d[a[i].t]=d[k]+1;
r++;
q[r]=a[i].t;
}
}
}
if(d[p]>=0)
return true;
return false;
}
inline int dfs(int k,int s)
{
if(k==p)
return s;
int t=s;
int i;
for(i=head[k];i!=0;i=a[i].next)
{
if(d[a[i].t]==d[k]+1&&a[i].f>0)
{
int xx=dfs(a[i].t,min(s,a[i].f));
a[i].f-=xx;
if(i%2==0)
a[i-1].f+=xx;
else
a[i+1].f+=xx;
s-=xx;
}
}
return t-s;
}
inline int maxflow()
{
int s=0;
while(bfs())
s+=dfs(0,2100000000);
return s;
}
inline int abs(int x)
{
if(x<0)
x=-x;
return x;
}
int map[21][21];
int main()
{
int r,c,d;
scanf("%d%d%d",&r,&c,&d);
p=r*c*2+1;
int pp=r*c;
int i,j;
char x;
for(i=1;i<=r;i++)
{
scanf("%c",&x);
for(j=1;j<=c;j++)
{
scanf("%c",&x);
map[i][j]=x-'0';
if(map[i][j]!=0)
{
edge++;
add((i-1)*c+j,(i-1)*c+j+pp,map[i][j]);
edge++;
add((i-1)*c+j+pp,(i-1)*c+j,0);
}
}
}
int ii,jj;
for(i=1;i<=r;i++)
{
for(j=1;j<=c;j++)
{
if(map[i][j]>0)
{
for(ii=1;ii<=r;ii++)
{
for(jj=1;jj<=c;jj++)
{
if(map[ii][jj]>0)
{
if((i-ii)*(i-ii)+(j-jj)*(j-jj)<=d*d)
{
edge++;
add((i-1)*c+j+pp,(ii-1)*c+jj,2100000000);
edge++;
add((ii-1)*c+jj,(i-1)*c+j+pp,0);
}
}
}
}
if(i<=d||r-i<d||j<=d||c-j<d)
{
edge++;
add((i-1)*c+j+pp,p,2100000000);
edge++;
add(p,2*((i-1)*c+j),0);
}
}
}
}
int sum=0;
for(i=1;i<=r;i++)
{
scanf("%c",&x);
for(j=1;j<=c;j++)
{
scanf("%c",&x);
if(x=='L')
{
sum++;
edge++;
add(0,(i-1)*c+j,1);
edge++;
add((i-1)*c+j,0,0);
}
}
}
printf("%d\n",sum-maxflow());
return 0;
}