1066: [SCOI2007]蜥蜴
Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 2410 Solved: 1179
[ 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..
........
........
00000000
02000000
00321100
02000000
00000000
........
........
..LLLL..
........
........
Sample Output
1
HINT
100%的数据满足:1<=r, c<=20, 1<=d<=4
最大流。
将每一个石柱拆成两个点i1,i2之间的边权为石柱的高度。
源点向有蜥蜴的所有石柱连边,边权为1
所有石柱若能够出边境则向汇点连边,边权inf
所有石柱i2向能到达的石柱i1连边,边权为inf;
最后求最大流。
有两个坑点:
1.石柱距离不是x^2+y^2
2.我一开始将head的初始值设为-1,但中间需要用到的时候判断为是否为0。于是超时很多次,怎么也查不出错。这就是不作死就不会死。。没事干为啥要赋值为-1呢。。哭瞎。。
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#define ll long long
#define inf 1e9
#define maxn 20050
using namespace std;
int tot=1;
struct Node{
int to,next,v;
}e[maxn];
int head[1010],r,c,d,a[21][21],s,t,h[maxn],cur[maxn],maxflow;
string ss;
queue<int> q;
int read()
{
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
void add(int x,int y,int v)
{
e[++tot]=(Node){y,head[x],v};head[x]=tot;
e[++tot]=(Node){x,head[y],0};head[y]=tot;
}
void makepic(int x,int y)
{
if(x<=d||y<=d||r-x<d||c-y<d)add((x-1)*c+y+r*c,r*c*2+5,inf);
for(int i=1;i<=r;i++)
for(int j=1;j<=c;j++)
{
if((i!=x||j!=y)&&(abs(i-x)+abs(j-y)<=d))
add(r*c+(x-1)*c+y,(i-1)*c+j,inf);
}
}
bool bfs()
{
for(int i=s;i<=t;i++)h[i]=-1;
q.push(s);h[s]=0;
while(!q.empty())
{
int x=q.front();q.pop();
for(int i=head[x];i;i=e[i].next)
if(e[i].v&&h[e[i].to]==-1)
{
h[e[i].to]=h[x]+1;q.push(e[i].to);
}
}
return h[t]!=-1;
}
int dfs(int x,int f)
{
if(x==t)return f;
int tmp,used=0;
for(int i=cur[x];i;i=e[i].next)
if(e[i].v&&h[e[i].to]==h[x]+1)
{
tmp=dfs(e[i].to,min(e[i].v,f-used));
e[i].v-=tmp;if(e[i].v)cur[x]=i;
e[i^1].v+=tmp;used+=tmp;
if(used==f)return f;
}
if(!used)h[x]=-1;
return used;
}
void dinic()
{
maxflow=0;
while(bfs())
{
for(int i=s;i<=t;i++)cur[i]=head[i];
maxflow+=dfs(s,inf);
}
}
int main()
{
r=read();c=read();d=read();
s=0;t=r*c*2+5;
int cnt=0;
memset(head,0,sizeof(head));
for(int i=1;i<=r;i++)
{
getline(cin,ss);
for(int j=0;j<c;j++)
{
a[i][j+1]=ss[j]-'0';
if(a[i][j+1]!=0)
{
int item=(i-1)*c+j+1;
add(item,r*c+item,a[i][j+1]);
}
}
}
for(int i=1;i<=r;i++)
{
getline(cin,ss);
for(int j=0;j<c;j++)
{
if(ss[j]=='L')cnt++,add(0,(i-1)*c+j+1,1);
if(a[i][j+1])makepic(i,j+1);
}
}
dinic();
cout<<cnt-maxflow<<endl;
return 0;
}