[BZOJ3585]mex(莫队+分块)

显然可以离线主席树，这里用莫队+分块做。分块的一个重要思想是实现修改与查询时间复杂度的均衡，这里莫队和分块互相弥补。

考虑暴力的分块做法，首先显然大于n的数直接忽略，于是将值域分成sqrt(n)份，每块记录块内的所有值是否在此当前区间内都已存在。

这样每次暴力从L到R分别放入这个表，最后从小到大询问每个块是否已满，若没有则在块内枚举第一个不存在的数。

注意到这样的总修改复杂度O(nq)，查询复杂度O(qsqrt(n))。

考虑莫队，将序列分成sqrt(n)份，使总修改复杂度变为O(nsqrt(n))。查询复杂度不变O(qsqrt(n))。

 #include<cstdio>
 #include<algorithm>
 #define rep(i,l,r) for (int i=(l); i<=(r); i++)
 using namespace std;

 const int N=,K=;
 int n,m,B,a[N],b[N],ans[N],cnt[K],s[K][K];
 struct P{ int l,r,id; }q[N];

 bool cmp(const P &x,const P &y){ return b[x.l]!=b[y.l] ? b[x.l]<b[y.l] : x.r<y.r; }

 void add(int x){
     if (x>n) return;
     int t=x%B; s[b[x]][t]++; if (s[b[x]][t]==) cnt[b[x]]++;
 }

 void del(int x){
     if (x>n) return;
     int t=x%B; s[b[x]][t]--; if (s[b[x]][t]==) cnt[b[x]]--;
 }

 int Que(){
     rep(i,,n/B+){
         int t=min(n,i*B-)-(i-)*B+;
         if (cnt[i]==t) continue;
         rep(j,(i-)*B,min(n,i*B-)) if (!s[i][j%B]) return j;
     }
     return n+;
 }

 int main(){
     freopen("bzoj3585.in","r",stdin);
     freopen("bzoj3585.out","w",stdout);
     scanf("%d%d",&n,&m); B=; b[]=;
     rep(i,,n) scanf("%d",&a[i]),b[i]=i/B+;
     rep(i,,m) scanf("%d%d",&q[i].l,&q[i].r),q[i].id=i;
     sort(q+,q+m+,cmp); int L=,R=;
     rep(i,,m){
         while (R<q[i].r) R++,add(a[R]);
         while (L>q[i].l) L--,add(a[L]);
         while (R>q[i].r) del(a[R]),R--;
         while (L<q[i].l) del(a[L]),L++;
         ans[q[i].id]=Que();
     }
     rep(i,,m) printf("%d\n",ans[i]);
     return ;
 }