通用代码- function [result]=myprim(a)//a为传入的每个点的距离矩阵) Z* t1 ~8 D' }& O, L
- result=[];//用result(3×n)矩阵来表示,第一行表示起点,第二行表示终点,第三行表示权值- h! t) w- k+ v; x3 F, }# J
- p=1;tb=2:length(a);1 N- l- {% A7 M3 Q- z8 q! |: L# o
- while size(result,2)~=length(a)-1
) W' r5 [8 @ s) V - temp=a(p,tb);temp=temp(:);
2 R, g0 o6 M/ R7 ` - d=min(temp);
9 I* h% Q; K3 g3 z$ N( r& [ - [jb,kb]=find(a(p,tb) == d,1);& {5 u* c& n\" P7 P
- j=p(jb);k=tb(kb);
/ J/ |- F\" A& H# b4 a* g - result=[result,[j;k;d]];p=[p,k];tb(find(tb==k))=[];
2 c5 w% n; o4 {\" v4 g: g3 b: f3 G - end; x L0 z+ N8 X% O0 R
复制代码 怎么看result我怕有些同学想用这个代码,看不懂result里面元素什么意思,在此提示# d _8 E+ F) J* I$ a( X# R1 U
result第一行表示起点,第二行表示终点,第三行表示权值) |1 s3 J1 w5 |7 l& t
例如% f+ }* o4 G. j' z+ k j4 Q5 X8 t
result = - 1 2 5 4 4 72 ~& y4 e) m' C: O
- 2 5 4 6 7 3
/ ?* E# P\" _3 k. K/ T' j - 50 40 10 30 42 45
/ R/ k/ q2 Z3 G) U8 W - 起点就是1,终点是2距离50
\" X/ N0 { ?. f8 [ - 然后从2,到5距离408 o! R# ]- Z( ?3 b$ E% j
- 从5到4距离10
4 N. E$ Z: O# j% W - 从4到6距离30+ a2 u; K o- Y( \4 [2 x% C! ]
- 从4到7距离42# G5 C0 W# h2 J( P3 ]9 H
- 从7到3距离45. e7 \/ `5 d: X8 Q6 e\" H5 ]. o
- 最后生成完一整棵树
, d- g( E. ^, W3 J8 \
复制代码
2 @0 x& a* q& r/ A/ l$ @ |