康拓展开和它的逆运算

康拓展开可以看做是一种hash函数，它可以把一个数列转化成一个数，这个数就是这个数列在它的全排列中处于第几位（按照字典序从小到大排列，编号从0开始）公式是 X=a[0]*(n-1)!+a[1]*(n-2)!+...+a[i]*(n-i-1)!+...+a[n-2]*1!+a[n-1]*0! 其中，a[i]为所给的数列中第 i （从0开始）个数在未出现的元素中是排在第几个（从0开始），如果存在重复的元素，那么把所有相同的元素看做一个。

康拓展开的逆运算是求出一个数列的全排列中的第 m （给定的）个排列。也就是说，已知 m = a[0]*(n-1)!+a[1]*(n-2)!+...+a[i]*(n-i-1)!+...+a[n-2]*1!+a[n-1]*0!，让确定a[0],a[1]……a[n-1]分别是多少。利用辗转相除法可以实现，比如当n=4,m=20的时候，有：

(其实应该是a[3]、a[2]、a[1]、a[0])

然后依次在数列中未出现过的元素中找地a[i]大的数就好。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

typedef long long ll;

ll fac[21];
int a[25];
int b[25];
bool f[25];

void init()
{
    fac[0] = 1;
    for (int i=1; i<25; i++)
        fac[i] = fac[i-1]*i;
}

ll contor(int *a, int n)  //康拓展开
{
    memset(f,false,sizeof(f));
    ll sum = 0;
    for (int i=0; i<n; i++)
    {
        memset(f,false,sizeof(f));
        int cnt = 0;
        for (int j=i+1; j<n; j++)
            if (a[j] < a[i] && !f[a[j]])
        {
            cnt++;
            f[a[j]] = true;
        }
        sum += cnt * fac[n-i-1];
    }

    return sum;
}

void inv_contor(int m, int n)  //康拓展开的逆运算
{
    memset(f,false,sizeof(f));
    memset(b,0,sizeof(b));
    for (int i=n-1; i>=0; i--)
    {
        int ai = m / fac[i];
        m = m % fac[i];
        int cnt = 0;
        int now = -1;
        for (int j=0; j<n; j++)
            if (!f[j] && a[j] > now)
        {
            if (cnt == ai)
            {
                b[n-i-1] = a[j];
                f[j] = true;
                break;
            }
            cnt++;
            now = a[j];
        }
    }

    for (int i=0; i<n; i++)
        printf("%d ",b[i]);
    printf("\n");
}

int main()
{
    int n, m;

    init();

    while (scanf("%d",&n)!=EOF)
    {
        for (int i=0; i<n; i++) scanf("%d",&a[i]);
        printf("%lld\n",contor(a,n));

        sort(a,a+n);

        scanf("%d",&m);
        inv_contor(m, n);
    }


    return 0;
}