BZOJ 2724: [Violet 6]蒲公英

补点题解 是不是太久没写blog了。。
在luogu月赛上看到这道题（原来莫队就好了嘛。那题是离线的。。而且只有次数）
找了一下 这道题原来就是经典的区间众数问题啊 bzoj这道是强制在线的

大家去膜陈老师（WJMZBMR）写的题解吧
带修改的我还不会呢。。。

主要思想是分块
由于数据较小 可以用 $O(n^{1.5}logn)$的复杂度搞过去
就是把不完整的块内的数进行二分求一下出现次数搞一下即可

$n \sqrt n$ 的我懒得打了 其实就是多一个预处理
PS:分块貌似分成 $\sqrt {n/logn}$会快一点 均值不等式？。

#include<bits/stdc++.h>
#define me(a,x) memset(a,x,sizeof a)
using namespace std;
const int N=40002,M=790;
inline int read(){
    char ch=getchar(); int x=0,f=1;
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0'; ch=getchar();}
    return x*f;
}
struct P{
    int x,y;
}f[M][M];
P max(P a,P b){
    if(a.x==b.x) return a.y<b.y?a:b;
    return a.x>b.x?a:b;
}
int bl[N],a[N],z[N],c[N],n,k;
vector<int>v[N];
void gao(int x){
    me(c,0); P u=(P){0,0};
    for(int i=x*k-k+1;i<=n;i++){
        c[a[i]]++; int t=bl[i];
        u=max(u,(P){c[a[i]],a[i]});
        f[x][t]=u;
    }
}
int get(int l,int r,int x){
    return upper_bound(v[x].begin(),v[x].end(),r)-lower_bound(v[x].begin(),v[x].end(),l);
}
int main()
{
    n=read(); int m=read(),i,j,u;
    k=sqrt(double(n)*log(2)/log(n));
    for(i=1;i<=n;i++)bl[i]=(i-1)/k+1;
    for(i=1;i<=n;i++)a[i]=z[i]=read();
    sort(z+1,z+1+n); u=unique(z+1,z+1+n)-z-1;
    for(i=1;i<=n;i++)
        a[i]=lower_bound(z+1,z+1+u,a[i])-z,v[a[i]].push_back(i);
    for(i=1;i<=bl[n];i++) gao(i);
    int last=0;
    while(m--){
        int l=(read()+last-1)%n+1,r=(read()+last-1)%n+1;
        if(l>r)l^=r^=l^=r; P ans=(P){0,0};
        if(bl[l]<bl[r]-1)ans=f[bl[l]+1][bl[r]-1];
        for(i=l;i<=r && bl[i]==bl[l];i++)ans=max(ans,(P){get(l,r,a[i]),a[i]});
        if(bl[l]!=bl[r])
            for(i=r;bl[i]==bl[r];i--)ans=max(ans,(P){get(l,r,a[i]),a[i]});
        printf("%d\n",last=z[ans.y]);
    }
    return 0;
}