djkstra最短路径算法

在解决有多种标尺的最短路径问题中,即当具有相同的最短路径要让总的花费最小等,有两种解决问题的模板代码,一种是直接用Dijkstra算法:只需要增加一个数组来存放新增的边权或点权或最短路径条数,然后在算法中修改优化d[v]的那个步骤即可,其他部分不需要改动;一种是使用Dijkstra+dfs:先在Dijkstra算法中记录所有的最短路径(只考虑距离),然后从这些最短路径中选出一条第二标尺最优的路径。
相关问题链接:travel plan
直接贴代码:
Dijkstra算法:

#include <cstdio>
#include <algorithm>
using namespace std;

const int MAXV = 501;
const int INF = 1000000000;
int n, m, s, D;
int g[MAXV][MAXV];
bool vis[MAXV] = {false};
int d[MAXV];
int cost[MAXV][MAXV] = {INF};
int c[MAXV];
int pre[MAXV];

void Dijkstra(int s)
{
    fill(d, d+MAXV, INF);
    d[s] = 0;
    fill(c, c+MAXV, INF);
    c[s] = 0;
    for(int i = 0; i < n; i++)
        pre[i] = i;
    for(int i = 0; i < n; i++){
        int u = -1, MIN = INF;
        for(int j = 0; j < n; j++){
            if(vis[j] == false && d[j] < MIN){
                u = j;
                MIN = d[j];
            }
        }
        if(u == -1) return;
        vis[u] = true;
        for(int v = 0; v < n; v++){
            if(vis[v] == false && g[u][v] != INF){
                if(g[u][v] + d[u] < d[v]){
                    d[v] = g[u][v] + d[u];
                    c[v] = c[u] + cost[u][v];
                    pre[v] = u;
                }
                else if(g[u][v] + d[u] == d[v]){
                    if(c[v] > c[u] + cost[u][v]){
                        c[v] = c[u] + cost[u][v];
                        pre[v] = u;
                    }
                }
            }
        }
    }
}

void dfs(int s, int v)
{
    if(s == v){
        printf("%d ", s);
        return;
    }
    dfs(s, pre[v]);
    printf("%d ", v);
}

int main()
{
    scanf("%d%d%d%d", &n,&m,&s,&D);
    for(int i = 0; i < n; i++)
        for(int j = 0; j < n; j++)
            g[i][j] = INF;
    for(int i = 0; i < m; i++){
        int pa1,pa2,dist,co;
        scanf("%d%d%d%d",&pa1,&pa2,&dist,&co);
        g[pa1][pa2] = dist;
        g[pa2][pa1] = dist;
        cost[pa1][pa2] = co;
        cost[pa2][pa1] = co;
    }
    Dijkstra(s);
    dfs(s, D);
    printf("%d %d\n", d[D], c[D]);
    return 0;
}

Dijkstra算法+深搜:

#include <cstdio>
#include <algorithm>
#include <vector>
using namespace std;

const int MAXV = 501;
const int INF = 1000000000;
int n, m, s, D;
int g[MAXV][MAXV];
bool vis[MAXV] = {false};
int d[MAXV];
int cost[MAXV][MAXV] = {INF};
int optValue = INF;
vector<int> pre[MAXV];
vector<int> tmpPath, path;

void Dijkstra(int s)
{
    fill(d, d+MAXV, INF);
    d[s] = 0;
    for(int i = 0; i < n; i++){
        int u = -1, MIN = INF;
        for(int j = 0; j < n; j++){
            if(vis[j] == false && d[j] < MIN){
                MIN = d[j];
                u = j;
            }
        }
        if(u == -1) return;
        vis[u] = true;
        for(int v = 0; v < n; v++){
            if(vis[v] == false && g[u][v] != INF){
                if(d[u]+g[u][v] < d[v]){
                    d[v] = d[u] + g[u][v];
                    pre[v].clear();
                    pre[v].push_back(u);
                }
                else if(d[u]+g[u][v]==d[v]){
                    pre[v].push_back(u);
                }
            }
        }
    }
}

