【bzoj 1066】蜥蜴（最大流）

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

思路

网络流经典题。建图时，把每根柱子拆成两个点，分别表示入和出。入点和出点之间用柱子的高度限流。把能互相到达的柱子连一条流量为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;
}