BZOJ 2724 分块统计

区间众数

先离散化，学到了lagoon的lower_bound+unique的离散化，比我写的简单多了

预处理分成sqrt(n)块，记录d[i][j]和p[i][j]分别表示从i块起始位置到j块终止位置的众数出现次数和这个数是谁

开一个数组，记录每个数的位置，使得同类的相邻，同类数的坐标升序排列。

对于询问：

对于同一块内或者相邻块内的，直接暴力。

对于不同块内的，统计出非整块的区间内的所有出现过的数字，然后对这些数字在二分存储位置的数组，统计块内这个数出现的次数。然后结合d[i][j]和p[i][j]就可以得出答案了

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cstdlib>
#include <cmath>

#define N 1000000
#define M 2000

using namespace std;

int sz,tot;
int lt[N],rt[N],ed[M],st[M],sum[N];
int q[N],cs[N],p[M][M],d[M][M];
int n,m;
int val[N],a[N];
int ef[N];
int bx,by;

inline void getblock()
{
    sz=sqrt(n); tot=n/sz;
    for(int i=1;i<=tot;i++)
    {
        st[i]=ed[i-1]+1;
        ed[i]=st[i]+sz-1;
    }
    if(ed[tot]!=n)
    {
        ++tot;
        st[tot]=ed[tot-1]+1;
        ed[tot]=n;
    }
    int cnum,csum;
    for(int i=1;i<=tot;i++)
    {
        csum=0;
        memset(cs,0,sizeof cs);
        for(int j=i;j<=tot;j++)
        {
            for(int k=st[j];k<=ed[j];k++)
            {
                cs[val[k]]++;
                if(cs[val[k]]>csum)
                {
                    csum=cs[val[k]];
                    cnum=val[k];
                }
                else if(cs[val[k]]==csum)
                {
                    cnum=min(cnum,val[k]);
                }
            }
            d[i][j]=csum;
            p[i][j]=cnum;
        }
    }
}

inline void read()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&val[i]);
        a[i]=val[i];
    }
    sort(a+1,a+1+n);
    int num=unique(a+1,a+1+n)-a-1;
    for(int i=1;i<=n;i++)
        val[i]=lower_bound(a+1,a+1+num,val[i])-a;
    for(int i=1;i<=n;i++) cs[val[i]]++;
    for(int i=1;i<=num;i++) sum[i]=sum[i-1]+cs[i];
    for(int i=1;i<=num;i++) lt[i]=sum[i-1]+1,rt[i]=sum[i-1];
    for(int i=1;i<=n;i++) ef[++rt[val[i]]]=i;
    
    getblock();
}

inline int force(int x,int y)
{
    for(int i=x;i<=y;i++) cs[val[i]]=0;
    int cnum,csum=0;
    for(int i=x;i<=y;i++)
    {
        cs[val[i]]++;
        if(cs[val[i]]>csum)
        {
            csum=cs[val[i]];
            cnum=val[i];
        }
        else if(cs[val[i]]==csum)
        {
            cnum=min(cnum,val[i]);
        }
    }
    return cnum;
}

inline int getsum(int ll,int rr,int x)
{
    return upper_bound(ef+rt[x-1]+1,ef+rt[x]+1,rr)-lower_bound(ef+rt[x-1]+1,ef+rt[x]+1,ll);
}


inline int query(int x,int y)
{
    bx=x/sz+(x%sz!=0);
    by=y/sz+(y%sz!=0);
    if(bx==by||bx+1==by) return force(x,y);
    
    int h=0;
    for(int i=x;i<=ed[bx];i++) cs[val[i]]=0;
    for(int i=st[by];i<=y;i++) cs[val[i]]=0;
    for(int i=x;i<=ed[bx];i++)
    {
        cs[val[i]]++;
        if(cs[val[i]]==1) q[++h]=val[i];
    }
    for(int i=st[by];i<=y;i++)
    {
        cs[val[i]]++;
        if(cs[val[i]]==1) q[++h]=val[i];
    }
    bx++; by--;
    int res=d[bx][by],ans=p[bx][by];
    for(int i=1;i<=h;i++)
    {
        int bnt=getsum(q[i]);
        if(bnt+cs[q[i]]>res)
        {
            ans=q[i];
            res=bnt+cs[q[i]];
        }
        else if(bnt+cs[q[i]]==res)
        {
            ans=min(ans,q[i]);
        }
    }
    return ans;
}

inline void go()
{
    int x,y,ans=0;
    while(m--)
    {
        scanf("%d%d",&x,&y);
        x=(x+ans-1)%n+1; y=(y+ans-1)%n+1;
        if(x>y) swap(x,y);
        ans=a[query(x,y)];
        printf("%d\n",ans); 
    }
}

int main()
{
    read();
    go();
    return 0;
}

再一次被二分干掉了。。。