[SCOI2007]蜥蜴

在一个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;
}