通用代码- function [result]=myprim(a)//a为传入的每个点的距离矩阵- [( }% A* Y$ e& \! y
- result=[];//用result(3×n)矩阵来表示,第一行表示起点,第二行表示终点,第三行表示权值
' c9 `: m* p9 p% v! D9 v7 a7 P - p=1;tb=2:length(a);# a2 o/ A6 [3 V: o8 s8 s+ X T
- while size(result,2)~=length(a)-19 n8 D7 r, Z$ m! ]
- temp=a(p,tb);temp=temp(:);
\" R0 ]6 N0 k$ z d\" J7 u - d=min(temp);
% N- f; k( E& `+ X\" J! w* b - [jb,kb]=find(a(p,tb) == d,1);
3 d: P( f7 ]$ O( y4 ?9 k - j=p(jb);k=tb(kb);
% v! T5 J* \( E; W' X; B - result=[result,[j;k;d]];p=[p,k];tb(find(tb==k))=[];9 |# V2 R( W0 ]5 y, g: x
- end9 h( C3 ^5 F6 b# P( E
复制代码 怎么看result我怕有些同学想用这个代码,看不懂result里面元素什么意思,在此提示
" A/ ^. u- q& E9 h3 z' xresult第一行表示起点,第二行表示终点,第三行表示权值
9 s- z/ E L/ D6 _例如5 i: x6 c+ d. f
result = - 1 2 5 4 4 7
# Y3 r- t& @$ z- X, N& [2 Z - 2 5 4 6 7 31 ]: z$ q: F- P1 ~2 Q) q% r
- 50 40 10 30 42 45# ]: Z; h: V7 P# s4 a
- 起点就是1,终点是2距离50
% w5 d8 ^( @! t - 然后从2,到5距离40. ^% C, g4 R, K+ Y9 s# d4 Q! X7 w
- 从5到4距离108 N/ v- x: K4 h3 W9 u
- 从4到6距离30+ B3 o: V- w4 p f
- 从4到7距离42
* N ?1 @6 s' W/ z5 R1 V - 从7到3距离45
N% z z/ q$ ~# S8 J6 ]' k - 最后生成完一整棵树
) \8 _& D! [: t. w' g/ e* ?6 b
复制代码
* M6 w: |3 z r. O( @5 z4 Z |