bzoj3301: [USACO2011 Feb] Cow Line

新姿势康托展开。。

------------------------------------

裸的康托展开&逆康托展开

　　康托展开就是一种特殊的hash，且是可逆的……

　　康托展开计算的是有多少种排列的字典序比这个小，所以编号应该+1；逆运算同理（-1）。

　　序列->序号：（康托展开）

　　　　对于每个数a[i]，数比它小的数有多少个在它之前没出现，记为b[i]，ans=1+∑b[i]∗(n−i)!ans=1+∑b[i]∗(n−i)!

　　序号->序列：（逆康托展开）

　　　　求第x个排列所对应的序列，先将x-1，然后对于a[i]，⌊x(n−i)!⌋⌊x(n−i)!⌋即为在它之后出现的比它小的数的个数，所以从小到大数一下有几个没出现的数，就知道a[i]是第几个数了。

------------------------------------

#include<bits/stdc++.h>
#define rep(i,l,r) for(int i=l;i<=r;++i)
using namespace std;
const int N=25;
int n,Q,a[N],f[N];
typedef long long ll;
char s[10];
ll ans,sum,fac[N],x;
int main(){
    scanf("%d%d",&n,&Q);
    fac[1]=1; fac[0]=1;
    rep(i,2,n) fac[i]=fac[i-1]*i;
    while(Q--){
        scanf("%s",s);
        if(s[0]=='P'){
            scanf("%lld",&x);--x;
            memset(f,0,sizeof f);
              rep(i,1,n){
                int t=x/fac[n-i],j,k;
                for(k=1,j=0;j<=t;k++) if (!f[k]) j++;
                   f[k-1]=1; a[i]=k-1;
                  x%=fac[n-i];
              }
              rep(i,1,n-1) printf("%d ",a[i]);
               printf("%d
",a[n]);
        }else {
            memset(f,0,sizeof f); ans=1; 
            rep(i,1,n) {
                scanf("%lld",&x);
                f[x]=1;
                sum=0;
                rep(j,1,x-1) sum+=(!f[j]);
                ans+=sum*fac[n-i];
            }
            printf("%lld
",ans);
        }
    }
}