通用代码- function [result]=myprim(a)//a为传入的每个点的距离矩阵
$ E/ k0 F, U; n: ~) F( [ - result=[];//用result(3×n)矩阵来表示,第一行表示起点,第二行表示终点,第三行表示权值) x @ b7 x6 \0 X, ~\" A
- p=1;tb=2:length(a); Z2 v. S( Q: f8 w5 z7 a3 a9 c
- while size(result,2)~=length(a)-11 @9 A; M+ ]$ [\" C
- temp=a(p,tb);temp=temp(:);
# H, |# a1 \+ a$ ?1 l - d=min(temp);0 |6 R g* n' T
- [jb,kb]=find(a(p,tb) == d,1);9 F& G, Q2 g# q2 }& ^
- j=p(jb);k=tb(kb);, S8 f: E9 l( e\" f+ r$ g3 F }
- result=[result,[j;k;d]];p=[p,k];tb(find(tb==k))=[];* |, X% I+ h0 y* e9 |3 l& Y
- end2 Z0 q$ ]. \3 V' T$ }
复制代码 怎么看result我怕有些同学想用这个代码,看不懂result里面元素什么意思,在此提示
! w6 K1 Z# I) Y, T" q0 Iresult第一行表示起点,第二行表示终点,第三行表示权值
8 D( w! h; e6 E" S* Y6 H例如. l0 p0 B$ H+ @, `# z3 t" ~
result = - 1 2 5 4 4 7
& u/ A+ a8 }$ V, _ - 2 5 4 6 7 3
' f6 D: ]& D. ~+ T\" M$ ?2 Q - 50 40 10 30 42 455 U+ D. }7 K, ^+ v9 ^7 d, Q8 V& Q
- 起点就是1,终点是2距离50# q' _* u( e/ i
- 然后从2,到5距离40
* f) w$ K* l, p: Y* S/ w - 从5到4距离10. q% A0 @: N; z3 c) F5 J2 \) K0 N
- 从4到6距离30- T8 m0 n0 k( l7 E. p; e; ~
- 从4到7距离42! r {! \4 Z: k3 e* H ^
- 从7到3距离453 a y! m# V- B& F+ k# f7 r
- 最后生成完一整棵树
* g; o( u. S/ D
复制代码 . r3 C7 M+ T0 M% |% e! `- {
|