康托展开

最新推荐文章于 2024-09-19 22:05:37 发布

原创最新推荐文章于 2024-09-19 22:05:37 发布 · 709 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#数论 #康托展开 #算法

数学专栏收录该内容

43 篇文章

订阅专栏

本文详细介绍了康托展开的概念及其应用，包括如何通过该方法将全排列转换为自然数，以及如何进行逆运算来求得特定位置的排列。此外，还提供了具体的计算实例帮助理解。

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

以下称第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;
}