通用代码- function [result]=myprim(a)//a为传入的每个点的距离矩阵
& V5 J& T, N; ` - result=[];//用result(3×n)矩阵来表示,第一行表示起点,第二行表示终点,第三行表示权值4 G4 p% }% u8 d
- p=1;tb=2:length(a);
* A\" R+ V% g. \\" L/ L# d - while size(result,2)~=length(a)-1
0 ^% p- f( z' O$ z. l7 L e) C( D+ e - temp=a(p,tb);temp=temp(:);. r+ ]1 S+ U4 a) i\" s$ U$ O: d/ ~- R' i
- d=min(temp);
/ }1 l# T, U& P* u: S7 G - [jb,kb]=find(a(p,tb) == d,1);$ w2 [% H; w8 M, A5 [% ]+ d
- j=p(jb);k=tb(kb);: _! F) @9 R7 D+ x* m8 j
- result=[result,[j;k;d]];p=[p,k];tb(find(tb==k))=[];
: Q7 O! H1 S1 ?& X9 j - end1 |* x\" j6 P M! j- u$ ]3 m
复制代码 怎么看result我怕有些同学想用这个代码,看不懂result里面元素什么意思,在此提示
2 ]& b7 f0 d3 E8 D+ I: t @7 H {- cresult第一行表示起点,第二行表示终点,第三行表示权值1 ^) T. M& Q+ v9 ^; T
例如
- }6 m' S* n" \result = - 1 2 5 4 4 7; z% U; F1 _' q# W2 U) ~$ ? |8 G2 n
- 2 5 4 6 7 3
6 S7 f* S( ^, C$ w' Z3 E - 50 40 10 30 42 45
, w* y, F& `+ `: ?( A - 起点就是1,终点是2距离50
3 b9 B. r' B. I( a! U2 @ - 然后从2,到5距离402 O, d7 o& Y7 o7 K
- 从5到4距离10
! D0 J\" Y5 F% I8 e\" Y# p) e# i/ I/ d - 从4到6距离30) l+ j0 v, P; I) I& ]: V
- 从4到7距离42
6 P& u# J! @; N5 k6 r, K3 P6 G - 从7到3距离459 Z0 a1 t; x& G' v2 Z7 _
- 最后生成完一整棵树0 f' X1 H0 i3 Y, R: U
复制代码
# r/ F+ _1 C9 C3 B- P( s# h |