void dfs(int v){
    if(v == s){
        tmpPath.push_back(s);
        int value = 0;
        //calculate the cost
        for(int i = tmpPath.size()-1; i > 0; i--){
            int id = tmpPath[i], nextId = tmpPath[i-1];
            value += cost[id][nextId];
        }
        if(value < optValue){
            optValue = value;
            path = tmpPath;
        }
        tmpPath.pop_back();
        return;
    }
    tmpPath.push_back(v);
    for(int i = 0; i < (int)pre[v].size(); i++){
        dfs(pre[v][i]);
    }
    tmpPath.pop_back();
}

int main()
{
    scanf("%d%d%d%d", &n,&m,&s,&D);
    for(int i = 0; i < n; i++)
        for(int j = 0; j < n; j++)
            g[i][j] = INF;
    for(int i = 0; i < m; i++){
        int pa1,pa2,dist,co;
        scanf("%d%d%d%d",&pa1,&pa2,&dist,&co);
        g[pa1][pa2] = dist;
        g[pa2][pa1] = dist;
        cost[pa1][pa2] = co;
        cost[pa2][pa1] = co;
    }
    Dijkstra(s);
    dfs(D);
    for(int i = path.size()-1; i >= 0; i--){
        printf("%d ", path[i]);
    }
    printf("%d %d\n", d[D], optValue);
    return 0;
}
  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 3
    评论
以下是C语言实现图的邻接矩阵存储和Djkstra算法的示例代码: ```c #include <stdio.h> #include <limits.h> // 用于定义INT_MAX #define MAX_VERTICES 100 // 最大顶点数 // 图的邻接矩阵存储结构 typedef struct { int n; // 顶点数 int weight[MAX_VERTICES][MAX_VERTICES]; // 邻接矩阵 } Graph; // Djkstra算法,求图G从顶点v出发到其他顶点的最短路径 void Djkstra(Graph G, int v, int *dist, int *prev) { int i, j, k, min; int visited[MAX_VERTICES]; // 标记顶点是否已经访问 // 初始化 for (i = 0; i < G.n; i++) { visited[i] = 0; dist[i] = G.weight[v][i]; if (dist[i] < INT_MAX) { prev[i] = v; } else { prev[i] = -1; } } visited[v] = 1; dist[v] = 0; // 主循环 for (i = 1; i < G.n; i++) { min = INT_MAX; // 找未访问顶点中距离最近的 for (j = 0; j < G.n; j++) { if (!visited[j] && dist[j] < min) { min = dist[j]; k = j; } } visited[k] = 1; // 更新距离和前驱顶点 for (j = 0; j < G.n; j++) { if (!visited[j] && G.weight[k][j] < INT_MAX && dist[k] + G.weight[k][j] < dist[j]) { dist[j] = dist[k] + G.weight[k][j]; prev[j] = k; } } } } int main() { Graph G; int i, j, v; int dist[MAX_VERTICES]; // 存放最短路径长度 int prev[MAX_VERTICES]; // 存放最短路径的前驱顶点 // 输入顶点数和邻接矩阵 printf("请输入顶点数:"); scanf("%d", &G.n); printf("请输入邻接矩阵:\n"); for (i = 0; i < G.n; i++) { for (j = 0; j < G.n; j++) { scanf("%d", &G.weight[i][j]); if (G.weight[i][j] == 0) { G.weight[i][j] = INT_MAX; } } } // 输入起始顶点 printf("请输入起始顶点:"); scanf("%d", &v); // 求最短路径 Djkstra(G, v, dist, prev); // 输出结果 for (i = 0; i < G.n; i++) { printf("从顶点%d到顶点%d的最短路径长度为%d,路径为:", v, i, dist[i]); j = i; while (j != v && prev[j] != -1) { printf("%d <- ", j); j = prev[j]; } if (j == v) { printf("%d\n", j); } else { printf("无\n"); } } return 0; } ``` 程序首先定义了一个Graph结构体,用于存储图的邻接矩阵;然后定义了Djkstra算法,求出从指定顶点出发到其他顶点的最短路径;最后在main函数中输入邻接矩阵和起始顶点,调用Djkstra函数求解最短路径,并输出结果。
评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值