通用代码- function [result]=myprim(a)//a为传入的每个点的距离矩阵
4 A, O A, f/ a) \7 O# |4 A: x - result=[];//用result(3×n)矩阵来表示,第一行表示起点,第二行表示终点,第三行表示权值
* y- D a* |% t% \, Q S& Y - p=1;tb=2:length(a);
\" y* g q\" b6 Y) N& G' s* M - while size(result,2)~=length(a)-1
; p Q- j# {. q5 W0 k! k - temp=a(p,tb);temp=temp(:);
. Q' \' r3 n3 f, u - d=min(temp);8 ^1 B3 ~4 k3 Q: A S% {! H% Y
- [jb,kb]=find(a(p,tb) == d,1);
9 q0 F+ q& k# L0 |$ Z8 u4 o3 k6 l2 h! V - j=p(jb);k=tb(kb);
* q5 ?: A' v3 {' m9 f - result=[result,[j;k;d]];p=[p,k];tb(find(tb==k))=[];\" V: z4 L9 G }# }9 F H
- end
/ G; c1 X! l; ]
复制代码 怎么看result我怕有些同学想用这个代码,看不懂result里面元素什么意思,在此提示
) G1 s' ?- l9 L! c3 I8 Uresult第一行表示起点,第二行表示终点,第三行表示权值
2 v; E9 Y: f! N' f3 D& Y1 G/ k0 B# l, K例如
2 a9 M5 b I6 @2 D/ V( N* Yresult = - 1 2 5 4 4 76 m; ]' U0 X- s; |
- 2 5 4 6 7 3
- Q7 p9 L! u4 ], a - 50 40 10 30 42 45
% r1 r* ?, T) N# k: U& A/ v, ? - 起点就是1,终点是2距离50
# Q: \' f: @0 L\" G A) \ - 然后从2,到5距离40
9 f1 J\" J! S, p2 g8 k. J. w- t' n - 从5到4距离10; v\" A, f3 Z$ c# \/ Q( C3 N
- 从4到6距离30
# Q9 q. F: q, F* o5 c( b - 从4到7距离42
, u9 q7 Y7 z! | h- v) t - 从7到3距离452 S2 u3 S/ P) y1 V5 N9 `; ~3 C# F
- 最后生成完一整棵树
0 k( Z+ v7 h8 C P\" t: u
复制代码
1 h2 u+ I4 s% Y |