Dijkstra算法介绍+正确性证明+性能分析

算法介绍

  • 源点s,数组d[u]表示s到u的最短距离,空集S,点集Q
  • 初始化:将源点s从点集中去掉,加入S,d[s]=0,∀v∈Q,d[v]=w[s][v]\forall v\in Q ,d[v]=w[s][v]vQ,d[v]=w[s][v]
  • 将Q中d[v]最小的点去掉加入S,并对u∈Q,w[v][u]<∞u\in Q,w[v][u]<\inftyuQ,w[v][u]<进行松弛操作:如果d[v]+w[v,u]<d[u],d[u]=d[v]+w[v,u]d[v]+w[v,u]<d[u],d[u]=d[v]+w[v,u]d[v]+w[v,u]<d[u],d[u]=d[v]+w[v,u]
  • 重复上述步骤直到Q为空集

正确性证明

Dijkstra算法条件:没有负边

设源点为s,d[v]表示源点到点v的举例,d[u,v]表示点u到点v的举例,w[u,v]表示从点u到点v的一条有向边的长度,δ(s,v)\delta(s,v)δ(s,v)表示源点到点v的最短路径长度。

最优子结构性质:任意两点之间的最短路径的子路径仍然是最短路径

证明:假设s到v的最短路径上有子路径u到w,那么δ(s,v)=d[s,u]+d[u,w]+d[w,v]\delta(s,v)=d[s,u]+d[u,w]+d[w,v]δ(s,v)=d[s,u]+d[u,w]+d[w,v]d[u,w]>δ(u,w)d[u,w]>\delta(u,w)d[u,w]>δ(u,w)

那么将d[u,w]d[u,w]d[u,w]替换为δ(u,w)\delta(u,w)δ(u,w)可以得到更好的路径,那么原路径就不知最短路径,这和δ(s,v)\delta(s,v)δ(s,v)的定义矛盾。因此该问题具有最优子结构性质。

引理1:初始化以后进行松弛操作前后∀v∈V,d[v]⩾δ(s,v)\forall v \in V, d[v]\geqslant \delta(s,v)vV,d[v]δ(s,v)

证明:由初始化条件,刚初始化以后∀v∈V,d[v]⩾δ(s,v)\forall v \in V, d[v]\geqslant \delta(s,v)vV,d[v]δ(s,v)

假设当对点v进行松弛操作时第一次导致上式不再成立,即d[v]=d[u]+w[u,v]<δ(s,v)d[v]=d[u]+w[u,v]<\delta(s,v)d[v]=d[u]+w[u,v]<δ(s,v)

因为点u在点v之前,所以点u满足上述不等式,d[u]⩾δ(s,u)d[u]\geqslant \delta(s,u)d[u]δ(s,u)

由定义容易得到w[u,v]⩾δ(u,v)w[u,v]\geqslant \delta(u,v)w[u,v]δ(u,v),因此d[v]=d[u]+w[u,v]⩾δ(s,u)+δ(u,v)d[v]=d[u]+w[u,v]\geqslant \delta(s,u)+\delta(u,v)d[v]=d[u]+w[u,v]δ(s,u)+δ(u,v)

