几种最短路径算法

Dijisktra

下面图片是从这个链接https://blog.csdn.net/qq_37796444/article/details/80663810找到,比较清楚
在这里插入图片描述

不能处理负边,时间复杂度O(n^2),堆优化O(nlogE)

基本思路:首先它是求单源点到其他边的最短路。方法就是创建一个集合来存放一些点,一开始放源点进去。然后每次取到集合外哪个点的距离最短,将其加入集合。(这个时候要更新,即通过该点再到其他点会不会比现在到其他点的距离短,会就更新)
在一开始输入图的信息后,肯定至少有一条边已经是最短边了,如果说能通过其他点再到目标点距离更短,那源点到这个其他点的距离不就更短了,不就是最短边了?

测试时,邻接矩阵太大了,爆掉了。
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
const int unreach=2147483647;
ll node[10010][10010];//邻接矩阵存图
ll dist[10010];//单源点到其他点的距离
int exist[10010];//该点是否已列入集合中
int n,m,s,F,G,W;
/*
    给出一个有向图,请输出从某一点出发到所有点的最短路径长度
    N、M、S,分别表示点的个数、有向边的个数、出发点的编号
    Fi、Gi、Wi,分别表示第i条有向边的出发点、目标点和长度
*/
void dijkstra(int s)
{
    for(int i=1;i<=n;i++){
        dist[i]=node[s][i];
    }
    exist[s]=1;
    ll newnode,length;
    for(int k=1;k<=n-1;k++){
        for(int i=1;i<=n;i++){
            if(!exist[i]){
                newnode=i;
                length=dist[i];
                break;
            }
        }
        for(int i=1;i<=n;i++){
            if(!exist[i]&&dist[i]<length){
                newnode=i;
                length=dist[i];
            }
        }
        exist[newnode]=1;
        for(int i=1;i<=n;i++){
            dist[i]=min(dist[i],dist[newnode]+node[newnode][i]);
        }
    }
}
int main()
{
    cin>>n>>m>>s;
    //初始化
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            node[i][j]=unreach;
        }
        node[i][i]=0;
    }
    for(int i=1;i<=n;i++)
        exist[i]=0;
    //
    for(int i=0;i<m;i++){
        scanf("%d%d%d",&F,&G,&W);
        if(W<node[F][G])node[F][G]=W;
    }
    dijkstra(s);
    for(int i=1;i<=n;i++)
        printf("%lld ",dist[i]);
    return 0;
}
所以用邻接表,这里用vector模拟
#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
using namespace std;
typedef long long ll;
const int unreach=2147483647;
struct dot
{
    int pos;
    int length;
};
vector<dot> node[10010];
ll dist[10010];//单源点到其他点的最短距离
ll edge[10010];//点到其他点的距离
int exist[10010];//该点是否已列入集合中
int n,m,s,F,G,W;
/*
    给出一个有向图,请输出从某一点出发到所有点的最短路径长度
    N、M、S,分别表示点的个数、有向边的个数、出发点的编号
    Fi、Gi、Wi,分别表示第i条有向边的出发点、目标点和长度
*/
void dijkstra(int s)
{
    for(int i=1;i<=n;i++){
        dist[i]=unreach;
    }
    for(int j=0;j<node[s].size();j++){
        dot t=node[s][j];
        if(t.length<dist[t.pos]){//两个点之间可能不止一条边
            dist[t.pos]=t.length;
        }
    }
    dist[s]=0;

    exist[s]=1;
    ll newnode,length;
    for(int k=1;k<=n-1;k++){
        for(int i=1;i<=n;i++){
            if(!exist[i]){
                newnode=i;
                length=dist[i];
                break;
            }
        }
        for(int i=1;i<=n;i++){
            if(!exist[i]&&dist[i]<length){
                newnode=i;
                length=dist[i];
            }
        }
        exist[newnode]=1;
        for(int i=1;i<=n;i++)
            edge[i]=unreach;
        for(int i=0;i<node[newnode].size();i++){
            dot t=node[newnode][i];
            //两个点之间可能不止一条边
            if(t.length<edge[t.pos]){
                edge[t.pos]=t.length;
            }
        }
        for(int i=1;i<=n;i++){
            dist[i]=min(dist[i],dist[newnode]+edge[i]);
        }
    }
}
int main()
{
    cin>>n>>m>>s;
    //
    for(int i=0;i<m;i++){
        scanf("%d%d%d",&F,&G,&W);
        dot t;
        t.pos=G,t.length=W;
        node[F].push_back(t);
    }
    dijkstra(s);
    for(int i=1;i<=n;i++)
        printf("%lld ",dist[i]);
    return 0;
}
最后用堆优化
#include<iostream>
#include<cstdio>
#include<vector>
#include<queue>
using namespace std;
const int maxn = 100000 + 10;
struct Node {
    int pos;
    int dis;
    Node() {}
    Node(int p, int d) {
        pos = p;
        dis = d;
    }
};
// 由于优先队列默认为大顶堆,所以重载小于号要用 a.dis > b.dis,来起到小顶堆的作用
bool operator<(const Node &a, const Node &b) {
    return a.dis > b.dis;
}
int n,m,s,u,v,w;
int dis[maxn];
int vis[maxn];
vector<Node> G[maxn];
priority_queue<Node> que;

void dij(int s) {
    // 初始化距离
    for(int i=1;i<=n;i++){
        dis[i]=1e9+10;
    }
    dis[s] = 0;
    que.push(Node(s, 0));
    while(!que.empty()) {
        Node tmp = que.top();
        que.pop();
        if(vis[tmp.pos])continue;
        vis[tmp.pos]=1;
        // 更新节点 tmp 的邻接节点,若邻接节点被更新,则将节点和 dis 加入小顶堆
        int len = G[tmp.pos].size();
        for(int i = 0; i < len; ++i) {
            int pos = G[tmp.pos][i].pos;
            int d = G[tmp.pos][i].dis;
            if(dis[pos] > tmp.dis + d) {
                dis[pos] = tmp.dis + d;
                que.push(Node(pos, dis[pos]));
            }
        }
    }
}
int main()
{
    cin>>n>>m>>s;
    for(int i=0;i<m;i++){
        scanf("%d%d%d",&u,&v,&w);
        Node node;
        node.pos=v,node.dis=w;
        G[u].push_back(node);
    }
    dij(s);
    for(int i=1;i<=n;i++)
        printf("%d ",dis[i]);
    return 0;
}

Floyd
/*
    |Floyd算法|
    |任意点对最短路算法|
    |求图中任意两点的最短距离的算法|
*/

for (int i = 0; i < n; i++) {   //  初始化为0  
    for (int j = 0; j < n; j++)  
        scanf("%lf", &dis[i][j]);  
}  
for (int k = 0; k < n; k++) {  
    for (int i = 0; i < n; i++) {  
        for (int j = 0; j < n; j++) {  
            dis[i][j] = min(dis[i][j], dis[i][k] + dis[k][j]);  
        }  
    }
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值