文章1——程序最短路.docVIP

最短路问题 0817010001聪 本文主要研究最短路问题，给出了Matlab的编程和winQSB在最短路问题的体现。 一、带负权的图的最短路问题 首先我们考虑的是对于适用性范围广的最短路问题，即允许赋权有向图存在着负权，但是我们考虑不应该存在负回路，因为存在负回路，它会不收敛。那么，我们首先考虑的这种方法蕴含的最重要的思想在于递推式，这其实的本质也就是逐次逼近算法思想的运用。而这里，我们知道这个最短路问题从起始点到任一点的最短路，用合适的算法最多只需要图的总的点数减去1次迭代必定收敛。那么，我们可以用编出程序： function D=short(T) n=length(T); D=T; m=1; while m=n for i=1:n for j=1:n if D(i,j)D(i,m)+D(m,j) D(i,j)=D(i,m)+D(m,j); end end end m=m+1; end D; 其中，T为图的矩阵存储。我们考虑课本266页的例子，我们在命令窗口键入 T=[0,-1,-2,3,inf,inf,inf,inf; 6,0,inf,inf,2,inf,inf,inf; inf,-3,0,-5,inf,1,inf,inf; inf,inf,inf,0,inf,inf,2,inf; inf,-1,inf,inf,0,inf,inf,inf; inf,inf,inf,inf,1,0,1,7; inf,inf,inf,-1,inf,inf,0,inf; inf,inf,inf,inf,-3,inf,-5,0]; D=short(T) 则我们运行后可以得到运行结果： D = 0 -5 -2 -7 -3 -1 -5 6 6 0 4 -1 2 5 1 12 3 -3 0 -5 -1 1 -3 8 Inf Inf Inf 0 Inf Inf 2 Inf 5 -1 3 -2 0 4 0 11 6 0 4 -1 1 0 1 7 Inf Inf Inf -1 Inf Inf 0 Inf 2 -4 0 -6 -3 1 -5 0 这里我们得到的结果不仅可以看出v1到其他点的最短路，就在矩阵第一行，而且我们还可以得到v2到其他点的最短路，其他点的情况依此类推。 二、Dijkstra算法编程方法1 对于Dijkstra算法，它的适用范围是=0,而此时，我们可以考虑求出它的起始点与其他点的最短路，我们一般用标号的方法去进行探究，有P标号、T标号、标号。下面，我们考虑求出两点间的最短路，并且求出P标号的顺序。 在编程前，我们首先引入行向量pp,Q,d分别来表示存放P标号的信息，P标号的顶点顺序，最短路的值。而我们令则我们可以得到程序： function [d Q]=shorta(T) pp(1:length(T))=0; pp(1)=1;Q=1; M=max(max(T)); d(1:length(T))=M;d(1)=0;K=1; while sum(pp)length(T) tt=find(pp==0); d(tt)=min(d(tt),d(K)+T(K,tt)); ttt=find(d(tt)==min(d(tt))); K=tt(ttt(1)); pp(K)=1; Q=[Q,K]; end d;Q; 我们考虑课本263页的例子，在命令窗口键入： T=[0,6,3,1,inf,inf,inf,inf,inf; inf,0,inf,inf,1,inf,inf,inf,inf; inf,2,0,2,inf,inf,inf,inf,inf; inf,inf,inf,0,inf,10,inf,inf,inf; inf,inf,inf,inf,0,4,3,6,inf; inf,inf,inf,inf,10,0,2,inf,inf; inf,inf,inf,inf,inf,inf,0,4,inf; inf,inf,inf,inf,inf,inf,inf,0,inf; inf,inf,in