康托展开

康托展开是一个全排列到一个自然数的双射,常用于构建哈希表时的空间压缩。 康托展开的实质是计算当前排列在所有由小到大全排列中的顺序,因此是可逆的。

以下称第x个全排列是都是指由小到大的顺序。

公式:

X=a[n]*(n-1)!+a[n-1]*(n-2)!+…+a[i]*(i-1)!+…+a[1]*0!

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

公式看起来可能不太明白,其实很简单。。

举个例子:

对abcdefghijkl(12个字符)进行字典序全排列(共有12!种排列),现在给你一个排列,让你求它是第几个排列(即展开)。比如给gfkedhjblcia,展开即为:X = 6*11!+ 5*10! + 8*9! + 4*8! + 3*7! + 3*6! + 4*5! + 1*4! + 3*3! + 1*2! + 1*1! + 0*0! = 260726925 (这答案是在排在他前面的序列数,他自己的位次应该是这个+1).

再看公式:X=a[n]*(n-1)!+a[n-1]*(n-2)!+…+a[i]*(i-1)!+…+a[1]*0!

a[n]是序列第n个数(从右往左数)在他之前(即在它右边)比它小的数的数量,像例子中a[12] = 6, 第12个数是g, 它右边比它小的数有a, b, c, d, e, f. 之后的a[11], a[10]....a[1]同理。

原理:

知道原理就很简单了,也不用死记硬背了,要计算这个序列是全排列的第几个,可以通过计算它前面有几个,再加一就是它是第几个。

从最高位开始算(最左边),这位比它小的,后面几位无论怎么排,序列肯定比它小,所以找出能放在这一位(在它前面不出现)且比它小的个数再乘上后面位数阶乘就行了。这一位相同的情况下,下一位同理。

还是没看懂的话看看这两个举例吧,

举例:

例如,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!  

例如,想知道321是{1,2,3}中第几个大的数可以这样考虑 :  
  第一位是3,当第一位的数小于3时,那排列数小于321 如 123、 213 ,小于3的数有1、2 。所以有2*2!个。再看小于第二位2的:小于2的数只有一个就是1 ,所以有1*1!=1 所以小于321的{1,2,3}排列数有2*2!+1*1!=5个  。所以321是第6个大的数。 2*2!+1*1!是康托展开。    
再举个例子:1324是{1,2,3,4}排列数中第几个大的数:
第一位是1小于1的数没有,是0个 0*3! 第二位是3小于3的数有1和2,但1已经在第一位了,所以只有一个数2 1*2! 。第三位是2小于2的数是1,但1在第一位,所以  有0个数 0*1! ,所以比1324小的排列有0*3!+1*2!+0*1!=2个,1324是第三个大数。


康托展开的逆运算:

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

就是上面的反过来计算,自己思考吧。

康托展开的代码:

题目地址:https://qduoj.com/problem/26/
/*
样例输入1 复制
abcdefghijkl
abcdefghiklj
gfkedhjblcia
样例输出1
1
4
260726926
*/
#include<iostream>
#include<cstdio>
using namespace std;
int main(void)
{
    long long fc[13] = {1};
    //计算阶乘
    for(int i = 1; i < 13; i++) fc[i] = fc[i-1] * i;
    char str[20];
    while(~scanf("%s", str))
    {
        long long ans = 1;
        for(int i = 0; i < 12; i++)
        {
            int count = 0;
            //找它后面比它小的数数量
            for(int j = i; j < 12; j++)
                if(str[j] < str[i]) count++;
            printf(" %d ", count);
            ans += count * fc[11-i];
        }
        printf("%lld\n", ans);
    }
    return 0;
}






  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值