通用代码- function [result]=myprim(a)//a为传入的每个点的距离矩阵 J: i) @+ F) n6 c
- result=[];//用result(3×n)矩阵来表示,第一行表示起点,第二行表示终点,第三行表示权值
1 l\" i1 J. u/ o - p=1;tb=2:length(a);
3 D$ F5 I\" m, k6 b! i& n - while size(result,2)~=length(a)-1
$ r! J2 a/ ?: }\" h- R) V - temp=a(p,tb);temp=temp(:);6 b$ U9 ]1 n c y, R5 G3 \
- d=min(temp);
7 w0 E. N; z( }, D( ^9 N/ v - [jb,kb]=find(a(p,tb) == d,1);
6 p2 I/ }, }7 W5 ?. H - j=p(jb);k=tb(kb);* V: P; D2 h+ [' u; Y/ d) X
- result=[result,[j;k;d]];p=[p,k];tb(find(tb==k))=[];# q4 Y\" E! Y* N. Q
- end2 W+ e9 k- W6 T3 C( U
复制代码 怎么看result我怕有些同学想用这个代码,看不懂result里面元素什么意思,在此提示
9 c$ h& {* A7 iresult第一行表示起点,第二行表示终点,第三行表示权值
' l3 X" k. I+ m例如- {3 p0 f/ }5 x8 V2 A; n4 A/ k9 Y5 y
result = - 1 2 5 4 4 7
\" I/ s- L. I\" I( W: l! `5 n - 2 5 4 6 7 3
! v6 t* Q0 D( E* `, m/ w - 50 40 10 30 42 45
& W- V\" r5 [ g# d# ~ - 起点就是1,终点是2距离50
8 l5 ]6 D: O5 l. H! j m& } - 然后从2,到5距离40* U. ]8 F1 Q) L8 H! H/ G. T. j6 |- V
- 从5到4距离10
\" C; _1 g, T3 F8 d) R. ] - 从4到6距离30
6 Z# g/ [1 P* L; G - 从4到7距离42% Y( z D: `8 `
- 从7到3距离45* c5 _' f, {; p' R1 X/ E; ]
- 最后生成完一整棵树 G; Q/ a3 a0 k2 k2 A
复制代码
+ X% W+ w; D3 [/ W) C |