matlab组合数表示方法,组合数的序号问题

附录---MATLAB代码

function [pmin,pvmin]=tspdp(D)

%TSPDP旅行销售员问题(TSP)的动态规划(Dynamic Programming)算法的递推实现

%

%[pmin, pvmin] = tspdp( D )

%输入参数：D --- 两两城市之间的距离矩阵，双精度二维数组

%输出参数：pmin --- 最短路线双精度，一维行数组

%pvmin --- 最短距离，双精度数

%

%Example: [pmin,pvmin]=tspdp(rand(5))

%

%注意：经分析，其时间复杂度为O(n*2^n)，也就是说n 每增大1，运行时间大致上变成

%原来的2 倍。因此，n不能太大，当取在12以内时，速度很快。若用C++等实现,在20以

%内都可以接受。

n=length(D);

pasmat=pascal(n+1);

fP=calfP(D,pasmat);

pvmin=inf;

for j=1:n-1

v=fP{n-1}(j,1)+D(j,n);

if pvmin>v

pvmin=v;

jmin=j;

end

end

pmin=zeros(1,n+1);

pmin(n:n+1)=[jmin, n];

C=1:n-1;

for j=n-1:-1:1

Csn=combineseqnum(C,n-1,j,pasmat);

I=find(pmin(j+1)==C);

pmin(j)=fP{j}(I,2,Csn);

C(I)=[];

end

%------------------------------------------------------------------------

function fP=calfP(D,pasmat)

n=length(D);

for k=1:n-1

fP{k}=zeros(k,2,pasmat(k+1,n-k));

end

for Csn=1:n-1

fP{1}(:,:,Csn)=[D(n,Csn),n];

end

for k=2:n-1

C=zeros(1,k);

Csn=0;

curi=1;

while curi

if C(curi)>=n-1-k+curi

curi=curi-1;

else

C(curi:end)=C(curi)+1:C(curi)+1+k-curi;

curi=k;

Csn=Csn+1;

S=C(2:k);

for I=1:k

CI=C(I); %后面要重复使用

Ssn=combineseqnum(S,n-1,k-1,pasmat);

vmin=inf;

for j=1:k-1

v=fP{k-1}(j,1,Ssn)+D(S(j),CI);

if vmin>v

vmin=v;

jmin=S(j);

end

end

fP{k}(I,:,Csn)=[vmin, jmin];

if I

S(I)=CI;

end

end

end

end

end

%-----------------------------------------------------------------------

function sn=combineseqnum(p,n,r,pasmat)

%COMBINESEQNUM求{1,2,...,n}中抽取r个数构成的组合在所有组合(升序)按字典序排

%序后所处的序号(从1开始编号)。

%

%sn = combineseqnum(p , n)

%输入参数：p --- {1,2,...,n}的一个组合，用数值向量表示，行或列均可，要求按

%升序排列

%n --- 元素的个数，即{1,2,...,n}中的 n

%r --- p中元素的个数

%pasmat --- 帕斯卡矩阵

%输出参数：sn --- p在所有组合中的按字典序排序后的编号

%

%Example: sn=combineseqnum([1 3],4)

%

if r>1

sn=pasmat(r+1,n-r+1);

for I=1:r

nn=n-p(I);

rr=r+1-I;

if nn>=rr

sn=sn-pasmat(rr+1,nn-rr+1);

end

end

else

sn=p;

end