数学建模——第一章 最短路问题.pptVIP

迭代次数 L(ui) U0 u1 u2 u3 u4 u5 u6 u7 1 2 3 4 5 6 7 8 0 ∞ ∞ ∞ ∞ ∞ ∞ ∞ 2 1 8 ∞ ∞ ∞ ∞ 2 8 ∞ ∞ 10 ∞ 8 3 ∞ 10 ∞ 8 6 10 12 7 10 12 9 12 12 1.2.3 每对顶点之间的最短路 求每对顶点之间最短路的算法是Floyd算法: 1.算法的基本思路: 直接在图的带权邻接矩阵中用插入顶点的方法依次构造出p个矩阵D(1),D(2), …,D(p),使最后得到的矩阵D(p)成为图的距离矩阵,同时也求出插入点矩阵以便得到两点间的最短路径. 2.算法原理: (1)求距离矩阵的方法: 把带权邻接矩阵W作为距离矩阵的初值,即: D(0)=(dij(0))p×p=W 10 D(1)=(dij(1)) p×p,其中dij(1)=min{dij(0),di1(0)+d1j(0)},dij(1)是从vi到vj的只允许以v1作为中间点的路径中最短路的长度. 20 D(2)=(dij(2)) p×p,其中dij(2)=min{dij(1),di2(1)+d1j(1)},dij(2),是从vi到vj的只允许以v1,v2作为中间点的路径中最短路的长度. … … … … P0 D(p)=(dij(p)) p×p,其中dij(p)=min{dij(p-1),dip(p-1)+dpj(p-1)},dij(p)是从vi到vj的只允许以v1,v2, …,vp作为中间点的路径中最短路的长度.即是从vi到vj中间可插入任何顶点的路径中最短路的长度,因此:D(P)即是距离矩阵. (2)求路由矩阵的方法: 在建立距离矩阵的同时可建立路由矩阵R,R=(rij) p×p,rij的含义是从vi到vj的最短路要经过点号为rij的点. R(0)=(rij(0)) p×p,rij(0)=j 每求得一个D(k)时,按下列方式产生相应的新的R(k): rij(k)= 即当vk被插入任何两点间的最短路径时,被记录在R(k)中,依次求D(p)时求得R(p),可由R(p)来查找任何点对之间最短路的路由. (3)查找最短路路由的方法: 若rij(p)=p1,则点p1是点i到j点的最短路的中间点,然后用同样的方法再分头查找,若: ①向点i追溯得: ②向点j追溯得: 则由点i到j的最短路的路由为:i,pk, …,p2,p1,q1,q2, …,qm,j 3.算法步骤: Floyd算法:求任意两点间的最短路: D(i,j):i到j的距离; R(i,j):i到j的插入点. 输入带权邻接矩阵W(i,j), (1)赋初值;对所有i,j: d(i,j)←W(i,j),r(i,j)←j,k←1 (2)更新d(i,j), r(i,j) :对所有i,j,若: d(i,k)+ d(k,j)d(i,j)则: d(i,j) ← d(i,k)+d(k,j), r(i,j) ← k (3)若k=p,停止;否则:k ← k+1,转(2). 下面给一个程序： Floyd-Warshall算法 Input:非负元素的n×n矩阵[cij] Output: n×n矩阵[dij] ，其中dij是在弧权为cij下，自i到j的最短路长度。 Begin for all i≠j do dij:=cij; for i=1, …,n do dij:=∞; for j=1, …,n do for i=1, …,n i≠j do for k=1, …,n k≠j do dik:=min{dik,dij+djk} end 例2:求下面加权图的任意两点间的距离与路径. 解: v1 v2 v3 v4 v5 9 3 2 2 4 7 插入v1得: 矩阵中带下划线的元素为经迭代比较以后有变化的元素,即需要引入中间点v1,从而R(1)中相应的位置换为1.