通用代码- function [result]=myprim(a)//a为传入的每个点的距离矩阵
- S4 M2 K; c( J7 W1 Q - result=[];//用result(3×n)矩阵来表示,第一行表示起点,第二行表示终点,第三行表示权值
& J4 |2 `; K t6 S) W( L, P - p=1;tb=2:length(a);
0 g- V3 N4 L0 i - while size(result,2)~=length(a)-1+ V. \8 Z$ E2 v. U/ Y* F
- temp=a(p,tb);temp=temp(:);6 l# Q3 B% Y' A6 Y
- d=min(temp);2 G3 w2 i7 q4 w
- [jb,kb]=find(a(p,tb) == d,1);7 }# t9 b* s' e\" y
- j=p(jb);k=tb(kb);
1 f$ F& J5 |7 Q4 Z6 b3 \1 b - result=[result,[j;k;d]];p=[p,k];tb(find(tb==k))=[];
# ^0 R$ b6 J% l% \. C - end) Z8 E1 s' Z. @\" ^. _
复制代码 怎么看result我怕有些同学想用这个代码,看不懂result里面元素什么意思,在此提示
4 V+ R; {, x1 F) [6 f- u/ R' Fresult第一行表示起点,第二行表示终点,第三行表示权值
+ | `6 l9 N& C4 ^) d4 i3 `# z例如9 l! u1 T: I* `! R& z/ E H" M
result = - 1 2 5 4 4 7) \2 _2 x# E! s! Q) Z) l3 U) o/ l
- 2 5 4 6 7 3 x. r; w3 s2 g8 a0 L) F( |
- 50 40 10 30 42 45/ F; f0 H- T9 s' N
- 起点就是1,终点是2距离50
+ s) Y5 M, w6 Q; q/ ^ - 然后从2,到5距离40
2 W$ n0 O! c* C0 i/ O- s2 O - 从5到4距离10
& _$ k$ H5 D2 O' \ - 从4到6距离302 ?3 R\" H# N% ^3 w$ X6 M
- 从4到7距离42; {) m. a* { c+ ^8 {
- 从7到3距离45
: |6 q- ^' e2 `/ `& Y - 最后生成完一整棵树+ ^4 }& n0 D# O9 M/ t6 g
复制代码
6 `9 F) P( W; Y; n1 a |