康托展开
原理:
其中, 为整数,并且
。ai表示原数的第i位在当前未出现的元素中是排在第几个
举个例子说明。
在 5个数的排列组合中,计算 34152的康托展开值。
首位是3,则小于3的数有两个,为1和2, ,则首位小于3的所有排列组合为
第二位是4,由于第一位小于4,1、2、3中一定会有1个充当第一位,所以排在4之下的只剩2个,所以其实计算的是在第二位之后小于4的个数。因此
第三位是1,则在其之后小于1的数有0个,所以
第四位是5,则在其之后小于5的数有1个,为2,所以
最后一位就不用计算啦,因为在它之后已经没有数了,所以 固定为0
根据公式:
所以比34152小的组合有61个,即34152是排第62。
//返回数组a中当下顺序的康拖映射
int cantor(int a[],int n)
{
int ans=0;
for(int i=0;i<n;i++)
{
int x=0;int c=1,m=1;//c记录后面的阶乘
for(int j=i+1;j<n;j++)
{
if(a[j]<a[i])x++;
m*=c;c++;
}
ans+=x*m;
}
return ans;
}
逆康托展开:
原理:
一开始已经提过了,康托展开是一个全排列到一个自然数的双射,因此是可逆的。即对于上述例子,在
给出61可以算出起排列组合为34152。由上述的计算过程可以容易的逆推回来,具体过程如下:
用 61 / 4! = 2余13,说明 ,说明比首位小的数有2个,所以首位为3。
用 13 / 3! = 2余1,说明 ,说明在第二位之后小于第二位的数有2个,所以第二位为4。
用 1 / 2! = 0余1,说明 ,说明在第三位之后没有小于第三位的数,所以第三位为1。
用 1 / 1! = 1余0,说明 ,说明在第二位之后小于第四位的数有1个,所以第四位为5。
最后一位自然就是剩下的数2。
通过以上分析,所求排列组合为 34152。
//康托展开逆运算
void decantor(int x, int n)
{
vector<int> v; // 存放当前可选数
vector<int> a; // 所求排列组合
for(int i=1;i<=n;i++)
v.push_back(i);
for(int i=m;i>=1;i--)
{
int r = x % FAC[i-1];
int t = x / FAC[i-1];
x = r;
sort(v.begin(),v.end());// 从小到大排序
a.push_back(v[t]); // 剩余数里第t+1个数为当前位
v.erase(v.begin()+t); // 移除选做当前位的数
}
}