极客时间已完结课程限时免费阅读

44 | 最短路径:地图软件是如何计算出最优出行路径的?

44 | 最短路径:地图软件是如何计算出最优出行路径的?-极客时间

44 | 最短路径:地图软件是如何计算出最优出行路径的?

讲述:冯永吉

时长14:34大小13.30M

基础篇的时候,我们学习了图的两种搜索算法,深度优先搜索和广度优先搜索。这两种算法主要是针对无权图的搜索算法。针对有权图,也就是图中的每条边都有一个权重,我们该如何计算两点之间的最短路径(经过的边的权重和最小)呢?今天,我就从地图软件的路线规划问题讲起,带你看看常用的最短路径算法(Shortest Path Algorithm)。
像 Google 地图、百度地图、高德地图这样的地图软件,我想你应该经常使用吧?如果想从家开车到公司,你只需要输入起始、结束地址,地图就会给你规划一条最优出行路线。这里的最优,有很多种定义,比如最短路线、最少用时路线、最少红绿灯路线等等。作为一名软件开发工程师,你是否思考过,地图软件的最优路线是如何计算出来的吗?底层依赖了什么算法呢?

算法解析

我们刚提到的最优问题包含三个:最短路线、最少用时和最少红绿灯。我们先解决最简单的,最短路线。
解决软件开发中的实际问题,最重要的一点就是建模,也就是将复杂的场景抽象成具体的数据结构。针对这个问题,我们该如何抽象成数据结构呢?
我们之前也提到过,图这种数据结构的表达能力很强,显然,把地图抽象成图最合适不过了。我们把每个岔路口看作一个顶点,岔路口与岔路口之间的路看作一条边,路的长度就是边的权重。如果路是单行道,我们就在两个顶点之间画一条有向边;如果路是双行道,我们就在两个顶点之间画两条方向不同的边。这样,整个地图就被抽象成一个有向有权图。
具体的代码实现,我放在下面了。于是,我们要求解的问题就转化为,在一个有向有权图中,求两个顶点间的最短路径。
public class Graph { // 有向有权图的邻接表表示
private LinkedList<Edge> adj[]; // 邻接表
private int v; // 顶点个数
public Graph(int v) {
this.v = v;
this.adj = new LinkedList[v];
for (int i = 0; i < v; ++i) {
this.adj[i] = new LinkedList<>();
}
}
public void addEdge(int s, int t, int w) { // 添加一条边
this.adj[s].add(new Edge(s, t, w));
}
private class Edge {
public int sid; // 边的起始顶点编号
public int tid; // 边的终止顶点编号
public int w; // 权重
public Edge(int sid, int tid, int w) {
this.sid = sid;
this.tid = tid;
this.w = w;
}
}
// 下面这个类是为了dijkstra实现用的
private class Vertex {
public int id; // 顶点编号ID
public int dist; // 从起始顶点到这个顶点的距离
public Vertex(int id, int dist) {
this.id = id;
this.dist = dist;
}
}
}
想要解决这个问题,有一个非常经典的算法,最短路径算法,更加准确地说,是单源最短路径算法(一个顶点到一个顶点)。提到最短路径算法,最出名的莫过于 Dijkstra 算法了。所以,我们现在来看,Dijkstra 算法是怎么工作的。
这个算法的原理稍微有点儿复杂,单纯的文字描述,不是很好懂。所以,我还是结合代码来讲解。
// 因为Java提供的优先级队列,没有暴露更新数据的接口,所以我们需要重新实现一个
private class PriorityQueue { // 根据vertex.dist构建小顶堆
private Vertex[] nodes;
private int count;
public PriorityQueue(int v) {
this.nodes = new Vertex[v+1];
this.count = v;
}
public Vertex poll() { // TODO: 留给读者实现... }
public void add(Vertex vertex) { // TODO: 留给读者实现...}
// 更新结点的值,并且从下往上堆化,重新符合堆的定义。时间复杂度O(logn)。
public void update(Vertex vertex) { // TODO: 留给读者实现...}
public boolean isEmpty() { // TODO: 留给读者实现...}
}
public void dijkstra(int s, int t) { // 从顶点s到顶点t的最短路径
int[] predecessor = new int[this.v]; // 用来还原最短路径
Vertex[] vertexes = new Vertex[this.v];
for (int i = 0; i < this.v; ++i) {
vertexes[i] = new Vertex(i, Integer.MAX_VALUE);
}
PriorityQueue queue = new PriorityQueue(this.v);// 小顶堆
boolean[] inqueue = new boolean[this.v]; // 标记是否进入过队列
vertexes[s].dist = 0;
queue.add(vertexes[s]);
inqueue[s] = true;
while (!queue.isEmpty()) {
Vertex minVertex= queue.poll(); // 取堆顶元素并删除
if (minVertex.id == t) break; // 最短路径产生了
for (int i = 0; i < adj[minVertex.id].size(); ++i) {
Edge e = adj[minVertex.id].get(i); // 取出一条minVetex相连的边
Vertex nextVertex = vertexes[e.tid]; // minVertex-->nextVertex
if (minVertex.dist + e.w < nextVertex.dist) { // 更新next的dist
nextVertex.dist = minVertex.dist + e.w;
predecessor[nextVertex.id] = minVertex.id;
if (inqueue[nextVertex.id] == true) {
queue.update(nextVertex); // 更新队列中的dist值
} else {
queue.add(nextVertex);
inqueue[nextVertex.id] = true;
}
}
}
}
// 输出最短路径
System.out.print(s);
print(s, t, predecessor);
}
private void print(int s, int t, int[] predecessor) {
if (s == t) return;
print(s, predecessor[t], predecessor);
System.out.print("->" + t);
}
我们用 vertexes 数组,记录从起始顶点到每个顶点的距离(dist)。起初,我们把所有顶点的 dist 都初始化为无穷大(也就是代码中的 Integer.MAX_VALUE)。我们把起始顶点的 dist 值初始化为 0,然后将其放到优先级队列中。
我们从优先级队列中取出 dist 最小的顶点 minVertex,然后考察这个顶点可达的所有顶点(代码中的 nextVertex)。如果 minVertex 的 dist 值加上 minVertex 与 nextVertex 之间边的权重 w 小于 nextVertex 当前的 dist 值,也就是说,存在另一条更短的路径,它经过 minVertex 到达 nextVertex。那我们就把 nextVertex 的 dist 更新为 minVertex 的 dist 值加上 w。然后,我们把 nextVertex 加入到优先级队列中。重复这个过程,直到找到终止顶点 t 或者队列为空。
以上就是 Dijkstra 算法的核心逻辑。除此之外,代码中还有两个额外的变量,predecessor 数组和 inqueue 数组。
predecessor 数组的作用是为了还原最短路径,它记录每个顶点的前驱顶点。最后,我们通过递归的方式,将这个路径打印出来。打印路径的 print 递归代码我就不详细讲了,这个跟我们在图的搜索中讲的打印路径方法一样。如果不理解的话,你可以回过头去看下那一节。
inqueue 数组是为了避免将一个顶点多次添加到优先级队列中。我们更新了某个顶点的 dist 值之后,如果这个顶点已经在优先级队列中了,就不要再将它重复添加进去了。
看完了代码和文字解释,你可能还是有点懵,那我就举个例子,再给你解释一下。
理解了 Dijkstra 的原理和代码实现,我们来看下,Dijkstra 算法的时间复杂度是多少?
在刚刚的代码实现中,最复杂就是 while 循环嵌套 for 循环那部分代码了。while 循环最多会执行 V 次(V 表示顶点的个数),而内部的 for 循环的执行次数不确定,跟每个顶点的相邻边的个数有关,我们分别记作 E0,E1,E2,……,E(V-1)。如果我们把这 V 个顶点的边都加起来,最大也不会超过图中所有边的个数 E(E 表示边的个数)。
for 循环内部的代码涉及从优先级队列取数据、往优先级队列中添加数据、更新优先级队列中的数据,这样三个主要的操作。我们知道,优先级队列是用堆来实现的,堆中的这几个操作,时间复杂度都是 O(logV)(堆中的元素个数不会超过顶点的个数 V)。
所以,综合这两部分,再利用乘法原则,整个代码的时间复杂度就是 O(E*logV)。
弄懂了 Dijkstra 算法,我们再来回答之前的问题,如何计算最优出行路线?
从理论上讲,用 Dijkstra 算法可以计算出两点之间的最短路径。但是,你有没有想过,对于一个超级大地图来说,岔路口、道路都非常多,对应到图这种数据结构上来说,就有非常多的顶点和边。如果为了计算两点之间的最短路径,在一个超级大图上动用 Dijkstra 算法,遍历所有的顶点和边,显然会非常耗时。那我们有没有什么优化的方法呢?
做工程不像做理论,一定要给出个最优解。理论上算法再好,如果执行效率太低,也无法应用到实际的工程中。对于软件开发工程师来说,我们经常要根据问题的实际背景,对解决方案权衡取舍。类似出行路线这种工程上的问题,我们没有必要非得求出个绝对最优解。很多时候,为了兼顾执行效率,我们只需要计算出一个可行的次优解就可以了
有了这个原则,你能想出刚刚那个问题的优化方案吗?
虽然地图很大,但是两点之间的最短路径或者说较好的出行路径,并不会很“发散”,只会出现在两点之间和两点附近的区块内。所以我们可以在整个大地图上,划出一个小的区块,这个小区块恰好可以覆盖住两个点,但又不会很大。我们只需要在这个小区块内部运行 Dijkstra 算法,这样就可以避免遍历整个大图,也就大大提高了执行效率。
不过你可能会说了,如果两点距离比较远,从北京海淀区某个地点,到上海黄浦区某个地点,那上面的这种处理方法,显然就不工作了,毕竟覆盖北京和上海的区块并不小。
我给你点提示,你可以现在打开地图 App,缩小放大一下地图,看下地图上的路线有什么变化,然后再思考,这个问题该怎么解决。
对于这样两点之间距离较远的路线规划,我们可以把北京海淀区或者北京看作一个顶点,把上海黄浦区或者上海看作一个顶点,先规划大的出行路线。比如,如何从北京到上海,必须要经过某几个顶点,或者某几条干道,然后再细化每个阶段的小路线。
这样,最短路径问题就解决了。我们再来看另外两个问题,最少时间和最少红绿灯。
前面讲最短路径的时候,每条边的权重是路的长度。在计算最少时间的时候,算法还是不变,我们只需要把边的权重,从路的长度变成经过这段路所需要的时间。不过,这个时间会根据拥堵情况时刻变化。如何计算车通过一段路的时间呢?这是一个蛮有意思的问题,你可以自己思考下。
每经过一条边,就要经过一个红绿灯。关于最少红绿灯的出行方案,实际上,我们只需要把每条边的权值改为 1 即可,算法还是不变,可以继续使用前面讲的 Dijkstra 算法。不过,边的权值为 1,也就相当于无权图了,我们还可以使用之前讲过的广度优先搜索算法。因为我们前面讲过,广度优先搜索算法计算出来的两点之间的路径,就是两点的最短路径。
不过,这里给出的所有方案都非常粗糙,只是为了给你展示,如何结合实际的场景,灵活地应用算法,让算法为我们所用,真实的地图软件的路径规划,要比这个复杂很多。而且,比起 Dijkstra 算法,地图软件用的更多的是类似 A* 的启发式搜索算法,不过也是在 Dijkstra 算法上的优化罢了,我们后面会讲到,这里暂且不展开。

总结引申

今天,我们学习了一种非常重要的图算法,Dijkstra 最短路径算法。实际上,最短路径算法还有很多,比如 Bellford 算法、Floyd 算法等等。如果感兴趣,你可以自己去研究。
关于 Dijkstra 算法,我只讲了原理和代码实现。对于正确性,我没有去证明。之所以这么做,是因为证明过程会涉及比较复杂的数学推导。这个并不是我们的重点,你只要掌握这个算法的思路就可以了。
这些算法实现思路非常经典,掌握了这些思路,我们可以拿来指导、解决其他问题。比如 Dijkstra 这个算法的核心思想,就可以拿来解决下面这个看似完全不相关的问题。这个问题是我之前工作中遇到的真实的问题,为了在较短的篇幅里把问题介绍清楚,我对背景做了一些简化。
我们有一个翻译系统,只能针对单个词来做翻译。如果要翻译一整个句子,我们需要将句子拆成一个一个的单词,再丢给翻译系统。针对每个单词,翻译系统会返回一组可选的翻译列表,并且针对每个翻译打一个分,表示这个翻译的可信程度。
针对每个单词,我们从可选列表中,选择其中一个翻译,组合起来就是整个句子的翻译。每个单词的翻译的得分之和,就是整个句子的翻译得分。随意搭配单词的翻译,会得到一个句子的不同翻译。针对整个句子,我们希望计算出得分最高的前 k 个翻译结果,你会怎么编程来实现呢?
当然,最简单的办法还是借助回溯算法,穷举所有的排列组合情况,然后选出得分最高的前 k 个翻译结果。但是,这样做的时间复杂度会比较高,是 O(m^n),其中,m 表示平均每个单词的可选翻译个数,n 表示一个句子中包含多少个单词。这个解决方案,你可以当作回溯算法的练习题,自己编程实现一下,我就不多说了。
实际上,这个问题可以借助 Dijkstra 算法的核心思想,非常高效地解决。每个单词的可选翻译是按照分数从大到小排列的,所以 肯定是得分最高组合结果。我们把 及得分作为一个对象,放入到优先级队列中。
我们每次从优先级队列中取出一个得分最高的组合,并基于这个组合进行扩展。扩展的策略是每个单词的翻译分别替换成下一个单词的翻译。比如 扩展后,会得到三个组合,。我们把扩展之后的组合,加到优先级队列中。重复这个过程,直到获取到 k 个翻译组合或者队列为空。
我们来看,这种实现思路的时间复杂度是多少?
假设句子包含 n 个单词,每个单词平均有 m 个可选的翻译,我们求得分最高的前 k 个组合结果。每次一个组合出队列,就对应着一个组合结果,我们希望得到 k 个,那就对应着 k 次出队操作。每次有一个组合出队列,就有 n 个组合入队列。优先级队列中出队和入队操作的时间复杂度都是 O(logX),X 表示队列中的组合个数。所以,总的时间复杂度就是 O(k*n*logX)。那 X 到底是多少呢?
k 次出入队列,队列中的总数据不会超过 k*n,也就是说,出队、入队操作的时间复杂度是 O(log(k*n))。所以,总的时间复杂度就是 O(k*n*log(k*n)),比之前的指数级时间复杂度降低了很多。

课后思考

在计算最短时间的出行路线中,如何获得通过某条路的时间呢?这个题目很有意思,我之前面试的时候也被问到过,你可以思考看看。
今天讲的出行路线问题,我假设的是开车出行,那如果是公交出行呢?如果混合地铁、公交、步行,又该如何规划路线呢?
欢迎留言和我分享,也欢迎点击“请朋友读”,把今天的内容分享给你的好友,和他一起讨论、学习。
分享给需要的人,Ta购买本课程,你将得20
生成海报并分享

赞 37

提建议

上一篇
43 | 拓扑排序:如何确定代码源文件的编译依赖关系?
下一篇
45 | 位图:如何实现网页爬虫中的URL去重功能?
unpreview
 写留言

精选留言(97)

  • 五岳寻仙
    2019-01-07
    课后思考题,自己能想到的。 1.获取通过某条路的时间:通过某条路的时间与①路长度②路况(是否平坦等)③拥堵情况④红绿灯个数等因素相关。获取这些因素后就可以建立一个回归模型(比如线性回归)来估算时间。其中①②④因素比较固定,容易获得。③是动态的,但也可以通过a.与交通部门合作获得路段拥堵情况;b.联合其他导航软件获得在该路段的在线人数;c.通过现在时间段正好在次路段的其他用户的真实情况等方式估算。 2.混合公交、地铁和步行时:地铁时刻表是固定的,容易估算。公交虽然没那么准时,大致时间是可以估计的,步行时间受路拥堵状况小,基本与道路长度成正比,也容易估算。总之,感觉公交、地铁、步行,时间估算会比开车更容易,也更准确些。
    展开
    共 7 条评论
    161
  • 李东勇
    2019-01-13
    有兴趣的可以看下LeetCode 上这道题: https://leetcode.com/problems/network-delay-time/ 用到的就是Dijkstra 算法
    共 2 条评论
    34
  • 徐凯
    2019-01-07
    @五岳寻仙的答案太棒了 👏 我感觉每条道路应该还有限速,这个因素也要考察。
    共 1 条评论
    31
  • Liam
    2019-01-08
    有2个疑问: 1 Dijkstra就是贪心算法吧? 2 它的解可能不是最优解

    作者回复: Dijkstra实际上可以看做动态规划:)

    共 8 条评论
    27
  • book尾汁
    2020-01-03
    感觉这个算法最不容易让人理解的地方就是优先级队列里的某个顶点B pop出去以后,会不会以后会有某个还没有入队列的节点C 使经过节点C到B到源点的距离更近,实际这是不会的,因为优先级队列每次都是pop出当前离源点距离最近的点,假如节点C经过B使B到源点的距离更近,那么C点在优先级队列一定会比B先pop出去,然后更新B到源点的距离,想明白这一点,这个算法就很好理解了。
    共 2 条评论
    20
  • 林大涛
    2019-02-13
    用小顶堆,就是为了确保每个阶段,堆顶的节点都是目前阶段的最短路径的节点。
    19
  • Shaohong
    2020-01-14
    geeksforgeeks上面对Dijkstra算法介绍很好懂。https://www.geeksforgeeks.org/dijkstras-shortest-path-algorithm-greedy-algo-7 根据这个算法,老师没有必要用inqueue数组。可以一开始就把所有Vertex都加入到PriorityQueue中,之后就只做取顶操作和更新操作。
    13
  • mrlay
    2019-07-20
    这是一个动态规划的问题,我最开始也以为是个贪心算法,会有这样的想法是每次都会去选择dist最小的那个顶点,殊不知这个顶点是在算完每个邻接顶点后(有可能是间接相邻的)选择最小的那个。 vertex 临时数组的作用就是为了临时记录从起点到该顶点dist, 用来更新最小优先队列。
    共 3 条评论
    11
  • yongxiang
    2019-01-12
    王老师,我输入代码运行后,实际出队列的顺序跟图中的不一样,实际(15,0)出队列 在(13,3)出队列前面。我看了代码,应该是修改(25,1)为 (13,3)的时候,小顶堆不会自动更新顺序。需要对22行进行如下修改,更新已经在队列中,又改了dist的Vertex的优先级: if (inQueue[nextVertex.id] == false){ queue.add(nextVertex); inQueue[nextVertex.id] = true; } else { // 更新已经在队列中,又改了dist的Vertex的优先级 queue.remove(nextVertex); queue.add(nextVertex); }
    展开

    作者回复: 嗯嗯 我更新下代码

    9
  • 王策
    2020-03-13
    文章翻译系统的例子与Leetcode上这一题类似:https://leetcode-cn.com/problems/find-k-pairs-with-smallest-sums/
    共 1 条评论
    8
  • Geek_94adb8
    2020-02-06
    使用内部类,代码的可读性很差!
    共 3 条评论
    7
  • xuery
    2019-03-22
    构建优先队列的update函数时,时间复杂度应该是O(n),因为小顶堆查找的时间复杂度是O(n),虽然查找之后向上堆化的时间复杂度时O(logn)
    共 1 条评论
    6
  • 社会你强哥
    2020-05-07
    不能只有代码,解释代码就完事了,能够从原理的角度来描述dijkstra是怎么做的么?
    5
  • zapup
    2020-04-04
    动态累计,谁最小谁行动
    4
  • book尾汁
    2020-01-02
    第十三天: 地图软件是如何计算出最短路径的 我们把每个岔路口看作一个顶点,岔路口与岔路口之间的路看作一条边,路的长度就是边的权重。如果路是单行道,我们就在两个顶点之间画一条有向边;如果路是双行道,我们就在两个顶点之间画两条方向不同的边,这样就构成了一幅有向有权图。这个问题就转化为求一个有向有权图中,两个顶点最短路径的问题。 上述问题可以使用dijkstra算法 原理:将每个顶点到起始顶点的距离初始为无穷,然后从起始点开始,将其加入一个优先级队列中,从优先级队列中取出到源点距离最小的顶点,然后比较其周围顶点离源点的距离是否大于其到源点的距离+其到周围顶点的距离,如果大于的话,更新周围顶点到源点的距离为较小的值以及其前序节点,并将其加入优先级队列中(如果已经加入过,就不需加入了),再取出优先级队列中的距离最小值,循环往复,直到取出终止顶点t,或者优先级队列为空。此时倒序输出终止顶点t的前序节点,前序节点的前序节点。。。,直到前序节点的前序节点为s,此时路径即为最短路径。 计算最少红绿灯依然可以采用上述的思路,构造一个有向无权图。 对于地点跨度范围比较大,可以分阶段来计算,找出哪些点是必经的点,然后拆分阶段。具体到每个区域,可以找个合适的区域将源点与中点覆盖进去,来减少顶点的数量。
    展开
    5
  • 三年过后
    2019-11-27
    优先队列代码: class PriorityQueue{//构建小顶堆 Vertex[] nodes; private int count;//队列个数 public PriorityQueue(int v){ nodes = new Vertex[v+1];//小顶堆,数组从小标1开始,好计算 this.count = 0;//初始0个元素 } public Vertex poll(){ Vertex v = nodes[1];//返回堆顶原始 nodes[1] = nodes[count];//将最后一个元素添加到堆顶,自上而下堆化 --count; heapifyUpToDown(1);//堆顶从上而下堆化 return v; } public void add(Vertex vertex){ nodes[++count] = vertex; vertex.i = count; heapifyDownToUp(count);//从下往上堆化 } public void update(Vertex vertex){ //查找,并更新 nodes[vertex.i].dist = vertex.dist; heapifyDownToUp(vertex.i);//从下往上堆化 } public boolean isEmpty(){ return this.count == 0 ? true : false; } //自上而下堆化 private void heapifyUpToDown(int i) { while(i<=count){ int maxPos = i; if((i*2)<=count && nodes[maxPos].dist > nodes[i*2].dist) maxPos = 2*i; else if((i*2+1)<=count && nodes[maxPos].dist > nodes[i*2+1].dist) maxPos = 2*i+1; else if(maxPos == i)break; swap(i,maxPos);//交换 i = maxPos; } } //从下往上堆化 private void heapifyDownToUp(int i) { while (i / 2 > 0 && nodes[i].dist < nodes[i / 2].dist) { swap(i,i/2);//交换 i = i / 2; } } /** * 数据交换 * @param i * @param maxPos */ private void swap(int i, int maxPos) { nodes[i].i = maxPos;//下标交换记录 nodes[maxPos].i = i; Vertex tmp = nodes[i]; nodes[i] = nodes[maxPos]; nodes[maxPos] = tmp; } }
    展开
    共 1 条评论
    4
  • Monday
    2019-08-30
    Dijkstra算法的最后到例子-那张图,各个节点的(0,无), (10,0)等等代表什么意思
    共 1 条评论
    3
  • mαnajay
    2019-05-18
    最开始以为是贪心算法, 再看了一遍,才发现优先级队列的作用(拥有类似回溯的功能),按照已添加队列中最短距离进行小顶堆, 每次 poll 的过程中,都有可能将之前已经计算过的顶点再拿出来, 遍历邻接表,重复到目的顶点
    3
  • 2021-05-12
    当权值变为1时, Dijkstra算法退化为BFS, prioQueue 退化为 queue, inqueue 数组退化为 visited 数组, predecessor 数组退化为 prev 数组
    2
  • Geek_6c75ea
    2021-02-13
    为啥我觉得dijkstra算法,就是带了权重的广度优先搜索?
    共 1 条评论
    2