function p=inversion1topermutation(I)
%INVERSION1TOPERMUTATION 把位置逆序表转换为对应的置换
% p = inversion1topermutation( I ) 输入位置逆序表I返回对应的置换p.
%算法的基本思想:等于0的最末逆序数所在位置就是当前最大的数所在的位置。
%$Author: WBC$ $Date: 2005/10/18$
n=length(I); %置换的长度
if any(I>(0:n-1))
error(['输入的位置逆序表 (' int2str(I(:)') ') 不合法']);
end
p=zeros(1,n);%初始化存储置换的向量
Izero=zeros(1,n);%存储等于0的逆序数在I中的位置
iz=find(I==0);
Izero(1:length(iz))=iz;%初始化Izero
iend=length(iz);%表示在I中等于0的逆序数的个数
%主循环,从大到小依次找到它们在置换中的位置
for i=n:-1:1
v=Izero(iend); %当前在I中的逆序数为0的最末位置
p(v)=i; %当前最大的数在位置v
iend=iend-1;%用过后iend少1
for l=v:n %对于v右边的逆序数都减1
I(l)=I(l)-1;
if I(l)==0 %记录下新的等于0的逆序数的位置
iend=iend+1;
Izero(iend)=l;
end
end
end