康托展开

最新推荐文章于 2020-02-20 16:04:52 发布

u011008379

最新推荐文章于 2020-02-20 16:04:52 发布

阅读量719

点赞数

分类专栏：相关概念

公式

其中,为整数,并且0<=<i,1<=i<=n。

的意义参见举例中的解释部分

举例

例如，3 5 7 4 1 2 9 6 8 展开为 98884。因为X=2*8!+3*7!+4*6!+2*5!+0*4!+0*3!+2*2!+0*1!+0*0!=98884.

解释：

排列的第一位是3，比3小的数有两个，以这样的数开始的排列有8!个，因此第一项为2*8!

排列的第二位是5，比5小的数有1、2、3、4，由于3已经出现，因此共有3个比5小的数，这样的排列有7!个，因此第二项为3*7!

以此类推，直至0*0!

用途

显然，n位（0~n-1）全排列后，其康托展开唯一且最大约为n!，因此可以由更小的空间来储存这些排列。由公式可将X逆推出唯一的一个排列。

康托展开的逆运算

既然康托展开是一个双射，那么一定可以通过康托展开值求出原排列，即可以求出n的全排列中第x大排列。

如n=5,x=96时：

首先用96-1得到95，说明x之前有95个排列.(将此数本身减去！)
用95去除4! 得到3余23，说明有3个数比第1位小，所以第一位是4.
用23去除3! 得到3余5，说明有3个数比第2位小，所以是4，但是4已出现过，因此是5.
用5去除2!得到2余1，类似地，这一位是3.
用1去除1!得到1余0，这一位是2.
最后一位只能是1.
所以这个数是45321.

按以上方法可以得出通用的算法。

下面转载自：http://baike.baidu.com/view/437641.htm?fr=aladdin。

下面的程序指示出了一个简洁的计算方法，和前文的运算思路略有不同，不需要检验某数码是否使用过，只需检查第(n+1-i)位之后比第(n+1-i)位小的位的数量，将这个数量作为公式中的a[i]。(1<=i<=n)

C++

 
        const int PermSize=12; 
       
        long long factory[PermSize]={0,1,2,6,24,120,720,5040,40320,362880,3628800,39916800}; 
       
        long long Cantor(stringbuf) 
       
        { 
       
        int i,j,counted; 
       
        long long result=0; 
       
        for 
        (i=0;i<PermSize;++i) 
       
        { 
       
        counted=0; 
       
        for 
        (j=i+1;j<PermSize;++j) 
       
        if 
        (buf[i]>buf[j]) 
       
        ++counted; 
       
        result=result+counted*factory[PermSize-i-1]; 
       
        } 
       
        return result; 
       
        }