2018.09.14 洛谷P3567 [POI2014]KUR-Couriers（主席树）

传送门

简单主席树啊。

但听说有随机算法可以秒掉%%%（本蒟蒻并不会）

直接维护值域内所有数的出现次数之和。

当这个值不大于区间总长度的一半时显然不存在合法的数。

这样在主席树上二分查值就行了。

代码：

#include<bits/stdc++.h>
#define N 500005
using namespace std;
inline int read(){
    int ans=0;
    char ch=getchar();
    while(!isdigit(ch))ch=getchar();
    while(isdigit(ch))ans=(ans<<3)+(ans<<1)+(ch^48),ch=getchar();
    return ans;
}
int n,m,tot=0,tmp,son[N*20][2],siz[N*20],rt[N],ql,qr;
inline int max(int a,int b){return a>b?a:b;}
inline void insert(int&p,int las,int l,int r,int v){
    son[p=++tot][0]=son[las][0],son[p][1]=son[las][1];
    if(l==r){siz[p]=siz[las]+1;return;}
    int mid=l+r>>1;
    if(v<=mid)insert(son[p][0],son[las][0],l,mid,v);
    else insert(son[p][1],son[las][1],mid+1,r,v);
    siz[p]=siz[son[p][0]]+siz[son[p][1]];
}
inline int query(int a,int b,int l,int r){
    if(l==r)return l;
    int mid=l+r>>1;
    if(2*(siz[son[b][0]]-siz[son[a][0]])>qr-ql+1)return query(son[a][0],son[b][0],l,mid);
    if(2*(siz[son[b][1]]-siz[son[a][1]])>qr-ql+1)return query(son[a][1],son[b][1],mid+1,r);
    return 0;
}
int main(){
    n=read(),m=read();
    for(int i=1;i<=n;++i)insert(rt[i],rt[i-1],1,n,(tmp=read()));
    while(m--){
        ql=read(),qr=read();
        printf("%d\n",query(rt[ql-1],rt[qr],1,n));
    }
    return 0;
}