通用代码- function [result]=myprim(a)//a为传入的每个点的距离矩阵* ]4 q% G; W7 ~. M# V6 z% p( `
- result=[];//用result(3×n)矩阵来表示,第一行表示起点,第二行表示终点,第三行表示权值
% `: Y, c' ^( \4 \* D5 F - p=1;tb=2:length(a);% i% ]/ e* g& s0 p2 r
- while size(result,2)~=length(a)-1
\" i) i, {) y/ Z7 ]/ g. {% M s/ e: b - temp=a(p,tb);temp=temp(:);
) {0 [8 z% N( r/ b' F* \( A5 z - d=min(temp);
' e' [7 v# T% v' {3 a5 U - [jb,kb]=find(a(p,tb) == d,1);$ Z) w. T3 V- k0 t( g& a$ k4 K
- j=p(jb);k=tb(kb);7 r( n' O& j4 A# ?4 c; d/ q
- result=[result,[j;k;d]];p=[p,k];tb(find(tb==k))=[];
; D) u) z) N( g1 s* M1 c - end& j0 n5 `0 k( I/ X
复制代码 怎么看result我怕有些同学想用这个代码,看不懂result里面元素什么意思,在此提示, @8 g$ W9 _. T' K3 }+ [
result第一行表示起点,第二行表示终点,第三行表示权值
. N* o& b1 y: X5 _) ?( w例如
5 \6 L. D7 |4 T1 Y1 ]1 ?" I: tresult = - 1 2 5 4 4 7
8 z$ f8 g0 f9 @2 Y& R: p1 n - 2 5 4 6 7 3
; M4 l5 K+ x& g- b- S2 b4 ^ - 50 40 10 30 42 459 C1 _\" E* ~& q+ e7 o) M# N
- 起点就是1,终点是2距离50
! e( _1 n% `$ L2 h: Z# V - 然后从2,到5距离40
! q) H) a$ ^: y. N5 ~$ | - 从5到4距离10
4 S& ~: z8 \( u: @! m. [ - 从4到6距离30
1 R% _: k5 o5 T- J3 l - 从4到7距离424 m% S3 @+ y$ G7 X/ ?; m\" n9 H- f
- 从7到3距离45
+ m+ G0 D/ S4 L5 I8 z - 最后生成完一整棵树 \/ ?\" {8 Z; }
复制代码
& D L6 H9 `/ P" j+ ?8 G |