【bzoj 2724】蒲公英（分块）

传送门biu~
分块，预处理f(i,j)为第i块到第j块的众数。每次查询区间众数时，可能作为答案的种类只有完整的块中的众数和不完整的块中的数最多$2*\sqrt n+1$种。二分求一下区间中出现的次数就可以了。

#include<bits/stdc++.h>
using namespace std;
int n,m,lastans,p;
int blocksiz,maxblock,block[40005];
int a[40005],b[40005],f[205][205],num[40005];
vector<int>color[40005];
inline int find(int x){return lower_bound(b+1,b+p+1,x)-b;}
inline int Query(int L,int R,int x){return upper_bound(color[x].begin(),color[x].end(),R)-lower_bound(color[x].begin(),color[x].end(),L);}
int main(){
    scanf("%d%d",&n,&m);
    blocksiz=sqrt(n);
    for(int i=1;i<=n;++i){
        scanf("%d",&a[i]);b[i]=a[i];
        block[i]=(i-1)/blocksiz+1;
    }
    maxblock=block[n];
    sort(b+1,b+n+1);
    p=unique(b+1,b+n+1)-b-1;
    for(int i=1;i<=n;++i){
        a[i]=find(a[i]);
        color[a[i]].push_back(i);
    }
    for(int i=1;i<=maxblock;++i){
        memset(num,0,sizeof(num));
        int now=0;
        for(int j=(i-1)*blocksiz+1;j<=n;++j){
            ++num[a[j]];
            if(num[a[j]]>num[now] || (num[a[j]]==num[now]&&a[j]<now) )  now=a[j];
            f[i][block[j]]=now;
        }
    }
    for(int i=1;i<=m;++i){
        int x,y;
        scanf("%d%d",&x,&y);
        x=(x+lastans-1)%n+1;y=(y+lastans-1)%n+1;
        if(x>y)     swap(x,y);
        int L=block[x]+1,R=block[y]-1;
        int ans=0,Max=-1;
        if(L<=R){
            ans=f[L][R];Max=Query(x,y,ans);
            for(int j=x;j<=(L-1)*blocksiz;++j){
                int tmp=Query(x,y,a[j]);
                if(tmp>Max || (tmp==Max&&ans>a[j]) )    ans=a[j],Max=tmp;
            }
            for(int j=R*blocksiz+1;j<=y;++j){
                int tmp=Query(x,y,a[j]);
                if(tmp>Max || (tmp==Max&&ans>a[j]) )    ans=a[j],Max=tmp;
            }
        }
        else{
            for(int j=x;j<=y;++j){
                int tmp=Query(x,y,a[j]);
                if(tmp>Max || (tmp==Max&&ans>a[j]) )    ans=a[j],Max=tmp;
            }
        }
        lastans=b[ans];
        printf("%d\n",lastans);
    }
    return 0;
}