康拓展开,康拓展开逆运算

康托展开

   康托展开的公式是 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/

代码是本人写的,可能写的不好讲得也不好,如果实在看不懂我很抱歉,请请忽略


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值