BZOJ2724: [Violet 6]蒲公英（分块）

传送门

• 题意

给一个数列，求区间众数（若有多个则输出较小值）。

• 题解

考虑分块，对于块i维护cnt[i][val]表示从第一块到i块中val出现了几次，维护ans[i][j]表示从第i块到第j块的众数（之后会用）。

因为没有修改操作，所以预处理cnt,ans;

对于ans，固定左端点，右端点暴力枚举。时间复杂度$O(n\sqrt{n})$。

对于cnt，暴力枚举，每扫描完一个区间就枚举val更新区间。时间复杂度$O(n\sqrt{n})$

对于查询操作，容易发现众数只会是所有完整块一起的众数或者不完整块中的数。（不完整块中的数可能使完整块中比众数次数小的数超过众数）前半部分已经预处理，个数计算可以通过前缀和实现，后半部分暴力枚举，加上完整块中的个数可以更新答案。时间复杂度$O(2\sqrt{n}·m)=O(m\sqrt{n})$

（md,特别坑的是bzoj居然不能cout,然后一直RE…）

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<string>
#include<cmath>
using namespace std;
const int Maxn=4e4+50;
const int S=200;

inline int read()
{
    char ch=getchar();int i=0,f=1;
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch)){i=(i<<1)+(i<<3)+ch-'0';ch=getchar();}
    return i*f;
}

int n,m,a[Maxn],c[Maxn],ans[S+50][S+50],cnt[S+50][Maxn],bg[S+50],ed[S+50],rk[Maxn],s;

struct node
{
    int val,id;
    friend inline bool operator <(const node &a,const node &b)
    {
        return a.val<b.val;
    }
}b[Maxn];

inline void l_b()
{
    sort(b+1,b+n+1);
    int t=0;
    for(int i=1;i<=n;i++)
    {
        if(i==1||b[i].val!=b[i-1].val)t++;
        a[b[i].id]=t;
        rk[t]=b[i].val;
    }
    return;
}

inline void pre()
{
    for(int i=1;i*S<=n;i++)
    {
        bg[i]=ed[i-1]+1;ed[i]=ed[i-1]+S;
        s=i;
    }
    if(s*S!=n)
    {
        ++s;
        bg[s]=ed[s-1]+1;
        ed[s]=n;
    }
    for(int i=1;i<=s;i++)
    {
        memset(c,0,sizeof(c));
        int maxx=0,A=0,h=i;
        for(int j=bg[i];j<=n;j++)
        {
            c[a[j]]++;
            if(c[a[j]]==maxx&&a[j]<A)A=a[j];
            else if(c[a[j]]>maxx)maxx=c[a[j]],A=a[j];
            if(j==ed[h])
            {
                ans[i][h]=A;
                h++;
            }
        }
    }
    memset(c,0,sizeof(c));
    int h=1;
    for(int j=1;j<=n;j++)
    {
        c[a[j]]++;
        if(j==ed[h])
        {
            for(int i=1;i<=n;i++)
            cnt[h][i]=c[i];
            h++;
        }
    }
    return;
}

int buf[50];
inline void W(int x)
{
    if(!x)putchar('0');
    if(x<0)putchar('-'),x=-x;
    while(x)buf[++buf[0]]=x%10,x/=10;
    while(buf[0])putchar('0'+buf[buf[0]--]);
    return;
}

int main()
{
    freopen("lx.in","r",stdin);
    n=read(),m=read();
    for(int i=1;i<=n;i++)b[i].val=read(),b[i].id=i;
    l_b();
    pre();
    int A=0;
    memset(c,0,sizeof(c));
    while(m--)
    {
        int l=((read()+rk[A]-1)%n)+1,r=((read()+rk[A]-1)%n)+1;
        if(l>r)swap(l,r);
        if(r-l+1>2*S)
        {
            int L=(l-1)/S+1,R=(r-1)/S+1;
            if(bg[L]!=l)L++;
            if(ed[R]!=r)R--;
            int maxx=0;
            A=0;
            A=ans[L][R],maxx=cnt[R][A]-cnt[L-1][A];
            for(int i=l;i<=bg[L]-1;i++)
            {
                c[a[i]]++;
                if((c[a[i]]+cnt[R][a[i]]-cnt[L-1][a[i]])==maxx)A=min(a[i],A);
                else if((c[a[i]]+cnt[R][a[i]]-cnt[L-1][a[i]])>maxx)A=a[i],maxx=(c[a[i]]+cnt[R][a[i]]-cnt[L-1][a[i]]);
            }
            for(int i=ed[R]+1;i<=r;i++)
            {
                c[a[i]]++;
                if((c[a[i]]+cnt[R][a[i]]-cnt[L-1][a[i]])==maxx)A=min(a[i],A);
                else if((c[a[i]]+cnt[R][a[i]]-cnt[L-1][a[i]])>maxx)A=a[i],maxx=(c[a[i]]+cnt[R][a[i]]-cnt[L-1][a[i]]);
            }
            for(int i=l;i<=bg[L]-1;i++)c[a[i]]--;
            for(int i=ed[R]+1;i<=r;i++)c[a[i]]--;
            W(rk[A]);
            putchar('\n');
        }
        else
        {
            int maxx=0;
            A=0;
            for(int i=l;i<=r;i++)
            {
                c[a[i]]++;
                if(c[a[i]]==maxx&&a[i]<A)A=a[i];
                else if(c[a[i]]>maxx)A=a[i],maxx=c[a[i]];
            }
            for(int i=l;i<=r;i++)
            {
                c[a[i]]--;
            }
            W(rk[A]);
            putchar('\n');
        }
    }
    return 0;
}