康拓展开+康拓展开的逆运算

康拓展开：

是一个全排列到一个自然数的双射，常用于构建哈希表时的空间压缩。康托展开的实质是计算当前排列在所有由小到大全排列中的顺序(名次)，因此是可逆的。

公式

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

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

举例

例如，一个数组s={1,2,3,4,5,6,7,8,9}，它的一个排列s1={3 ,5 ,7 ,4 ,1 ,2, 9 ,6 ,8} 展开为 98884。因为X=2*8!+3*7!+4*6!+2*5!+0*4!+0*3!+2*2!+0*1!+0*0!=98884.

解释：

排列的第一位是3，比3小的数有两个，以这样的数开始的排列有8!个，因此第一项为2*8!

排列的第二位是5，比5小的数有1、2、3、4，由于3已经出现，因此共有3个比5小的数，这样的排列有7!个，因此第二项为3*7!

以此类推，直至0*0!

用途

显然，n位（0~n-1）全排列后，其康托展开唯一且最大约为n!，因此可以由更小的空间来储存这些排列。由公式可将X逆推出对应的全排列。

康托展开的逆运算：

根据某个排列在总的排列中的名次来确定这个排列。

[既然康托展开是一个双射，那么一定可以通过康托展开值求出原排列，即可以求出n的全排列中第x大排列。]

如n=5,x=96时：最开始最小的序列s[5]={1,2,3,4,5} 首先用96-1得到95，说明x之前有95个排列.(将此数本身减去！) 用95去除4! 得到3余23，说明有3个数比第1位小，所以第一位是s[3]=4. 用23去除3! 得到3余5，说明有3个数比第2位小，所以是4，但是4已出现过，因此是5.(除去出现过的4,s[3]=5) 用5去除2!得到2余1，类似地，这一位是3.(除去了出现过的4和5，s[2]=3) 用1去除1!得到1余0，这一位是2.(除去出现过的4和5和3，s[1]=2) 最后一位只能是1.(除去出现过的4和5和3和2) 所以这个数是45321.

 

//求康拓展开   以及  康拓展开的逆运算
#include<iostream>
#include<cstring>
#include<string.h>
#include<algorithm>
#include<cmath>
using namespace std;

#define MAX 1000000

int factorial(int n) //求n的阶乘 
{
	int product=1;
	for(int i=1;i<=n;i++)
		product*=i;
	return product;
}

//康拓展开 
int kangtuo(int a[],int len)
{
	int ans=0;
	for(int i=0;i<len;i++){
		int lessnum=0;
		for(int j=i+1;j<len;j++){  //求未出现的元素中，a[i]排第几，即有几个比a[i]小 
			if(a[j]<a[i]) lessnum++; 
		}
		ans+=lessnum*factorial(len-1-i); 
		// X=a[n]*(n-1)!+a[n-1]*(n-2)!+…+a[i]*(i-1)!+…+a[1]*0!
	}
	return ans+1; //返回该排列在全排列中的位置，ans从0开始 ,所以返回ans+1 
}

//康拓展开的逆运算 (只能用于求‘连续数字的数组’)
/* 
void ktinverse(int len,int pos, int* &ans)
{
	bool visit[len];
	memset(visit,0,sizeof(visit));
	pos--;
	int t,j;
	for(int i=0;i<len;i++){
		t=pos/factorial(len-1-i);
		for( j=0;j<len;j++){
			if(!visit[j]){
				if(t==0) break;
				t--;
			}
		}
		ans[i] = j+1;
		visit[j] = true;
		pos%=factorial(len-1-i);
	}
}
*/

//康拓展开的逆运算 
void ktinverse(int len,int pos, int* &ans,int a[])  
{  //len:排列长度 ， pos：该排列在全排列中的位置，a[]全排列中最小的排列 
	pos--;
	int t;
	for(int i=0;i<len;i++){
		t=pos/factorial(len-1-i); 
		ans[i] = a[t];
		a[t]=MAX;
		sort(a,a+len);
		pos%=factorial(len-1-i);
	}
}



int main()
{
	int a[5]={4,5,3,2,1};
	cout << kangtuo(a,5)<< endl;
	int *ans = new int[5];
	sort(a,a+5);
	ktinverse(5,96,ans,a);
	for(int i=0;i<5;i++)
	cout << ans[i] << endl;
	
	
} 

前面的知识点转载自：

http://www.acmerblog.com/cantor-expansion-3647.html