想认真的写篇题解。html
题目来自:https://www.luogu.org/problemnew/show/P1126ios
机器人移动学会(RMI
)如今正尝试用机器人搬运物品。机器人的形状是一个直径$1.6米的球。在试验阶段,机器人被用于在一个储藏室中搬运货物。储藏室是一个N×M的网格,有些格子为不可移动的障碍。机器人的中心老是在格点上,固然,机器人必须在最短的时间内把物品搬运到指定的地方。机器人接受的指令有:向前移动1步(Creep
);向前移动2步(Walk
);向前移动3步(Run
);向左转(Left
);向右转(Right
)。每一个指令所须要的时间为1秒。请你计算一下机器人完成任务所需的最少时间。c++
第一行为两个正整数N,M(N,M≤50),下面N行是储藏室的构造,0表示无障碍,1表示有障碍,数字之间用一个空格隔开。接着一行有4个整数和1个大写字母,分别为起始点和目标点左上角网格的行与列,起始时的面对方向(东E,南S,西W,北N),数与数,数与字母之间均用一个空格隔开。终点的面向方向是任意的。算法
一个整数,表示机器人完成任务所需的最少时间。若是没法到达,输出-1。数组
之前在刷刘汝佳的紫书《算法竞赛入门经典》时作过这道题,如今又再次遇到,幸不汝(辱)命,一次就过了。数据结构
咱们都知道,BFS 具备从起点到目标节点(或状态)路径最短的特性,可是使用 BSF 这一特性时须要注意,只有当全部的边权重相同(通常为1)时,它才具备此性质,边权不等时不具备。BFS 每次从一个节点只经历一次转移,当求单纯的求距离时能够认为每次转移的边权为一,转移次数最少的路径必定是距离最短的。咱们能够用两张图来直观的表现这个特性:测试
在图上:this
初始节点为 0,每次从一个节点向四周节点扩散,访问全部距离为一(相邻,通过一次状态转换)的节点。spa
第一次扩散访问了全部蓝色节点,并无找到目标节点,继续扩散。第二次扩散访问了两个粉色节点,虚线节点并无被访问。由于在扩散到左边的粉色节点时,咱们已经找到了目标节点。那么不去搜索第三个粉色节点(虚线)节点不会丢解吗?在图中咱们也确实能看到,虚线节点在经历一次扩散后,到达紫色节点,紫色节点再扩散一次,也可到达目标节点。不过,这个选择不管是从上述图片仍是从个人描述文字来看,他的距离都不短。那么为何会这样呢?我说一下我我的的理解。看下面一幅图:调试
在树上:
在树上的 BFS 被称为层次遍历。他的工做原理就如其名,每次访问一层的节点,同一层的节点有一个特色就是,他们的层数是相同的,也即到根节点的距离。当扩散到第三层时,(从左向右)第三个粉色节点做为目标节点被发现,此时咱们能够对比三种状况:
由此,能够得出若是只有一个解,那么第一次被发现的目标节点 即第三个粉色节点必定是距离根节点最近的。
在图上也能够类比层的概念,得出相同的结论。
BFS 的性质讨论完,再来具体考虑这道题。这道题应该能明显看出来,是在图上寻找最短路的问题。不过首先须要对数据进行一些分析处理,才能更好的应用 BFS。此题惟一有些麻烦的就是,机器人具备半径,将每一个格子单独处理不方便。由给出的图例能够发现,机器人必定会占据四个格子的空间,而一旦一个格子障碍物出现,那么这个障碍物格子所在的四个四方格上的中心格点机器人都不能走(机器人只走格点)。
因此在读完输入键一个图后,能够再建立一个简化版的图,mini map。它将原图中每四个格点当作一个点,而机器人只走这四个格子的中心,一旦一个四方格中有一个障碍,那么这个四方格认为不可走。这样一来机器人移动、机器人半径、单个方格障碍物问题就简化成了最多见的状况:在一个图上的点避开障碍物到达另外一个点的最短路。(这中简化和机器人在原图上的移动是等价的,能够模拟一下)。
如下是这部分处理代码:
for(int i = 0; i < n; ++i) { for(int j = 0; j < m; ++j) { cin >> graph[i][j]; } } int p = 0,q = 0; for(int i = 0; i < n-1; ++i) { q = 0; for(int j = 0; j < m-1; ++j) { if(graph[i][j] || graph[i][j+1] || graph[i+1][j] || graph[i+1][j+1]) mini[p][q] = 1;//标记为障碍,不然为 0 表示可达 q++; } p++; }
建出了图,如今考虑状态的转移。此题中,每一个节点能够有三种转移状况,向左转,向右转,移动。
每种状况耗时为 1(能够认为是距离,权重)。关于转向问题,能够定义一个方向数组,按顺时针顺序给出北东南西,而后不管朝向那个方位,向左转就是数组索引减一,向右转就是加一,检索方向数组得到新的方位。完成一次转向后,就完成了一次状态转移,将新的节点入队列,这个新的节点在下一次被取出考虑进行转移状态时,他的当前方向就是移动的方向。
宏和数据结构:
#define MAX 50 #define DIR 4 #define N 0 #define E 1 #define S 2 #define W 3 struct Node { int x,y; int step = 0; int d; Node (int x,int y,int step,int dir) { this->x = x,this->y = y,this->step = step,this->d = dir; } // Node () {} // int pre = 0; // int loc = 0; }; int graph[MAX + 1][MAX + 1]; int mini[MAX + 1][MAX + 1]; //mini graph int vis[MAX + 1][MAX + 1][DIR]; // N E S W int dx[] = {-1,0, 1, 0}; int dy[] = {0, 1, 0,-1}; int n,m;
BFS 代码:
Node bfs(int sx,int sy,int tx,int ty,int dir) { queue<Node> q; Node start = Node(sx,sy,0,dir); q.push(start); // path[pi++] = start; vis[sx][sy][dir] = 1; while(!q.empty()) { Node u = q.front(); q.pop(); if(u.x == tx && u.y == ty) return u; //尝试向左转 if(!vis[u.x][u.y][(u.d-1+DIR)%DIR]) { vis[u.x][u.y][(u.d-1+DIR)%DIR] = 1; Node nu = Node(u.x,u.y,u.step+1,(u.d-1+DIR)%DIR); // nu.pre = u.loc; //记录路径信息 // nu.loc = pi; // path[pi++] = nu; q.push(nu); } //尝试向右转 if(!vis[u.x][u.y][(u.d+1)%DIR]) { vis[u.x][u.y][(u.d+1)%DIR] = 1; Node nu = Node(u.x,u.y,u.step+1,(u.d+1)%DIR); // nu.pre = u.loc; // nu.loc = pi; // path[pi++] = nu; q.push(nu); } //尝试移动 creep,walk,run for(int i = 1; i <= 3; ++i) { Node nu = Node(u.x,u.y,u.step+1,u.d); //判断是是否有障碍物。只要有一个,就没必要移动了。 bool isok = true; for(int j = 0; j < i; ++j) { if(mini[nu.x+dx[u.d]*j][nu.y+dy[u.d]*j]) { isok = false; break; } } if(isok == false) break; nu.x += dx[u.d] * i; nu.y += dy[u.d] * i; if(inrange(nu) && !vis[nu.x][nu.y][nu.d]) { vis[nu.x][nu.y][nu.d] = 1; // nu.pre = u.loc; // nu.loc = pi; // path[pi++] = nu; q.push(nu); } } } return Node(-1,-1,-1,-1); }
其中注释代码是用来记录最短路的路径信息。在这里是我用来测试程序的。
在调试的时候遇到一个 bug 花了我几个小时,机器人移动方式有 creep,walk,run,分别移动 1 、2 、3 步,在移动步数大于 1 时,不能只判断目标节点是不是障碍物,而要判断每一步是否有障碍物。
完整程序:
#include <iostream> #include <algorithm> #include <cstring> #include <stdlib.h> #include <memory.h> #include <queue> using namespace std; #define MAX 50 #define DIR 4 #define N 0 #define E 1 #define S 2 #define W 3 struct Node { int x,y; int step = 0; int d; Node (int x,int y,int step,int dir) { this->x = x,this->y = y,this->step = step,this->d = dir; } // Node () {} // int pre = 0; // int loc = 0; }; int graph[MAX + 1][MAX + 1]; int mini[MAX + 1][MAX + 1]; //mini graph int vis[MAX + 1][MAX + 1][DIR]; // N E S W int dx[] = {-1,0, 1, 0}; int dy[] = {0, 1, 0,-1}; int n,m; // Node path[MAX * MAX + 1]; // int pi = 0; bool inrange(Node nd) { return nd.x >= 0 && nd.x <= n-2 && nd.y >= 0 && nd.y <= m-2; } Node bfs(int sx,int sy,int tx,int ty,int dir) { queue<Node> q; Node start = Node(sx,sy,0,dir); q.push(start); // path[pi++] = start; vis[sx][sy][dir] = 1; while(!q.empty()) { Node u = q.front(); q.pop(); if(u.x == tx && u.y == ty) return u; //尝试向左转 if(!vis[u.x][u.y][(u.d-1+DIR)%DIR]) { vis[u.x][u.y][(u.d-1+DIR)%DIR] = 1; Node nu = Node(u.x,u.y,u.step+1,(u.d-1+DIR)%DIR); // nu.pre = u.loc; //记录路径信息 // nu.loc = pi; // path[pi++] = nu; q.push(nu); } //尝试向右转 if(!vis[u.x][u.y][(u.d+1)%DIR]) { vis[u.x][u.y][(u.d+1)%DIR] = 1; Node nu = Node(u.x,u.y,u.step+1,(u.d+1)%DIR); // nu.pre = u.loc; // nu.loc = pi; // path[pi++] = nu; q.push(nu); } //尝试移动 creep,walk,run for(int i = 1; i <= 3; ++i) { Node nu = Node(u.x,u.y,u.step+1,u.d); ////判断是是否有障碍物。只要有一个,就没必要移动了。 bool isok = true; for(int j = 1; j <= i; ++j) { if(mini[nu.x+dx[u.d]*j][nu.y+dy[u.d]*j]) { isok = false; break; } } if(isok == false) break; nu.x += dx[u.d] * i; nu.y += dy[u.d] * i; if(inrange(nu) && !vis[nu.x][nu.y][nu.d]) { vis[nu.x][nu.y][nu.d] = 1; // nu.pre = u.loc; // nu.loc = pi; // path[pi++] = nu; q.push(nu); } } } return Node(-1,-1,-1,-1); } int main(int argc, char const *argv[]) { freopen("/home/skipper/Documents/code/刷题/洛谷 OJ/重启/in.txt","r",stdin); int sx,sy,tx,ty; char dir; cin >> n >> m; for(int i = 0; i < n; ++i) { for(int j = 0; j < m; ++j) { cin >> graph[i][j]; } } int p = 0,q = 0; for(int i = 0; i < n-1; ++i) { q = 0; for(int j = 0; j < m-1; ++j) { if(graph[i][j] || graph[i][j+1] || graph[i+1][j] || graph[i+1][j+1]) mini[p][q] = 1; q++; } p++; } int d; cin >> sx >> sy >> tx >> ty >> dir; switch(dir) { case 'S': d = S;break; case 'N': d = N;break; case 'E': d = E;break; case 'W': d = W;break; }; //test 查看 mini map // for(int i = 0; i < p; ++i) { // for(int j = 0; j < q; ++j) { // cout << mini[i][j] << " "; // } // cout << endl; // } Node res = bfs(sx-1,sy-1,tx-1,ty-1,d); //输出路径 // int pp = res.loc; // do { // cout << path[pp].x << "," << path[pp].y << endl; // pp = path[pp].pre; // }while(pp); cout << res.step << endl; return 0; }
若有错误,欢迎指正评论。
<br><br> 做者:<a href='https://www.cnblogs.com/backwords/'>Skipper</a> 出处:https://www.cnblogs.com/backwords/p/10542486.html 本博客中未标明转载的文章归做者 <a href='https://www.cnblogs.com/backwords/'>Skipper</a> 和博客园共有,欢迎转载,但未经做者赞成必须保留此段声明,且在文章页面明显位置给出原文链接,不然保留追究法律责任的权利。