mex（主席树）

题意

有一个长度为$n$的数组 { a 1 , a 2 , … , a n } \{a_1,a_2,\dots, a_n\} {a1​,a2​,…,an​}。$m$次询问，每次询问一个区间内最小没有出现过的自然数。

数据范围

$1\le n,m\le 2\times 10^5$
$0\le a_i\le 10^9$

思路

很经典的一道应用主席树的题目。每个区间记录这个区间所有数最近出现位置的最小值，如果说某个数小于等于这个最小值，那么这个区间所有数都在那个数之后出现。
对于查询区间$[l,r]$，我们搜索$r$代表的那个权值线段树，如果左子树存储的最小值小于$l$，那么答案一定在左子树中；反之，一定在右子树中。
这里需要注意一点，就是离散化的时候应该把$a_i+1$以及$0$也放进去。原因是保证区间的连续性。即如果一个区间没有出现$0$，那么结果一定是$0$；如果$0\sim a_i$都出现了，但是$a_i+1$未出现，那么答案一定是$a_i+1$；

代码

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>

using namespace std;

const int N = 200010 * 2;

int n, m;
int a[N];
vector<int> nums;

int root[N], idx;

struct Node
{
    int l, r;
    int val;
}tr[N * 22];

int find(int x)
{
    return lower_bound(nums.begin(), nums.end(), x) - nums.begin();
}

int build(int l, int r)
{
    int p = ++ idx;
    if(l == r) return p;
    int mid = l + r >> 1;
    tr[p].l = build(l, mid), tr[p].r = build(mid + 1, r);
    return p;
}

int insert(int p, int l, int r, int x, int id)
{
    int q = ++ idx;
    tr[q] = tr[p];
    if(l == r) {
        tr[q].val = id;
        return q;
    }
    int mid = l + r >> 1;
    if(x <= mid) tr[q].l = insert(tr[p].l, l, mid, x, id);
    else tr[q].r = insert(tr[p].r, mid + 1, r, x, id);
    tr[q].val = min(tr[tr[q].l].val, tr[tr[q].r].val);
    return q;
}

int query(int p, int l, int r, int x)
{
    if(l == r) return r;
    int mid = l + r >> 1;
    if(tr[tr[p].l].val < x) return query(tr[p].l, l, mid, x);
    else return query(tr[p].r, mid + 1, r, x);
}

int main()
{
    scanf("%d%d", &n, &m);
    for(int i = 1; i <= n; i ++) {
        scanf("%d", &a[i]);
        nums.push_back(a[i]);
        nums.push_back(a[i] + 1);
    }
    nums.push_back(0);
    sort(nums.begin(), nums.end());
    nums.erase(unique(nums.begin(), nums.end()), nums.end());
    root[0] = build(0, nums.size() - 1);
    for(int i = 1; i <= n; i ++) {
        root[i] = insert(root[i - 1], 0, nums.size() - 1, find(a[i]), i);
    }
    while(m --) {
        int l, r;
        scanf("%d%d", &l, &r);
        printf("%d\n", nums[query(root[r], 0, nums.size() - 1, l)]);
    }
    return 0;
}