1066: [SCOI2007]蜥蜴

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

 

 

问题转化成最多有多少只蜥蜴可以逃出去。。就变成了最大流。。。石柱的高度只需，把石柱拆成两个点，中间连一条流量为高度的边。。。然后建图即可。。。

#include<iostream>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<string>
#include<map>
#include<queue>
#include<vector>
#include<set>
#define inf 1000000000
#define maxn 805
#define maxm 500000+5
#define eps 1e-10
#define ll long long
#define for0(i,n) for(int i=0;i<=(n);i++)
#define for1(i,n) for(int i=1;i<=(n);i++)
#define for2(i,x,y) for(int i=(x);i<=(y);i++)
#define for3(i,x,y) for(int i=(x);i>=(y);i--)
#define for4(i,x) for(int i=head[x],y=e[i].go;i;i=e[i].next,y=e[i].go)
using namespace std;
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=10*x+ch-'0';ch=getchar();}
    return x*f;
}
int  n,m,s,t,maxflow,tot=1,head[maxn],cur[maxn],h[maxn],cnt,sum;
int mp[25][25],mp2[25][25],id[25][25],r,c,d;
queue<int>q;
struct edge{int go,next,v;}e[maxm];
bool check(int x1,int y1,int x2,int y2){
    if((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)<=d*d)return 1;
    return 0;
}
void insert(int x,int y,int v)
{
    e[++tot]=(edge){y,head[x],v};head[x]=tot;
    e[++tot]=(edge){x,head[y],0};head[y]=tot;
}
bool bfs()
{
    for(int i=s;i<=t;i++)h[i]=-1;
    q.push(s);h[s]=0;
    while(!q.empty())
    {
        int x=q.front();q.pop();
        for(int i=head[x];i;i=e[i].next)
         if(e[i].v&&h[e[i].go]==-1)
         {
            h[e[i].go]=h[x]+1;q.push(e[i].go);
         }
    }
    return h[t]!=-1;
}
int dfs(int x,int f)
{
    if(x==t) return f;
    int tmp,used=0;
    for(int i=cur[x];i;i=e[i].next)
     if(e[i].v&&h[e[i].go]==h[x]+1)
    {
        tmp=dfs(e[i].go,min(e[i].v,f-used));
        e[i].v-=tmp;if(e[i].v)cur[x]=i;
        e[i^1].v+=tmp;used+=tmp;
        if(used==f)return f;       
    }
    if(!used) h[x]=-1;
    return used;
}
void dinic()
{
    maxflow=0;
    while(bfs())
    {
        for (int i=s;i<=t;i++)cur[i]=head[i];maxflow+=dfs(s,inf);
    }
}

int main(){
    //freopen("input.txt","r",stdin);
    //freopen("output.txt","w",stdout);
    r=read();c=read();d=read();s=0;t=801;
    for1(i,r)
        for1(j,c){
            sum++;id[i][j]=sum;
            char ch=getchar();
            while(ch<'0'||ch>'9')ch=getchar();
            mp[i][j]=ch-'0';
        }
    for1(i,r)
        for1(j,c){
            char ch=getchar();
            while(ch!='.'&&ch!='L')ch=getchar();
            if(ch=='L'){insert(0,id[i][j],1);cnt++;}
        }
    for1(i,r)
        for1(j,d){
            insert(id[i][j]+400,801,inf);
            insert(id[i][c-j+1]+400,801,inf);
        }
    for1(i,c)
        for1(j,d){
            insert(id[j][i]+400,801,inf);
            insert(id[r-j+1][i]+400,801,inf);
        }
    for(int x1=1;x1<=r;x1++)
        for(int y1=1;y1<=c;y1++)
            for(int x2=x1-d;x2<=x1+d;x2++)
                for(int y2=y1-d;y2<=y1+d;y2++)
                     if(check(x1,y1,x2,y2)&&(x1!=x2||y1!=y2))insert(id[x1][y1]+400,id[x2][y2],inf); 
    for1(i,r)
        for1(j,c)
            if(mp[i][j])insert(id[i][j],id[i][j]+400,mp[i][j]);
    dinic();
    printf("%d",cnt-maxflow);
    return 0;
}