康托展开求全排列
解析
把一个整数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;
}