HDU多校第六场 1002 Nonsense Time —— LIS删点

HDU多校第六场 1002 Nonsense Time —— LIS删点

题目链接：点我啊╭(╯^╰)╮

题目大意：

    长度为 n n n 的排列 p p p，一开始全部冻结
    每次永久释放 p k i p_{k_{i}} pki​​
    求每次释放之后的 L I S LIS LIS

解题思路：

    要查找 x x x 是否在 L I S LIS LIS 中，考虑用树状数组维护
    树状数组 t [ i ] t[i] t[i] 记录的是每个点 L I S LIS LIS 的位置，而不是值

核心：期望 L I S LIS LIS，树状数组维护位置

#include<bits/stdc++.h>
#define rint register int
#define deb(x) cerr<<#x<<" = "<<(x)<<'n';
using namespace std;
typedef long long ll;
const int maxn = 4e5 + 5;
int T, n, a[maxn], k[maxn], ans[maxn];
int l[maxn], r[maxn], used[maxn];
int t[maxn], f[maxn], pre[maxn]; inline int lb(int x){return x & -x;
}inline void update(int x, int p){for(int i=x; i<=n; i+=lb(i))if(f[t[i]] < f[p]) t[i] = p;
}inline int query(int x){int ret = 0;for(int i=x; i; i-=lb(i))if(f[t[i]] > f[ret]) ret = t[i];return ret;
}inline void lis(){for(int i=0; i<=n+1; i++) t[i] = f[i] = used[i] = 0;for(int i=r[0]; i<=n+1; i=r[i]){int q = query(a[i]);f[i] = f[q] + 1;pre[i] = q;update(a[i], i);}for(int i=n+1; i; i=pre[i]) used[i] = 1;
}int main() {scanf("%d", &T);while(T--){scanf("%d", &n);a[n+1] = n + 1;for(int i=1; i<=n; i++) scanf("%d", a+i);for(int i=1; i<=n; i++) scanf("%d", k+i);for(int i=0; i<=n+1; i++) l[i] = i-1, r[i] = i+1;lis();for(int i=n; i; i--){ans[i] = f[n+1] - 1;l[r[k[i]]] = l[k[i]];r[l[k[i]]] = r[k[i]];if(used[k[i]]) lis();}for(int i=1; i<=n; i++) printf("%d%c", ans[i], i<n ? ' ' : 'n');}
}