康托展开的公式是 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 等于啥?
(1)a4 = "D" 这个元素在子数组 ["D", "B", "A", "C"] 中是第几大的元素。"A"是第0大的元素,"B"是第1大的元素,"C" 是第2大的元素,"D"是第3大的元素,所以 a4 = 3。
(2)a3 = "B" 这个元素在子数组 ["B", "A", "C"] 中是第几大的元素。"A"是第0大的元素,"B"是第1大的元素,"C" 是第2大的元素,所以 a3 = 1。
(3)a2 = "A" 这个元素在子数组 ["A", "C"] 中是第几大的元素。"A"是第0大的元素,"C"是第1大的元素,所以 a2 = 0。
(4)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
康拓展开代码:
#include<stdio.h>
#include<string.h>
typedef long long LL;
LL f[20]={0,1};
void init()//得到0~19阶乘值
{
int i=2;
for(;i<20;i++)
{
f[i]=i*f[i-1];
}
}
LL getNum(char str[])//康拓展开
{
int len=strlen(str);
int i,j;
LL ans=0;
for(i=0;i<len;i++)
{
int tmp=0;
for(j=i+1;j<len;j++)//找到i后面有几个比它小的,就知道他排第几了
{
if(str[i]>str[j])tmp++;
}
ans=ans+tmp*f[len-1-i];
}
return ans;
}
int main()
{
init();
char str[1000]="bckfqlajhemgiodnp";
printf("%lld\n",getNum(str));
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 个排列。
康拓展开的逆运算代码:
这里我用vector代替过程中用到的一些数组,这样会简便一点
vector :http://www.cnblogs.com/scandy-yuan/archive/2013/01/07/2849735.html
#include<stdio.h>
#include<vector> //这里我用vector代替数组,会简便一点
#include<string.h>
using namespace std;
typedef long long LL;
const char str[200]="abcdefghijklmnopq"; //第0个字符串 即未排列的初始字符串
LL f[20]={0,1};//阶乘结果
void init()//得到0~19内的阶乘结果,存到f[]中
{
int i;
for(i=2;i<=20;i++)
{
f[i]=i*f[i-1];
}
}
void getStr(LL n)
{
int len=strlen(str); //字符串的长度
int i,j=0;
int cnt=0; //根据n计算出cnt,表示一共进行了几个数的全排列
//比如cnt=0表示第n个序号的排列,前面进行了一位的全排列
//cnt=1表示第n个序号的排列,进行了两位的全排列
//即 cnt 表示要到第n个序号的排列,需要对cnt+1位进行全排列
for(i=0;i<20;i++)
{
if(n<f[i]){ //比f[i]的阶乘小,说明没有到f[i]这么大
cnt=i-1; //则一共全排列了cnt+1个字符
break;
}
}
//那么原字符串就可以分为两个部分,一个部分是不用全排保持原状的字符串
//另一部分就是要进行全排的字符串
vector<char>cv(str+len-cnt-1,str+len); //1.截取字符串全排列的字符串
vector<char>rv(str,str+len-cnt-1);//2.存放结果的字符串,前面没有参与全排列的不变照搬进来
vector<int>nv;//存放要取的字符在 1.截取字符串中的 位置
vector<char>::iterator it;
while(cnt>0) //利用辗转相除法得到每一个字符的位置
{
nv.push_back(n/f[cnt]);
n=n%f[cnt];
cnt--;
}
nv.push_back(0);
i=0;
while(cv.size()>0) //从参与全排列的字符串中 取出指定位置的字符
{
int pos=nv[i];
rv.push_back(cv[pos]);
it=cv.begin();
while(it!=cv.end())
{
if(*it==cv[pos]){
break;
}
it++;
}
i++;
cv.erase(it);
}
it=rv.begin(); //输出结果
while(it!=rv.end())
{
printf("%c",*it);
it++;
}
}
int main()
{
init(); //初始化得到0~19的阶乘值
LL n=22952601027516;
getStr(n);
return 0;
}
以上文字讲解内容转自:http://archive.cnblogs.com/a/2026276/
代码是本人写的,可能写的不好讲得也不好,如果实在看不懂我很抱歉,请请忽略