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。
对于地图内任意两个石柱,如果间距小于d,就将其中一根石柱的顶部与另一根石柱的底部相连,其连线容量为INF。
构图完成,剩下就是跑一遍最大流,然后用蜥蜴数量减去最大流就是最终结果。
by hzwer.com
代码
#include<cstdio>
#include<iostream>
#include<cstring>
#define inf 1000000000
int ret[20000],Next[20000],len[20000],gap[800],dis[8000];
int num,S,T,r,c,d,tot,ans,mp[21][21],mark[21][21],q[802],h[802],Head[8002];
using namespace std;
int insert(int u,int v,int l)
{
tot++;
ret[tot]=v;
len[tot]=l;
Next[tot]=Head[u];
Head[u]=tot;
}
bool judge(int x1,int y1,int x2,int y2)
{
if(((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2))<=(d*d)&&mp[x1][y1]&&mp[x2][y2])return 1;
return 0;
}
void build()
{
for(int x1=1;x1<=r;x1++)
for(int y1=1;y1<=c;y1++)
for(int x2=x1-d;x2<=x1+d;x2++)
for(int y2=y1-d;y2<=y1+d;y2++)
if(judge(x1,y1,x2,y2)&&(x1!=x2||y1!=y2))insert(mark[x1][y1]+400,mark[x2][y2],inf),insert(mark[x2][y2],mark[x1][y1]+400,0);
for(int i=1;i<=r;i++)
for(int j=1;j<=c;j++)
if(mp[i][j])insert(mark[i][j],mark[i][j]+400,mp[i][j]),insert(mark[i][j]+400,mark[i][j],0);
}
int isap(int u,int aug)
{
if (u==T) return aug;
if (aug==0) return 0;
int min_d=num-1,flow=0;
for (int i=Head[u];i;i=Next[i])
{
int v=ret[i];
if (len[i]>0)
{
if (dis[u]==dis[v]+1)
{
int t=isap(v,min(len[i],aug-flow));
len[i]-=t;
len[i^1]+=t;
flow+=t;
}
if (flow==aug) return flow;
if (dis[S]==num) return flow;
min_d=min(min_d,dis[v]);
}
}
if (flow==0)
{
gap[dis[u]]--;
if (gap[dis[u]]==0) dis[S]=num;
dis[u]=min_d+1;
gap[dis[u]]++;
}
return flow;
}
int main()
{
scanf("%d%d%d",&r,&c,&d);
char ch[21];
for (int i=1;i<=r;i++)
{
scanf("%s",ch);
for (int j=1;j<=c;j++)
mp[i][j]=ch[j-1]-'0';
}
int cnt=0;
for (int i=1;i<=r;i++)
{
for (int j=1;j<=c;j++)
{
cnt++;
mark[i][j]=cnt;
}
}
tot=1;S=r*c*2+400;
int sum=0;
for (int i=1;i<=r;i++)
{
scanf("%s",ch);
for (int j=1;j<=c;j++)
if (ch[j-1]=='L')
{
sum++;
insert(S,mark[i][j],1);
insert(mark[i][j],S,0);
}
}
T=r*c*2+401;
for (int i=1;i<=d;i++)
{
for (int j=d+1;j<=r-d;j++)
{
insert(mark[j][i]+400,T,inf);
insert(T,mark[j][i]+400,0);
insert(mark[j][c-i+1]+400,T,inf);
insert(T,mark[j][c-i+1]+400,0);
}
}
for (int i=1;i<=d;i++)
{
for (int j=1;j<=c;j++)
{
insert(mark[i][j]+400,T,inf);
insert(T,mark[i][j]+400,0);
insert(mark[r-i+1][j]+400,T,inf);
insert(T,mark[r-i+1][j]+400,0);
}
}
build();
num=r*c*2+2;
while (dis[S]<num)
{
ans+=isap(S,inf);
}
cout<<sum-ans<<endl;
}