康托展开

康托展开

　　康托展开的公式是 X=an*(n-1)!+an-1*(n-2)!+...+ai*(i-1)!+...+a2*1!+a1*0! 其中，ai为当前未出现的元素中是排在第几个（从0开始）。
　　这个公式可能看着让人头大，最好举个例子来说明一下。例如，有一个数组 s = ["A", "B", "C", "D"]，它的一个排列 s1 = ["D", "B", "A", "C"]，现在要把 s1 映射成 X。n 指的是数组的长度，也就是4，所以
X(s1) = a4*3! + a3*2! + a2*1! + a1*0!
关键问题是 a4、a3、a2 和 a1 等于啥？
a4 = "D" 这个元素在子数组 ["D", "B", "A", "C"] 中是第几大的元素。"A"是第0大的元素，"B"是第1大的元素，"C" 是第2大的元素，"D"是第3大的元素，所以 a4 = 3。
a3 = "B" 这个元素在子数组 ["B", "A", "C"] 中是第几大的元素。"A"是第0大的元素，"B"是第1大的元素，"C" 是第2大的元素，所以 a3 = 1。
a2 = "A" 这个元素在子数组 ["A", "C"] 中是第几大的元素。"A"是第0大的元素，"C"是第1大的元素，所以 a2 = 0。
a1 = "C" 这个元素在子数组 ["C"] 中是第几大的元素。"C" 是第0大的元素，所以 a1 = 0。（因为子数组只有1个元素，所以a1总是为0）
所以，X(s1) = 3*3! + 1*2! + 0*1! + 0*0! = 20

A B C | 0
A C B | 1
B A C | 2
B C A | 3
C A B | 4
C B A | 5

康托展开表示的是当前排列在n个不同元素的全排列中的名次。比如213在这3个数所有排列中排第3。

 

那么，对于n个数的排列，康托展开为：

 

其中表示第i个元素在未出现的元素中排列第几。举个简单的例子：

 

对于排列4213来说，4在4213中排第3，注意从0开始，2在213中排第1，1在13中排第0，3在3中排第0，即：

 

，这样得到4213在所有排列中排第ans=20

 

代码实现：(从0开始计数)

#include <iostream>
#include <cstring>
#include <vector>
#include <algorithm>
using namespace std;
typedef long long LL;
LL fact[20];
void getFact(int n)
{
	fact[0]=1;
	for(int i=1;i<n;i++)
	{
		fact[i]=fact[i-1]*i;
	}
}
//康托展开,给出排列方式,求序号(从0开始) 
LL Cantor(int a[],int n)
{
	LL ans=0;
	for(int i=0;i<n;i++)
	{
		int cnt=0;
		for(int j=i+1;j<n;j++)
			if(a[j]<a[i])
				cnt++;
		ans+=cnt*fact[n-i-1];
	}
	return ans;
}

int main()
{
	getFact(20);
	int a[]={4,2,1,3};
	cout<<Cantor(a,4)<<endl;
	
	return 0;
}

通过康托逆展开生成全排列

　　如果已知 s = ["A", "B", "C", "D"]，X(s1) = 20，能否推出 s1 = ["D", "B", "A", "C"] 呢？
　　因为已知 X(s1) = a4*3! + a3*2! + a2*1! + a1*0! = 20，所以问题变成由 20 能否唯一地映射出一组 a4、a3、a2、a1？如果不考虑 ai 的取值范围，有
3*3! + 1*2! + 0*1! + 0*0! = 20
2*3! + 4*2! + 0*1! + 0*0! = 20
1*3! + 7*2! + 0*1! + 0*0! = 20
0*3! + 10*2! + 0*1! + 0*0! = 20
0*3! + 0*2! + 20*1! + 0*0! = 20
等等。但是满足 0 <= ai <= n-1 的只有第一组。可以使用辗转相除的方法得到 ai，如下图所示：

知道了a4、a3、a2、a1的值，就可以知道s1[0] 是子数组["A", "B", "C", "D"]中第3大的元素 "D"，s1[1] 是子数组 ["A", "B", "C"] 中第1大的元素"B"，s1[2] 是子数组 ["A", "C"] 中第0大的元素"A"，s[3] 是子数组 ["C"] 中第0大的元素"C"，所以s1 = ["D", "B", "A", "C"]。
这样我们就能写出一个函数 Permutation3()，它可以返回  s 的第 m 个排列。

康托展开的逆运算：就是根据某个排列的在总的排列中的名次来确定这个排列。比如：

 

求1234所有排列中排第20的是啥，那么就利用辗转相除法确定康托展开中的系数，然后每次输出当前未出现过的第个元素。

 

代码实现康托展开逆运算：

#include <iostream>
#include <cstring>
#include <vector>
#include <algorithm>
using namespace std;
typedef long long LL;
LL fact[20];
void getFact(int n)
{
	fact[0]=1;
	for(int i=1;i<n;i++)
	{
		fact[i]=fact[i-1]*i;
	}
}

//1-m组成的第n个排列, 
void CntorInver(LL n,LL m)
{
	
	//n--;//从1开始计数则-- 
	vector<int>v;
	vector<int>a;
	for(int i=1;i<=m;i++)
		v.push_back(i);
	for(int i=m;i>0;i--)
	{
		LL r = n % fact[i-1];
		LL t = n / fact[i-1];
		n = r;
		sort(v.begin(),v.end());
		a.push_back(v[t]);
		v.erase(v.begin()+t);
	}
	vector<int>::iterator it;
	for(it = a.begin();it!=a.end();it++)
		cout<<*it;
	cout<<endl;
}
int main()
{
	getFact(20);
	int a[]={4,2,1,3};
	cout<<Cantor(a,4)<<endl;
	CntorInver(20,4);
	return 0;
}

参考博客:

http://blog.csdn.net/zhongkeli/article/details/6966805

http://blog.csdn.net/acdreamers/article/details/7982067