|
|
1.1 ! root 1: C = < ! 2: 1 2 15 . . . ! 3: 2 1 3 . . . ! 4: 3 2 4 11 . . ! 5: 4 3 5 . . . ! 6: 5 4 6 7 . . ! 7: 6 5 8 . . . ! 8: 7 5 9 30 . . ! 9: 8 6 9 10 11 . ! 10: 9 7 8 30 . . ! 11: 10 8 12 30 31 34 ! 12: 11 3 8 12 13 . ! 13: 12 10 11 34 36 . ! 14: 13 11 14 . . . ! 15: 14 13 15 16 38 . ! 16: 15 1 14 . . . ! 17: 16 14 17 20 35 37 ! 18: 17 16 18 . . . ! 19: 18 17 19 . . . ! 20: 19 18 20 . . . ! 21: 20 16 19 21 . . ! 22: 21 20 22 . . . ! 23: 22 21 23 . . . ! 24: 23 22 24 35 . . ! 25: 24 23 25 39 . . ! 26: 25 24 . . . . ! 27: 26 27 33 39 . . ! 28: 27 26 32 . . . ! 29: 28 29 32 . . . ! 30: 29 28 30 . . . ! 31: 30 7 9 10 29 . ! 32: 31 10 32 . . . ! 33: 32 27 28 31 34 . ! 34: 33 26 34 . . . ! 35: 34 10 12 32 33 35 ! 36: 35 16 23 34 36 . ! 37: 36 12 35 38 . . ! 38: 37 16 38 . . . ! 39: 38 14 36 37 . . ! 40: 39 24 26 . . . ! 41: >; ! 42: ! 43: <n, m> = size(C); ! 44: A = 0*ones(n,n); ! 45: for i=1:n, for j=2:m, k=c(i,j); if k>0, a(i,k)=1; ! 46: check = norm(A-A',1), if check > 0, return ! 47: <X,D> = eig(A+eye); ! 48: D = diag(D); D = D(n:-1:1) ! 49: X = X(:,n:-1:1); ! 50: <x(:,1)/sum(x(:,1)) x(:,2) x(:,3) x(:,19)>
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.