HDU 2732 Leapin' Lizards（最大流）

Leapin' Lizards

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=2732

题目意思看起来比较难懂。

题目是说一个n*m的迷宫中，有每个格子有柱子。柱子高度为0~3，高度为0的柱子是不能站的（高度为0就是没有柱子）

在一些有柱子的格子上有一些蜥蜴，一次最多跳距离d,相邻格子的距离是1，只要跳出迷宫就是安全的。

这个距离是曼哈顿距离（好像是的）。

蜥蜴一次最多跳距离d,但是起跳的地方的柱子高度会减一，一个柱子同一时间只能有一个蜥蜴

要求最少几个不能逃出迷宫。

具体意思看题目吧，很多描述不清呢~~还是英语好啊

做法就是最大流。

要拆点，两个点的容量就是柱子的高度。

加一个源点和汇点。

//============================================================================
// Name        : HDU.cpp
// Author      : 
// Version     :
// Copyright   : Your copyright notice
// Description : Hello World in C++, Ansi-style
//============================================================================

#include <iostream>
#include <string.h>
#include <stdio.h>
#include <algorithm>
using namespace std;

const int MAXN=2200;
const int MAXM=200020;
const int INF=0x3f3f3f3f;
struct Node
{
    int to,next,cap;
}edge[MAXM];//注意是MAXM
int tol;
int head[MAXN];
int gap[MAXN],dis[MAXN],pre[MAXN],cur[MAXN];

void init()
{
    tol=0;
    memset(head,-1,sizeof(head));
}
void addedge(int u,int v,int w,int rw=0)
{
    edge[tol].to=v;edge[tol].cap=w;edge[tol].next=head[u];head[u]=tol++;
    edge[tol].to=u;edge[tol].cap=rw;edge[tol].next=head[v];head[v]=tol++;
}
int sap(int start,int end,int nodenum)
{
    memset(dis,0,sizeof(dis));
    memset(gap,0,sizeof(gap));
    memcpy(cur,head,sizeof(head));
    int u=pre[start]=start,maxflow=0,aug=-1;
    gap[0]=nodenum;
    while(dis[start]<nodenum)
    {
        loop:
        for(int &i=cur[u];i!=-1;i=edge[i].next)
        {
            int v=edge[i].to;
            if(edge[i].cap&&dis[u]==dis[v]+1)
            {
                if(aug==-1||aug>edge[i].cap)
                    aug=edge[i].cap;
                pre[v]=u;
                u=v;
                if(v==end)
                {
                    maxflow+=aug;
                    for(u=pre[u];v!=start;v=u,u=pre[u])
                    {
                        edge[cur[u]].cap-=aug;
                        edge[cur[u]^1].cap+=aug;
                    }
                    aug=-1;
                }
                goto loop;
            }
        }
        int mindis=nodenum;
        for(int i=head[u];i!=-1;i=edge[i].next)
        {
            int v=edge[i].to;
            if(edge[i].cap&&mindis>dis[v])
            {
                cur[u]=i;
                mindis=dis[v];
            }
        }
        if((--gap[dis[u]])==0)break;
        gap[dis[u]=mindis+1]++;
        u=pre[u];
    }
    return maxflow;
}

char g1[30][30];
char g2[30][30];
int mat[30][30];

int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    int T;
    int n,d;
    int iCase=0;
    scanf("%d",&T);
    while(T--)
    {
        iCase++;
        scanf("%d%d",&n,&d);
        init();
        int tol=0;
        for(int i=0;i<n;i++)
            scanf("%s",&g1[i]);
        int m=strlen(g1[0]);
        memset(mat,0,sizeof(mat));
        for(int i=0;i<n;i++)
            for(int j=0;j<m;j++)
                if(g1[i][j]>'0')
                {
                    mat[i][j]=++tol;
                    addedge(2*tol-1,2*tol,g1[i][j]-'0');
                }

        int start=0,end=2*tol+1,nodenum=2*tol+2;
        //进行拆点，加上源点和汇点，共2*tol+2个点。
        int sum=0;//总数
        for(int i=0;i<n;i++)
        {
            scanf("%s",&g2[i]);
            for(int j=0;j<m;j++)
                if(g2[i][j]=='L')
                {
                    sum++;
                    addedge(start,2*mat[i][j]-1,1);
                }
        }
        for(int i=0;i<n;i++)
            for(int j=0;j<m;j++)
                if(mat[i][j])
                {
                    for(int x=-d;x<=d;x++)
                        for(int y=abs(x)-d;y<=d-abs(x);y++)
                        {
                            int newi=i+x;
                            int newj=j+y;
                            if(newi<0||newi>=n||newj<0||newj>=m)continue;
                            if(mat[newi][newj]==0)continue;
                            if(newi==i&&newj==j)continue;
                            addedge(2*mat[i][j],2*mat[newi][newj]-1,INF);
                        }
                    if(i<d||j<d||n-i<=d||m-j<=d)
                        addedge(2*mat[i][j],end,INF);
                }
        int ans=sum-sap(start,end,nodenum);
        if(ans==0)printf("Case #%d: no lizard was left behind.\n",iCase);
        else if(ans==1)printf("Case #%d: 1 lizard was left behind.\n",iCase);
        else printf("Case #%d: %d lizards were left behind.\n",iCase,ans);
    }
    return 0;
}