实际问题引入
7 j( `5 I( y/ f" X
0 H( e7 f' K6 g) C( m
实这道题的答案,其实就是找到这个图的最小生成树。' h% F1 X& _# ^* Q- P1 e* _
% G) |) ^6 O8 u8 a* Q6 v) M) \
Kruskal算法
, D5 |; D! \- q5 O2 h Y; [, `此算法可以称为“加边法”,初始最小生成树的边数为 0,每迭代一次就选择一条满足条件的最小代价的边,加入到最小生成树边的集合里面。
& E0 x4 M$ I" H& H: [! n0 |% E+ \其实核心思想就是贪心思想:通过局部最优达到整体最优
, O1 l% O( o- v( b% {3 X6 j, r
, c0 A3 f6 R; ]6 G: Q将所有的边权进行排序
) m" i; v5 S9 K/ Z/ d, r不断迭代选择权最小的边,直到所有的点被连起来(边数=节点数-1)。. [ c. f3 d. W. Q; r" ~( a- G
在迭代期间,如果边构成了环,就要丢弃该边,因为树中是不存在环的!8 P3 n) ~2 ?2 s) |0 F
整体代码展示$ C$ r1 J- I2 ]1 P8 }
在matlab中,最小生成树的生成直接用minspantree()函数就行。- s=[1,1,1,1,2,2,3,3,4,4,5,5,6];5 k5 B7 `7 ~9 m7 ]5 o
- t=[2,3,4,5,3,6,5,7,5,6,6,7,7];
8 y8 b* [4 T* F: u- N - w=[35,24,10,25,25,20,15,11,12,30,15,25,18];
, R, R# d, G- ~9 X- O# ], ^ - names={'1','2','3','4','5','6','7'};, v4 n( @' a4 N9 D, v' a
- G=graph(s,t,w,names);
( D# y4 @) {) f' }& I& V/ K - p=plot(G,"EdgeLabel",G.Edges.Weight);
! _7 A; X3 V c1 v; f\" ]- F: L - % 求解最小生成树; K- Z# p. V' L: Z( h$ e
- T=minspantree(G,"Method","sparse");9 ~\" K. {- ?# y0 y) |( [
- % sparse代表的是Kruskal算法
4 y$ O8 `; {/ ]' O/ j2 |* ~0 ~* f - % dense代表的是Prim算法2 C! a# L- _' l, ~# P
. x4 ]0 U7 T! J5 z- y$ D2 ~$ V- % sparse:Kruskal算法
' D# g: I {0 d* x - % 算法按权重对所有的边排序,然后将不构成循环的边添加到树中
7 U& {1 Y& d, j6 w2 r1 X - p=plot(G,"EdgeLabel",G.Edges.Weight);9 |0 j# ]. N6 }. U; ?
- highlight(p,T,"NodeColor","red","EdgeColor","red"); * v* |! b9 F: i p9 \
- % 将最小生成树的边设置为红色!5 O# }$ \6 u' p7 ^
复制代码
4 |- D+ o. ?* V. b9 \+ o. @( ?9 F6 J+ z. h9 R4 _
生成的最小生成树:. h( R- r$ ]9 J
. J2 ~' V5 F; N% l. W1 m3 N
我们也可以把最小生成树的边和节点打印出来,也可以把整段路的权加起来看看:
G5 f1 \# z$ p: u/ J( y
尾声看到这里,相信我们已经学会Kruskal算法寻找最小生成树的过程了,当然,这离数学建模的要求,离我们的目标还非常遥远,博主在不断学习的过程中,也希望可以通过分享学习日记的方式带动大家! # o2 ?: a5 D( x( c1 \+ Y- H" ~" i
1 H. k" e1 f" Q5 k |