算法竞赛进阶指南0x61 最短路

news/2024/11/30 20:46:02/

        对于一张有向图,我们一般有邻接矩阵和邻接表两种存储方式。对于无向图,可以把无向边看作两条方向相反的有向边,从而采用与有向图一样的存储方式。 

$$ 邻接矩阵的空间复杂度为 O(n^2),因此我们一般不采用这种方式 $$

        我们用数组模拟链表:长度为n的表头数组head记录了从每个节点出发的第一条边在 ver 和 edge 数组中的存储位置,长度为 m 的边集数组 ver 和 edge 记录了每条边的终点和边权。长度为 m 的数组 next 模拟了链表指针,表示从相同节点出发的下一条边在 ver 和 edge数组中的存储位置。

$$ 邻接表的空间复杂度为 O(n + m) $$

void add(int x, int y, int z) {ver[++tot] = y, edge[tot] = z, Next[tot] = head[x], head[x] = tot;
}//访问从 x 出发的所有边
for (int i = head[x]; i; i = Next[i]) {int y = ver[i], z = edge[i];//找到了一条有向边(x, y),权值为 z
}

单源最短路径

        单源最短路径问题(Single Source Shortest Psth,SSSP问题) 是说,给定一张有向图 G = (V,E),V是点集,E是边集,|V| = n,|E| = m,节点以 [1,n] 之间的连续整数编号,(x,y,z) 描述一条从 x 出发,到达 y,长度为 z 的有向边。设 1 号点为起点,求长度为 n 的数组 dist,其中dist[i],表示从起点 1 到 节点 i 的最短路径的长度。

Dijkstra算法

        $$ Dijkstra算法的流程如下: $$ 

        $$ 1.初始化 dist[1] = 0,其余节点的 dist 值为正无穷大 $$

        $$ 2.找出一个未被标记的、dist[x] 最小的节点 x,然后标记节点 x $$

        $$ 3.扫描节点 x 的所有出边 (x,y,z),若 dist[y] > dist[x]  + z,则使用 dist[x] +  z 更新 dist[y] $$

        $$ 4.重复上述 2 ~ 3 两个步骤,直到所有节点都被标记 $$

        Dijkstra 算法基于贪心思想,它只适用于所有边的长度都是非负数的图。当边长都是非负数时,全局最小值不可能再被其他节点更新,故在第一步中选出的节点 x 必然满足:dist[x] 已经是起点到 x 的最短路径。我们不断选择全局最小值进行标记和扩展,最终可得到起点 1 到每个节点的最短路径的长度。

