求最短路——Dijstra算法（朴素版）

一、相关笔记

二、算法模板（含注释）

#include<iostream>
#include<cstring>
#include<algorithm>

using namespace std;

const int N = 510;

int g[N][N];  // 稠密图利用邻接矩阵来存储
int dist[N];  // 存储每个点到源点的最近距离
bool visited[N];  // 存储每个点是否已经确定其到源点的最小距离

int n, m; // 点数，边数

int Dijstra(){
    memset(dist, 0x3f3f3f3f, sizeof(dist)); // 初始化，将所有点到源点的距离初始化为无穷
    dist[1] = 0;  // 将源点到自己的距离改为0

    for(int i = 0; i < n; ++ i){  // 进行n次迭代，每次都可以求出一个点到源点的最短距离
        //在未确定到源点的距离的点中寻找当前距离源点最近的点， 将此距离确定为该点距离源点的最近距离
        int t = -1;  // 用来存放当前被访问的点
        for(int j = 1; j <= n; ++ j){
            if(!visited[j] &&(t == -1 || dist[j] < dist[t])){
                t = j;
            }
        }
        visited[t] = true;

        //用该点更新其他（未确定的）点的最近距离
        for(int j = 1; j <= n; ++ j){
        if(!visited[j]){
            dist[j] = min(dist[j], dist[t]+g[t][j]);
        }
      }
    }

    if(dist[n] == 0x3f3f3f3f) return -1;
    return dist[n];
}

int main(){
    // 建立图
    scanf("%d%d", &n, &m);
    memset(g, 0x3f3f3f3f, sizeof(g));

    while(m--){
       int x, y, z;
       scanf("%d%d%d", &x, &y, &z);
       g[x][y] = min(g[x][y], z);  // 注意点：如果有重边，则取权值较小者
    }

    // 求最短路
    int t = Dijstra();
    printf("%d", t);
}

 

 

相关模板题

https://www.acwing.com/activity/content/code/content/7371478/#:~:text=%23include%3Ciostream,d%22%2C%20t)%3B%0A%7D

三、相关题目

1.（PTA作业18-7-1）城市间紧急救援

题目链接

PTA | 程序设计类实验辅助教学平台 (pintia.cn)

（1）题目内容

7-1 城市间紧急救援

分数 25

全屏浏览题目

作者 陈越

单位 浙江大学

作为一个城市的应急救援队伍的负责人，你有一张特殊的全国地图。在地图上显示有多个分散的城市和一些连接城市的快速道路。每个城市的救援队数量和每一条连接两个城市的快速道路长度都标在地图上。当其他城市有紧急求助电话给你的时候，你的任务是带领你的救援队尽快赶往事发地，同时，一路上召集尽可能多的救援队。

输入格式：

输入第一行给出4个正整数N、M、S、D，其中N（2≤N≤500）是城市的个数，顺便假设城市的编号为0 ~ (N−1)；M是快速道路的条数；S是出发地的城市编号；D是目的地的城市编号。

第二行给出N个正整数，其中第i个数是第i个城市的救援队的数目，数字间以空格分隔。随后的M行中，每行给出一条快速道路的信息，分别是：城市1、城市2、快速道路的长度，中间用空格分开，数字均为整数且不超过500。输入保证救援可行且最优解唯一。

输出格式：

第一行输出最短路径的条数和能够召集的最多的救援队数量。第二行输出从S到D的路径中经过的城市编号。数字间以空格分隔，输出结尾不能有多余空格。

输入样例：

4 5 0 3
20 30 40 10
0 1 1
1 3 2
0 3 3
0 2 2
2 3 2

输出样例： 

2 60
0 1 3

代码长度限制

16 KB

时间限制

400 ms

内存限制

64 MB

（2）思路分析

（3）代码实现（含注释）

方法1：

#include<iostream>
#include<cstring>
#include<algorithm>
#include<vector>

using namespace std;

const int N = 510;

int g[N][N];
int dist[N];
int visited[N];
int tnum[N];

int n, m, s, d;
int cnt, maxTNum;

vector<int>ans,p;

