关于康托展开式的问题我已经在第一篇文章提到了,需要的朋友可以点进去看看。
ACM题目之排列序数
康托逆展开式就是给你一个数组,然后让你求出其全排列第n大的序列.
这里借用百度百科给的大家解释下
例1 {1,2,3,4,5}的全排列,并且已经从小到大排序完毕
(1)找出第96个数
首先用96-1得到95
用95去除4! 得到3余23
有3个数比它小的数是4
所以第一位是4
用23去除3! 得到3余5
有3个数比它小的数是4但4已经在之前出现过了所以第二位是5(4在之前出现过,所以实际比5小的数是3个)
用5去除2!得到2余1
有2个数比它小的数是3,第三位是3
用1去除1!得到1余0
有1个数比它小的数是2,第二位是2
最后一个数只能是1
所以这个数是45321
(2)找出第16个数
首先用16-1得到15
用15去除4!得到0余15
用15去除3!得到2余3
用3去除2!得到1余1
用1去除1!得到1余0
有0个数比它小的数是1
有2个数比它小的数是3 但由于1已经在之前出现过了所以是4(因为1在之前出现过了所以实际比4小的数是2)
有1个数比它小的数是2 但由于1已经在之前出现过了所以是3(因为1在之前出现过了所以实际比3小的数是1)
有1个数比它小得数是2 但由于1,3,4已经在之前出现过了所以是5(因为1,3,4在之前出现过了所以实际比5小的数是1)
最后一个数只能是2
所以这个数是14352。
code:
//康托逆展开式
//求n个数的全排列中第m大的组合:
//n>20就需要用高精度来做
#include <iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int fact[10]= {1};
int res[10];
int n,m;
int main()
{
for(int i=1; i<=10; ++i)
fact[i]=fact[i-1]*i; //阶乘数组
while(scanf("%d%d",&n,&m)!=EOF)
{
bool *vis=new bool[n];
memset(res,0,sizeof(res));
memset(vis,false,sizeof(vis));
//先判断下m是否比n!还要大
if(m>fact[n]){cout<<"impossible"<<endl;continue;}
int i,j,x,t;
m--;
for(i=1; i<=n; ++i)
{
x=m/fact[n-i];
m%=fact[n-i];
for(j=1,t=0; t<=x; j++)
if(!vis[j])t++;
j--;
res[i-1]=j;
vis[j]=true;
}
for(i=0; i<n; ++i)
printf("%d",res[i]);
printf("\n");
}
// cout << "Hello world!" << endl;
return 0;
}