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

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

妈妈我再也不会把数组开小了Ｏ（≧口≦）Ｏ

ＣＯＤＥ：

#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;   }