数学建模社区-数学中国
标题:
常用模型&算法总结—图&网络模型应用—树:连通性、最小生成树
[打印本页]
作者:
浅夏110
时间:
2020-5-20 09:49
标题:
常用模型&算法总结—图&网络模型应用—树:连通性、最小生成树
1 基本概念
" q c) D& N1 e w2 _, e" }1 N
连通的无圈图叫做树,记之为T 。若图G 满足V(G) =V(T ) , E(T ) ⊂ E(G) , 则称T 是G 的生成树。图G 连通的充分必要条件为G 有生成树。一个连通图的生成树 的个数很多,用τ (G) 表示G 的生成树的个数,则有公式
/ J: \' ?( V2 v& a$ F( v- D" X7 l E
- \+ d- ~- a u ~
/ X H K: q1 k3 e, g% O: s1 E
树有下面常用的五个充要条件。
1 G! m+ u) q5 w! V* k( d
9 ~; z% S3 @4 h, l7 l8 x
定理 1 (i)G 是树当且仅当G 中任二顶点之间有且仅有一条轨道。
2 c% D& B! z7 b% I. ~, B1 l9 ~
; g9 j! R5 }. _0 Z; t
(ii)G 是树当且仅当G 无圈,且ε =ν −1。
2 P: {( k& D) s: D$ O3 C! z; q. l
8 L/ Q/ l- R, Q _0 y0 Q# i
(iii)G 是树当且仅当G 连通,且ε =ν −1。
0 R- |4 z6 L E( J$ e" E
5 V7 k- w( a* l4 f6 I5 Z
(iv)G 是树当且仅当G 连通,且∀e∈ E(G) ,G − e 不连通。
2 Z7 l _4 l+ [
3 C) U( j# K/ @3 R' h; K
(v)G 是树当且仅当G 无圈,∀e∉ E(G) ,G + e 恰有一个圈。
, x9 m8 X* g* L2 T3 U
9 z3 Z8 Y: f; b
2 应用—连线问题
( U! a/ U- R! m
欲修筑连接 n 个城市的铁路,已知i 城与 j 城之间的铁路造价为Cij ,设计一个线 路图,使总造价最低。
! _' g1 p$ [2 w
$ E9 J4 a+ z r& @/ l* d
连线问题的数学模型是在连通赋权图上求权最小的生成树。赋权图的具最小权的生 成树叫做最小生成树。
4 U4 T V7 V+ B3 B
" R# r p& W) k+ T6 v7 F) i* w2 w! j
下面介绍构造最小生成树的两种常用算法。
7 T& e7 F- `* {: K
6 p; N! ^" f( L4 }5 @0 p7 Y
2.1 prim 算法构造最小生成树
$ ` D" a$ s& T& G! M8 i
设置两个集合 P 和Q ,其中 P 用于存放G 的最小生成树中的顶点,集合Q 存放G的最小生成树中的边。令集合 P 的初值为 { P = v1 }(假设构造最小生成树时,从顶点 v1 出发),集合Q 的初值为Q = Φ 。
8 V) N. k4 o: I) W, }' m
5 K1 L5 Z" [8 j! }: ?% N
prim 算法的思想是,从所有 p ∈ P ,v ∈V − P 的边 中,选取具有最小权值的边 pv ,将顶点 v 加入集合 P 中,将边 pv 加入集合Q 中,如 此不断重复,直到 P =V 时,最小生成树构造完毕,这时集合Q 中包含了最小生成树 的所有边。
$ s8 Z) W8 n1 Y `4 E2 w' J
/ ^- ^ j9 u& m4 m: q2 ^. O' j( K& N
! Z' g/ r1 K o6 N$ C1 D
+ R' f3 _0 a" f1 Q
例 13 用 prim 算法求图 5 的最小生成树。 我们用 的第一、二、三行分别表示生成树边的起点、终点、权集合。Matlab 程序如下:
- X/ ^2 _( ^5 v4 s/ ~
+ \ y" H+ c1 r4 E8 \
! t5 G5 U4 n% |# s, _3 i$ o
clc;clear;
: [8 _) F/ b7 s" Y) f) `; q* o" v' X
a=zeros(7);
/ p8 E' X# j5 M+ j" P
a(1,2)=50; a(1,3)=60;
/ L+ |; T+ g' S7 |
a(2,4)=65; a(2,5)=40;
$ |# y6 n/ H" ~; |) t1 o
a(3,4)=52;a(3,7)=45;
+ O. W4 x3 K! P b' F( `9 X2 f
a(4,5)=50; a(4,6)=30;a(4,7)=42;
4 q1 e0 t, `6 {
a(5,6)=70;
0 E I/ M8 U ^) ?8 q2 f* C
a=a+a';a(find(a==0))=inf;
. H1 ]; B3 k6 t F
result=[];p=1;tb=2:length(a);
$ D6 g4 _' \+ R0 e0 I; t ^- R" d
while length(result)~=length(a)-1
0 O$ T& Q' p" N/ a) c: R! A+ W
temp=a(p,tb);temp=temp(
;
! ~8 y) g0 o2 n
d=min(temp);
8 }9 { E# \) L5 L
[jb,kb]=find(a(p,tb)==d);
, F# O; C3 a$ K$ D& T* e
j=p(jb(1));k=tb(kb(1));
/ j4 |1 Y8 j5 w% D
result=[result,[j;k;d]];p=[p,k];tb(find(tb==k))=[];
4 B; C8 l' X' I" g( A' @
end
. G2 i5 c3 y, |# L7 c) I6 {3 A
result
6 c, z. |* `+ X" v3 k7 G
' i p1 A/ G( Z3 n
2.1 Kruskal 算法构造最小生成树
. r1 f8 h* i5 \2 ~
科茹斯克尔(Kruskal)算法是一个好算法。Kruskal 算法如下:
6 J! G3 ^; I( K' h. i
% g# n- R- F5 Z$ b8 N
, ] L' a$ k$ y8 E" P* g
/ g" z# y A" H( g
例 14 用 Kruskal 算法构造例 3 的最小生成树。 我们用 存放各边端点的信息,当选中某一边之后,就将此边对应的顶点序 号中较大序号u 改记为此边的另一序号v ,同时把后面边中所有序号为u 的改记为v 。
5 ]. p" P7 Z/ a8 @
8 I7 X E' r' s0 V
此方法的几何意义是:将序号u 的这个顶点收缩到v 顶点,u 顶点不复存在。后面继续 寻查时,发现某边的两个顶点序号相同时,认为已被收缩掉,失去了被选取的资格。 Matlab 程序如下:
/ ?! f" g8 e. D5 }
& z9 m2 y; P8 v5 H7 g& X2 [' ~/ l
clc;clear;
" w! L4 n8 B( ^
a(1,2)=50; a(1,3)=60; a(2,4)=65; a(2,5)=40;
1 e: x0 W9 V: t% t8 X! I+ W/ i5 G2 a
a(3,4)=52;a(3,7)=45; a(4,5)=50; a(4,6)=30;
* M# V: d6 n. {+ n9 `
a(4,7)=42; a(5,6)=70;
2 j3 Q% m$ q; k% ]
[i,j,b]=find(a);
. v* o& O! t- R" O8 b( Z
data=[i';j';b'];index=data(1:2,
;
1 X* v% a8 Y! e
loop=max(size(a))-1;
' @, G' c* T! J0 H
result=[];
( S9 C. S# B- A# R( D# B
while length(result)<loop
6 V: [% U! t" o. }; y. A
temp=min(data(3,
);
3 N- G4 H( _9 l. g9 b2 W
flag=find(data(3,
==temp);
/ _5 y, u5 Y. C; p& ?
flag=flag(1);
) U+ r ?( _0 M* h) M
v1=data(1,flag);v2=data(2,flag);
4 w- P" ^$ A# e' a
if index(1,flag)~=index(2,flag)
) f: u0 M- F5 u3 H1 u! {/ M
result=[result,data(:,flag)];
9 W; m) C- ]& F( X/ V' X* B
end
( ]7 S; n+ l+ A
index(find(index==v2))=v1;
$ N+ B# g- d) B4 L
data(:,flag)=[];
- Y7 ~2 ]8 \3 b) B& P" J
index(:,flag)=[];
5 }- u6 C( `) z) O6 O$ |% H
end
3 M6 r& s3 T9 I) s
result
7 r3 c" a4 P6 E0 H, S
0 D; ]: k+ P9 [% B N5 r% ^4 ?! t
4 Q6 F8 u$ y2 T; @6 M. ~; z
/ F5 f& m2 }% Y l! R
————————————————
) h! f$ D. c' ]* _) x' `3 R1 S, Z7 w
版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
s7 {) A% s9 y/ d# N! B9 J
原文链接:https://blog.csdn.net/qq_29831163/article/details/89785681
) k" r* @* }& }, N
5 ^6 z) E" I' o
; g; X/ Y; L' ~3 [& R# m3 f2 J
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5