BZOJ-1066 蜥蜴 最大流+拆点+超级源超级汇

1066: [SCOI2007]蜥蜴
Time Limit: 1 Sec Memory Limit: 162 MB
Submit: 2582 Solved: 1272
[Submit][Status][Discuss]

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

Source
Pku 2711 Leapin’ Lizards

这道题啊,我想啊。。(模仿自谢老板)
一眼就看出是网络流了,不过需要建图,一般网络流题目难点就在建图上,所以在建图上费了心
建立超级源,与各初始有蜥蜴的柱子连边,边权为1
建立超级汇,找可以跳出去的柱子与之连边,边权为INF
拆点:将每个柱子拆成两个点(入点出点)入点到出点边权为柱子高度
暴力将能够互相跳到的点连边(其中一点的出点连另一个点的入点),边权为INF
Dinic一遍得到最多能跳出的蜥蜴数,总数减掉即可。。。

妈妈我再也不会把数组开小了O(≧口≦)O

CODE:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int r,c,d;
int q[200000],h,t;
int dis[2000];
struct data{
    int to,next,v;
}edge[500001];
int cnt=1,head[2000]={0};
int n;
int zh[25][25]={0};// 记录柱子高度 
int xys[25][25]={0};//蜥蜴数 
int bh[25][25]={0}; //编号(建图用) 
int tot=0;//总蜥蜴数 

void add(int u,int v,int w)
{
    cnt++;
    edge[cnt].next=head[u];
    head[u]=cnt;
    edge[cnt].to=v;
    edge[cnt].v=w;
}

void init()
{
    scanf("%d%d%d",&r,&c,&d);
    for (int i=1; i<=r; i++)
        {
            char x[30];
            scanf("%s",&x);
            for (int j=1; j<=c; j++)
                zh[i][j]=x[j-1]-48;
        }
    for (int i=1; i<=r; i++)
        {
            char x[30];
            scanf("%s",&x);
            for (int j=1; j<=c; j++)
                if (x[j-1]=='L')
                    {
                        tot++;
                        xys[i][j]=1;
                    }
        }
}//读入 

void make()//超级源点为1,超级汇为n(最后的num) 
{
    int num=2;
    for (int i=1; i<=r; i++)
        for (int j=1; j<=c; j++)
            {
                if (zh[i][j]>0)
                    {
                        bh[i][j]=num;//bh【】记录这个柱子的入点的编号(出点编号为入点编号+1) 
                        add(num,num+1,zh[i][j]);
                        add(num+1,num,0);
                        num+=2;
                    }
            }//拆点 
    for (int i=1; i<=r; i++)
        for (int j=1; j<=c; j++)
            {
                if (xys[i][j]>0)
                    {
                        add(1,bh[i][j],1);
                        add(bh[i][j],1,0);
                    }
            }//把初始有蜥蜴的与超级源连边 
    for (int i=1; i<=r; i++)
        for (int j=1; j<=c; j++)
            {
                if (zh[i][j]>0)
                    {
                        for (int a=max(1,i-d); a<=min(r,i+d); a++)
                            for (int b=max(1,j-d); b<=min(c,j+d); b++)
                                {
                                    if (zh[a][b]>0 && (a!=i || b!=j))
                                        if ((a-i)*(a-i)+(b-j)*(b-j)<=d*d)
                                            {
                                                add(bh[i][j]+1,bh[a][b],0x7fffffff);
                                                add(bh[a][b],bh[i][j]+1,0);
                                            }
                                }                                               
                    }
            }//暴力把能调到的两个柱子连边 
    for (int i=1; i<=r; i++)
        for (int j=1; j<=c; j++)
            if (zh[i][j]>0 && (i-d<=0 || i+d>r || j-d<=0 || j+d>c))
                {
                    add(bh[i][j]+1,num,0x7fffffff);
                    add(num,bh[i][j]+1,0);
                }//把能跳出去的柱子与超级汇连边 
    n=num;
}//建图 

bool bfs()
{
    memset(dis,-1,sizeof(dis));
    q[1]=1;dis[1]=1;
    h=0; t=1;
    while (h<t)
        {
            int j=q[++h],i=head[j];
            while (i)
                {
                    if (dis[edge[i].to]<0 && edge[i].v>0)
                        {
                            dis[edge[i].to]=dis[j]+1;
                            q[++t]=edge[i].to;
                        }
                    i=edge[i].next;
                }
        }
    if (dis[n]>0)
        return true;
    else
        return false;
}

int dfs(int loc,int low)
{
    int ans=0;
    if (loc==n) return low;
    int i=head[loc];
    while (i)
        {
            if (edge[i].v>0 && dis[edge[i].to]==dis[loc]+1 && (ans=dfs(edge[i].to,min(low,edge[i].v))))
                {
                    edge[i].v-=ans;
                    edge[i^1].v+=ans;
                    return ans;
                }
            i=edge[i].next;
        }
    return 0;
}

int main()
{
    init();
    make();
    int ans=0;
    while (bfs())
        {
            int now;
            while ((now=dfs(1,0x7fffffff)))
                ans+=now;
        }
    printf("%d",tot-ans);
    return 0;   
}
posted @ 2016-01-02 17:37  DaD3zZ  阅读(274)  评论(0编辑  收藏  举报