通用代码- function [result]=myprim(a)//a为传入的每个点的距离矩阵+ m) y8 J) `3 L2 B F! P2 A% |( u
- result=[];//用result(3×n)矩阵来表示,第一行表示起点,第二行表示终点,第三行表示权值- D7 i$ D' v& p8 H; v; n, q8 i
- p=1;tb=2:length(a);
: f\" n( @( F1 k - while size(result,2)~=length(a)-16 S( n7 w0 L% W G) b- x
- temp=a(p,tb);temp=temp(:);6 H, p& r7 s* R ?: S: ?
- d=min(temp);* @6 l0 V- f) l! P: g0 D
- [jb,kb]=find(a(p,tb) == d,1);
: [7 ]4 C0 o- c; m( X5 g. h5 u0 K+ S - j=p(jb);k=tb(kb);
3 p3 c7 p+ P [\" Z. S - result=[result,[j;k;d]];p=[p,k];tb(find(tb==k))=[];- C6 w1 @, }& X3 C
- end
# I; u& l5 q; e0 \- W, X7 p. E
复制代码 怎么看result我怕有些同学想用这个代码,看不懂result里面元素什么意思,在此提示
( s' u0 A- _& A7 z& {" ?$ gresult第一行表示起点,第二行表示终点,第三行表示权值
- S/ }; n5 d( x例如- }! @( V7 \# F2 F4 J* S6 C# E
result = - 1 2 5 4 4 7
) O$ M9 g% E7 } - 2 5 4 6 7 3
# p6 k0 D9 x' y: w% z - 50 40 10 30 42 45) y x\" P9 I* `: w
- 起点就是1,终点是2距离50
# y V5 z) ~' I7 n6 O - 然后从2,到5距离40+ s8 H$ ]- w+ G6 N! b0 H
- 从5到4距离10% t2 V2 j( R\" K7 O3 f6 _
- 从4到6距离30
5 b: H7 o' ]' c! J6 b/ _' x- O - 从4到7距离42
j1 E9 c# Q, D% V, N - 从7到3距离45. R$ G7 b% o. H g
- 最后生成完一整棵树' ?2 j( Y6 Q( U2 @, d4 `& s
复制代码
7 e* K1 p' Q3 Q K) L4 C |