两种经典最短路径问题_Dijkstra和Floyd算法
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
4
最短路径问题
定义：设P(u,v)是加权图G中从u到v的路径,则该路
径上的边权之和称为该路径的权,记为w(P). 从u到v 的路径中权最小者 P*(u,v)称为u到v的最短路径.
3 5
2
8
3
5
6 1
4
2 7
12 10
1
8
7
7
9
6
3
5
2
11
8
9
9
2
10
5
最短路径算法
1
Dijkstra算法
2
8 7
4
6 1
D便是最廉价的航费表， 要求飞行路线，由path矩 阵可以得到，比如2到5的 路线：path(2,5)=4, path(4,5)=5,因此，应为 2→4 →5
17
1
2
4
4
1
6
3 5
2
8 5 7 8
3
6 1
4
2 7
12 10
1
7
9
6
3
5
2 2
11
8
9
9
10
3
引例2：最廉价航费表的制定
某公司在六个城市C1,C2,C3,C4,C5,C6都有分公司， 公司成员经常往来于它们之间，已知从Ci到Cj的直达航 班票价由下述矩阵的第 i 行，第 j 列元素给出（ 表示无 直达航班），该公司想算出一张任意两个城市之间的最 廉价路线航费表。 0 50 40 25 10 50 0 15 20 25 15 0 10 20 40 20 10 0 10 25 25 20 10 0 55 10 25 25 55 0
8
3 5
3 6
3
5
4
1
6 2
7 2
12 10
7 8
9 9
5
10
11
使用范围:
1)
2
9
寻求从一固定顶点到其余各点的最短路径;
2)
3)ቤተ መጻሕፍቲ ባይዱ
有向图、无向图和混合图;
权非负.
算法思路： 采用标号作业法,每次迭代产生一个永久标号, 从而生长一颗以v0为根的最短路树,在这颗树上每个
顶点与根节点之间的路径皆为最短路径.
11
MATLAB程序（Floyd算法）
function [D,path,min1,path1]=floyd(a,start,terminal) D=a;n=size(D,1);path=zeros(n,n); for i=1:n if nargin==3 for j=1:n min1=D(start,terminal); if D(i,j)~=inf m(1)=start; path(i,j)=j; i=1; end, end, end path1=[ ]; for k=1:n while path(m(i),terminal)~=terminal for i=1:n k=i+1; for j=1:n m(k)=path(m(i),terminal); if D(i,k)+D(k,j)<D(i,j) i=i+1; D(i,j)=D(i,k)+D(k,j); end path(i,j)=path(i,k); m(i+1)=terminal; end, end, end,end path1=m; end
13
引例1的Matlab求解
edge= [ 2,3,1,3,3,5,4, 4,1,7,6,6,5, 5,11, 1,8,6,9,10,8,9, 9,10;...
3,4,2,7,5,3,5,11,7,6,7,5,6,11, 5, 8,1,9,5,11,9,8,10,9;... 3,5,8,5,6,6,1,12,7,9,9,2,2,10,10,8,8,3,7, 2, 9,9, 2, 2]; n=11; weight=inf*ones(n, n); for i=1:n weight(i, i)=0; end for i=1:size(edge,2) end 8 1 7 8
[D, path]=floyd(a)
运行便可输出结果。
0 50 40 25 10
50 40 25 10 0 15 20 25 15 0 10 20 20 10 0 10 25 20 10 0 55 25 25 55 0
16
运行输出结果： D =
6
Dijkstra算法——算法步骤
S: 具有永久标号的顶点集; l(v): v的标记; f(v):v的父顶点,用以确定最短路径; 输入加权图的带权邻接矩阵w=[w(vi,vj)]nxm.
1) 2)
初始化
令l(v0)=0,S=; vv0 ,l(v)=;
更新l(v), f(v) 寻找不在S中的顶点u,使l(u)为最小.把u加入到S中, 然后对所有不在S中的顶点v,如l(v)>l(u)+w(u,v),则 更新l(v),f(v), 即 l(v)l(u)+w(u,v),f(v)u;
9
最短路径算法
1
2
8 7
8
3 5
3 6
3
5
4
1
6 2
7 2
12 10
7 8
9 9
5
10
11
Floyd算法
使用范围:
1) 2)
2
9
求每对顶点的最短路径; 有向图、无向图和混合图;
算法思想: 直接在图的带权邻接矩阵中用插入顶点的方法依 次递推地构造出n个矩阵D(1), D(2), …, D(n), D(n) 是图的距离矩阵, 同时引入一个后继点矩阵记录两 点间的最短路径.
10
Floyd算法——算法步骤
d(i,j) : i到j的距离;
path(i,j): i到j的路径上i的后继点;
输入带权邻接矩阵a(i,j). 1）赋初值 对所有i,j, d(i,j)a(i,j) , path(i,j)j,k=l. 2）更新d(i,j) , path(i,j) 对所有i,j, 若d(i,k)+d(k,j)<d(i,j),则 d(i,j)d(i,k)+d(k,j) , path(i,j)path(i,k) , k k+1 3）重复2)直到k=n+1
因此顶点1到顶点11的最短路径为1→8 →9 →10 →11, 其长度为21。
15
引例2的Matlab求解
建立脚本m文件如下：
a= [ 0,50,inf,40,25,10;50,0,15,20,inf,25;inf,15,0,10,20,inf;… 40,20,10,0,10,25;25,inf,20,10,0,55;10,25,inf,25,55,0];
8
最短路径算法
Dijkstra算法程序的使用说明：
调用格式为
[min,path]=dijkstra(w,start,terminal),
其中输入变量w为所求图的带权邻接矩阵，start, terminal分别为路径的起点和终点的号码。返回start 到terminal的最短路径path及其长度min. 注意：顶点的编号从1开始连续编号。
最短路径问题
Dijkstra算法和Floyd算法
1
主要内容
引例1：最短运输路线问题
引例2：最廉价航费表的制定
Dijkstra算法
Floyd算法 两个例子的求解
2
引例1：最短运输路线问题
如图的交通网络，每条弧上的数字代表车辆在该路段行 驶所需的时间，有向边表示单行道，无向边表示可双向 行驶。若有一批货物要从 1号顶点运往 11号顶点，问运 货车应沿哪条线路行驶，才能最快地到达目的地？
12
最短路径算法
Floyd算法程序的使用说明：
1. [D, path]=floyd(a), 返回矩阵D, path 。其中a是所求 图的带权邻接矩阵，D(i,j)表示i到j的最短距离; path(i,j)表示i与j之间的最短路径上顶点i的后继点.
2. [D, path, min1, path1]= floyd(a,i,j) 返回矩阵D, path; 并返回i与j之间的最短距离min1和最短路径path1.
7
3)
重复步骤2), 直到所有顶点都在S中为止.
MATLAB程序（Dijkstra算法）
function [min,path]=dijkstra(w,start,terminal) n=size(w,1); label(start)=0; f(start)=start; for i=1:n min=label(terminal); if i~=start path(1)=terminal; label(i)=inf; i=1; end, end s(1)=start; u=start; while path(i)~=start path(i+1)=f(path(i)); while length(s)<n i=i+1 ; for i=1:n ③ end ① ins=0; for j=1:length(s) path(i)=start; L=length(path); if i==s(j) path=path(L:-1:1); ins=1; end, end if ins==0 v=i; if label(v)>(label(u)+w(u,v)) label(v)=(label(u)+w(u,v)); f(v)=u; end, end, end v1=0; k=inf; ② for i=1:n ins=0; for j=1:length(s) if i==s(j) ins=1; end, end if ins==0 v=i; if k>label(v) k=label(v); v1=v; end, end, end s(length(s)+1)=v1; u=v1; end
0 35 45 35 25 10 path = 1 6 6 2 5 3 5 4 5 4 6 6 35 0 15 20 30 25 45 15 0 10 20 35 35 20 10 0 10 25 25 30 20 10 0 35 10 25 35 25 35 0
5
5 1
2
2 4
3
3 3
4
4 4
5
5 5
2
7
3
5 9
3
5
6
4 1
5
10
9 8 weight(edge(1, i), edge(2, i))=edge(3, i);
6 2 3 7 2 9
1 12 0 2
11
[dis, path]=dijkstra(weight, 1, 11)
14
引例1的求解
运行上页程序输出： dis = 21 path = 1 8 9 10 11
最短路径问题
定义：设P(u,v)是加权图G中从u到v的路径,则该路
径上的边权之和称为该路径的权,记为w(P). 从u到v 的路径中权最小者 P*(u,v)称为u到v的最短路径.
3 5
2
8
3
5
6 1
4
2 7
12 10
1
8
7
7
9
6
3
5
2
11
8
9
9
2
10
5
最短路径算法
1
Dijkstra算法
2
8 7
4
6 1
D便是最廉价的航费表， 要求飞行路线，由path矩 阵可以得到，比如2到5的 路线：path(2,5)=4, path(4,5)=5,因此，应为 2→4 →5
17
1
2
4
4
1
6
3 5
2
8 5 7 8
3
6 1
4
2 7
12 10
1
7
9
6
3
5
2 2
11
8
9
9
10
3
引例2：最廉价航费表的制定
某公司在六个城市C1,C2,C3,C4,C5,C6都有分公司， 公司成员经常往来于它们之间，已知从Ci到Cj的直达航 班票价由下述矩阵的第 i 行，第 j 列元素给出（ 表示无 直达航班），该公司想算出一张任意两个城市之间的最 廉价路线航费表。 0 50 40 25 10 50 0 15 20 25 15 0 10 20 40 20 10 0 10 25 25 20 10 0 55 10 25 25 55 0
8
3 5
3 6
3
5
4
1
6 2
7 2
12 10
7 8
9 9
5
10
11
使用范围:
1)
2
9
寻求从一固定顶点到其余各点的最短路径;
2)
3)ቤተ መጻሕፍቲ ባይዱ
有向图、无向图和混合图;
权非负.
算法思路： 采用标号作业法,每次迭代产生一个永久标号, 从而生长一颗以v0为根的最短路树,在这颗树上每个
顶点与根节点之间的路径皆为最短路径.
11
MATLAB程序（Floyd算法）
function [D,path,min1,path1]=floyd(a,start,terminal) D=a;n=size(D,1);path=zeros(n,n); for i=1:n if nargin==3 for j=1:n min1=D(start,terminal); if D(i,j)~=inf m(1)=start; path(i,j)=j; i=1; end, end, end path1=[ ]; for k=1:n while path(m(i),terminal)~=terminal for i=1:n k=i+1; for j=1:n m(k)=path(m(i),terminal); if D(i,k)+D(k,j)<D(i,j) i=i+1; D(i,j)=D(i,k)+D(k,j); end path(i,j)=path(i,k); m(i+1)=terminal; end, end, end,end path1=m; end
13
引例1的Matlab求解
edge= [ 2,3,1,3,3,5,4, 4,1,7,6,6,5, 5,11, 1,8,6,9,10,8,9, 9,10;...
3,4,2,7,5,3,5,11,7,6,7,5,6,11, 5, 8,1,9,5,11,9,8,10,9;... 3,5,8,5,6,6,1,12,7,9,9,2,2,10,10,8,8,3,7, 2, 9,9, 2, 2]; n=11; weight=inf*ones(n, n); for i=1:n weight(i, i)=0; end for i=1:size(edge,2) end 8 1 7 8
[D, path]=floyd(a)
运行便可输出结果。
0 50 40 25 10
50 40 25 10 0 15 20 25 15 0 10 20 20 10 0 10 25 20 10 0 55 25 25 55 0
16
运行输出结果： D =
6
Dijkstra算法——算法步骤
S: 具有永久标号的顶点集; l(v): v的标记; f(v):v的父顶点,用以确定最短路径; 输入加权图的带权邻接矩阵w=[w(vi,vj)]nxm.
1) 2)
初始化
令l(v0)=0,S=; vv0 ,l(v)=;
更新l(v), f(v) 寻找不在S中的顶点u,使l(u)为最小.把u加入到S中, 然后对所有不在S中的顶点v,如l(v)>l(u)+w(u,v),则 更新l(v),f(v), 即 l(v)l(u)+w(u,v),f(v)u;
9
最短路径算法
1
2
8 7
8
3 5
3 6
3
5
4
1
6 2
7 2
12 10
7 8
9 9
5
10
11
Floyd算法
使用范围:
1) 2)
2
9
求每对顶点的最短路径; 有向图、无向图和混合图;
算法思想: 直接在图的带权邻接矩阵中用插入顶点的方法依 次递推地构造出n个矩阵D(1), D(2), …, D(n), D(n) 是图的距离矩阵, 同时引入一个后继点矩阵记录两 点间的最短路径.
10
Floyd算法——算法步骤
d(i,j) : i到j的距离;
path(i,j): i到j的路径上i的后继点;
输入带权邻接矩阵a(i,j). 1）赋初值 对所有i,j, d(i,j)a(i,j) , path(i,j)j,k=l. 2）更新d(i,j) , path(i,j) 对所有i,j, 若d(i,k)+d(k,j)<d(i,j),则 d(i,j)d(i,k)+d(k,j) , path(i,j)path(i,k) , k k+1 3）重复2)直到k=n+1
因此顶点1到顶点11的最短路径为1→8 →9 →10 →11, 其长度为21。
15
引例2的Matlab求解
建立脚本m文件如下：
a= [ 0,50,inf,40,25,10;50,0,15,20,inf,25;inf,15,0,10,20,inf;… 40,20,10,0,10,25;25,inf,20,10,0,55;10,25,inf,25,55,0];
8
最短路径算法
Dijkstra算法程序的使用说明：
调用格式为
[min,path]=dijkstra(w,start,terminal),
其中输入变量w为所求图的带权邻接矩阵，start, terminal分别为路径的起点和终点的号码。返回start 到terminal的最短路径path及其长度min. 注意：顶点的编号从1开始连续编号。
最短路径问题
Dijkstra算法和Floyd算法
1
主要内容
引例1：最短运输路线问题
引例2：最廉价航费表的制定
Dijkstra算法
Floyd算法 两个例子的求解
2
引例1：最短运输路线问题
如图的交通网络，每条弧上的数字代表车辆在该路段行 驶所需的时间，有向边表示单行道，无向边表示可双向 行驶。若有一批货物要从 1号顶点运往 11号顶点，问运 货车应沿哪条线路行驶，才能最快地到达目的地？
12
最短路径算法
Floyd算法程序的使用说明：
1. [D, path]=floyd(a), 返回矩阵D, path 。其中a是所求 图的带权邻接矩阵，D(i,j)表示i到j的最短距离; path(i,j)表示i与j之间的最短路径上顶点i的后继点.
2. [D, path, min1, path1]= floyd(a,i,j) 返回矩阵D, path; 并返回i与j之间的最短距离min1和最短路径path1.
7
3)
重复步骤2), 直到所有顶点都在S中为止.
MATLAB程序（Dijkstra算法）
function [min,path]=dijkstra(w,start,terminal) n=size(w,1); label(start)=0; f(start)=start; for i=1:n min=label(terminal); if i~=start path(1)=terminal; label(i)=inf; i=1; end, end s(1)=start; u=start; while path(i)~=start path(i+1)=f(path(i)); while length(s)<n i=i+1 ; for i=1:n ③ end ① ins=0; for j=1:length(s) path(i)=start; L=length(path); if i==s(j) path=path(L:-1:1); ins=1; end, end if ins==0 v=i; if label(v)>(label(u)+w(u,v)) label(v)=(label(u)+w(u,v)); f(v)=u; end, end, end v1=0; k=inf; ② for i=1:n ins=0; for j=1:length(s) if i==s(j) ins=1; end, end if ins==0 v=i; if k>label(v) k=label(v); v1=v; end, end, end s(length(s)+1)=v1; u=v1; end
0 35 45 35 25 10 path = 1 6 6 2 5 3 5 4 5 4 6 6 35 0 15 20 30 25 45 15 0 10 20 35 35 20 10 0 10 25 25 30 20 10 0 35 10 25 35 25 35 0
5
5 1
2
2 4
3
3 3
4
4 4
5
5 5
2
7
3
5 9
3
5
6
4 1
5
10
9 8 weight(edge(1, i), edge(2, i))=edge(3, i);
6 2 3 7 2 9
1 12 0 2
11
[dis, path]=dijkstra(weight, 1, 11)
14
引例1的求解
运行上页程序输出： dis = 21 path = 1 8 9 10 11