void Dijstra(){
    memset(dist, 0x3f3f3f3f, sizeof(dist));
    dist[s] = 0;

    for(int i = 0; i < n; ++ i){

        int t = -1;
        for(int j = 0; j < n; ++ j){
             if(!visited[j] && (t == -1 || dist[j] < dist[t])){
                t = j;
             }
        }
        visited[t] = 1;

        for(int j = 0; j < n; ++ j){
            if(!visited[j]){
                dist[j] = min(dist[j], dist[t] + g[t][j]);
            }
        }
    }
}

void DFS(int s, int tsum){
    if(s == d){
        if(tsum > maxTNum){
            maxTNum = tsum;
            ans.clear();
            ans = p;
        }
      ++ cnt;
      return ;
    }

    for(int i = 0; i < n; ++ i){
        if(!visited[i] && dist[i] == dist[s] + g[s][i]){
            p.push_back(i);
            visited[i] = 1;
            DFS(i, tsum + tnum[i]);
            p.pop_back();
            visited[i] = 0;
        }
    }
}

int main(){
    scanf("%d%d%d%d", &n, &m, &s, &d);
    for(int i = 0; i < n; ++ i){
        scanf("%d", &tnum[i]);
    }
    memset(g, 0x3f3f3f3f, sizeof(g));

    while(m--){
        int x, y, z;
        scanf("%d%d%d", &x, &y, &z);
        g[x][y] = min(g[x][y], z);
        g[y][x] = g[x][y];
    }

    Dijstra();

    memset(visited, 0, sizeof(visited));
    p.push_back(s);
    visited[s] = 1;
    DFS(s, tnum[s]);

    printf("%d %d\n", cnt, maxTNum);
    for(int i = 0; i < ans.size(); ++ i){
        if(i == 0) printf("%d", ans[i]);
        else printf(" %d", ans[i]);
    }
    return 0;
}

方法2：

#include<iostream>
#include<cstring>
#include<algorithm>

using namespace std;

const int N = 510;

int g[N][N];
int dist[N];
int visited[N];

int tnum[N];
int tsum[N];
int path[N];
int cntpath[N];

int n, m, s, d;

void Dijstra(){
    //初始化
    memset(dist, 0x3f3f3f3f, sizeof(dist));//初始化距离数组
    //对源点进行特殊处理
    dist[s] = 0;
    tsum[s] = tnum[s];
    cntpath[s] = 1;
    path[s] = -1;
    
    // 求出n个点的答案数组（本题的第一个难点：需要维护多个答案数组）  包括：到源点的最短距离，该点的前一点（后面通过递归可得出完整路径），到该点的最短路数目，到该点的救援队总数量
    for(int i = 0; i < n; ++ i){
        int t = -1;
        for(int j = 0; j < n; ++ j){
            if(!visited[j] && (t == -1 || dist[j] < dist[t])){
                t = j;
            }
        }
        visited[t] = 1;

        for(int j = 0; j < n; ++ j){// 本题的第二个难点：可能有多个最短路径
            if(!visited[j] && dist[j] > dist[t] + g[t][j]){
                dist[j] = dist[t] + g[t][j];
                path[j] = t;
                tsum[j] = tnum[j] + tsum[t];
                cntpath[j] = cntpath[t];
            }else if(!visited[j] && dist[j] == dist[t] + g[t][j]){
                cntpath[j] += cntpath[t];
                if(tsum[j] < tnum[j] + tsum[t]){ //如果该条路径能得到更多救援队则替换为该路径
                    tsum[j] = tnum[j] + tsum[t];
                    path[j] = t;
                }
            }
        }
    }
}

void PrintPath(int d){  // 递归打印路径
    if(path[d] == -1) return;
    PrintPath(path[d]);
    printf("%d ", path[d]);
}

int main(){
    //建立图
    scanf("%d%d%d%d", &n, &m, &s, &d);
    for(int i = 0; i < n; ++ i){
        scanf("%d", &tnum[i]);
    }
    memset(g, 0x3f3f3f3f, sizeof(g));
    while(m--){
        int x, y, z;
        scanf("%d%d%d", &x, &y, &z);
        g[x][y] = min(g[x][y], z);
        g[y][x] = g[x][y];
    }
    // 求最短路，维护答案数组
    Dijstra();
    printf("%d %d\n", cntpath[d], tsum[d]);
    // 打印路径
    PrintPath(d);
    printf("%d", d);
    return 0;
}