通用代码- function [result]=myprim(a)//a为传入的每个点的距离矩阵
4 S( M! f; F' D4 ^0 P8 r - result=[];//用result(3×n)矩阵来表示,第一行表示起点,第二行表示终点,第三行表示权值/ u: v) U9 v- ]1 o) {4 R/ b4 E& `
- p=1;tb=2:length(a);
) B# A G, z\" W& Y' I - while size(result,2)~=length(a)-12 [) F+ w1 ~3 J
- temp=a(p,tb);temp=temp(:);
* y- s. O! a( _; p\" u$ ~' E$ f - d=min(temp);
& h4 @4 s/ l0 l5 S( q - [jb,kb]=find(a(p,tb) == d,1);6 k. C. g\" W* d. p8 d% N A
- j=p(jb);k=tb(kb);3 g4 H6 Z/ q! ]9 J5 k: c. n
- result=[result,[j;k;d]];p=[p,k];tb(find(tb==k))=[];6 P( |0 D8 z% X8 }8 e# Z- {1 ^
- end
2 W7 F1 w$ s- f3 n+ n3 y% o
复制代码 怎么看result我怕有些同学想用这个代码,看不懂result里面元素什么意思,在此提示
4 U3 }* E0 l8 c2 Cresult第一行表示起点,第二行表示终点,第三行表示权值
0 J5 R& `7 T& C6 @( t% g例如
) Q8 `6 M! d) B' q c6 Bresult = - 1 2 5 4 4 74 { ^9 y$ ]7 J
- 2 5 4 6 7 3
) S, r b9 T& {- u6 X* T( y3 L - 50 40 10 30 42 45
6 `) q# s# [- s/ K% G6 N - 起点就是1,终点是2距离50
& p$ n; }. l+ Y' U5 | - 然后从2,到5距离40
[1 s. i5 R; W, U' I3 a4 ]3 ] - 从5到4距离10
. Q! Q- a& ]# {/ s( P/ N3 w - 从4到6距离30 p3 e1 ?. y% J$ I% ~$ W
- 从4到7距离42
9 J- M' a2 x0 ]1 R - 从7到3距离45
' u; d7 l: u* X' X8 v - 最后生成完一整棵树1 g9 @( M' h% `\" {4 I
复制代码 " S `" `, K- i. `) \0 a3 G
|