康托展开和康托展开的逆运算

八数码问题不用康托展开判断重复8s,用康托展开判断重复30MS。康托展开最大最明显的作用就是在判断状态是否重复方面了,其实属于hash的一个技巧。

一、康托展开
【问题背景】对于一个有n个不同元素的集合{1,2,3,4,...,n}的从小到大排序(从大到小 同理)的全排列 显然它有n!项。如n=4,那么就有4!=4×3×2×1=24项。
与自然数1,2,3,4,-----n!与之一一对应。比如 1~4四个数的全排列按字典序如下:
1234:第1个
1243:第2个
1324:第3个
1342:第4个
1423:第5个
1432:第6个
2134:第7个
2143:第8个
2314: 第9个
2341:第10个
2413:第11个
2431:第12个
3124:第13个
3142:第14个
3214:第15个
3241:第16个
3412:第17个
3421:第18个
4123:第19个
4132:第20个
4213:第21个
4231:第22个
4312:第23个
4321:第24个

【主要问题】
例1:求4132是第几个排列? 看上面就知道答案就是:20。 那么是怎么算的呢?
解:总共4个数,所以n=4.ans:=0;
第一个数是4,研究比4小的并且还没有出现过的数有3个:1,2,3。那么ans:=ans+ 3*(n-1)!
所以 ans:= ans+ 3* 3*2*1 =18
第二个数是1,研究比1小的并且还没有出现过的数为 0个。那么ans:=ans+ 0 * (n-2)!,那么ans不变。
第三个数是3,研究比3小的并且还没有出现过的数为1个:2。那么ans:=ans+ 1* (n-3)!,那么ans:=18+1* 1=19
第四个数是2,研究比2小的并且还没有出现过的数为0个。那么ans不变。其实最后一个可以不研究了,比它大和比它小的全都出现过了。 最后ans怎么等于19啊??代表它前面有19个排列嘛,那么4132自己就是第20个罗( 最后ans:=ans+1)
例2:问45231是第几个排列?
4 5 2 3 1
ans:= 3*4! + 3*3! + 1*2! + 1*1! + 0*0! + 1 =94

#include<cstdio>
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll n;
ll a[17],d[17];
 
ll kangtuo()
{
    bool bo[17];memset(bo,false,sizeof(bo));
    ll sum=0;
    for(ll i=1;i<=n-1;i++)
    {
        ll k=0;
        for(ll j=i+1;j<=n;j++)   //看第i个数字,它后面有多少数字比它小
     		if(a[i]>a[j])k++;
        sum=sum+k*d[n-i];
      //  cout<<i<<"   "<<k<<endl;
     //   bo[a[i]]=true;
    }
    sum++;
    return sum;
}
int main(){
    scanf("%lld",&n);
    for(int i=1;i<=n;i++)scanf("%lld",&a[i]);
    d[0]=1;for(int i=1;i<=15;i++)d[i]=d[i-1]*i;
    printf("%lld\n",kangtuo());
    return 0;
}

  当数据范围较大时,可以用树状数组进行优化。。洛谷P5367

#include<iostream>
#include<cstdio>
using namespace std;
typedef long long ll;
const int maxn=1e6+10;
const ll p=998244353;
ll n,a[maxn],ans;
ll f[maxn],c[maxn];
void update(int x,int k)
{
	for(int i=x;i<=n;i+=i&-i)
		c[i]+=k;
}
int deal(int x)
{
	int t=0;
	for(int i=x;i>0;i-=i&-i)
		t+=c[i];
	return t;
}
int main()
{
	cin>>n;
	f[0]=1;
	for(int i=1;i<=n;i++)
	f[i]=f[i-1]*i%p;
	for(int i=1;i<=n;i++) update(i,1);
	for(int i=1;i<=n;i++)
	{
		cin>>a[i];
		ans=(ans+(deal(a[i])-1)*f[n-i]%p+p)%p;
		//查找比a[i]小的数字有多少个(还没有出现过的数字中) 
		update(a[i],-1);
	}
	cout<<ans+1;
	return 0;
}

  

二、康托展开的逆运算
我把之前的div改为/,div是pascal语言的运算符号,/是C++的,意思都一样就是求两个整数除法运算后的商(不理会余数)
例3:1~5从小到大全排列中,找出第96个排列?
解:首先设x1x2x3x4x5, (x1等于?不知道),用96-1得到95,表示x1x2x3x4x5前面有95个序列。
第一个数x1,假设x1目前有k个比x1小的并且还没有出现过的数,那么
k:= 95 / (n-1)! = 95 / 24=3, 也就是有3个比x1小并且没有出现过的数,那么x1=4.
95变成95-3×24=23
第二个数x2,假设x2目前有k个比x2小的并且还没有出现过的数,那么
k:= 23 / (n-2)! = 23 / 6 = 3, 也就是有3个比x2小并且没有出现过的数,那么x2=5.(有3个数比它小的数是4,但4已经在之前出现过了,所以是5)
23变成 23 – 3 * 6 = 5
第三个数用5去除2! 得到2余1,那第三个数就是还没有出现过的第3小的数3(5变成1)
第四个数用1去除1! 得到1余0,那第四个数就是还没有出现过的第2小的数2
第五个数就是最后还没有出现的那个。
所以这个数是45321

#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
ll n;
ll a[17],d[17];
void _kangtuo(ll sum)
{
	sum--;
	bool bo[17];memset(bo,true,sizeof(bo));
	for(ll i=1;i<=n;i++)
	{
		ll k=sum/d[n-i];
		sum-=k*d[n-i];
		ll j;
		for(j=1;j<=n;j++)
			if(bo[j]==true)
			{
				if(k==0)break;
				k--;
			}
		bo[j]=false;
		a[i]=j;
	}
	for(ll i=1;i<n;i++)printf("%lld ",a[i]);
	printf("%lld\n",a[n]);
}
int main(){
	scanf("%lld",&n);
	d[0]=1;for(int i=1;i<=15;i++)d[i]=d[i-1]*i;
	ll k;scanf("%lld",&k);
	_kangtuo(k);
	return 0;
}

  

posted @ 2020-08-23 17:34  我微笑不代表我快乐  阅读(151)  评论(0编辑  收藏  举报