• 算法-11| 最短路径| Dijkstra算法


    有向有权图

    图的两种搜索算法,深度优先搜索和广度优先搜索。这两种算法主要是针对无权图的搜索算法。

    针对有权图,也就是图中的每条边都有一个权重,该如何计算两点之间的最短路径(经过的边的权重和最小)呢?常用的最短路径算法(Shortest Path Algorithm)。

    地图软件的最优路线是如何计算出来?底层依赖什么算法?这里的最优,比如最短路线、最少用时路线、最少红绿灯路线等等。

    最短路线:解决软件开发中的实际问题最重要的一点就是建模,也就是将复杂的场景抽象成具体的数据结构。

    把地图抽象成图,把每个岔路口看作一个顶点,岔路口与岔路口之间的路看作一条边,路的长度就是边的权重。 如果路是单行道,我们就在两个顶点之间画一条有向边;如果路是双行道,我们就在两个顶点之 间画两条方向不同的边。这样,整个地图就被抽象成一个有向有

    权图。于是,要求解的问题就转化为,在一个有向有权图中,求两个顶点间的最短路径。

    import java.util.LinkedList;
    // 有向有权图的邻接表表示
    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 implements Comparable<Vertex> {
            public int id; // 顶点编号 ID
            public int dist; // 从起始顶点到这个顶点的距离
            public Vertex(int id, int dist) {
                this.id = id;
                this.dist = dist;
            }
            @Override
            public int compareTo(Vertex o) { // 按照 dist 从小到大排序
                if (o.dist > this.dist) return -1;
                else return 1;
            }
        }
    }

    Dijkstra 算法

    解决这个问题,有一个非常经典的算法,最短路径算法,单源最短路径算 法(一个顶点到一个顶点),最出名的莫过于 Dijkstra 算法。Dijkstra 算法的工作原理如下:

        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 < v; ++i) { // 初始化 dist 为无穷大
                vertexes[i] = new Vertex(i, Integer.MAX_VALUE);
            }
            PriorityQueue<Vertex> queue = new PriorityQueue<>(); // 小顶堆
            boolean[] inQueue = new boolean[this.v]; // 标记是否进入过队列
            queue.add(vertexes[s]); // 先把起始顶点放到队列中
            vertexes[s].dist = 0;
            inqueue[s] = true;
            while (!queue.isEmpty()) {
                Vertex minVertex= queue.poll(); // 取 dist 最小的顶点
                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
                    // 找到一条到 nextVertex 更短的路径
                    if (minVertex.dist + e.w < nextVertex.dist) {
                        nextVertex.dist = minVertex.dist + e.w; // 更新 dist
                        predecessor[nextVertex.id] = minVertex.id; // 更新前驱顶点
                        if (inqueue[nextVertex.id] == false) { // 如果没有在队列中
                            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 数组的作用是为了还原最短路径,它记录每个顶点的前驱顶点。最后,我们通过 递归的方式,将这个路径打印出来。这个跟图的搜索中打印路径方法一样。

    inqueue 数组是为了避免将一个顶点多次添加到优先级队列中。更新了某个顶点的 dist 值 之后,如果这个顶点已经在优先级队列中了,就不要再将它重复添加进去了。

    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 算

    法, 这样就可以避免遍历整个大图,也就大大提高了执行效率。

    如果两点距离比较远,从北京海淀区某个地点,到上海黄浦区某个地点,那 上面的这种处理方法,显然就不工作了,毕竟覆盖北京和上海的区块并不小。

    对于这样两点之间距离较远的路线规划,可以把北京海淀区或者北京看作一个顶点,把上海黄浦区或者上海看作一个顶点,先规划大的出行路线。比如如何从北京到上海,必须要经过某几个顶点,或者某几条干道,再细化每个阶段的小路线。 这样,最短路径问题就解决了

    最少时间和最少红绿灯

    最短路径问题时,每条边的权重是路的长度。在计算最少时间的时候,算法还是不变,只需要把边的权重,从路的长度变成经过这段路所需要的时间。不过,这个时间会根据拥堵情况时刻变化。如何计算车通过一段路的时间呢?

    每经过一条边,就要经过一个红绿灯。关于最少红绿灯的出行方案,实际上,我们只需要把每条边的权值改为 1 即可,算法还是不变,继续使用Dijkstra算法。不过,边的权值为1,也就相当于无权图了,还可以使用广度优先搜索算法。

    广度优先搜索算法计算出来的两点之间的路径,就是两点的最短路径。

    真实的地图软件的路径规划,要比这个复杂很多;比起Dijkstra 算法,地图软件用的更多的是类似 A* 的启发式搜索算法,不过也是在 Dijkstra 算法上的优化罢了。

    总结

    一种非常重要的图算法,Dijkstra 最短路径算法。实际上,最短路径算法还有很多,比如 Bellford 算法、Floyd 算法等等。

     Dkijstra 算法不是贪心算法,它是动态规划算法,求得的解全局最优解。

    广度优先算法遍历无权图中一点s到另 一点t的路径,就是它的最短路径。

    在计算最短时间的出行路线中,如何获得通过某条路的时间呢?这个题目很有意思;

    获取通过某条路的时间:通过某条路的时间与①路长度②路况(是否平坦等)③拥堵情况④红 绿灯个数等因素相关。获取这些因素后就可以建立一个回归模型(比如线性回归)来估算时间。

    其中①②④因素比较固定,容易获得。③是动态的,但也可以通过

    • a.与交通部门合作获得路段 拥堵情况;
    • b.联合其他导航软件获得在该路段的在线人数;
    • c.通过现在时间段正好在次路段的 其他用户的真实情况等方式估算。
  • 相关阅读:
    C++笔记
    mongodb文件损坏的恢复--无可恢复数据
    wiredtiger--初学数据恢复
    mongodb异常恢复
    守护进程,互斥锁, IPC ,Queue队列,生产消费着模型
    进程,多进程,进程与程序的区别,程序运行的三种状态,multiprocessing模块中的Process功能,和join函数,和其他属性,僵尸与孤儿进程
    socket之UDP协议,并发编程介绍,操作系统发展史
    半连接数,粘包问题,自定义报头
    socket模块
    网络编程介绍,C/S 架构,网络通讯协议,osi七层
  • 原文地址:https://www.cnblogs.com/shengyang17/p/15186190.html
Copyright © 2020-2023  润新知