【分块】bzoj2724 [Violet 6]蒲公英

分块，离散化，预处理出：

①前i块中x出现的次数(差分)；

②第i块到第j块中的众数是谁，出现了多少次。

询问的时候，对于整块的部分直接获得答案；对于零散的部分，暴力统计每个数出现的次数，加上差分的结果，尝试更新ans。

#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
using namespace std;
int n,m,sum,sz,num[40001],l[220],r[220],plv[40001][220],mode[220][220],mplv[220][220];
int a[40001],en,Time[40001],x,y,ma[40001],ans;
struct Point{int v,p;}b[40001];
bool operator < (const Point &a,const Point &b){return a.v<b.v;}
int Res,Num;char C,CH[12];
inline int G()
{
    Res=0;C='*';
    while(C<'0'||C>'9')C=getchar();
    while(C>='0'&&C<='9'){Res=Res*10+(C-'0');C=getchar();}
    return Res;
}
inline void P(int x)
{
    Num=0;if(!x){putchar('0');puts("");return;}
    while(x>0)CH[++Num]=x%10,x/=10;
    while(Num)putchar(CH[Num--]+48);
    puts("");
}
void makeblock()
{
    sz=(int)sqrt((double)n); if(!sz) sz=1;
    for(sum=1;sum*sz<n;sum++)
      {
        l[sum]=r[sum-1]+1;
        r[sum]=sum*sz;
        for(int i=l[sum];i<=r[sum];i++) num[i]=sum;
      }
    l[sum]=r[sum-1]+1;
    r[sum]=n;
    for(int i=l[sum];i<=r[sum];i++) num[i]=sum;
}
void LiSan()
{
    sort(b+1,b+n+1);
    for(int i=1;i<=n;i++)
      {
        if(b[i].v!=b[i-1].v) en++;
        ma[a[b[i].p]=en]=b[i].v;
      }
}
void makeplv()
{
    for(int i=1;i<=n;i++)
      for(int j=num[i];j<=sum;j++)
        plv[a[i]][j]++;
}
void makemode()
{
    for(int i=1;i<=sum;i++)
      {
        memset(Time,0,sizeof(Time));
        int modenow,modeplv=0;
        for(int j=i;j<=sum;j++)
          {
            for(int k=l[j];k<=r[j];k++)
              {
                Time[a[k]]++;
                if(Time[a[k]]>modeplv||(Time[a[k]]==modeplv&&a[k]<modenow))
                  {
                    modenow=a[k];
                    modeplv=Time[a[k]];
                  }
              }
            mode[i][j]=modenow;
            mplv[i][j]=modeplv;
          }
      } memset(Time,0,sizeof(Time));
}
int Getplv(const int &v,const int &L,const int &R){return plv[v][R]-plv[v][L-1];}
int main()
{
    n=G(); m=G();
    for(int i=1;i<=n;i++) {b[i].v=G(); b[i].p=i;}
    makeblock(); LiSan(); makeplv(); makemode();
    for(int i=1;i<=m;i++)
      {
        x=G(); y=G(); x=(x+ans-1)%n+1; y=(y+ans-1)%n+1;
        if(x>y) swap(x,y);
        int modenow,modeplv=0;
        if(num[x]+1>=num[y])
          {
            for(int j=x;j<=y;j++)
              {
                Time[a[j]]++;
                if(Time[a[j]]>modeplv||(Time[a[j]]==modeplv&&a[j]<modenow))
                  {
                    modenow=a[j];
                    modeplv=Time[a[j]];
                  }
              }
            for(int j=x;j<=y;j++) Time[a[j]]--;
          }
        else
          {
            modenow=mode[num[x]+1][num[y]-1];
            modeplv=mplv[num[x]+1][num[y]-1];
            for(int j=x;j<=r[num[x]];j++)
              {
                Time[a[j]]++; int t=Time[a[j]]+Getplv(a[j],num[x]+1,num[y]-1);
                if(t>modeplv||(t==modeplv&&a[j]<modenow))
                  {
                    modenow=a[j];
                    modeplv=t;
                  }
              }
            for(int j=l[num[y]];j<=y;j++)
              {
                Time[a[j]]++; int t=Time[a[j]]+Getplv(a[j],num[x]+1,num[y]-1);
                if(t>modeplv||(t==modeplv&&a[j]<modenow))
                  {
                    modenow=a[j];
                    modeplv=t;
                  }
              }
            for(int j=x;j<=r[num[x]];j++) Time[a[j]]--;
            for(int j=l[num[y]];j<=y;j++) Time[a[j]]--;
          }
        P(ans=ma[modenow]);
      }
    return 0;
}