【C++算法】Dijkstra求最短路(朴素算法 + 堆优化)

题目描述

样例

3 3
1 2 2
2 3 1
1 3 4
结果 3

---

算法1

(朴素算法)

适用于 n^2 ~ m的稠密图，二维数组g表示从i到j的距离，dist表示距离数组，设置dist[1] = 0,表示从第一个点开始，从第一个点更新相邻的边的距离，用 t 搜素出全场最短的路（定理：最短路的子路也是最短），并且标记t，因为最短，否则会一直选这个t，用最短的t更新t的邻边，进入下一个循环，再找到离1最短的点，再度更新…

时间复杂度 $O(n^2)$

C++ 代码

#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int N = 510, M = 1e5+10;
int dist[N], g[N][N];
int n,m;
bool st[N];
int dijkstra()
{
    memset(dist,0x3f,sizeof dist);
    dist[1] = 0;
    
    for(int i = 0; i < n; i++)
    {
        int t = -1;
        for(int j = 1; j <= n;j++)
        {
            if(!st[j] && (t == -1 || dist[t] > dist[j]))
                t = j;
        }
        
        for(int j = 1; j <= n;j++)
            dist[j] = min(dist[j], dist[t] + g[t][j]);
    
        st[t] = true;
    }
    return dist[n];
}
int main()
{
    memset(g,0x3f,sizeof g);
    cin>>n>>m;
    for(int i = 0; i < m; i++)
    {
        int a,b,c;  cin>>a>>b>>c;
        g[a][b] = min(g[a][b],c);
    }
    if(dijkstra() == 0x3f3f3f3f)
        cout<<-1<<endl;
        else cout<<dist[n]<<endl;
    
}

---

算法2

(堆优化)

适用于 n ~ m的稀疏图，利用C++STL的priority_queue，特性为底层为堆，优先队列的队头为最小。

就代替了朴素算法的t寻找全场最短距离，用邻接表存储所有的边。

时间复杂度 $O(m\ast log(n))$

C++ 代码

#include<iostream>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
typedef pair<int,int>PII;
const int N = 2e5;
int e[N],ne[N],h[N],w[N],idx;
int dist[N];
bool st[N];
int n,m;
void add(int a,int b,int c)
{
    e[idx] = b,ne[idx] = h[a],w[idx] = c, h[a] = idx++;
}

void dijkstra()
{
    memset(dist,0x3f,sizeof dist);
    dist[1] = 0;
    priority_queue<PII,vector<PII>,greater<PII>>q;
    q.push({0,1});
    
    while(q.size())
    {
        auto t = q.top(); q.pop();
        int ver = t.second;
        if(st[ver])
            continue;
        st[ver] = true;
        for(int i = h[ver]; ~i ; i = ne[i])
        {
            int j = e[i];
            if(dist[j] > dist[ver] + w[i])
            {
                dist[j] = dist[ver] + w[i];
                q.push({dist[j],j});
            }
        }
    }
}

int main()
{
    cin>>n>>m;
    
    memset(h,-1,sizeof h);
    for(int i = 0; i < m; i++)
    {
        int a,b,c;  cin>>a>>b>>c;
        add(a,b,c);
    }
    dijkstra();
    if(dist[n] == 0x3f3f3f3f) cout<<-1<<endl;
    else cout<<dist[n]<<endl;
}