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