BZOJ1066[SCOI2007]蜥蜴

1066: [SCOI2007]蜥蜴

Time Limit: 1 Sec   Memory Limit: 162 MB
Submit: 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..
........
........

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