tem.m

2005年电力部数学建模大赛题目： 共有61个比赛项目

clear;
imp=importdata('sj.txt');
[row,col]=size(imp);
rel=zeros(col);
for i=1:col
    for k=i+1:col
        for j=1:row
            if imp(j,i)==1 & imp(j,k)==1
                rel(i,k)=rel(i,k)+1;
                rel(k,i)=rel(i,k);
            end
        end
    end
end
x=zeros(1,col);y=zeros(1,col);t=zeros(col,col);ny=zeros(1,col);cy=zeros(1,col);count=zeros(1,col);result=zeros(1,col)
for i=1:col
    y(i)=1;
end
%****************************************
x(1)=1;y(1)=0;%cy(1)=100;
for i=2:col
    ny(i)=1;cy(i)=rel(1,i);%ny表示与
end
for k=1:col-1
    for i=2:col
        if y(i)==1 small=cy(i);
            break;
        end
        
    end
    for j=1+i:col
        if y(j)==1 & cy(j)<cy(i)
           small=cy(j);
           i=j;
        end
    end
    temp=i;

    t(temp,ny(temp))=1;
    count(temp)=count(temp)+1;
    count(ny(temp))=count(ny(temp))+1;
    x(temp)=1;
    y(temp)=0;
    for i=1:col
        if y(i)==1
            w=i;
            if rel(temp,w)<=cy(w)
                ny(w)=temp;
                cy(w)=rel(temp,w);
            elseif rel(temp,w)>cy(w) & count(ny(w))>=2
                    for i=1:col
                       if x(i)==1 & count(i)<2
                          ny(w)=i;
                          cy(w)=rel(i,w);
                          break;
                       end
                    end
                   for j=1+i:col
                       if x(j)==1 & count(j)<2 & cy(j)<cy(i)
                          ny(w)=j;
                          cy(w)=rel(j,w);
                          i=j;
                       end
                   end
               end
       
 
            end
        end
    end
end
%********************************************
index=1;
for i=1:col
    if count(i)==1 
       result(index)=i;
      
       break
   end
end
for i=1:col
%while index<=col
  for j=1:col
    if t(result(index),j)==1
       result(index)=j;
       t(result(index),j)=10;
       index=index+1;
       %elseif t(j,result(index))==1
          
       %   result(index)=j;
        %  t(j,result(index))=10
         % index=index+1;
    end
  end
end
disp(result);