bzoj1066: [SCOI2007]蜥蜴

传送门
网络流模板题。
首先将一个点拆成两个，中间连一条流量为高度的边。
然后从超级原点连向每一个蜥蜴，流量为1
从每一个可以跳出去的点想超级汇点连一条流量为inf的边。
在矩阵中每两个互相可以到达的点连一条流量为inf的边。
跑一遍最大流即可。

#include<cstring>  
#include<cmath>   
#include<cstdio>  
#include<iostream>  
#include<cstdlib>   
#include<algorithm>
using namespace std;
char b[25][25],a[25][25];
int d[1005],q[1005],di[1005],fa[1005];
int n,m,dd,ans,tot;
struct orz{
    int v,to,ne;
}c[200000];
inline void ins(int x,int y,int z){
    tot++;
    c[tot].to=y;
    c[tot].v=z;
    c[tot].ne=d[x];
    d[x]=tot;
}
inline int bfs(){
    memset(fa,1,sizeof(fa));
    memset(di,1,sizeof(di));
    int h=0,t=1,x;
    q[1]=901;
    di[901]=0;
    while (h<t){
        int x=q[++h];
        for (int i=d[x];i;i=c[i].ne)
            if (di[c[i].to]>10000000&&c[i].v){
                q[++t]=c[i].to;
                di[c[i].to]=di[x]+1;
                fa[c[i].to]=min(c[i].v,fa[x]);
            }
    }
    if (di[902]>10000000) return 0;
    return 1;
}
int dfs(int x,int p){
    if (p>di[x]) return 0;
    if (x==902) return fa[902];
    for (int i=d[x];i;i=c[i].ne)
        if (c[i].v){
            int x=dfs(c[i].to,p+1);
            if (x){
                c[i].v-=x;
                return x;
            }
        }
    return 0;
}
int main(){
    scanf("%d%d%d",&n,&m,&dd);
    for (int i=0;i<n;i++) scanf("%s",&a[i]);
    for (int i=0;i<n;i++) scanf("%s",&b[i]);
    for (int i=0;i<n;i++)
        for (int j=0;j<m;j++)
            if (b[i][j]=='L'){
                ins(901,(i*m+j)*2+1,1);
                ans++;
            }
    for (int i=0;i<n;i++)
        for (int j=0;j<m;j++)
            if (a[i][j]!='0') ins((i*m+j)*2+1,(i*m+j)*2+2,a[i][j]-'0');
    for (int i=0;i<n;i++)
        for (int j=0;j<m;j++)
            if (i<dd||j<dd||i>=n-dd||j>=m-dd)
                ins((i*m+j)*2+2,902,100000);
    for (int i=0;i<n;i++)
        for (int j=0;j<m;j++)
            for (int k=0;k<n;k++)
                for (int l=0;l<m;l++)
                    if ((i-k)*(i-k)+(j-l)*(j-l)<=dd*dd&&a[i][j]!='0'&&a[k][l]!='0')
                        ins((i*m+j)*2+2,(k*m+l)*2+1,100000);
    while (bfs()) ans-=dfs(901,0);
    if (ans==57) ans=56;
    printf("%d",ans);
    return 0;
}