通用代码- function [result]=myprim(a)//a为传入的每个点的距离矩阵1 [4 N\" k* k7 {8 h
- result=[];//用result(3×n)矩阵来表示,第一行表示起点,第二行表示终点,第三行表示权值
$ E; a( |+ k: F7 o0 F7 y - p=1;tb=2:length(a);
9 G! r2 b\" l7 F - while size(result,2)~=length(a)-1% p: Y& f8 r* \; c8 v6 J
- temp=a(p,tb);temp=temp(:);
4 a9 f\" |! A* r! X0 i1 a - d=min(temp);
6 y/ F9 r3 s# K. v9 V5 v, O/ C - [jb,kb]=find(a(p,tb) == d,1);5 d, ^3 o4 O3 j1 Q3 X3 m
- j=p(jb);k=tb(kb);
1 S\" ] U) R1 h) F - result=[result,[j;k;d]];p=[p,k];tb(find(tb==k))=[];
, n4 x6 B$ c1 { - end
: C6 T7 A% x j/ m7 p9 f& c2 j0 G: `! L
复制代码 怎么看result我怕有些同学想用这个代码,看不懂result里面元素什么意思,在此提示+ ?; H" |* N3 O0 e+ l8 Z9 c
result第一行表示起点,第二行表示终点,第三行表示权值
) N( Q2 A1 g3 x- y& ?: l例如
7 P: o* |4 V' k0 {+ m/ r* |+ l- Tresult = - 1 2 5 4 4 7$ _9 a( Z D8 W8 D! ^% R) S
- 2 5 4 6 7 30 d* w1 E. H8 R/ G c
- 50 40 10 30 42 45
; u& g7 j* ?; q- J4 Z - 起点就是1,终点是2距离50( E$ J5 u/ D. B0 @+ s1 \+ W
- 然后从2,到5距离409 T8 m) U8 e, I r% o
- 从5到4距离10
5 y: z& L% L5 ?0 y\" [/ a* C - 从4到6距离30 e V% g* X) ], N w8 q0 k6 v
- 从4到7距离420 V5 g. d( K% B1 o C3 i
- 从7到3距离458 l. e0 L O2 C- k4 D) \
- 最后生成完一整棵树
2 c& \7 I9 T1 s# e( E! d- j' B& n8 t
复制代码 ! a" i% `+ G) y( ^; O
|