通用代码- function [result]=myprim(a)//a为传入的每个点的距离矩阵, Q) h0 V$ _; a: @
- result=[];//用result(3×n)矩阵来表示,第一行表示起点,第二行表示终点,第三行表示权值5 U# ^& a- R, L% H; }1 z/ f
- p=1;tb=2:length(a);0 ]2 W% }0 t0 c. ?9 S/ F
- while size(result,2)~=length(a)-1; C6 j- |7 v3 e8 R9 g
- temp=a(p,tb);temp=temp(:);
) N\" `% `& z2 K8 k! Y+ U# ~% h - d=min(temp);
& O5 D/ ~7 j3 P5 W9 D8 o' d. z0 h' x - [jb,kb]=find(a(p,tb) == d,1);, ?! p. t' V9 @
- j=p(jb);k=tb(kb);\" c! T6 G6 q9 [& t, W# X6 F
- result=[result,[j;k;d]];p=[p,k];tb(find(tb==k))=[];* I+ d* ~- W& X- d; ~1 o1 v5 J3 w8 i( s
- end
) e! f3 g( s\" Y) l. h5 C+ O1 S2 v+ w
复制代码 怎么看result我怕有些同学想用这个代码,看不懂result里面元素什么意思,在此提示$ k) [' w& ^' `
result第一行表示起点,第二行表示终点,第三行表示权值
; }3 P* _$ O6 h; f J5 u% p例如
) f9 W3 s8 z8 Y# B7 l1 {# D! jresult = - 1 2 5 4 4 72 C; r2 e. u' C& s
- 2 5 4 6 7 3
6 x1 ?2 w7 g1 { - 50 40 10 30 42 45
4 r* C8 O* T' V$ ?3 i - 起点就是1,终点是2距离50\" p+ T( R/ H* p* @# j$ [0 ]! f
- 然后从2,到5距离40: m: K! q3 ?; q8 \3 |
- 从5到4距离105 I7 u( s\" M4 r1 t7 P
- 从4到6距离30
/ b+ L7 l: P2 D* Z) u/ A& \ i - 从4到7距离42; |7 t; k; [# n3 ~- r
- 从7到3距离45
\" W# F, {$ m1 b$ ^6 @ - 最后生成完一整棵树
. L$ u* q7 A$ U3 I* E$ H5 i( h
复制代码 1 y0 @! |& L H" S5 b m" ]1 d
|