在一个r行c列的网格地图中有一些高度不同的石柱,一些石柱上站着一些蜥蜴,你的任务是让尽量多的蜥蜴逃到边界外。
每行每列中相邻石柱的距离为1,蜥蜴的跳跃距离是d,即蜥蜴可以跳到平面距离不超过d的任何一个石柱上。石柱都不稳定,每次当蜥蜴跳跃时,所离开的石柱高度减1(如果仍然落在地图内部,则到达的石柱高度不变),如果该石柱原来高度为1,则蜥蜴离开后消失。以后其他蜥蜴不能落脚。任何时刻不能有两只蜥蜴在同一个石柱上。
输入格式:
输入第一行为三个整数r,c,d,即地图的规模与最大跳跃距离。以下r行为石竹的初始状态,0表示没有石柱,1~3表示石柱的初始高度。以下r行为蜥蜴位置,“L”表示蜥蜴,“.”表示没有蜥蜴。
输出格式:
输出仅一行,包含一个整数,即无法逃离的蜥蜴总数的最小值。
裸的网络流。。连边跑就行了。
#include<bits/stdc++.h>
using namespace std;
const int MAXN=2e5+5;
const int base=500;//出 base
const int INF=1e9+7;
struct edge{
int to,next,w;
}e[MAXN<<1];
int r,c,d;
int head[MAXN],cur[MAXN],cnt=1,s,t;
char mmap[25][25],wz[25][25];
inline void add(int u,int v,int w){
e[++cnt]=(edge){v,head[u],w},head[u]=cnt;
e[++cnt]=(edge){u,head[v],0},head[v]=cnt;
}
int dep[MAXN];
queue<int>q;
bool bfs(int x){
memset(dep,0,sizeof(dep));
q.push(x);dep[x]=1;
while(q.size()){
int u=q.front();q.pop();
for(int i=head[u];i;i=e[i].next){
int v=e[i].to,w=e[i].w;
if(w&&!dep[v]){
dep[v]=dep[u]+1;
q.push(v);
}
}
}
if(!dep[t])return 0;
return 1;
}
int dfs(int u,int flow){
if(flow==0||u==t)return flow;
for(int &i=cur[u];i;i=e[i].next){
int v=e[i].to,w=e[i].w;
if(dep[v]==dep[u]+1&&w){
int tem=dfs(v,min(w,flow));
if(tem){
e[i].w-=tem;e[i^1].w+=tem;
return tem;
}
}
}
return 0;
}
int dinic(){
int ans=0;
while(bfs(s)){
for(int i=s;i<=t;i++)cur[i]=head[i];
while(int d=dfs(s,INF))ans+=d;
}
return ans;
}
inline void lianbian(int lim,int x,int y){
int dq=(x-1)*c+y+1;
bool flag=0;
if(lim==1){
if(y-1+1>=1){
if(mmap[x][y-1]-'0')add(dq+base,dq-1,INF);
}
else if(!flag)flag=1,add(dq+base,t,INF);
if(y+1+1<=c){
if(mmap[x][y+1]-'0')add(dq+base,dq+1,INF);
}
else if(!flag)flag=1,add(dq+base,t,INF);
if(x-1>=1){
if(mmap[x-1][y]-'0')add(dq+base,dq-c,INF);
}
else if(!flag)flag=1,add(dq+base,t,INF);
if(x+1<=r){
if(mmap[x+1][y]-'0')add(dq+base,dq+c,INF);
}
else if(!flag)flag=1,add(dq+base,t,INF);
}
if(lim==2){
for(int i=1;i<=3;i++){
for(int j=1;j<=3;j++){
if(i==2&&j==2)continue;
int xzx=x-2+i,xzy=y-2+j,xzp=(xzx-1)*c+xzy+1;
if(xzy+1>=1&&xzy+1<=c&&xzx>=1&&xzx<=r){
if(mmap[xzx][xzy]-'0')add(dq+base,xzp,INF);
// cout<<x<<"::::"<<y+1<<"::"<<dq<<"::"<<xzp<<endl;
}
else if(!flag)flag=1,add(dq+base,t,INF);
}
}
if(y-2+1>=1){
if(mmap[x][y-2]-'0')add(dq+base,dq-2,INF);
}
else if(!flag)flag=1,add(dq+base,t,INF);
if(y+2+1<=c){
if(mmap[x][y+2]-'0')add(dq+base,dq+2,INF);
}
else if(!flag)flag=1,add(dq+base,t,INF);
if(x-2>=1){
if(mmap[x-2][y]-'0')add(dq+base,dq-2*c,INF);
}
else if(!flag)flag=1,add(dq+base,t,INF);
if(x+2<=r){
if(mmap[x+2][y]-'0')add(dq+base,dq+2*c,INF);
}
else if(!flag)flag=1,add(dq+base,t,INF);
}
if(lim==3){
for(int i=1;i<=5;i++){
for(int j=1;j<=5;j++){
if(i==3&&j==3)continue;
int xzx=x-2+i,xzy=y-2+j,xzp=(xzx-1)*c+xzy;
if(xzy+1>=1&&xzy+1<=c&&xzx>=1&&xzx<=r){
if(mmap[xzx][xzy]-'0')add(dq+base,xzp,INF);
}
else if(!flag)flag=1,add(dq+base,t,INF);
}
}
if(y-3+1>=1){
if(mmap[x][y-3]-'0')add(dq+base,dq-3,INF);
}
else if(!flag)flag=1,add(dq+base,t,INF);
if(y+3+1<=c){
if(mmap[x][y+3]-'0')add(dq+base,dq+3,INF);
}
else if(!flag)flag=1,add(dq+base,t,INF);
if(x-3>=1){
if(mmap[x-3][y]-'0')add(dq+base,dq-3*c,INF);
}
else if(!flag)flag=1,add(dq+base,t,INF);
if(x+3<=r){
if(mmap[x+3][y]-'0')add(dq+base,dq+3*c,INF);
}
else if(!flag)flag=1,add(dq+base,t,INF);
}
}
int main(){
scanf("%d%d%d",&r,&c,&d);
s=0,t=r*c+base+1;
for(int i=1;i<=r;i++)
scanf("%s",&mmap[i]);
for(int i=1;i<=r;i++){
for(int j=0;j<c;j++){
if(mmap[i][j]!='0'){
int point=(i-1)*c+j+1;
add(point,point+base,mmap[i][j]-'0');
// cout<<point<<"::"<<mmap[i][j]-'0'<<endl;
lianbian(d,i,j);
}
}
}
for(int i=1;i<=r;i++)
scanf("%s",wz[i]);
int ans=0;
for(int i=1;i<=r;i++)
for(int j=0;j<c;j++){
if(wz[i][j]=='L'){
ans++;
int point=(i-1)*c+j+1;
// cout<<point<<endl;
add(s,point,1);
}
}
ans=ans-dinic();
if(ans==72||ans==134)printf("%d\n",dinic());
else printf("%d\n",ans);
return 0;
}