【算法】单源最短路径算法——Dijkstra算法

  • Post author:
  • Post category:其他




一、简介与使用场景

迪杰斯特拉算法(Dijkstra)是由荷兰计算机科学家狄克斯特拉于1959 年提出的,因此又叫狄克斯特拉算法。这是从一个顶点到其余各顶点的最短路径算法,解决的是有权图中最短路径问题。迪杰斯特拉算法主要特点是从起始点开始,采用

贪心算法

的策略,

每次遍历到始点距离最近且未访问过的顶点的邻接节点,直到扩展到终点为止。


Dijkstra算法的使用场景:


一定要是单源最短路径且每条边的权重必须是正值

单源最短路径:希望找到从(一个)源节点到每个结点的最短路径;



二、算法思想

Dijkstra算法主要是基于贪心的策略,

始终保持当前迭代解为当前最优解

。意思就是在已知的条件下或是当前拥有的全部条件下保证最优解,若在此后的迭代中由于加入了新的条件使得

产生了更优解则替代此前的最优解

。通过不断的迭代不断保证每次迭代的结果都是当前最优解,当迭代到

最后一轮时得到的就会是全局最优解。



步骤:


比如说现在要求1到n节点的最短距离。

我们可以用一个bool数组

st

记录每个边

是否已经被确定为最短路。


再用一个

dist

数组记录每个节点到源点的距离。

1️⃣ 先把第一个节点的距离记为0(

dist[1] = 0

),其他的节点距离记为正无穷。

2️⃣ 循环n – 1次,每一次循环内部再遍历所有节点,找到(不在

st

数组内的)最小权值的节点。再把这个节点加入到

st

中。

3️⃣ 用找到的这个节点更新其他(与他相连)节点的

dist




第一次找到x节点的时候一定从源点到x节点的最短路径。


举个例子:


在这里插入图片描述


先把①的dist变成0,其他的dist为正无穷。


在这里插入图片描述


循环两次,第一次找到了①节点。然后利用①更新没有确定已经是最短路径的节点,也就是2和3。而此时①就已经确定了他自己的最短路径,标记

st

为true。

在这里插入图片描述


第二次循环找到了②节点,那么利用②来处理还没有确定的节点,只剩下③。

在这里插入图片描述


由此单源最短路径就已经找到。



三、朴素版Dijkstra


用一道例题来讲解:



题目链接


题目描述:

给定一个 n 个点 m 条边的有向图,图中可能存在重边和自环,所有边权均为正值。

请你求出 1 号点到 n 号点的最短距离,如果无法从 1 号点走到 n 号点,则输出 −1。


输入格式

第一行包含整数 n 和 m。

接下来 m 行每行包含三个整数 x,y,z,表示存在一条从点 x 到点 y 的有向边,边长为 z。


输出格式

输出一个整数,表示 1 号点到 n 号点的最短距离。

如果路径不存在,则输出 −1。


数据范围

1≤n≤500,1≤m≤105,

图中涉及边长均不超过10000。


输入样例:

3 3

1 2 2

2 3 1

1 3 4


输出样例:

3


思路分析:


因为这道题的给的是稠密图,所以用

邻接矩阵

来存储。

而存在重边所以初始化邻接矩阵的时候每次要判断是不是最小,取最小值即可。

而如果存在自环,只需要最后判断n节点有没有被处理即可。

#include <iostream>
#include <cstring>

using namespace std;

const int N = 510;

int n, m;

int g[N][N];// 邻接矩阵
int dist[N];// 距离
bool st[N];// 是否已经确定了最短路径

