[OpenJudge] 2.5基本算法之搜索 仙岛求药

一、原题

2727:仙岛求药

总时间限制: 1000ms 内存限制: 65536kB

描述

少年李逍遥的婶婶病了，王小虎介绍他去一趟仙灵岛，向仙女姐姐要仙丹救婶婶。叛逆但孝顺的李逍遥闯进了仙灵岛，克服了千险万难来到岛的中心，发现仙药摆在了迷阵的深处。迷阵由M×N个方格组成，有的方格内有可以瞬秒李逍遥的怪物，而有的方格内则是安全。现在李逍遥想尽快找到仙药，显然他应避开有怪物的方格，并经过最少的方格，而且那里会有神秘人物等待着他。现在要求你来帮助他实现这个目标。

下图显示了一个迷阵的样例及李逍遥找到仙药的路线.

输入

输入有多组测试数据. 每组测试数据以两个非零整数 M 和 N 开始，两者均不大于20。M 表示迷阵行数, N 表示迷阵列数。接下来有 M 行, 每行包含N个字符,不同字符分别代表不同含义:
1) ‘@’：少年李逍遥所在的位置；
2) ‘.’：可以安全通行的方格；
3) ‘#’：有怪物的方格；
4) ‘*’：仙药所在位置。
当在一行中读入的是两个零时，表示输入结束。

输出

对于每组测试数据，分别输出一行，该行包含李逍遥找到仙药需要穿过的最少的方格数目(计数包括初始位置的方块)。如果他不可能找到仙药, 则输出 -1。

样例输入

8 8
.@##...#
#....#.#
#.#.##..
..#.###.
#.#...#.
..###.#.
...#.*..
.#...###
6 5
.*.#.
.#...
..##.
.....
.#...
....@
9 6
.#..#. 
.#.*.# 
.####. 
..#... 
..#... 
..#... 
..#... 
#.@.## 
.#..#. 
0 0

样例输出

10
8
-1

二、分析

这道题要有广搜（广搜是什么，点击这里）来做，比如说之前我用深搜（深搜是什么，点击这里）来做就Time Limit Exceeded（超时）了，这是我的错误程序：

#include<cstdio>
#include<iostream>
using namespace std;
int ud[5]={0,1,-1,0,0},lr[5]={0,0,0,1,-1},ans=100000,tmp,sx,sy,tx,ty,m,n;//ud数组控制上下，lr数组控制作用
char a[25][25];

void dfs(int x,int y,int tmp){//深搜
	if(x==tx&&y==ty)
		if(tmp<ans)
			ans=tmp;
	if(x<1||y<1||x>m||y>n||a[x][y]=='#')//越界或撞墙退出
		return;
	for(int i=1;i<=4;i++){//往上下左右四个方向探索
		a[x][y]='#';//当前坐标堵住
		dfs(x+ud[i],y+lr[i],tmp+1);//探索
		a[x][y]='.';//当前坐标设为可以通过的路
	}
}

int main()
{
	while(1){
		cin>>m>>n;
		if(m==0&&n==0)
			return 0;
		for(int i=1;i<=m;i++)
			for(int j=1;j<=n;j++){
				cin>>a[i][j];
				if(a[i][j]=='@'){sx=i;sy=j;}//找到开始点
				if(a[i][j]=='*'){tx=i;ty=j;}//找到结束点
			}
		dfs(sx,sy,0);//深搜
		if(ans!=100000){//输出
			cout<<ans<<endl;
			ans=100000;
		}
		else
			cout<<"-1"<<endl;
	}
}

然而，我就测了一个很大的数据，程序竟然过了很久都没有输出结果，甚至还爆掉了，用freopen也是一样，下面就是我的超大数据：

20 20
@...................
....................
....................
....................
....................
....................
....................
....................
....................
....................
....................
....................
....................
....................
....................
....................
....................
....................
....................
...................*

所以，这道题经过我的探究，必须用广搜来做，对于不懂广搜的人，这里是一张广搜的概念图。

三、源代码

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
char a[201][201];
int b[40401],c[40401];
int e[40401];
int x[5]={0,1,0,-1,0};
int y[5]={0,0,1,0,-1};//<span style="font-family: Arial, Helvetica, sans-serif;">方向坐标移动，方便后面的for循环</span>
int tot=0;
void print(int d)
{
	tot++;
	if(e[d]!=0)
		print(e[d]);
	return;
}
int main()
{
	int R,C;
	while(1){
		memset(b,0,sizeof(b));
		memset(c,0,sizeof(c));
		memset(e,0,sizeof(e));//清零
		int head=0,tail=1,xz=0,yz=0,f=0;
		tot=0;
		scanf("%d%d",&R,&C);
		if(!R&&!C)
			return 0;
		for(int i=1;i<=R;i++)//广搜
			for(int r=1;r<=C;r++){
				cin>>a[i][r];
				if(a[i][r]=='@'){
					a[i][r]='#';
					b[tail]=i;
					c[tail]=r;
					continue;
				}
				if(a[i][r]=='*'){
					xz=i;
					yz=r;
					continue;
				}
			}
		do{
			head++;
			for(int i=1;i<=4;i++){
				if(b[head]+x[i]>=1&&b[head]+x[i]<=R&&c[head]+y[i]>=1&&c[head]+y[i]<=C&&a[b[head]+x[i]][c[head]+y[i]]!='#'){
					tail++;
				e[tail]=head;
				b[tail]=b[head]+x[i];
				c[tail]=c[head]+y[i];
				if(b[tail]==xz&&c[tail]==yz){
					f=1;
					print(tail);
					printf("%d\n",tot-1);
					break;
				}
				a[b[tail]][c[tail]]='#';
				}
			}
			if(f==1)
				break;
			else
				continue;
		}while(head!=tail);
		if(f!=1)
			printf("-1\n");
		else
			continue;
	}
}