BZOJ3173 TJOI2013最长上升子序列(Treap+ZKW线段树)

传送门

Description

给定一个序列，初始为空。现在我们将1到N的数字插入到序列中，每次将一个数字插入到一个特定的位置。每插入一个数字，我们都想知道此时最长上升子序列长度是多少？

Input

第一行一个整数N，表示我们要将1到N插入序列中，接下是N个数字，第k个数字Xk，表示我们将k插入到位置Xk（0<=Xk<=k-1,1<=k<=N）

Output

N行，第i行表示i插入Xi位置后序列的最长上升子序列的长度是多少。

Sample Input

3

0 0 2

Sample Output

1

1

2

HINT

100%的数据 n<=100000

因为每一次加进来的都是最大的值，所以是不会更新其他的答案的，所以我们可以先把序列搞出来，离线乱搞搞就好了，没事就试了一下zkw线段树，发现还挺快的，不过还是要比树状数组慢不少。

/**************************************************************
Problem: 3173
User: geng4512
Language: C++
Result: Accepted
Time:644 ms
Memory:4712 kb
****************************************************************/

#include<cstdio>
#define MAXN 100005
struct node { int c[2], sz, rnd; } t[MAXN];
unsigned sd = 2333;
int Sz, a[MAXN], cnt, n, rt, mx[MAXN<<2], M, ans[MAXN];
inline int Max(int a, int b) { return a > b ? a : b; }
inline void GET(int &n) {
n = 0; char c;
do c = getchar(); while('0' > c || c > '9');
do n = n * 10 + c - '0', c = getchar(); while('0' <= c && c <= '9');
}
inline unsigned Ran() { return sd = sd * sd + 12580; }
inline void Upd(int k) { t[k].sz = t[t[k].c[0]].sz + t[t[k].c[1]].sz + 1; }
inline void Rot(int &k, bool f) {
int tmp = t[k].c[f]; t[k].c[f] = t[tmp].c[!f]; t[tmp].c[!f] = k;
Upd(k); Upd(tmp); k = tmp;
}
void Insert(int &k, int sz) {
if(!k) { k = ++ Sz; t[k].sz = 1; t[k].rnd = Ran(); return; }
++ t[k].sz;
bool f = t[t[k].c[0]].sz < sz;
Insert(t[k].c[f], sz - (t[t[k].c[0]].sz+1)*f);
if(t[k].rnd > t[t[k].c[f]].rnd) Rot(k, f);
}
void dfs(int u) {
if(!u) return;
dfs(t[u].c[0]);
a[++ cnt] = u;
dfs(t[u].c[1]);
}
namespace Seg {
void Insert(int x, int v) {
for(mx[x += M] = v, x >>= 1; x; x >>= 1)
mx[x] = Max(mx[x<<1], mx[x<<1|1]);
}
int Query(int r) {
int ans = 0;
for(r += M + 1; r ^ 1; r >>= 1)
if(r & 1) ans = Max(ans, mx[r^1]);
return ans;
}
}
int main() {
GET(n); int t;
for(int i = 1; i <= n; ++ i) {
GET(t);
Insert(rt, t);
}
dfs(rt);
for(M = 1; M < n + 2; M <<= 1);
for(int i = 1, tmp; i <= n; ++ i) {
tmp = Seg::Query(a[i]) + 1;
Seg::Insert(a[i], tmp);
ans[a[i]] = tmp;
}
for(int i = 1; i <= n; ++ i) {
ans[i] = Max(ans[i], ans[i-1]);
printf("%d\n", ans[i]);
}
return 0;
}