康托展开求全排列

康托展开求全排列

解析

把一个整数X展开成如下形式：

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

其中a[i]为当前未出现的元素中是排在第几个（从0开始），并且0<=a[i]

1.求序列的次序

带入公式即可求出

2.求给定次序的序列

由康托展开可知，给定cnt是由n-1阶到0阶共n项a[m] * (m-1)!组成(m∈[1,n])。可以从n-1阶开始，求出最大的a[m]，满足a[m] * (m-1)!

例题可见vijos 全排列

测试样例：20 999999999

结果：1 2 3 4 5 6 7 10 9 8 16 18 15 17 20 12 14 13 11 19

注：这里的a[i]并没有在代码中用数组表示，是在1-n的遍历中求出最大的满足条件的i赋值给A[cur]

#include<iostream>
#include<cstdio>
using namespace std;
long long  count = 0;
int A[1000];
unsigned long long  n,m;
//阶乘
unsigned long long  fact(int s)
{
    unsigned long long  sum = 1;
    for(int i = 1;i<=s;i++)
        sum*=i;
    return sum;
}
//康托展开求第cnt个展开式的第cur个元素，cur表示数组长度
void contor(int * A,int cur,unsigned long long cnt)
{
    if(cur == n) return ;//长度等于n回溯
    for(int i = 1;i<=n;i++)//从1-n中取合适的值向A[cur]赋值
    {
        int flag = 1,level=0;
        for(int j=0;j<cur;j++)
        {
            if(A[j] == i) flag = 0;//如果前面出现i，则不赋值i
            else if(A[j] < i) level ++;
             //看数组前面出现了几个比i小的的值，i-level即未出现的值中所排次序
        }
        //找到i所能赋的最大值，然后递减这个值
        if(cnt <= (i-level) * (unsigned long long )fact(n - cur -1)&&flag == 1)
        {
            A[cur] = i;
            cnt-=(i-level-1)* fact(n - cur - 1);
            contor(A,cur+1,cnt);//再进行康托展开
            break;
        }
    }
}
int main()
{
    cin>>n>>m;
    contor(A,0,m);
    for(int i=0;i<n;i++)
        cout<<A[i]<<" ";
    cout<<endl;
}