小R的树(权限题)

解：考场上爆0了......

回想怎么求两个排列的最长公共子序列。

回想怎么求1~n每个数恰出现两次的两个序列的最长公共子序列。就是每个数替换为它在另一个序列里的出现位置，降序。

所以我们可以把这每个空位都倒序填入m个数，然后暴力，最后输出方案。

考虑优化。发现在每个空位的时候，这m个数都是单降的。直接拿指针在单调栈上扫，可以O(top + m)转移。

关于记录方案，每个位置记录以它结尾的lis中的前一个非-1位置，以及在这之间有多少个-1。

然后就可以把空间优化到O(n)。

#include <bits/stdc++.h>

typedef long long LL;
const int N = 200010;
const LL INF = 0x3f3f3f3f3f3f3f3fll;;

struct Node {
    int a, b;
    Node(int A = 0, int B = 0) {
        a = A;
        b = B;
    }
}frp[N], fr[N];

LL a[N], b[N], p[N];
int top, n, m;
bool vis[N];

int main() {
    scanf("%d", &n);
    for(int i = 1; i <= n; i++) {
        scanf("%lld", &a[i]);
    }
    scanf("%d", &m);
    for(int i = 1; i <= m; i++) {
        scanf("%lld", &b[i]);
    }
    std::sort(b + 1, b + m + 1);
    std::reverse(b + 1, b + m + 1);
    a[n + 1] = INF;
    for(int i = 1; i <= n + 1; i++) {
        if(a[i] != -1) {
            int l = 0, r = top;
            while(l < r) {
                int mid = (l + r + 1) >> 1;
                if(p[mid] < a[i]) l = mid;
                else r = mid - 1;
            }
            fr[i] = frp[r];
            if(r == top) {
                p[++top] = a[i];
                frp[top] = Node(i, 0);
            }
            else if(p[r + 1] > a[i]) {
                p[r + 1] = a[i];
                frp[r + 1] = Node(i, 0);
            }
        }
        else {
            int p1 = top;
            for(int j = 1; j <= m; j++) {
                while(p1 && p[p1] >= b[j]) {
                    p1--;
                }
                Node temp = frp[p1];
                temp.b++;
                if(p1 == top) {
                    p[++top] = b[j];
                    frp[top] = temp;
                }
                else if(p[p1 + 1] > b[j]) {
                    p[p1 + 1] = b[j];
                    frp[p1 + 1] = temp;
                }
            }
        }
    }

    int now = frp[top].a, p1 = n, p2 = 1;
    while(now) {
        if(fr[now].b) {
            while(p1 > now) p1--;
            while(b[p2] >= a[now]) p2++;
            for(int i = 1; i <= fr[now].b; i++) {
                while(a[p1] != -1) p1--;
                a[p1] = b[p2];
                vis[p2] = 1;
                p2++;
            }
        }
        now = fr[now].a;
    }

    p2 = 1;
    for(int i = 1; i <= n; i++) {
        if(a[i] == -1) {
            while(vis[p2]) p2++;
            a[i] = b[p2];
            p2++;
        }
    }

    for(int i = 1; i <= n; i++) {
        printf("%lld ", a[i]);
    }
    return 0;
}