hdu多校第五场1005 （hdu6628） permutation 1 排列/康托展开/暴力

题意：

定义一个排列的差分为后一项减前一项之差构成的数列，求对于n个数的排列，差分的字典序第k小的那个，n<=20，k<=1e4。

题解：

暴力打表找一遍规律，会发现，对于n个数的排列，如果想找到差分的字典序第k小的，如果k<=(n-1)!，那么对应的那个排列就是把第一位赋值为n，后面的是1~n-1的元素本身排列字典序第k小的。

比如，4个元素的排列的差分字典序最小的前6个分别是

4，1，2，3

4，1，3，2

4，2，1，3

4，2，3，1

4，3，1，2

4，3，2，1

当n为10或更多的时候，(n-1)!>1e4，便可用康托逆展开直接计算。

当n为9以下时，原先的想法是在本地暴力排序，对于每个n都取前1e4个元素交表。后来发现hdu限制提交代码大小，分析了一下，9!约等于3e5，暴力打表，过了。

#include<iostream>
#include<vector>
#include<algorithm>
#define LL long long
using namespace std;
int MAXN;
LL fact[];
int tmp[];
struct Node{
    int px[];//排列
    int cha[];//排列的差分
    int indexpx;
    int indexcha;
    friend bool operator > (const Node &a,const Node &b){
        for(int i=;i<MAXN;i++){
            if(a.cha[i]!=b.cha[i])return a.cha[i]>b.cha[i];
        }
    }
    friend bool operator < (const Node &a,const Node &b){
        for(int i=;i<MAXN;i++){
            if(a.cha[i]!=b.cha[i])return a.cha[i]<b.cha[i];
        }
    }
    //比较差分的字典序
}node[];
int ans[][];
void make_ans(){
    for(int i=;i<=MAXN;i++){
        node[].px[i]=i;
        node[].cha[i-]=;
    }
    node[].indexpx=;
    int i,n;
    for(i=;;i++){
        for(int j=;j<=MAXN;j++){
            node[i].px[j]=node[i-].px[j];
        }
        if(!next_permutation(&node[i].px[],&node[i].px[MAXN+])){
            n=i-;
            break;
        }

        for(int j=;j<MAXN;j++){
            node[i].cha[j]=node[i].px[j+]-node[i].px[j];
        }

        node[i].indexpx=i;
    }
    sort(node+,node++n);
    for(int i=;i<=min(n,);i++){
        ans[MAXN][i]=node[i].indexpx;
    }
    //暴力打表预处理
}
void make_fact(){
    fact[]=;
    for(int i=;i<=;i++){
        fact[i]=fact[i-]*i;
    }
}
void Cantor_invexp(int *p,int len,LL rank){
    //康托逆展开
    int temp[len];
    for(int i=;i<len;i++){
        temp[i]=i+;
    }
    for(int i=;i<=len;i++){
        int a=rank/fact[len-i];
        rank%=fact[len-i];
        for(int j=;j<len;j++){
            if(a== && temp[j]>){
                p[i]=temp[j];
                temp[j]=;
                break;
            }else if(temp[j]>){
                a--;
            }
        }
    }
    return ;
}

int main(){
    make_fact();
    for(int i=;i<=;i++){
        MAXN=i;
        make_ans();
    }

//    return 0;
    int t;
    scanf("%d",&t);
    while(t--){
        int n,k;
        scanf("%d %d",&n,&k);
        if(n>=){
            printf("%d",n);
            Cantor_invexp(tmp,n-,1LL*k-);
            for(int i=;i<=n-;i++){
                printf(" %d",tmp[i]);
            }
            printf("\n");
        }else{
            Cantor_invexp(tmp,n,1LL*ans[n][k]-);
            for(int i=;i<=n;i++){
                printf("%d",tmp[i]);
                if(i<n)printf(" ");
                else printf("\n");
            }
        }
    }
    return ;
}