BZOJ1066 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

HINT

100%的数据满足：1<=r, c<=20, 1<=d<=4

一道比较裸的（套路）最大流板子

把一个石柱点分成进和出两个部分，中间连一条容量为高度的边

然后就可以做完了

//yangkai

#include<bits/stdc++.h>

using namespace std;

#define INF 0x3f3f3f3f

#define N 10010

struct Edge{int u,v,cap,flow;};

struct Dinic{

    int s,t;

    int vis[N],d[N];

    vector<Edge> E;

    vector<int> G[N];

    void add(int u,int v,int cap){

        E.push_back((Edge){u,v,cap,0});

        E.push_back((Edge){v,u,0,0});

        int cnt=E.size();

        G[u].push_back(cnt-2);

        G[v].push_back(cnt-1);

    }

    bool BFS(){

        memset(vis,0,sizeof(vis));

        memset(d,0,sizeof(d));

        queue<int> Q;

        Q.push(s);

        vis[s]=1;

        while(!Q.empty()){

            int u=Q.front();Q.pop();

            for(int i=0;i<G[u].size();i++){

                Edge &e=E[G[u][i]];

                if(!vis[e.v]&&e.cap>e.flow){

                    vis[e.v]=1;

                    d[e.v]=d[u]+1;

                    Q.push(e.v);

                }

            }

        }

        return vis[t];

    }

    int DFS(int u,int a){

        if(u==t||!a)return a;

        int flow=0;

        for(int i=0;i<G[u].size();i++){

            Edge &e=E[G[u][i]];

            if(d[e.v]!=d[u]+1)continue;

            int f=DFS(e.v,min(a,e.cap-e.flow));

            e.flow+=f;

            E[G[u][i]^1].flow-=f;

            flow+=f;

            a-=f;

            if(!a)break;

        }

        return flow;

    }

    int Maxflow(){

        int flow=0;

        while(BFS())

            flow+=DFS(s,INF);

        return flow;

    }

}dinic;

int w,h,d,sum=0;

int g[30][30],p[30][30];

char op[30];

bool judge_in(int x1,int x2,int y1,int y2){

    return (x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)<=d*d;

}

int getid(int x,int y){

    return (x-1)*h+y;

}

void build(){

    dinic.s=0;dinic.t=w*h*2+1;

    for(int i=1;i<=w;i++)

        for(int j=1;j<=h;j++)

            if(p[i][j])dinic.add(dinic.s,getid(i,j)*2-1,1);

    for(int i=1;i<=w;i++)

        for(int j=1;j<=h;j++)

            if(i<=d||i>=(w-d+1)||j<=d||j>=(h-d+1))

                dinic.add(getid(i,j)*2,dinic.t,INF);

    for(int i=1;i<=w;i++)

        for(int j=1;j<=h;j++)

            for(int ii=1;ii<=w;ii++)

                for(int jj=1;jj<=h;jj++)

                    if(judge_in(i,ii,j,jj))

                        dinic.add(getid(i,j)*2,getid(ii,jj)*2-1,INF);

    for(int i=1;i<=w;i++)

        for(int j=1;j<=h;j++)

            if(g[i][j])dinic.add(getid(i,j)*2-1,getid(i,j)*2,g[i][j]);

}

int main(){

    cin>>w>>h>>d;

    for(int i=1;i<=w;i++){

        scanf("%s",op+1);

        for(int j=1;j<=h;j++)g[i][j]=op[j]-'0';

    }

    for(int i=1;i<=w;i++){

        scanf("%s",op+1);

        for(int j=1;j<=h;j++){

            p[i][j]=(op[j]=='L')?1:0;

            sum+=p[i][j];

        }

    }

    build();

    cout<<sum-dinic.Maxflow();

    return 0;

}