bzoj 1066 蜥蜴 最大流

套路拆点，把石柱高度作为出点入点之间流量，最大流。

#include<cstdio>
#include<iostream>
#include<queue>
#include<algorithm>
using namespace std;
struct bian{
    int to,cap,nxt;
}b[64005];
int fst[1005],tot=1;
int cur[1005];
void build(int f,int t,int d)
{
    b[++tot].to=t;
    b[tot].cap=d;
    b[tot].nxt=fst[f];
    fst[f]=tot;
}
int r,c,d;
void insert(int f,int t,int d)
{
    build(f,t,d);
    build(t,f,d);
}
char s[25];
int mp[250][250];
int p(int i,int j)
{return (i-1)*c+j;}
int inf=1e8;
int S=0,T;
bool check(int i,int j)
{
    if(i<=d||j<=d)
        return true;
    if(i>r-d||j>c-d)
        return true;
    return false;
}
queue<int> q;
int dis[1005];
bool bfs()
{
    for(int i=0;i<=T;i++)
        dis[i]=0;
    dis[S]=1;
    q.push(S);
    int u,v;
    while(!q.empty())
    {
        u=q.front();q.pop();
        for(int i=fst[u];i;i=b[i].nxt)
        {
            v=b[i].to;
            if(b[i].cap&&!dis[v])
            {
                dis[v]=dis[u]+1;
                q.push(v);
            }
        }
    }
    if(dis[T])
        return true;
    return false;
}
int dfs(int x,int a)
{
    if(x==T||a==0) return a;
    int v,ans=0;
    for(int &i=cur[x];i;i=b[i].nxt)
    {
        v=b[i].to;
        if(dis[v]==dis[x]+1&&b[i].cap)
        {
            int f=dfs(v,min(a,b[i].cap));
            ans+=f;
            a-=f;
            b[i].cap-=f;
            b[i^1].cap+=f;
        }
    }
    if(!ans) dis[x]=inf;
    return ans;
}
int maxflow()
{
    int ans=0;
    while(bfs())
    {
        for(int i=0;i<=T;i++)
            cur[i]=fst[i];
        ans+=dfs(S,inf);
    }
    return ans;
}
bool judge(int i,int j)
{   return (i>0&&i<=r&&j>0&&j<=c);}
int G[405];
int main()
{
    scanf("%d%d%d",&r,&c,&d);
    T=r*c*2+1;
    for(int i=1;i<=r;i++)
    {
        scanf("%s",s+1);
        for(int j=1;j<=c;j++)
            mp[i][j]=s[j]-'0';
    }
    for(int i=1;i<=r;i++)
    for(int j=1;j<=c;j++)
    {
        if(mp[i][j])
        {
            int q=p(i,j);
            build(q,q+r*c,mp[i][j]);
            build(q+r*c,q,0);
            for(int I=-d;I<=d;I++)
            for(int J=-d;J<=d;J++)
            {
                if(abs(I)+abs(J)>d)continue;
                if(!judge(i+I,j+J))continue;
                if(I==0&&J==0)continue;
                if(!mp[i+I][j+J])continue;
                build(r*c+p(i,j),p(i+I,j+J),inf);
                build(p(i+I,j+J),r*c+p(i,j),0);
            }
        }
    }
    int tot=0;
    for(int i=1;i<=r;i++)
    {
        scanf("%s",s+1);
        for(int j=1;j<=c;j++)
            if(s[j]=='L')
            {
                build(S,p(i,j),1);
                build(p(i,j),S,0);
                tot++;
            }
    }
    for(int i=1;i<=r;i++)
    for(int j=1;j<=c;j++)
    {
        if(check(i,j))
        {
            build(p(i,j)+r*c,T,inf);
            build(T,p(i,j)+r*c,0);
        }
    }
    printf("%d",tot-maxflow());
    return 0;
}