较优H圈matlab实现

大家好，我是小鸭酱，博客地址为：http://www.cnblogs.com/xiaoyajiang

%解决完备图中的较优H圈

clc

clear

w = [ inf 6 1 8 3 1 ;...        %权赋矩阵

6 inf 1 1 6 9 ;...

1 1 inf 3 6 7 ;...

8 1 3 inf 1 2 ;...

3 6 6 1 inf 1 ;...

1 9 7 2 1 inf ] ;

v = [ 11 12 13 14 15 16 ;...   %两点之间的路径，如25表示从点2到点5

21 22 23 24 25 26 ;...

31 32 33 34 35 36 ;...

41 42 43 44 45 46 ;...

51 52 53 54 55 56 ;...

61 62 63 64 65 66 ] ;

p = [ 12 23 34 45 56 61 ]    %保存某个H圈的每两点间的路径，H圈便是由这些小的路径构成的

W = zeros(1,6) ;

S=0;

for h = 1 : 6

d= rem(p(h),10) ;

c= (p(h) - d)/10 ;

W(h) = w(c,d) ;

S= S + W(h) ;

end

S

for i=1:1000

x =  randperm(6) ;

m = x(1) ;

n = x(2) ;

j = rem(p(m),10) ;

i = (p(m) - j)/10 ;

z = rem(p(n),10) ;

y = (p(n) - z)/10 ;

if ( w(i,j) + w(y,z) > w(i,y) + w(j,z) )

s = 0;

for b = 1 : 6

if ( p(b) == v(i,y) | p(b) == v(j,z) | p(b) == v(y,i)| p(b) == v(z,j) )

s = 1 ;

break ;

end

end

if s

i = i + 1 ;

break ;

end

p(m) = v(i,y) ;

p(n) = v(j,z) ;

end

end

W = zeros(1,6) ;

S=0;

for h = 1 : 6

d= rem(p(h),10) ;

c= (p(h) - d)/10 ;

W(h) = w(c,d) ;

S= S + W(h) ;

end

p

S