zuhe.m

来自「一个十分简短却非常有用的产生组合数的程序。」· M 代码 · 共 34 行

function [y,zhi,yy]=zuhe(B)
N=1;h=1;n=1;k=1;
BHL=size(B);BH=BHL(1);
for h=1:BH
    Bb=B(h,:);BB=Bb(Bb~=0);BL(h)=sum(size(BB))-1;
    N=N*BL(h);y(1,h)=B(h,1);m(h)=1;
end
for n=1:N
for k=1:BH
    if m(k)>BL(k)
        m(k)=1;
        if k~=BH
            m(k+1)=m(k+1)+1;
        end
    end
      y(n,k)=B(k,m(k));
end
     m(1)=m(1)+1;
end
s=sum(y');
[zhi,num]=min(s);
yy=y(num,:);
%y为组合后的矩阵
%zhi为最少值
%yy为取得最少值时各行所选取的数