首先什么是康托展开?
X=a[n]*(n-1)!+a[n-1]*(n-2)!+...+a[i]*(i-1)!+...+a[1]*0! ,其中a[i]为当前未出现的元素中是排在第几个(从0开始)。这就是康托展开;
给出一个1~n的的n个数的一种排列方式, 用康托展开可以求出该排列方式是第几种排列方式;
如:123, 132, 213, 231, 312, 321 中
3 1 2 是第五种;
康托展开为 2*2!+0*1!+0*0! = 4;表示有四种情况比3 1 2 小, 它是第五种排列;
看第一位数:剩下的数中比3小的有2个,后两位 有2!种排列情况;(3*2!)
当第一位为3时, 看第二位:剩下的数中比1小的有0个,后一位 有1!种排列情况;(0*1!)
当第二位为1时, 看第三位:剩下的数中比2小的有0个,后零位 有0!种排列情况;(0*0!)
加起来就是4;
同理知道一个排列的康托展开是几 ,可以求出这个排列;
假设n=5, 康托展开是23;
23/(4!)=0;说明有0个数比第一位小,第一位是1; 23%(4!)=23;
23/(3!)=3; 说明剩下的数中有3个数比第二位小,第二位是5; 23%(3!)=5;
5/(2!)=2; 说明剩下的数中有2个数比第三位小,第三位是4; 5%(2!)=1;
1/(1!)=1; 说明剩下的数中有1个数比第四位小,第四位是3; 1%(1!)=0;
0/(0!)=0; 说明剩下的数中有0个数比第五位小,第五位是2;
所以排列为1 5 4 3 2;
现在用康拓展开验证一下:0*4!+3*3!+2*2!+1*1!+0*0!=23;
最后是代码实现:
#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <algorithm>
using namespace std;
int fact[]={1, 1, 2, 6, 24, 120, 720, 5040, 40320, 362880};
int n;
//康托展开:
int canto_ex(int a[]){
int sum=0;
for(int i=0; i<n; i++){
int tmp=0;
for(int j=i+1; j<n; j++)
if(a[j]<a[i]) tmp++;
sum+=tmp*fact[n-1-i];
}
return sum;
}
//康托逆展开:
void re_canto_ex(int sum, int a[]){
int k=sum, x;
int vis[1000];
for(int i=1; i<n; i++) vis[i]=0;
for(int i=0; i<n; i++){
x=k/fact[n-1-i];
k=k%fact[n-1-i];
for(int j=1; j<=n; j++){
if(vis[j]) continue;
if(x==0) a[i]=j, vis[j]=1;
x--;
}
}
return;
}