蓝桥杯：试题 算法训练 星际交流 康托展开

题目

资源限制
时间限制：1.0s 内存限制：256.0MB

问题描述
　　人类终于登上了火星的土地并且见到了神秘的火星人。人类和火星人都无法理解对方的语言，但是我们的科学家发明了一种用数字交流的方法。这种交流方法是这样 的，首先，火星人把一个非常大的数字告诉人类科学家，科学家破解这个数字的含义后，再把一个很小的数字加到这个大数上面，把结果告诉火星人，作为人类的回 答。
　　火星人用一种非常简单的方式来表示数字——掰手指。火星人只有一只手，但这只手上有成千上万的手指，这些手指排成一列，分别编号为1，2，3……。火星人的任意两根手指都能随意交换位置，他们就是通过这方法计数的。
　　一个火星人用一个人类的手演示了如何用手指计数。如果把五根手指——拇指、食指、中指、无名指和小指分别编号为1，2，3，4和5，当它们按正常顺序排列 时，形成了5位数12345，当你交换无名指和小指的位置时，会形成5位数12354，当你把五个手指的顺序完全颠倒时，会形成54321，在所有能够形 成的120个5位数中，12345最小，它表示1；12354第二小，它表示2；54321最大，它表示120。下表展示了只有3根手指时能够形成的6个 3位数和它们代表的数字：

三进制数
123
132
213
231
312
321
代表的数字
1
2
3
4
5
6

现在你有幸成为了第一个和火星人交流的地球人。一个火星人会让你看他的手指，科学家会告诉你要加上去的很小的数。你的任务是，把火星人用手指表示的数与科 学家告诉你的数相加，并根据相加的结果改变火星人手指的排列顺序。输入数据保证这个结果不会超出火星人手指能表示的范围。
　　
输入格式
　　包括三行，第一行有一个正整数N，表示火星人手指的数目（1 <= N <= 10000）。第二行是一个正整数M，表示要加上去的小整数（1 <= M <= 100）。下一行是1到N这N个整数的一个排列，用空格隔开，表示火星人手指的排列顺序。
　　
输出格式
　　只有一行，这一行含有N个整数，表示改变后的火星人手指的排列顺序。每两个相邻的数中间用一个空格分开，不能有多余的空格。
　　
样例输入

5
3
1 2 3 4 5

样例输出

1 2 4 5 3

数据规模和约定
对于30%的数据，N<=15；
对于60%的数据，N<=50；
对于全部的数据，N<=10000；

思路

可以AC的思路：【蓝桥杯：试题 算法训练 星际交流 全排】

其实也可以通过计算康托展开，推导出到底是第几个排列，然后加上m，在通过这个数字求出对应的序列

具体实现见【康托展开式与其逆运算的实现】

由序列求排列号，就是：

• 计算每一项的系数，第x位系数是当前位x后面，有多少个数大于当前数字
• 然后数字就是每一个位的系数 * 对应的阶乘

逆运算：
由数字求排列，同样

• 通过短除法确定位置 i 的系数 x
• 选取 第x个未被访问的数字作为当前位 i 的数字

代码

（不能AC因为阶乘到21 longloong 就溢出 了）

#include <bits/stdc++.h>

using namespace std;

#define maxlen 19
typedef long long ll;
int n, bias;
int a[maxlen], vis[maxlen], fac[maxlen];

ll fac_(int x)
{
	ll s=1;
	for(int i=2; i<=x; i++) s*=i;
	return s;
}

// 计算当前序列是第几个排列 
int kth(int* begin, int* end)
{
	int len=end-begin, sum=0;
	for(int i=0; i<len; i++)
	{
		int cnt = 0;
		// 当前位系数是：有多少个数比当前数小 
		for(int j=i+1; j<len; j++)
			if(begin[i]>begin[j]) cnt++;
		sum += cnt*fac[len-i-1];
	}
	return sum;
}

// 通过数字x得到对应序号的排列 
void getseq(int x)
{
	memset(vis, 0, sizeof(vis));
	int nums[maxlen];
	// 计算系数 
	for(int i=0; i<n; i++)
		nums[i]=x/fac[n-i-1], x%=fac[n-i-1];
	for(int i=0; i<n; i++)
	{
		int j;
		// 找第nums[i]+1个没访问的数字 
		for(j=0; j<n&&nums[i]>=0; j++)
			if(!vis[j]) nums[i]--; 
		nums[i]=j, vis[j-1]=1;
	}
	// 输出 
	for(int i=0; i<n; i++) cout<<nums[i]<<" ";
	cout<<endl; 
}

int main()
{
	cin>>n>>bias;
	for(int i=0; i<n; i++) 
		cin>>a[i], fac[i]=fac_(i);
	
	int k = kth(a, a+n);
	k += bias;
	getseq(k);
	
	return 0;
}
/*
3 0
1 2 3 
*/