【bzoj 1066】蜥蜴（最大流）

传送门biu~
网络流经典题。建图时，把每根柱子拆成两个点，分别表示入和出。入点和出点之间用柱子的高度限流。把能互相到达的柱子连一条流量为INF的弧，再把能让蜥蜴跳出地图的柱子连一条到汇点流量为INF的弧，把源点连向所有一开始有蜥蜴的柱子，流量为INF。跑一遍最大流就是能跳出地图的蜥蜴的最大数目，用蜥蜴总数减去即可。

#include<bits/stdc++.h>
#define S (2*n*m+1)
#define T (2*n*m+2)
using namespace std;
const int INF=1e9;
int n,m,d,num;
char a[25][25];
int head[1005],nex[50005],to[50005],cap[50005],tp=1;
int fir[1005],dep[1005];
void add(int x,int y,int z){
    nex[++tp]=head[x];
    head[x]=tp;
    to[tp]=y;
    cap[tp]=z;
    nex[++tp]=head[y];
    head[y]=tp;
    to[tp]=x;
    cap[tp]=0;
}
int dfs(int x,int now){
    if(x==T || !now)    return now;
    int c=0;
    for(int &i=fir[x];i;i=nex[i]){
        if(dep[to[i]]!=dep[x]+1)    continue;
        int f=dfs(to[i],min(now,cap[i]));
        now-=f;
        cap[i]-=f;
        cap[i^1]+=f;
        c+=f;
        if(now==0)  break;
    }
    return c;
}
int bfs(){
    memset(dep,0,sizeof(dep));
    queue<int>q;
    q.push(S);dep[S]=1;
    while(!q.empty()){
        int x=q.front();q.pop();
        for(int i=head[x];i;i=nex[i]){
            if(!dep[to[i]] && cap[i]){
                dep[to[i]]=dep[x]+1;
                q.push(to[i]);
            }
        }
    }
    return dep[T];
}
int Dinic(){
    int c=0;
    while(bfs()){
        for(int i=1;i<=2*n*m+2;++i)     fir[i]=head[i];
        c+=dfs(S,INF);
    }
    return c;
}
int main(){
    scanf("%d%d%d",&n,&m,&d);
    for(int i=1;i<=n;++i)   scanf("%s",a[i]+1);
    for(int i=1;i<=n;++i){
        for(int j=1;j<=m;++j){
            if(a[i][j]!='0'){
                add((i-1)*m+j,n*m+(i-1)*m+j,a[i][j]-'0');
                if(min(min(i,n+1-i),min(j,m+1-j))<=d)   add(n*m+(i-1)*m+j,T,INF);
                for(int p=-d;p<=d;++p){
                    for(int q=-d;q<=d;++q){
                        if(abs(p)+abs(q)>d) continue;
                        int nx=i+p,ny=j+q;
                        if(nx<1 || nx>n || ny<1 || ny>m || a[nx][ny]=='0')  continue;
                        add(n*m+(i-1)*m+j,(nx-1)*m+ny,INF);
                    }
                }
            }
        }
    }
    for(int i=1;i<=n;++i)   scanf("%s",a[i]+1);
    for(int i=1;i<=n;++i){
        for(int j=1;j<=m;++j){
            if(a[i][j]=='L'){
                add(S,(i-1)*m+j,1);
                ++num;
            }
        }
    }
    printf("%d",num-Dinic());
    return 0;
}