一眼看status吓得我还以为是分块。。。其实是水题。
考虑给定一个集合,如何得到FRBSUM。
按某种顺序从集合中取数,设之前取的数可以取到的最大值为
mx
(雾)且
[1,mx]
都可以取到,然后现在新加了一个数
x
,如果
至此我们可以得到一个算法,设当前查询集合为S,先从1开始,判断S内小于1的数的和,然后再用这个和+1扔进S内查询一下,直至查询得到的结果不变。
于是就可以很愉快地用主席树乱搞了。
由于每次倍增最多增加一倍,于是查询次数是
O(logn)
的,每次查询的时间复杂度是
O(logn)
。总时间复杂度是
O(nlog2n)
,空间复杂度
O(nlogSmax)
。
#include <bits/stdc++.h>
#define rep(i,a,b) for(int i = a , _ = b ; i <= _ ; i ++)
#define per(i,a,b) for(int i = a , _ = b ; i >= _ ; i --)
#define maxn 100007
#define maxs 4000007
inline int rd() {
char c = getchar();
while (!isdigit(c)) c = getchar() ; int x = c - '0';
while (isdigit(c = getchar())) x = x * 10 + c - '0';
return x;
}
inline void upmax(int&a , int b) { if (a < b) a = b ; }
typedef int arr[maxn];
typedef int seg[maxs];
arr rt , a;
seg lc , rc , sum;
int tot , n , m , lim;
void update(int pr , int&nr , int l , int r , int p) {
if (!nr) nr = ++ tot;
sum[nr] = sum[pr] + p;
if (l == r) return;
int m = l + r >> 1;
if (p <= m)
rc[nr] = rc[pr] , update(lc[pr] , lc[nr] , l , m , p);
else
lc[nr] = lc[pr] , update(rc[pr] , rc[nr] , m + 1 , r , p);
}
void input() {
n = rd();
rep (i , 1 , n) a[i] = rd() , upmax(lim , a[i]);
rep (i , 1 , n) update(rt[i - 1] , rt[i] , 1 , lim , a[i]);
}
int query(int pr , int nr , int l , int r , int p) {
if (l > p) return 0;
if (!(sum[nr] - sum[pr])) return 0;
if (r <= p) return sum[nr] - sum[pr];
int m = l + r >> 1;
int ret = query(lc[pr] , lc[nr] , l , m , p);
if (p > m) ret += query(rc[pr] , rc[nr] , m + 1 , r , p);
return ret;
}
inline int query(int l , int r) {
for (int t = 1 , x;;t = x + 1)
if ((x = query(rt[l] , rt[r] , 1 , lim , t)) < t)
return t;
}
void solve() {
m = rd();
rep (i , 1 , m) {
int l = rd() , r = rd();
int ans = query(l - 1 , r);
printf("%d\n" , ans);
}
}
int main() {
#ifndef ONLINE_JUDGE
freopen("data.txt" , "r" , stdin);
#endif
input();
solve();
return 0;
}