2020-11-17

11.16 做题记录

康顿展开与逆康顿展开

学习这个东西的原因是因为codeforces一道题：
https://codeforces.com/contest/1443/problem/E
大概的意思是这样的：对于一个排列，有两种操作，第个是求一个范围内的和，另一个操作是将排列顺序顺延若干次。
啪，很快啊 我就把题读完后一点点思路都没有，因为根据传统思路，全排列的数量级都是大的离谱，所以我年轻人不讲武德，直接翻到了下面的算法tag， brute force？？？ two pointers？？？
这道题数量级虽然是很大很大，但是就是因为特别大，所以有时候会形成一种惯性思维就是不敢往这方面想。仔细想，14位全排列数字达到了8e10的大小，而在这题里，q范围是1e5，x范围是1e5，相乘也不过1e10，所以无论n多大，在这个题的范围里，最多只改变了14位。
所以我们只需要找到全排列与序号的双射函数，就可以很轻松的解决这个问题。
而这个函数就是康拓展开和逆康托展开。
直接给出来函数吧。。。
X=a[n]*(n-1)!+a[n-1]*(n-2)!+…+a[i]*(i-1)!+…+a[1]*0!
其中，a[i]为整数，代表当前未出现的元素中排第几个。
所以逆康托展开也非常简单了，下面直接给出ac代码吧

AC代码

#include<bits/stdc++.h>
using namespace std;
#define rep(i, a, b) for(int i = a; i <= b; i++)
#define dow(i, a, b) for(int i = a; i >= b; i--)
#define maxn (int)2e5 + 5
#define pb push_back
#define ll long long
ll f[maxn];
ll z[20], vis[20], a[20];
void init(int n){
    z[0] = 1;
    rep(i, 1, 15) z[i] = z[i - 1] * i;
    rep(i, 1, n) f[i] = f[i - 1] + i;
}
void solve(){
    int n, q;
    cin >> n >> q;
    init(n);
    ll id = 0;
    rep(i, 1, q){
        int t;
        cin >> t;
        if(t == 1){
            int l, r;
            cin >> l >> r;
            cout << f[r] - f[l - 1] << endl;
        }
        else{
            int x;
            cin >> x;
            int len = min(n, 15);
            id += x;
            for (int j = 1; j <= len; j++) vis[j] = 0;
            int64_t tmp = id;
            for (int j = n - len + 1; j <= n; j++) {
                int k = tmp / z[n - j] + 1;
                tmp %= z[n-j];
                int p = 1;
                while (true) {
                    k -= !vis[p];
                    if (k == 0) break;
                    p++;
                }
                vis[p] = 1;
                a[j - (n - len + 1)] = n - len + p;
            }
            for (int j = n - len + 1; j <= n; j++) {
                f[j] = f[j-1] + a[j - (n - len + 1)];
            }
//            vector<int> r;
//            rep(j, 1, len) r.pb(j);
//            id += (ll)x;
//            ll w = id;
//            vector<int> aws;
//            dow(j, len, 1){
//                int s = w / z[j - 1];
//                w = w % z[j - 1];
//                aws.pb(r[s]);
//                r.erase(r.begin() + s);
//            }
//            rep(j, n - len + 1, n) f[j] = f[j - 1] + aws[j - (n - len + 1)];
//            if(n >= 14){
//                dow(j, 13, 0){
//                    f[n - j] = f[n - j - 1] + aws[13 - j];
//                }
//            }
//            else{
//                dow(j, n - 1, 0){
//                    f[n - j] = f[n - j - 1] + aws[n - 1 - j];
//                }
//            }
        }
    }

}
int main(){
    //int t;
    // cin >> t;
    // while(t--){
    solve();
    // }
}