图的最短路

最新推荐文章于 2024-09-27 11:17:47 发布

小柒5437

最新推荐文章于 2024-09-27 11:17:47 发布

阅读量138

点赞数

文章标签：图论算法数据结构

本文链接：https://blog.csdn.net/m0_73079611/article/details/128508842

版权

这篇博客介绍了图的最短路问题，讲解了Floyd算法和Dijkstra算法的原理及优化。Floyd通过三重循环实现任意两点间最短路径，Dijkstra则采用堆优化解决单源最短路径问题，确保在大规模数据下依然高效。文中还给出了两种算法的代码实现，并提供了样例分析。

摘要由CSDN通过智能技术生成

图的最短路

知识点回顾

Floyd算法

三重大循环暴力解决,再加上就地滚动,使空间优化

for(int i=1;i<=n;i++)
{
  for(int j=1;j<=n;j++)
  {
    for(int z=1;z<=n;z++)
    {
      f[j][z]=min(f[j][z],f[j][i]+f[i][z]);
    }
  }
}

Dijkstra算法

原始代码

memset(dis,0x3f,sizeof(dis));
dis[1]=0;
for(int i=1;i<=n;i++)
{
  int v;
  for(int j=1;j<=n;j++)//找到当前未找到的点的dis的值最小的点
  {
    if(v==0 || dis[j]<dis[v] && !vis[j]) v=j;
  }
  vis[k]=1;
  for(int i=1;i<=G[k].size();i++)
  {
    if(vis[G[k][i-1]]==1)continue;
    if(dis[G[k][i-1]]>dis[k]+f[k][G[k][i-1]])
    {
      dis[G[k][i-1]]=dis[k]+f[k][G[k][i-1]];
    }
  }
}

堆优化和链式前向星

在例题中有,就不赘述了

例题

【模板】floyd

题目背景

模板题，无背景

题目描述

给出n个点，m条边的无向图，求每个点到其他点的距离之和%998244354的值

输入格式

第一行两个数n，m含义如上
从第二行开始，共m行，每行三个数x，y，l，代表从x到y点的长度为l

输出格式

n行，每行一个数，第i行代表点i到其他点的距离之和

样例 #1

样例输入 #1

2 1
1 2 4

样例输出 #1

4
4

样例 #2

样例输入 #2

样例输出 #2

提示

模板题，保证图联通$ n\le500 $ , $ m\le10000 $ , $ 1 \le x,y \le n $ , $ l \le 10^9 $

思路

三个暴力循环解决

代码实现

#include<bits/stdc++.h>
using namespace std;
const int maxn=550;
int f[maxn][maxn]={};
int ret[maxn]={};
int main()
{
    int n,m;
    cin>>n>>m;
    memset(ret,0,sizeof(ret));
    memset(f,0x3f,sizeof(f));
    for(int i=1;i<=n;i++) f[i][i]=0;
    for(int i=1;i<=m;i++)
    {
        int x,y;
        int t;
        cin>>x>>y;
        if(f[x][y]!=0x3f3f3f3f) 
        {
            cin>>t;
            f[x][y]=min(f[x][y],t);
        }
        else cin>>f[x][y];
        
        f[y][x]=f[x][y];
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
            for(int z=1;z<=n;z++)
            {
                f[j][z]=min(f[j][z],f[j][i]+f[i][z]);
            }
        }
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
            ret[i]=(ret[i]+f[i][j])%998244354;
        }
    }
    for(int i=1;i<=n;i++)
    {
        cout<<ret[i]<<endl;
    }
    return 0;
}

【模板】单源最短路径（标准版）

题目背景

2018 年 7 月 19 日，某位同学在 NOI Day 1 T1 归程一题里非常熟练地使用了一个广为人知的算法求最短路。

然后呢？

$\rightarrow 60$ ；

$\text{Ag} \rightarrow \text{Cu}$ ；

最终，他因此没能与理想的大学达成契约。

小 F 衷心祝愿大家不再重蹈覆辙。

题目描述

给定一个 $n$ 个点， $m$ 条有向边的带非负权图，请你计算从 $s$ 出发，到每个点的距离。

数据保证你能从 $s$ 出发到任意点。

输入格式

第一行为三个正整数 $n, m, s$ 。
第二行起 $m$ 行，每行三个非负整数 $u_i, v_i, w_i$ ，表示从 $u_i$ 到 $v_i$ 有一条权值为 $w_i$ 的有向边。

输出格式

输出一行 $n$ 个空格分隔的非负整数，表示 $s$ 到每个点的距离。

样例 #1

样例输入 #1

样例输出 #1

0 2 4 3

提示

样例解释请参考数据随机的模板题。

$\leq n \leq 10^5$ ；

$\leq m \leq 2\times 10^5$ ；

$s = 1$ ；

$\leq u_i, v_i\leq n$ ；

$\leq w_i \leq 10 ^ 9$ ,

$\leq \sum w_i \leq 10 ^ 9$ 。

本题数据可能会持续更新，但不会重测，望周知。

2018.09.04 数据更新 from @zzq

思路

考虑到数据在 $ 10^5 $ 左右,单单地用原始版本的Dijkstra $ O(n^2) $ ,会爆炸.

所以要用到堆优化的Dijkstra,复杂度在 $ mlog(n) $ 刚好,主要是用于优化查找还没找到的dis的值最小的点.

同时可以在输入的时候用链式前向星优化.

代码实现

#include<bits/stdc++.h>
using namespace std;
int head[150005]={},dis[150005]={},cnt;
bool vis[150005]={};
int n,m,s;
int u,v,w;
struct node
{
    int dis;
    int pos;
    bool operator<(const node &n) const//用于后续priority_queue的排序,注意,只能重载<,这里和sort是相反的.
    {
        return dis>n.dis; //设置为升序
    }
};
struct edge //设一个边的类,用于链式前向星的输入
{
    int to ,dis ,next;
}e[500005];
void add(int u ,int v,int d)//加边
{
    cnt++;
    e[cnt].dis=d;
    e[cnt].to=v;
    e[cnt].next=head[u];
    head[u]=cnt;
}

priority_queue<node> q;
void dijkstra()
{
    dis[s]=0;
    q.push((node){0,s});
    while(!q.empty())
    {
        node tmp =q.top();
        q.pop();
        int x=tmp.pos,d=tmp.dis;
        if(vis[x]) continue;
        vis[x]=1;
        for(int i=head[x];i;i=e[i].next)//对x的出边进行遍历
        {
            int y=e[i].to;
            if(dis[y]>dis[x]+e[i].dis)
            {
                dis[y]=dis[x]+e[i].dis;
                if(!vis[y])//把未遍历过的点都弄到队列里面去,这里一定包含我们需要的点,又由于dis小的点会自动排到最前面,所以我们需要的答案也一定排到了最前面.
                {
                    q.push((node){dis[y],y});
                }
            }
        }
    }
}

int main()
{
    cin>>n>>m>>s;
    for(int i=1;i<=n;i++) dis[i]=0x7fffffff;
    for(int i=1;i<=m;i++)
    {
        cin>>u>>v>>w;
        add(u,v,w);
    }
    dijkstra();
    for(int i=1;i<=n;i++)
    {
        cout<<dis[i];
        if(i==n)cout<<endl;
        else cout<<" ";
    }
    return 0;
}