第4篇1：图论模型(最优连线问题、最短路问题)

1, vi 是e j的起点 mij 1, vi 是e j的终点 0, vi与e j 不关联
例1
v1 e1 v2
e2 e5
v4
1 e4 1 M 0 0
v3
1 1 0 1 0 0 1 1 0 1 1 0 1 0 0 0
e3
② 邻接矩阵
rij ( n ) p1 ，则点p1是点i到点j最短路的中间点， 若 R ( k ) (rij ( k ) ) nn , (k 1, 2,..., n) rij ( k ) 的意思是从ui到uj的最短路要经过点号为 rij ( k ) 的
然后用同样的方法再分头查找：
rip1 ( n ) p2 , rip2 ( n ) p3 ,..., ripk ( n ) pk ①向点i追溯得： rp1 j ( n ) q1 , rq1 j ( n ) q2 ,,..., rqm j ( n ) j ②向点j追溯得：
以u1作为可行中间点的路径中最短路的长度。得距离矩 阵 ； D (1) (dij (1) ) nn
dij (2) min{dij (1) , di 2 (1) d 2 j (1)} , D (2) (dij (2) ) nn , 第2步： d ij ( 2) 是从ui到uj的只允许以u1、u2作为可行中间点
Floyd算法程序：
%floyd.m（Floyd程序 for k=1:n 需要调用子函数） for i=1:n function [D,R]=floyd(a) for j=1:n n=size(a,1); if D(i,k)+D(k,j)<D(i,j) D=a; D(i,j)=D(i,k)+D(k,j); for i=1:n R(i,j)=R(i,k); for j=1:n end R(i,j)=j; end end end end k R D R end
%road2.m （floyd算法程序） a=[0 9 inf 3 inf; 9 0 2 inf 7; inf 2 0 2 4; 3 inf 2 0 inf; inf 7 4 inf 0]; [D,R]=floyd(a)
最后输出的矩阵D中元素aij即为ui和uj之间的最 短距离。
查找最短路径的方法: Floyd算法程序在上述各步中还输出一个路径矩 阵 点。
0
8
u8
1
u4
1
10
u7
第三步：在其余点中找从u1直接到达或从u1经u4 到达或从u1经u2到达的的距离最近的点u5， min{8,11,11,3,9}。
u2
2 7 5 3 9 1
u5
2
3
u1
u3
u6
1
5 10 2 5 6
0
8
u8
1
u4
1
10
u7
第四步：min{8,11,11,9,8,6,12}，u6。
2
1
u5
2
7 5
3
3 9
u1
u3
u6
0
8
8
10
1
2
6
5 5 6
u8
1
u4
1
10
10
u7
第七步：min{12,11,15}，u8。
u2
2
1
u5
2
7 5
3
3 9
u1
u3
u6
0
8
8
10
1
2
6
5 5 6
11 u8
1
u4
1
10
10
u7
%road1.m (dijkstra算法程序)
w=[0 2 1 8 inf inf inf inf; 2 0 inf 6 1 inf inf inf; 1 inf 0 7 inf inf 9 inf; 8 6 7 0 5 1 2 inf; inf 1 inf 5 0 3 inf 9; inf inf inf 1 3 0 4 6; inf inf 9 2 inf 4 0 3; inf inf inf inf 9 6 3 0]; n=size(w,1); w1=w(1,:); for i=1:n l(i)=w1(i); z(i)=1; end s=[];
无向图G，邻接矩阵A=（aij） 1, vi与v j 相邻（即vi与v j 之间有边） aij 0, vi与v j 不相邻 有向图G，邻接矩阵A=（aij）
1, 从vi到v j 有有向边 aij 0, 从vi到v j 无有向边 有向赋权图G，邻接矩阵A=（aij） wij , 从vi到v j 有有向边,且wij 为其权 aij 0, 若i j , 从v 到v 无有向边 i j
ch8 图论模型
图论是离散数学的重要分支，在物理学、化学、 系统工程、电力通讯、编码理论、可靠性理论、科 学管理、电子计算机等各个领域又具有极其广泛的 应用。 图论的历史可以追朔到1736年，这一年29岁的 瑞士大数学家Euler发表了图论的第一篇论文，解决 了著名的哥尼斯堡七桥问题。
现实生活中的公路交通网、铁路交通网、灌溉 网、自来水（石油、天然气）管道网、电话线网计 算机通讯网、输电线网等，都可以用上述图的方式 来描述和分析解决问题。
lv; v; s(k+1)=v; k=k+1; u=s(k); end l z 输出结果为： l=0 2 1 7 3 6 9 12 z=1 1 1 6 2 5 4 5
注：l输出的是u1到u1、u2、…、u8各个顶点的 最短路径距离。 z输出的是最短路径中u1、u2、…、u8的父节点。
%求从u0到uj0的最短路径 disp('起点为u1.'); j=input('输入终点u'); disp('下面求从起点u1到终点'); j, disp('的最短路径。'); lj=[]; while j~=1 lj=[[j],lj]; j=z(j); end lj=[[1],lj]; lj 例如求u1到u8的最短路径，程序执行后输出为：1 2 5 6 8 各位有兴趣还可以考虑将图可视化，点击屏幕输入终点以及 在图形上输出显示最短路径。
对于有向图的顶点的度数，还可分为出度 d (v) 和 入度 d (v) 。 定理：