const int N = 1e3 + 10;
int a[N][N], d[N], n, m, s;
bool v[N];
void dijkstra(int s) {std::memset(d, 0x3f, sizeof d);d[s] = 0;for (int i = 1; i < n; i++) {//重复进行 n - 1次int x = 0;//找到未标记节点中 dist 最小的for (int j = 1; j <= n; j++) {if (!v[i] && (x == 0 || d[j] < d[x])) {x = j;}}v[x] = 1;//用全局最小值点 x 更新其他节点for (int y = 1; y <= n; y++) {d[y] = std::min(d[y], d[x] + a[x][y]);}}
}int main() {std::cin >> n >> m >> s;//构建邻接矩阵std::memset(a, 0x3f, sizeof a);for (int i = 1; i <= n; i++) {a[i][i] = 0;}for (int i = 1; i <= m; i++) {int u, v, w;std::cin >> u >> v >> w;a[u][v] = std::min(a[u][v], w);}//求单源最短路径dijkstra(s);for (int i = 1; i <= n; i++) {std::cout << d[i] << " \n"[i == n];}return 0;
}

        $$ 上面程序的时间复杂度为O(n^2)$$

        $$ 主要瓶颈在于第一步的寻找全局最小值的过程 $$

$$ 可以用二叉堆(C++ STL priority_queue) 对 dist 数组进行维护 $$

$$ 用 O(logn) 的时间获取最小值并从堆中删除 $$

$$ 用O(logn) 的时间执行一条边的扩展和更新 $$

$$ 最终可在 O((m + n)logn) 的时间内实现 Dijkstra 算法 $$

const int N = 1e5 + 10, M = 1e6 + 10;
int head[N], ver[M], edge[M], next[M], d[N];
bool v[N];
int n, m, s, tot;
std::priority_queue<std::pair<int, int> > q;void add(int u, int v, int w) {ver[++tot] = v, edge[tot] = w, next[tot] = head[u], head[u] = tot;
}void dijkstra(int s) {std::memset(d, 0x3f, sizeof d);d[s] = 0;q.push(std::make_pair(0, s));while (q.size()) {//取出栈顶int x = q.top().second;q.pop();if (v[x]) continue;v[x] = true;//扫描所有出边for (int i = head[x]; i; i = next[i]) {int y = ver[i], z = edge[i];if (d[y] > d[x] + z) {//更新,把新的二元组插入堆d[y] = d[x] + z;q.push(std::make_pair(-d[y], y));}}}
}int main() {std::cin >> n >> m >> s;//构建邻接表for (int i = 1; i <= m; i++) {int u, v, w;add(u, v, w);}dijkstra(s);for (int i = 1; i <= n; i++) {std::cout << d[i] << " \n"[i == n];}return 0;
}

Bellman - Ford 算法和 SPFA 算法

       $$ 给定一张有向图,若对于图中的某一条边 (x, y, z),有 dist[y] \le dist[x] + z 成立 $$

$$ 则称该边满足三角形不等式。 $$

$$ 若所有边都满足三角形不等式,则 dist 数组就是所求最短路 $$

$$ 下面介绍SPFA算法 $$

const int N = 1e5 + 10, M = 1e6 + 10;
int head[N], ver[M], edge[M], next[M], d[N];
int n, m, tot, s;
std::queue<int> q;
bool v[N];void add(int u, int v, int w) {ver[++tot] = v, edge[tot] = w, next[tot] = head[u], head[u] = tot;
}void spfa(int s) {std::memset(d, 0x3f, sizeof d);d[s] = 0, v[s] = true;q.push(s);while (q.size()) {//取出对头int x = q.front();q.pop();v[x] = false;//扫描所有出边for (int i = head[x]; i; i = next[i]) {int y = ver[i], z = edge[i];if (d[y] > d[x] + z) {//更新,把新的二元组插入堆d[y] = d[x] + z;if (!v[y]) {q.push(y), v[y] = true;}}}}
}

http://www.ppmy.cn/news/39070.html

相关文章

Node-CJS和ESM混合开发

Node — CJS和ESM混合开发 在JavaScript中&#xff0c;有两种主要的模块系统&#xff1a;CommonJS和ES Modules&#xff08;ESM&#xff09;。 CJS通过require引入&#xff0c;module.exports导出 ESM通过import引入&#xff0c;export导出 判断JS使用CJS还是ESM .cjs后缀的…

超级进化吧switch case in java

switch case的进化历程&#x1f690;简介区别Java 6及以下 | 整数当道Java 7 | String来贺Java 8 | 常量不够, 表达式来凑Java 12 | ->替: yield能返回Java 14 | ->凑多值,转型在心中转型Java 15 | 哎呦喂,都能定义变量了Java 17 | 整数不够, 长整型来凑总结因为长情是古…

【消息队列】聊一下如何避免消息的重复消费

什么是重复消费 一条消息在传输过程中&#xff0c;为了保证消息的不丢失&#xff0c;可能会多少量的消息进行重试&#xff0c;这样就可能导致Broker接受到的消息出现重复&#xff0c;如果说下游系统没有针对业务上的处理&#xff0c;那么可能导致同一笔借款或者支付订单出现重…

【网口交换机:交换机KSZ9897学习-笔记-资料汇总-记录】

【网口交换机&#xff1a;交换机KSZ9897学习-笔记-资料汇总-记录】1、概述2、 自己的学习与摸索之路第一阶段&#xff1a;随意在网上查找相关资料第二阶段&#xff1a;针对性在网上资料第三阶段&#xff1a;测试并且使用开发板第四阶段&#xff1a;针对性使用工具进行测试。2、…

操作系统-内存管理

一、总论 1.1 硬件术语 ​ 为了不让读者懵逼&#xff08;主要是我自己也懵逼&#xff09;&#xff0c;所以特地整理一下在后面会用到术语。 ​ 我们电脑上有个东西叫做内存&#xff0c;他的大小比较小&#xff0c;像我的电脑就是 16 GB 的。它是由 ROM 和 RAM 组成的&#x…

(5)(5.9) 推力损失和偏航不平衡警告

文章目录 前言 1 潜在的推力损失 2 偏航不平衡 前言 如果你看到推力损失或偏航不平衡的警告&#xff0c;这个页面概述了一些应该做的检查和修改来解决这个问题。在大多数情况下&#xff0c;这些警告是由于错误的硬件选择或设置造成的。 这些警告是为了检测推进系统的硬件故…

【C语言】预处理和程序环境

目录 程序的环境 运行环境 翻译环境 编译的过程 预编译阶段 编译阶段 汇编阶段 链接阶段&#xff08;不属于编译阶段&#xff09; 预处理详解 预定义符号 #define #define定义标识符 #define定义宏 #define的替换规则 #和##的使用 带副作用的宏参数 宏和…

二叉树的5个性质【要点:完全二叉树的性质】

只讲不会的 普通二叉树就要讲排列顺序了&#xff01;&#xff01;&#xff01; 预备&#xff1a;满二叉树&#xff1a;1.前提是它必须是二叉树 2.每个结点&#xff08;除了终端结点外&#xff09;都是2个子女。 要点1&#xff1a;关于普通的树的结点的计算&#xff0…