由三角不等式d[v⩾δ(s,u)+δ(u,v)⩾δ(s,v)d[v\geqslant \delta(s,u)+\delta(u,v)\geqslant \delta(s,v)d[vδ(s,u)+δ(u,v)δ(s,v),与上面假设矛盾

因此初始化以后进行松弛操作前后∀v∈V,d[v]⩾δ(s,v)\forall v \in V, d[v]\geqslant \delta(s,v)vV,d[v]δ(s,v),引理1成立

引理2:假设存在s到v的最短路径:s->…->u->v,d[u]=δ(s,u)d[u]=\delta(s,u)d[u]=δ(s,u),那么在点u对v进行松弛操作d[v]=d[u]+w[u,v]d[v]=d[u]+w[u,v]d[v]=d[u]+w[u,v]之后,d[v]=δ(s,v)d[v]=\delta(s,v)d[v]=δ(s,v)

证明:由最优子结构性质,s到v的路径中s到u为点u的最短路径,δ(s,v)=δ(s,u)+w[u,v]\delta(s,v)=\delta(s,u)+w[u,v]δ(s,v)=δ(s,u)+w[u,v]

d[u]+w[u,v]=δ(s,u)+w[u,v]=δ(s,v)d[u]+w[u,v]=\delta(s,u)+w[u,v]=\delta(s,v)d[u]+w[u,v]=δ(s,u)+w[u,v]=δ(s,v),由引理1,在松弛操作以前d[v]⩾δ(s,v)d[v]\geqslant \delta(s,v)d[v]δ(s,v),因此松弛操作会保证d[v]=δ(s,v)d[v]=\delta(s,v)d[v]=δ(s,v),等证。

当算法结束后,d[v]=δ(s,v)d[v]=\delta(s,v)d[v]=δ(s,v)

证明:因为算法中将点加入确定的集合后就不再进行修改,所以相当于证明加入确定集合时d[v]=δ(s,v)d[v]=\delta(s,v)d[v]=δ(s,v)

由引理1,假设加入确定集合时第一个出现d[v]>δ(s,v)d[v]>\delta(s,v)d[v]>δ(s,v)为节点v(开始使用反证法)

假设s到v的最短路径为s->…->x->y->…->v,x是该路径在确定集合中的最后一个元素,y是该路径不在确定集合中的第一个元素。

因为v是第一个出现的点,而x在v之前加入,所以d[x]=δ(s,x)d[x]=\delta(s,x)d[x]=δ(s,x)

由最优子结构,s->…->x->y是s到y的一条最短路径,此时x已经进行了松弛操作

由引理2,d[y]=δ(s,y)d[y]=\delta(s,y)d[y]=δ(s,y)

因为y是s到v路径上的一点,所以δ(s,v)⩾δ(s,y)\delta(s,v)\geqslant \delta(s,y)δ(s,v)δ(s,y)

加入确定集合时d[v]d[v]d[v]是所有未加入点中最小的,所以d[v]⩽d[y]=δ(s,y)⩽δ(s,v)d[v]\leqslant d[y]=\delta(s,y)\leqslant\delta(s,v)d[v]d[y]=δ(s,y)δ(s,v),与假设矛盾。证毕。

复杂度分析

初始化的复杂度为Θ(∣V∣)\Theta(|V|)Θ(V)

总共要循环|V|次寻找点集中d[]最小的元素,复杂度为O(∣V∣T寻找最小元素)O(|V|T_{寻找最小元素})O(VT)

对于每个点都需要对所指向的点进行降低键值操作,由握手定理,总共需要降低出度的总和即|E|,复杂度为O(∣E∣T降低键值)O(|E|T_{降低键值})O(ET)

因此总共的复杂度为O(∣V∣T寻找最小元素+∣E∣T降低键值)O(|V|T_{寻找最小元素}+|E|T_{降低键值})O(VT+ET),通过使用不同的数据结构我们可以得到不同的复杂度。特殊的,如果我们使用数组,复杂度为O(∣V∣2)O(|V|^2)O(V2),如果我们使用二叉堆,复杂度为O((V+E)logV)O((V+E)logV)O((V+E)logV),如果我们使用斐波那契堆,可以得到最优秀的复杂度O(E+VlogV)O(E+VlogV)O(E+VlogV)

和BFS的关系

特殊地,如果边的权值全部为1(或者相等),那么BFS就是Dijkstra算法,复杂度为O(V+E)O(V+E)O(V+E)

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/383647.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Linux C 实现一个简单的线程池

https://www.cnblogs.com/GyForever1004/p/9185240.html 线程池的定义 线程池是一种多线程处理形式&#xff0c;处理过程中将任务添加到队列&#xff0c;然后在创建线程后自动启动这些任务。线程池线程都是后台线程。每个线程都使用默认的堆栈大小&#xff0c;以默认的优先级…

斐波那契数列求解+尾递归

1.普通递归 这里观察f[4]的递归树代替f[10]的递归树&#xff08;后者比较大&#xff0c;画不下&#xff09;。 使用递归求解的时候复杂度为T(n)T(n−1)T(n−2)T(n)T(n-1)T(n-2)T(n)T(n−1)T(n−2)&#xff0c;观察递归树&#xff0c;发现降速最快的是最右边每次减2&#xff0c…

循环服务器,并发服务器模型以及I/O多路转接模型

https://blog.csdn.net/xinianbuxiu/article/details/53455784 一、基于TCP/IP协议的基本循环服务器 tcp_server.c #include <stdio.h> #include <stdlib.h> #include <string.h> #include <sys/types.h> #include <sys/socket.h> #incl…

c++继承父类的子类,如何调用父类的同名函数?

https://blog.csdn.net/qq_26399665/article/details/52080215 子类调用父类的同名函数&#xff1a; 子类和父类返回值参数相同&#xff0c;函数名相同&#xff0c;有virtual关键字&#xff0c;则由对象的类型决定调用哪个函数。 子类和父类只要函数名相同&#xff0c;没有vi…

LCS最长公共子串

问题介绍 LCS问题(longest common subsequence problem)指的是求解两个字符串最长公共子序列问题。这里的子序列是可以不连续的。LCS问题广泛地出现在计算生物学中&#xff08;DNA序列、系统生成树等等&#xff09;。这里介绍如何解决LCS问题&#xff0c;以及算法的正确性证明…

将字符串中的空格用%20替换

如果不需要原地操作&#xff0c;则一遍遍历&#xff0c;将非空串复制&#xff0c;遇到空格加上%20&#xff0c;如果需要原地操作&#xff0c;首先进行遍历出空格的个数x,然后扩容2x,从后往前遍历实现。如果非空格字符串比空格字符串多的多的时候而且字符串非常长的时候使用原地…

12步轻松搞定python装饰器

http://python.jobbole.com/81683/ 呵呵&#xff01;作为一名教python的老师&#xff0c;我发现学生们基本上一开始很难搞定python的装饰器&#xff0c;也许因为装饰器确实很难懂。搞定装饰器需要你了解一些函数式编程的概念&#xff0c;当然还有理解在python中定义和调用函数…

操作系统【六】虚拟内存

传统存储管理方式的不足 一次性&#xff1a;作业必须一次性全部装入内存后才能开始运行。这会造成&#xff1a;当作也很大时不能全部装入内存&#xff1b;当大量作业要求运行时&#xff0c;由于内存无法容纳所有作业&#xff0c;因此只有少量作业能够运行&#xff0c;导致多道…

python装饰器详解

https://blog.csdn.net/xiangxianghehe/article/details/77170585 你会Python嘛&#xff1f; 我会&#xff01; 那你给我讲下Python装饰器吧&#xff01; Python装饰器啊&#xff1f;我没用过哎 以上是我一个哥们面试时候发生的真实对白。 ———————————————-分…

SQL Server【一】简介和基本概念和命令

数据结构和数据库的区别 数据库是应用软件级别研究数据的存储和操作&#xff08;主要针对磁盘上的数据&#xff09; 数据结构是在系统软件级别研究数据的存储和操作&#xff08;主要是针对内存中的数据&#xff09; 对硬盘数操作是数据库的强项&#xff0c;是数据库研究的核心…

SQL Server【二】单表查询

查询 计算列 select * from emp; -- *通配符&#xff0c;表示所有的字段 -- from emp 从emp表查询select empno, ename from emp; select ename as "员工姓名", sal*12 as "年薪" from emp;-- as可以省略&#xff0c;用于设置字段名 -- 注意用双引号将字…

SQL Server【三】连接查询

将两个表或者两个以上的表以一定的连接条件连接起来&#xff0c;从中检索出满足条件的数据。 内连接 使用inner join&#xff0c;inner可以省略 -- 查询员工的姓名和部门名称 select "E".ename as "员工姓名", "D".dname as "部门名称&q…

Linux下网络socket编程——实现服务器(select)与多个客户端通信

一、关于socket通信 服务器端工作流程&#xff1a; 调用 socket() 函数创建套接字 用 bind() 函数将创建的套接字与服务端IP地址绑定调用listen()函数监听socket() 函数创建的套接字&#xff0c;等待客户端连接 当客户端请求到来之后调用 accept()函数接受连接请求&#xff0c…

SQL Server【四】

identity 主键自动增长&#xff0c;用户不需要为identity修饰的主键赋值 create table student (std_id int primary key identity(10,5),--(10,5)可以省略&#xff0c;默认为(1,1)std_name nvarchar(200) not null ) select * from student insert into student values (张三…

TCP服务器/客户端实例(C/C )

1.1、Linux下的TCP服务器&#xff1a; #include <stdio.h> #include <stdlib.h> #include <string.h> #include <unistd.h> #include <arpa/inet.h> #include <sys/types.h> #include <sys/socket.h>void error_handling(char *mess…

pip代理解决pip下载失败问题

在用pip下载各种库的时候发现速度实在是太慢了&#xff0c;还会有各种奇奇怪怪的问题&#xff0c;动不动就玄学失败。 在网上找来找去找到知乎上一位大佬的回答&#xff1a;传送门&#xff0c;用了豆瓣的代理。哇咔咔&#xff0c;妈妈再也不用担心我下载失败了。 代理&#x…

实现Linux select IO复用C/S服务器代码

服务器端#include<stdio.h> #include<unistd.h> #include<stdlib.h> #include<string.h> #include<sys/socket.h> #include<sys/stat.h> #include<arpa/inet.h> #include <sys/select.h>#define MAXBUF 256 #define MAXLISTEN…

Bellman-Ford算法和SPFA算法

Belloman-Ford算法 算法介绍 Dijkstra可以解决单源无负边最短路径问题。但是当遇到含有负边的单源最短路径问题就需要使用Bellman-Ford算法来解决。Bellman-Ford算法还可以检测出负环。 算法步骤 源点s,数组d[u]d[u]d[u]表示s到u的最短距离初始化&#xff1a;d[s]0d[s]0d[s…

C语言实现单链表操作

SLIST_H #ifndef __SLIST_H__ #define __SLIST_H__ #include<cstdio> #include<malloc.h> #include<assert.h> typedef int ElemType; typedef struct Node { //定义单链表中的结点信息 ElemType data; //结点的数据域 struct Node *next; //结点的指针…

计算机网络【4】传输层

概述 传输层是只有主机才有的层次 传输层的功能&#xff1a; 传输层提供进程和进程之间的逻辑通信&#xff08;网络层提供主机与主机之间的逻辑通信&#xff09;复用和分用传输层对收到的报文进行差错检测 传输层有两个协议&#xff1a; 面向连接的传输层控制协议TCP&…