d (v ) 2 m
i 1 i
n
4 图的矩阵表示 ① 关联矩阵 无向图G，关联矩阵M=（mij）
1, vi与e j 关联 mij 0, vi与e j 不关联
有向图G，关联矩阵M=（mij）
2 Floyd算法
Dijkstra算法的缺点是只能求起点到其它个顶点的 最短路径，不能求出任意两点之间的路径, Floyd算法可 解决该问题。 Floyd算法简介：
D (0) (dij (0) ) nn W 将加权矩阵W作为距离矩阵的初值 d ij (0) 和ui经过u1到 第1步：比较ui和uj之间直达的距离 di1(0) d1 j (0) ，最小值作为第一步的 达uj的路径距离 (1) (0) (0) (0) d ij (1) 值dij min{dij , di1 d1 j } ， 是从ui到uj的只允许
七桥问题
著名数学家欧拉
无向图 1 定义：由顶点和边组成的图形称为图。 有向图 赋权图 2 边e与顶点u、v相关联。顶点u与v相邻。
e u v
图的基本概念
边e1与e2相邻。 u=v时，边e称为环。
e1
e2
3度
定义：与顶点v关联的边的数目称为顶点v的度数， 记为d(v)。（注：环算2度。）
的路径中最短路的长度；
…… dij ( n ) min{dij ( n 1) , di 2 ( n 1) d 2 j (n 1) } , D (n ) (dij (n ) )nn , 第n步： d ij ( n )是从ui到uj的只允许以u1、u2、……、un作为 可行中间点的路径中最短路的长度； 则最后输出的矩阵D中元素aij即为ui和uj之间的最短 距离。
u6
5 5
u8
6
2
u4
10
u7
解：算法原理为蚂蚁算法（探索算法），每次新连 接一个点。所有新到一个点最短路程中最短的那个 店，作为新增点。 第一步：找从u1出发到达的距离最近的点u4， min{2,8,1}，将这个距离写进u4的圈中。将从u1到 u4的边描成红色，u4成为永久标记的点；
u2
2 7 5 3 9 1
注：inf指 无穷大。
s(1)=1; u=s(1); k=1; l; z; while k<n for i=1:n for j=1:k if i~=s(j) if l(i)>l(u)+w(u,i) l(i)=l(u)+w(u,i); z(i)=u; end end end end l; z;
ll=l; for i=1:n for j=1:k if i~=s(j) ll(i)=ll(i); else ll(i)=inf; end end end lv=inf; for i=1:n if ll(i)<lv lv=ll(i); v=i; end end
（其思想是：在剩余点集中找连接到U中顶点的最小权重的边，添加到生成树中。
（显然不会产生回路）。仍然是以局部的最优谋求全局的最优。）
上例中，采用Prim算法最小生成树的生 长过程：
A 8 5 B
7 E
6 水厂
1
3
10
D
9
C
例：如何设计海底管道网。（ Prim 算法）
1 2 6 8 7 1 9 1 3 7 6 4 4 2 9 2
算法如下：
（1）选择边 e1 ，使得 w(e1 ) 最小； （2）若已经选定边 e1 , e2 ,..., ei ，则从剩余边集