ZH奶酪:【数据结构与算法】搜索之BFS
1、目标
通过本文,希望可以达到以下目标,当遇到任意问题时,可以:
1、很快建立状态空间;
2、提出一个合理算法;
3、简单估计时空性能;
2、搜索分类
2.1、盲目搜索
按照预定的控制策略进行搜索,在搜索过程中获得的中间信息不用来改进控制策略;
常见算法:
1、广度优先搜索(Breadth First Search);
2、深度优先搜索(Depth First Search);
3、纯随机搜索、重复式搜索、迭代加深搜索、迭代加宽搜索、柱形搜索;
2.2、启发式搜索
在搜索中加入了与问题有关的启发性信息,用以指导搜索朝着最有希望的方向发展,加速问题的求解过程并达到最优解。
常见算法:
1、A*算法;
2、IDA*算法;
3、状态空间
状态
问题在某一时刻进展状况的数学描述;
状态转移
问题从一种状态转移到另一种或几种状态的操作;
状态空间
一个“图”,图结点对应于状态,点之间的边对应于状态转移;
搜索
寻找一种可行的操作序列,从起始状态经过一系列状态转移,达到目标状态;
4、过河问题
某人带一条狗、一只鸡、一筐米过河,但小船除了需要人划动外,最多可以带一个物品,而人不在场时,狗要吃鸡,鸡要吃米。问此人应该如何过河?
4.1、问题分析
状态:建立四元组(人,狗,鸡,米)。0表示在左岸,1表示在右岸。
起始状态:(0,0,0,0)
终止状态:(1,1,1,1)
状态转移规则(操作,用1-x来表示过一次河,无论是到左岸还是到右岸都可以):
(a,b,c,d)
→(1-a,1-b,c,d)(当a=b)【人带狗过河】
→(1-a,b,1-c,d)(当a=c)【人带鸡过河】
→(1-a,b,c,1-d)(当a=d)【人带米过河】
→(1-a,b,c,d) 【人自己过河】
约束条件:
(a,b,c,d)中,a≠b时b≠c,a≠c时c≠d
4.2、直观的搜索过程
首先从初始状态开始:
然后再对两个可能状态进行进一步搜索:
搜索在“图”中进行,但是图不需要事先建立(“隐式图”);
搜索过程就是对图的遍历过程,可以得到一课“搜索树”;
搜索树的结点个数、分支数、深度,决定着搜索的效率;
存储和算法
普通状态可以用4个整数表示;
压缩状态可以用4个bit表示;
用BFS或DFS;
5、BFS基本结构
BFS搜索过程没有回溯,是一种牺牲空间换取时间的方法。时间复杂度:O(V+E)
1 定义一个队列; 2 起始点加入队列; 3 while(队列不空){ 4 取出队头结点; 5 若它是所求的目标状态,跳出循环; 6 否则,从它扩展出子节点,全都添加到队尾; 7 } 8 若循环中找到目标,输出结果; 9 否则输出无解;
6、BFS例题
6.1 HDOJ1242 Rescue
Description:
Angel was caught by the MOLIGPY! He was put in prison by Moligpy. The prison is described as a N * M (N, M <= 200) matrix. There are WALLs, ROADs, and GUARDs in the prison.
Angel's friends want to save
Angel. Their task is: approach Angel. We assume that "approach Angel" is
to get to the position where Angel stays. When there's a guard in the
grid, we must kill him (or her?) to move into the grid. We assume that
we moving up, down, right, left takes us 1 unit time, and killing a
guard takes 1 unit time, too. And we are strong enough to kill all the
guards.
You have to calculate the minimal time to approach Angel.
(We can move only UP, DOWN, LEFT and RIGHT, to the neighbor grid within
bound, of course.)
Input:
First line contains two integers stand for N and M.
Then N lines
follows, every line has M characters. "." stands for road, "a" stands
for Angel, and "r" stands for each of Angel's friend.
Process to the end of the file.
Output:
For each test case, your program should output a single integer, standing for the minimal time needed. If such a number does no exist, you should output a line containing "Poor ANGEL has to stay in the prison all his life."
Sample Input:
1 7 8 2 #.#####. 3 #.a#..r. 4 #..#x... 5 ..#..#.# 6 #...##.. 7 .#...... 8 ........
Sample Output:
1 13
Code 1:
1 //这是一个比较标准的bfs,没有经过任何优化,但是思路比较清晰,容易看懂。 2 #include <iostream> 3 #include <cstdio> 4 #include <queue> 5 using namespace std; 6 //node结构体 7 typedef struct 8 { 9 int x; 10 int y; 11 int len; 12 }node; 13 //全局变量定义 14 #define M 202 15 char Map[M][M];//地图 16 int mask[M][M];//访问标志 17 queue<node> q;//队列,只在bfs中用到 18 int bx,by,ex,ey,w,h;//起点、终点、宽、高 19 int step[4][2] = {//四个方向 20 //0-up 21 0,-1, 22 //1-right 23 1,0, 24 //2-down 25 0,1, 26 //3-left 27 -1,0 28 }; 29 30 void readMap(int m,int n);//读取地图 31 void bfs();//bfs 32 int tryXY(int x,int y);//尝试x、y点 33 34 void main() 35 { 36 while (scanf("%d %d",&h,&w)==2) 37 { 38 readMap(h,w); 39 bfs(); 40 cout<<mask[ex][ey]<<endl; 41 } 42 } 43 44 void readMap(int m,int n)//m-h,n-w 45 { 46 int i,j; 47 for (i=0;i<m;i++) 48 { 49 for (j=0;j<n;j++) 50 { 51 cin>>Map[i][j]; 52 mask[i][j] = -1;//标志为未访问 53 if (Map[i][j] == 'r') 54 {//friend为起点,且化为road 55 Map[i][j] = '.'; 56 bx = i; by = j; 57 } 58 if (Map[i][j] == 'a') 59 {//angel为终点,且化为road 60 Map[i][j] = '.'; 61 ex = i; ey = j; 62 } 63 } 64 } 65 } 66 67 void bfs() 68 { 69 node n,m;//m为队头,n为m衍生的测试点 70 int i,ret; 71 //起点 72 m.x = bx; 73 m.y = by; 74 m.len = 0;//len 75 mask[bx][by] = 0;//ask 76 q.push(m);//push 77 //处理 78 while (q.size()) 79 { 80 //get front 81 m = q.front(); 82 q.pop(); 83 //test node m 84 for (i=0 ; i<4 ; i++) 85 { 86 n.x = m.x+step[i][0]; n.y = m.y+step[i][1]; n.len = m.len+1; 87 ret = tryXY(n.x,n.y); 88 switch(ret) 89 { 90 case -1:// 91 break; 92 case 0: 93 case 1: 94 n.len += ret; 95 //如果未访问,保存花销;如果已经访问,保存花销最少。 96 if (mask[n.x][n.y] == -1 || n.len<mask[n.x][n.y]) 97 { 98 mask[n.x][n.y] = n.len;//已经访问且保存的是最小花销 99 q.push(n); 100 } 101 break; 102 } 103 } 104 105 } 106 } 107 int tryXY(int x,int y) 108 { 109 int ret; 110 //越界或遇到墙 111 if (!(x>=0 && x<w && y>=0 && y<h) || (Map[x][y] == '#')) 112 ret = -1; 113 //road or angel 114 if (Map[x][y] == '.') 115 ret = 0; 116 //guard 117 if (Map[x][y] == 'x') 118 ret = 1; 119 return ret; 120 }
然而这个代码还有一些不足~~
1、82行代码之后应该判断是否到了目标点,达到目标则return,此时就是花销最少的路径,而不需要全部遍历结束;
2、96行的代码中“<”判断是多余的,因为如果之前已经访问过了,那么现在的一定比之前的路径长;
Code 2:
1 #include<iostream> 2 #include<cstdio> 3 #include<queue> 4 using namespace std; 5 6 //定义结点结构体 7 typedef struct{ 8 int x; 9 int y; 10 int len; 11 }node; 12 //定义全局变量 13 #define M 202 14 char Map[M][M]; //地图 15 int visited[M][M]; //地图访问标识 16 queue<node> q; //队列,在bfs中用到 17 int bx,by,ex,ey,w,h;//起点、终点、地图宽、地图高 18 int step[4][2] = { 19 0,-1, //move up 20 1,0, //move right 21 0,1, //move down 22 -1,0 //move left 23 }; 24 25 void readMap(int m,int n); //读取地图 26 void bfs(); //bfs搜索路径 27 int costXY(int x,int y); //尝试x,y是否可行 28 29 void main(){ 30 while(scanf("%d %d",&w,&h) == 2){ 31 readMap(w,h); 32 bfs(); 33 if(visited[ex][ey]!=-1){ 34 cout<<visited[ex][ey]<<endl; 35 }else{ 36 cout<<"Poor ANGEL has to stay in the prison all his life."<<endl; 37 } 38 } 39 } 40 void readMap(int m,int n){ 41 int i,j; 42 for(i=0;i<m;i++){ 43 for(j=0;j<n;j++){ 44 cin>>Map[i][j]; 45 visited[i][j] = -1;//-1表示未访问 46 if(Map[i][j] == 'r'){//store the start point 47 Map[i][j] = '.'; 48 bx = i; 49 by = j; 50 } 51 if(Map[i][j] == 'a'){//store the end point 52 Map[i][j] = '.'; 53 ex = i; 54 ey = j; 55 } 56 } 57 } 58 } 59 60 void bfs(){ 61 node n,m;//m is queue head, n is a copied m for test; 62 int ret; 63 //save start point into the queue 64 m.x = bx; 65 m.y = by; 66 m.len = 0; 67 visited[bx][by] = 0; 68 q.push(m); 69 //bfs 70 while(q.size()){ 71 //get the head point 72 m = q.front(); 73 q.pop(); 74 //if arrive end point return 75 if(m.x == ex && m.y == ey){ 76 visited[ex][ey] = m.len; 77 return ; 78 } 79 //test node m 80 for(int i=0;i<4;i++){ 81 n.x = m.x + step[i][0]; 82 n.y = m.y + step[i][1]; 83 n.len = m.len + 1; 84 ret = costXY(n.x,n.y); 85 switch(ret){ 86 case -1: 87 break; 88 case 0: 89 case 1: 90 n.len += ret; 91 //if first time visit, save the cost 92 //if has been visited, save the minimum cost 93 if(visited[n.x][n.y] == -1){ 94 visited[n.x][n.y] = n.len; 95 q.push(n); 96 } 97 break; 98 } 99 } 100 } 101 } 102 103 int costXY(int x,int y){ 104 int ret; 105 //out of map or a wall 106 if(!(x>=0 && x<=w && y>=0 && y<=h) || Map[x][y] == '#'){ 107 ret = -1; 108 } 109 //road or angel 110 if(Map[x][y] == '.'){ 111 ret = 0; 112 } 113 //guard, need one more unit time 114 if(Map[x][y] == 'x'){ 115 ret = 1; 116 } 117 return ret; 118 }
6.2 HDOJ1372 Knight Moves
Description:
A friend of you is doing research on the Traveling Knight Problem (TKP) where you are to find the shortest closed tour of knight moves that visits each square of a given set of n squares on a chessboard exactly once. He thinks that the most difficult part of the problem is determining the smallest number of knight moves between two given squares and that, once you have accomplished this, finding the tour would be easy.
Of course you know that it is vice versa. So you offer him to write a program that solves the "difficult" part.
Your
job is to write a program that takes two squares a and b as input and
then determines the number of knight moves on a shortest route from a to
b.
Input:
The input file will contain one or more test cases. Each test case consists of one line containing two squares separated by one space. A square is a string consisting of a letter (a-h) representing the column and a digit (1-8) representing the row on the chessboard.
Output:
For each test case, print one line saying "To get from xx to yy takes n knight moves.".
Sample Input:
1 e2 e4 2 a1 b2 3 b2 c3 4 a1 h8 5 a1 h7 6 h8 a1 7 b1 c3 8 f6 f6
Sample Output:
1 To get from e2 to e4 takes 2 knight moves. 2 To get from a1 to b2 takes 4 knight moves. 3 To get from b2 to c3 takes 2 knight moves. 4 To get from a1 to h8 takes 6 knight moves. 5 To get from a1 to h7 takes 5 knight moves. 6 To get from h8 to a1 takes 6 knight moves. 7 To get from b1 to c3 takes 1 knight moves. 8 To get from f6 to f6 takes 0 knight moves.
Code:
1 #include<iostream> 2 #include<cstdio> 3 #include<queue> 4 using namespace std; 5 6 //定义node结构体 7 typedef struct{ 8 int x,y,len; 9 }node; 10 //八个方向 11 int step[8][2] = { 12 -2,-1, //left2-up1 13 -1,-2, //left1-up2 14 1,-2, //right1-up2 15 2,-1, //right2-up1 16 2,1, //right2-down1 17 1,2, //right2-down2 18 -1,2, //left1-down2 19 -2,1 //left2-down1 20 }; 21 #define M 10 22 int Map[M][M];//chessboard 23 int bx,by,ex,ey;//begin point,end point 24 char a[3],b[3]; 25 26 void bfs(int x,int y); 27 28 void main(){ 29 while(scanf("%s%s",a,b)!=EOF){ 30 bx = a[0]-'a'+1; //the range of index is from 1 to 8 31 by = a[1]-'0'; //the range of index is from 1 to 8 32 ex = b[0]-'a'+1; //the range of index is from 1 to 8 33 ey = b[1]-'0'; //the range of index is from 1 to 8 34 bfs(bx,by); 35 } 36 return 0; 37 } 38 39 void bfs(int x,int y){ 40 node t,p; 41 queue<node> q; 42 //BFS Step1:push first node t into queue q 43 t.x = x; 44 t.y = y; 45 t.len = 0; 46 q.push(t); 47 //while(!q.empty()) 48 while(q.size()){ 49 //BFS Step2:get the head of queue 50 t = q.front(); 51 q.pop(); 52 //BFS Step3:arrive the end point,return. first time arriving is the minimum. 53 if(t.x == ex && t.y == ey){ 54 printf("To get from %s to %s takes %d knight moves.\n",a,b,t.num); 55 return; 56 } 57 //BFS Step4:push the next steps into queue q 58 for(int i=0;i<8;i++){ 59 node next; 60 next.x = t.x + step[i][0]; 61 next.y = t.y + step[i][1]; 62 //no out of range and not been visited 63 if(!(next.x>0 && next.x<9 && next.y>0 && next.y<9) && !Map[next.x][next.y]){ 64 next.len = t.len+1; 65 q.push(next); 66 } 67 } 68 } 69 }
![字节跳动内推](https://img2020.cnblogs.com/blog/323808/202004/323808-20200415172728571-1992447742.jpg)
找我内推: 字节跳动各种岗位
作者:
ZH奶酪(张贺)
邮箱:
cheesezh@qq.com
出处:
http://www.cnblogs.com/CheeseZH/
*
本文版权归作者和博客园共有,欢迎转载,但未经作者同意必须保留此段声明,且在文章页面明显位置给出原文连接,否则保留追究法律责任的权利。