bzoj1066 [SCOI2007]蜥蜴（建图最大流）

把每个点拆成两个点，一个表示进入这个点，记作a，一个表示离开这个点，记作b。a->b，边权为高度，表示最多只有这么多的蜥蜴可以通过这个点。然后任意两个可以互达的石柱x,y，建边bx->ay,边权为inf，建立超级源点，向所有一开始有蜥蜴的石柱建边，边权为inf，建立超级汇点，所有可以跳出地图的点向超级汇点连边，边权为inf。然后跑dinic就好啦。

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define N 1000
#define inf 0x3f3f3f3f
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,m,d,mp[21][21],h[N],num=1,tot=0,lev[N];
struct edge{
    int to,next,val;
}data[170010];
inline void add(int x,int y,int val){
    data[++num].to=y;data[num].next=h[x];h[x]=num;data[num].val=val;
    data[++num].to=x;data[num].next=h[y];h[y]=num;data[num].val=0;
}
inline bool bfs(){
    queue<int>q;memset(lev,0,sizeof(lev));
    q.push(0);lev[0]=1;
    while(!q.empty()){
        int x=q.front();q.pop();
        for(int i=h[x];i;i=data[i].next){
            int y=data[i].to;if(lev[y]||!data[i].val) continue;
            lev[y]=lev[x]+1;q.push(y);
        }
    }return lev[900];
}
int dinic(int x,int low){
    if(x==900) return low;int tmp=low;
    for(int i=h[x];i;i=data[i].next){
        int y=data[i].to;if(lev[y]!=lev[x]+1||!data[i].val) continue;
        int res=dinic(y,min(low,data[i].val));
        if(!res) lev[y]=0;tmp-=res;data[i].val-=res;data[i^1].val+=res;
        if(!tmp) return low;
    }return low-tmp;
}
inline int maxflow(){
    int res=0;
    while(bfs()) res+=dinic(0,inf);return res;
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();d=read();
    for(int i=1;i<=n;++i){
        char s[30];scanf("%s",s+1);
        for(int j=1;j<=m;++j) mp[i][j]=s[j]-'0';
    }for(int x=1;x<=n;++x)
        for(int y=1;y<=m;++y){
            if(!mp[x][y]) continue;
            if(x-d<1||y-d<1||x+d>n||y+d>m) add((x-1)*m+y,900,inf);
        }
    for(int x1=1;x1<=n;++x1)
        for(int y1=1;y1<=m;++y1){
            if(!mp[x1][y1]) continue;add((x1-1)*m+y1+400,(x1-1)*m+y1,mp[x1][y1]);
            for(int x2=1;x2<=n;++x2)
                for(int y2=1;y2<=m;++y2){
                    if(!mp[x2][y2]||(x1==x2&&y1==y2)||abs(x2-x1)+abs(y2-y1)>d) continue;
                    add((x1-1)*m+y1,(x2-1)*m+y2+400,inf);
                }
        }
    for(int i=1;i<=n;++i){
        char s[30];scanf("%s",s+1);
        for(int j=1;j<=m;++j) if(s[j]=='L') add(0,(i-1)*m+j+400,1),++tot;
    }printf("%d\n",tot-maxflow());
    return 0;
}