[F][5]bzoj 2724: [Violet 6]蒲公英

题目就不粘了…
题意：在线求区间众数

---

【分析】
在线区间众数的分块做法比较多，这里提供一个思路：

首先离散化一下比较方便。

最初可能会有一个想法，是不是众数只可能是完整的块的众数，或者不完整的块出现的数呢？显然很容易得出反例。

应该是完整的所有块的众数，和不完整块中出现的数。

所以我们可以预处理f(i,j)表示第 i 块到第 j 块的众数（枚举 i 开个桶扫一遍）。

桶的话可以用map！超好用！不MLE！棒极了！

那么只要能快速得出一个数在某个区间内出现次数即可，每次只要比较至多2√n+1个元素的出现次数，这题就解决了。

由于没有修改，只要离散化以后，给每个数 x 开个vector，按顺序存下 x 出现的位置，每次询问 x 时把区间的左右端点放进对应 vector 二分一下即可。

根据均值不等式，可以算出分块大小大概是√(n/logn)

---

【代码】

//bzoj 2724 [Violet 6]蒲公英#include<algorithm>#include<iostream>#include<cstring>#include<cstdio>#include<vector>#include<cmath>#include<map>#define M(a) memset(a,0,sizeof a)#define fo(i,j,k) for(i=j;i<=k;i++)using namespace std;const int mxn=50005;vector <int> f[mxn];int n,m,id,x;map <int,int> mp;int a[mxn],b[mxn],be[mxn],res[500][500];int lef[500],rig[500],vis[mxn];inline int read(){     int x=0,f=1;char ch=getchar();    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}    return x*f;}inline void yu(int s){    int i,j,k,sc,now,mx=0,tmp=0;    M(vis);    fo(i,lef[s],n)    {        int t=be[i];        vis[a[i]]++;        if(vis[a[i]]>mx || (vis[a[i]]==mx && b[a[i]]<b[tmp]))          mx=vis[a[i]],tmp=a[i];        if(i==rig[t])          res[s][t]=tmp;    }}inline int calc(int l,int r,int i){    return upper_bound(f[i].begin(),f[i].end(),r)-lower_bound(f[i].begin(),f[i].end(),l);}inline int query(int l,int r){    int i,j,L=be[l],R=be[r];    if(L==R || R==L+1)    {        M(vis);        int mx=0,tmp=0;        fo(i,l,r)        {          vis[a[i]]++;          if(vis[a[i]]>mx || (vis[a[i]]==mx && b[a[i]]<tmp))            tmp=b[a[i]],mx=vis[a[i]];      }      return tmp;    }    int k=res[L+1][R-1],mx=0,tmp=0,rlt,czy;    mx=calc(l,r,k),tmp=b[k];    fo(i,l,rig[L])    {      czy=calc(l,r,a[i]);      if(czy>mx || (czy==mx && b[a[i]]<tmp))        mx=czy,tmp=b[a[i]];    }    fo(i,lef[R],r)    {      czy=calc(l,r,a[i]);      if(czy>mx || (czy==mx && b[a[i]]<tmp))        mx=czy,tmp=b[a[i]];    }    return tmp;}int main(){    int i,j,k,l,r,ans=0;    n=read(),m=read();    fo(i,1,n)    {        a[i]=read();        if(!mp[a[i]])        {            mp[a[i]]=++id;            b[id]=a[i];        }        a[i]=mp[a[i]];        f[a[i]].push_back(i);           }    x=sqrt(n);    fo(i,1,x)    {        lef[i]=rig[i-1]+1;        rig[i]=i*x;    }    rig[x]=n;    fo(i,1,x)      fo(j,lef[i],rig[i])        be[j]=i;    fo(i,1,x) yu(i);    while(m--)    {        l=read(),r=read();        l=(l+ans-1)%n+1,r=(r+ans-1)%n+1;        if(l>r) swap(l,r);        ans=query(l,r);        printf("%d\n",ans);    }    return 0;}