bzoj3585: mex

链接

　　http://www.lydsy.com/JudgeOnline/problem.php?id=3585

题解

　　达成成就：无调试一遍过样例一遍$AC$。
　　这也说明了这道题多么$SB$。
　　肯定得先上莫队。下面考虑怎样维护$mex$，有个性质很显然，就是我们的$mex$不会大于$N$。
　　因为最极端的情况就是$N$个数分别是从$0$到$N-1$，只有这个时候$mex$才会取到$N$，其它时候一定都是小于$N$的。对权值进行分块，用一个权值数组$cnt[]$表示每个数字目前的出现次数，每个块则记录这个块内多少个$cnt$是非零的。
　　查询直接暴力$O(\sqrt N)$查询就好了。

代码

//莫队+分块
#include <cstdio>
#include <algorithm>
#define maxn 200010
using namespace std;
int N, M, num[maxn], q[maxn][3], size, block[maxn], cnt[maxn], a[maxn];
inline int read(int x=0)
{
    char c=getchar();
    while(c<48 or c>57)c=getchar();
    while(c>=48 and c<=57)x=(x<<1)+(x<<3)+c-48, c=getchar();
    return x;
}
inline void add(int &x){if(x>=N)return;if(++cnt[x]==1)block[x>>size]++;}
inline void del(int &x){if(x>=N)return;if(--cnt[x]==0)block[x>>size]--;}
inline int quiry()
{
    int i, j;
    for(i=0;i<=N>>size;i++)
        if(block[i]!=1<<size)for(j=i<<size;;j++)if(!cnt[j])return j;
}
inline bool cmp(int a, int b)
{return q[a][0]>>size==q[b][0]>>size?q[a][1]<q[b][1]:q[a][0]>>size<q[b][0]>>size;}
void init()
{
    int i;
    N=read(), M=read();
    for(i=1;i<=N;i++)a[i]=read();for(size=1;(1<<size)*(1<<size)<N;size++);
    for(i=1;i<=M;i++)q[i][0]=read(), q[i][1]=read(), num[i]=i;
    sort(num+1,num+N+1,cmp);
}
void work()
{
    int l=0, r=2, i, x, y;
    add(a[1]);
    for(i=1;i<=M;i++)
    {
        x=q[num[i]][0], y=q[num[i]][1];
        while(l>=x)add(a[l]),l--;
        while(r<=y)add(a[r]),r++;
        while(l<x-1)l++,del(a[l]);
        while(r>y+1)r--,del(a[r]);
        q[num[i]][2]=quiry();
    }
    for(i=1;i<=M;i++)printf("%d\n",q[i][2]);
}
int main()
{
    init();
    work();
    return 0;
}