深入解析Dijkstra's Algorithm —— 高效解决有向图中的单点出发最短路径问题

什么是Dijkstra算法?

Dijkstra算法是用来寻找最短路径最著名的算法之一。具体来说,Dijkstra算法主要用来寻找一个边的权值不为负的有向图中的任意一点到其他任意结点(在两点相互联通的情况下)之间的最小路径。如果利用Dijkstra算法找出从一点出发,到图中其他所有点的最短路径,事实上我们就构造出了一个最短路径树(shortest-path tree)。

Dijkstra最短路径算法因其简单高效的特征而被广泛应用于互联网的路由协议中,例如IS-IS(Intermediate System to Intermediate System)和OSPF(Open Shortest Path First)协议。

Dijkstra算法的历史

1956年,当Edsger W. Dijkstra在阿姆斯特丹数学中心(Mathematical Center in Amsterdam)工作室,为了说明一台名为ARMAC的新电脑的计算能力,发明了最小路径问题。Dijkstra期初的目的是提出一个能被非计算机背景的人所能理解的问题,同时这一问题能被计算机有效解决。因此他对荷兰64个城市简化过的交通图(之所以是64个城市,是为了使得6比特的内存能够完全编码这些城市)计算最短路径,同时实现了Dijkstra算法的最初版本。

算法

Dijkstra算法实际上是一个贪婪算法(Greedy algorithm)。因为该算法总是试图优先访问每一步循环中距离起始点最近的下一个结点。Dijkstra算法的过程如下图所示。

Dijkstra Animation

下面我们分别给出算法的步骤及其正确性的证明。

初始化

  1. 给定图中的一个结点s作为起始点。
  2. 给定一个数组dist[]存储图中所有结点到s的距离。将dist[s]初始化为0。对于图中的其他结点v,初始化dist[v]为无穷大。初始化为无穷大的意义在于我们假设其余所有结点在当前情况下尚未与s联通。随着算法的执行,dist[v]会保存图中从sv的最短路径的距离。
  3. 给定一个minimum Heap,记为Q。堆顶为当前情况下距离s最近的结点及相应的距离。将(s, 0)放入堆中。
  4. 给定一个Set,记为S,保存所有已经访问过的结点。Set初始为空。基于Dijkstra算法的性质,我们总是以最短的路径遍历每一个结点,因此对于任一结点,一旦我们已经访问过,就代表着我们已经得到了从s到达这一结点的最短路径。

计算最短路径

  1. Q不为空的情况下,取出堆顶的元素(v, [dist[v]) —— 也就是当前距离s最近的结点v,及其距离dist[v]
  2. 如果vS中,则代表我们已经访问过v的最短路径。那么跳过当前v,重复步骤1。
  3. 否则,将v放在S中。
  4. 对于每一个与v相邻的结点t
    • 如果dist[v] + weight(v, t) < dist[t],则更新dist[t] = dist[v] + weight(v, t)。同时将(t, dist[t])放进Q中。
    • 否则,不做任何处理。

当算法结束后,dist[]中保存图中每一个除s之外的结点到s的最短路径的权重值(或长度)。如果从sv不存在联通的路径,则dist[v] = ∞

证明算法正确性

假设对于每个已经访问过的结点vdist[v]存储从起始点sv的最短路径。

当算法初始化时,dist[]中只包含dist[s] = 0,其正确性显而易见。

对于其余n-1个结点,假设u已经被访问且v尚未被访问,同时uv之间存在一条边u -> v,其权重为weight(u,v),那么一定有dist[v] = dist[u] + weight(u, v)。否则的话,假设存在另一条更短的路径dist[t]满足dist[t] + weight(t, v),则根据上述算法,t一定先于u被访问,则与我们当前的假设产生了矛盾。该论断对于余下的所有结点都成立。

因此Dijkstra算法一定能给出从出发点到其余所有结点(在可以到达的情况下)的最短路径。

复杂度分析

设图中总计有E条边,N个结点。

时间复杂度:O(ElogE)。因为所使用的最小堆最大可达O(E)大小,同时我们从其中将每个元素取出来一次。

空间复杂度:O(N+E)。其中O(N)为存储dist所用空间。O(E)为存储图的邻接链表及最小堆所用空间。

Java实现

class DijkstraShortestPath {
    /**
     * Given a list of (source, target, weight) edge pairs, return the shortest distance from s to any
     * other nodes in the graph. Any unreachable node has a distance of Integer.MAX_VALUE. 
     * @param edges List of tuple representation of edges containing [source, target weight].
     * @param N The graph contains nodes from 1 to N. 
     * @param s Start node of the shortest path tree
     * @return Shortest path from s to other nodes in the graph.
     */
    public Map<Integer, Integer> dijkstraShortestPath(int[][] edges, int N, int s) {
        Map<Integer, List<int[]>> adjList = new HashMap<>();
        for (int[] edge : edges) {
            adjList.computeIfAbsent(edge[0], k -> new ArrayList<>()).add(new int[]{edge[1], edge[2]});
        }
        
        Map<Integer, Integer> dist = new HashMap<>();
        PriorityQueue<int[]> minHeap = new PriorityQueue<>((edge1, edge2) -> {
            return edge1[1] - edge2[1];
        });
        minHeap.offer(new int[]{s, 0});
        
        while (!minHeap.isEmpty()) {
            int[] curr = minHeap.poll();
            if (dist.containsKey(curr[0])) continue;
            dist.put(curr[0], curr[1]);
            if (adjList.containsKey(curr[0])) {
                for (int[] edge : adjList.get(curr[0])) {
                    minHeap.offer(new int[]{edge[0], edge[1] + curr[1]});
                }
            }
        }
        
        for (int i = 1; i <= N; i++) {
            if (!dist.containsKey(i)) {
                dist.put(i, Integer.MAX_VALUE);
            }
        }
        
        return dist;
    }
}
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 204,445评论 6 478
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,889评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 151,047评论 0 337
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,760评论 1 276
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,745评论 5 367
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,638评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,011评论 3 398
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,669评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,923评论 1 299
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,655评论 2 321
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,740评论 1 330
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,406评论 4 320
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,995评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,961评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,197评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 45,023评论 2 350
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,483评论 2 342

推荐阅读更多精彩内容

  • 参见贪心算法——最短路径Dijkstra算法参见动态规划 目录 0.最短路径问题0.1 最短路径问题描述 0.1....
    王侦阅读 4,766评论 1 9
  • dijkstra算法介绍:即迪杰斯特拉算法,是从一个顶点到其余各顶点的最短路径算法,解决的是有向图中最短路径问题。...
    俊爷拒做学渣阅读 24,608评论 6 17
  • 1 概述 最短路径是图中的常见问题,最典型的应用是:当我们用百度地图或高德地图引导我们去某个地方时,它通常会给出一...
    CodingTech阅读 1,490评论 4 9
  • 目录 1.流网络及最大流问题1.1 流网络1.2 最大流问题1.3 最大流问题的三种解法对比 2.Ford-Ful...
    王侦阅读 4,574评论 0 3
  • https://zh.visualgo.net/graphds 浅谈图形结构https://zh.visualgo...
    狼之独步阅读 4,113评论 0 0