康托展开

康托展开

康托展开是建立$f:A_n\iff {\mathbb{N}}_n$，即$n$的一个排列到自然数的子集的的一个双射函数。

康托展开

一个排列的康托展开对应的自然数的含义是，这个排列之前有多少排列小于它。因此，我们可以建立如下过程：

例如：$54321$

1. 令$ans=0$。
2. 首先我们看到的数字是$5$，在$5$之后的序列$4321$中有4个比他小的数，$1234$均可以在$5$之前打头，因此令$ans+=4\times 4!$，意思是$5$的前面有四个可能的前序分别生成$4!$的排列。
3. 我们看到的数字是$4$，在$4$之后的序列$321$中有3个比他小的数，$123$均可以在$3$之前打头，因此令$ans+=3\times 3!$，意思是$3$的前面有三个可能的前序分别生成$3!$的排列。
4. 以此类推得到$ans$的值。

逆康托展开

例如：$119$

1. 首先令集合 S = { 1 , 2 , 3 , 4 , 5 } S=\{1,2,3,4,5\} S={1,2,3,4,5}。
2. $1$号位上的位权为$4!$，因此$119/4!=4$,$119\mathrm{m}\mathrm{o}\mathrm{d}4!=23$，因此说明第一个位置有4个比他小的数字，正好是$S_4=5$，因此第一位上的数字是$5$，并删除$S_4$。（除法是整除，$S$的下标从$0$开始）。
3. 以此类推，直到数字均被填满。

例题

P2525

#include <bits/stdc++.h>

using namespace std;

#define FR freopen("in.txt","r",stdin)

typedef long long ll;
int f[9];
int arr[10];
int n;

int cantor()
{
    int val = 0;
    for(int i = 1;i<=n;i++)
    {
        int cnt = 0;
        for(int j = i+1;j<=n;j++) if(arr[j] < arr[i]) cnt++;
        val += f[n - i] * cnt;
    }
    return val;
}

void invcantor(int seq)
{
    vector<int> can;

    for(int i = 1;i<=n;i++)
        can.push_back(i);

    for(int i = 1;i<=n;i++)
    {
        int fact = seq / f[n - i];
        seq %= f[n - i];

        arr[i] = can[fact];
        can.erase(can.begin() + fact);
    }
}

int main()
{
    cin >> n;
    f[0] = f[1] = 1;
    for(int i = 2; i<=8; i++)
        f[i] = i * f[i-1];


    for(int i = 1; i<=n; i++)
        cin >> arr[i];

    int k = cantor();

    if(k == 0)
    {
        cout << "ERROR";
    }
    else
    {
        invcantor(k - 1);
        copy(arr+1,arr+1+n,ostream_iterator<int>(cout, " "));
    }
    return 0;
}