把位置逆序表转换为对应的置换

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