bzoj3585: mex

题意

求区间的mex

题解

好题啊啊啊啊啊！！
一直以为主席树不能做区间mex的题。。
然而今天发现了十分神的做法

以权值为下标，
minn[i]=j，表示i所代表的权值为[l,r]的区间，最早出现的位置为j
例：0 2 1 0 1 3 2
对于root[7]来说，
[0,3]=1, 0最早出现在第1个位置
[1,3]=2, 2最早出现在第2个位置
查询[l,r]时，在root[r]里查询
如果左子区间的minn>=l，说明权值小的左边一半都最早出现在l以后，就往右孩子找
否则，说明权值小的左边一半有没有出现在l之后的，就往左孩子找
因为是在root[r]里，所以保证不超过右边界

感觉这个做法真的牛逼

#include<cstdio>
#include<cstring>
const int N=200005*25;
int n,m;
int s1[N],s2[N],minn[N];
int root[N];
int num=0;
int mymin (int x,int y){return x<y?x:y;}
void change (int &rt1,int rt2,int L,int R,int x,int y)
{
    if (rt1==0) rt1=++num;
    if (L==R)
    {
        minn[rt1]=y;
        return ;
    }
    int mid=(L+R)>>1;
    if (x<=mid) 
    {
        change(s1[rt1],s1[rt2],L,mid,x,y);
        s2[rt1]=s2[rt2];
    }
    else
    {
        change(s2[rt1],s2[rt2],mid+1,R,x,y);
        s1[rt1]=s1[rt2];
    }
    minn[rt1]=mymin(minn[s1[rt1]],minn[s2[rt1]]);
}
int find (int now,int L,int R,int x)
{
    if (L==R) return L;
    int mid=(L+R)>>1;
    if (minn[s1[now]]>=x) return find(s2[now],mid+1,R,x);
    else return find(s1[now],L,mid,x);
}
int main()
{
    scanf("%d%d",&n,&m);
    for (int u=1;u<=n;u++)
    {
        int x;
        scanf("%d",&x);
        if (x>n) x=n;
        change(root[u],root[u-1],0,n,x,u);
    }
    for (int u=1;u<=m;u++)
    {
        int l,r;
        scanf("%d%d",&l,&r);
        printf("%d\n",find(root[r],0,n,l));
    }
    return 0;
}