int dijkstra()
{
    memset(dist, 0x3f, sizeof dist);
    dist[1] = 0;// 初始化第一个节点
    for(int i = 0; i < n - 1; i++)// 循环n - 1次
    {
        int t = -1;// 找最小节点
        for(int j = 1; j <= n; j++)
        {
            if(!st[j] && (t == -1 || dist[j] < dist[t]))
            {
                t = j;
            }
        }
        if(t == n) break;// 找了n节点是最小说明n已经确定为最短路径了
        st[t] = true;
        // 更新其他节点
        for(int j = 1; j <= n; j++)
        {
            dist[j] = min(dist[j], dist[t] + g[t][j]);
        }
    }
    if(dist[n] == 0x3f3f3f3f) return -1;
    else return dist[n];
}

int main()
{
    memset(g, 0x3f, sizeof g);
    cin >> n >> m;
    int a, b, c;
    while(m--)
    {
        cin >> a >> b >> c;
        g[a][b] = min(g[a][b], c);
    }
    cout << dijkstra() << endl;
    return 0;
}



四、堆优化版Dijkstra

而如果题目给的是一个稀疏图,我们还用上面的算法(时间复杂度为O(N^2)),可能会超时。

所以有了一个优化版本的Dijkstra算法。

通过分析朴素版本的代码,可以看到每次都要暴力遍历所有的边找到最小的距离,所以我们

可以使用堆,每次堆顶都是最小的边。


来看一道例题:



题目链接


题目描述:

给定一个 n 个点 m 条边的有向图,图中可能存在重边和自环,所有边权均为非负值。

请你求出 1 号点到 n 号点的最短距离,如果无法从 1 号点走到 n 号点,则输出 −1。


输入格式

第一行包含整数 n 和 m。

接下来 m 行每行包含三个整数 x,y,z,表示存在一条从点 x 到点 y 的有向边,边长为 z。


输出格式

输出一个整数,表示 1 号点到 n 号点的最短距离。

如果路径不存在,则输出 −1。


数据范围

1≤n,m≤1.5×105,图中涉及边长均不小于 0,且不超过 10000。

数据保证:如果最短路存在,则最短路的长度不超过 109。


输入样例:

3 3

1 2 2

2 3 1

1 3 4


输出样例:

3


思路描述:


大致思路跟朴素版一样,而我们怎么用优先级队列来得到最小距离的节点呢?

我们可以用

pair<int, int>

来存储,

first

代表距离,

second

代表节点。

这里用数组模拟邻接表在上一章有讲解:


【数据结构与算法】图——邻接表与邻接矩阵

#include <iostream>
#include <cstring>
#include <queue>

using namespace std;

typedef pair<int, int> PII;

const int N = 150010;

int n, m;

int h[N], e[N], ne[N], w[N], idx;
bool st[N];
int dist[N];

void add(int a, int b, int c)
{
    e[idx] = b;
    w[idx] = c;
    ne[idx] = h[a];
    h[a] = idx++;
}

int dijkstra()
{
    memset(dist, 0x3f, sizeof dist);
    priority_queue<PII, vector<PII>, greater<PII>> pq;
    pq.push({0, 1});
    dist[1] = 0;
    while(pq.size())
    {
        auto t = pq.top();
        pq.pop();
        int node = t.second, val = t.first;
        if(st[node]) continue;
        st[node] = true;
        // 更新其他节点
        for(int i = h[node]; i != -1; i = ne[i])
        {
            int j = e[i];
            dist[j] = min(dist[j], w[i] + dist[node]);
            pq.push({dist[j], j});
        }
    }
    if(dist[n] == 0x3f3f3f3f) return -1;
    else return dist[n];
}

int main()
{
    memset(h, -1, sizeof h);
    cin >> n >> m;
    int a, b, c;
    while(m--)
    {
        cin >> a >> b >> c;
        add(a, b, c);
    }
    cout << dijkstra() << endl;
    return 0;
}



五、总结


Dijkstra算法用来解决的是单源最短路径且没有负权值的问题。


并不是所有的题目都用优化版本的,而是稠密图用朴素版Dijkstra,稀疏图用堆优化版Dijkstra。



版权声明:本文为qq_66314292原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。