codeforces840D Destiny -- 可持久化线段树

%%%Vectorxj
对于询问$(l,r,k)$，令$p={l+r\over k}$。假如将这个区间排序，那么对于这个区间所有出现次数大于 $p$ 的数，它一定在 $p,2p\cdots kp$ 出现过。于是只要用主席树求出区间第 $p,2p\cdots kp$ 大，再判断一下出现次数是否大于 p <script type="math/tex" id="MathJax-Element-285">p</script> 就可以了。

代码

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

using namespace std;

const int N=300010;

int n,m,cnt;
int a[N],rt[N],tot[N*80],ls[N*80],rs[N*80];

void Build(int &g,int l,int r){
    g=++cnt; tot[g]=0;
    if(l==r) return ;
    int mid=l+r>>1;
    Build(ls[g],l,mid); Build(rs[g],mid+1,r);
}

void Add(int &g1,int g2,int x,int l,int r){
    g1=++cnt; ls[g1]=ls[g2]; rs[g1]=rs[g2]; tot[g1]=tot[g2]+1;
    if(l==r) return ;
    int mid=l+r>>1;
    if(x<=mid) Add(ls[g1],ls[g2],x,l,mid);
    else Add(rs[g1],rs[g2],x,mid+1,r);
}

int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&a[i]);
    Build(rt[0],1,n);
    for(int i=1;i<=n;i++)
        rt[i]=rt[i-1],Add(rt[i],rt[i],a[i],1,n);
    while(m--){
        int l,r,k,ans=-1;
        scanf("%d%d%d",&l,&r,&k);
        for(int i=1;i<=k;i++){
            int L=1,R=n,g1=rt[r],g2=rt[l-1],now=i*(r-l+1)/k;
            while(L<R){
                int mid=L+R>>1;
                if(tot[ls[g1]]-tot[ls[g2]]>=now)
                    g1=ls[g1],g2=ls[g2],R=mid;
                else
                    g1=rs[g1],g2=rs[g2],L=mid+1,now-=tot[ls[g1]]-tot[ls[g2]];
            }
            if((tot[g1]-tot[g2])*k>(r-l+1) && (L<ans || ans==-1)) ans=L;
        }
        printf("%d\n",ans);
    }
    return 0;
}