lock5.m

全国大学生数学建模竞赛 锁具装箱程序

for i=1:5880
    c2(i,:)=c(i,:);
end
r1=zeros(1,500);
for i=1:500
    num3=0;temp2=5880;
    for j=1:120
        temp2=temp2-1;
        temp1=fix(rand(1)*temp2)+1;    
        c3(j,:)=c2(temp1,:);
        for k=1:temp1-1
            c2(k,:)=c(k,:);
        end
        for k=temp1:temp2
            c2(k,:)=c(k+1,:);
        end

    end
    for j=1:120-1
        for k=j:120
            num2=0;flag1=0;
            for l=1:5
                if c3(j,l)==c3(k,l) num2=num2+1;end
            end
            if num2==4 flag1=1;end
            if (abs(sum(c3(j,:),2)-sum(c3(k,:),2))==1)&&(flag1==1) num3=num3+1;end
        end
    end
    r1(num3+1)=r1(num3+1)+1;    
end