2724: [Violet 6]蒲公英

题目链接

题目大意：区间众数，强制在线

题解：维护f[x][y]表示第i块到第j块的答案
只要能快速得出一个数在某个区间内出现次数即可
这里有两种方法，第一种是用vector二分，见hzwer

还有一种方法是维护sum[i][j]表示数字i在第1块到第j块出现了多少次
然后一阵操作就可以做了……

虽然理论上少一个log，但是我写的两种速度差不多……

块大小约为$\sqrt{n/\log n}$

我的收获：分块套路·……

#include <cstdio>
#include <vector>
#include <map>
#include <cmath>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
using namespace std;

const int N=50005;

int n,m,blo,ans,tot;
int a[N],z[N],pos[N],cnt[N],val[N];
int f[505][505],sum[N][505];

void pre(int x)
{   
    int ret=0;
    for(int i=(x-1)*blo+1;i<=n;i++){       
        cnt[a[i]]++;        
        if(cnt[a[i]]>cnt[ret]||(cnt[a[i]]==cnt[ret]&&a[i]<ret)) ret=a[i];
        f[x][pos[i]]=ret;
    }
    memset(cnt,0,sizeof(cnt));
}

int query(int x,int y)
{
    if(pos[x]==pos[y]){
        int now=0;
        for(int i=x;i<=y;i++){
            cnt[a[i]]++;
            if(cnt[a[i]]>cnt[now]||cnt[a[i]]==cnt[now]&&a[i]<now) now=a[i];
        }
        for(int i=x;i<=y;i++) cnt[a[i]]--;
        return now;
    }
    int now=f[pos[x]+1][pos[y]-1];
    cnt[now]+=sum[now][pos[y]-1]-sum[now][pos[x]];
    for(int i=x;i<=pos[x]*blo;i++)
    {
        if(!cnt[a[i]]) cnt[a[i]]+=sum[a[i]][pos[y]-1]-sum[a[i]][pos[x]];
        cnt[a[i]]++;
        if(cnt[a[i]]>cnt[now]||cnt[a[i]]==cnt[now]&&a[i]<now) now=a[i];
    }
    for(int i=(pos[y]-1)*blo+1;i<=y;i++)
    {
        if(!cnt[a[i]]) cnt[a[i]]+=sum[a[i]][pos[y]-1]-sum[a[i]][pos[x]];
        cnt[a[i]]++;
        if(cnt[a[i]]>cnt[now]||cnt[a[i]]==cnt[now]&&a[i]<now) now=a[i];
    }
    cnt[f[pos[x]+1][pos[y]-1]]=0;
    for(int i=x;i<=pos[x]*blo;i++) cnt[a[i]]=0;
    for(int i=(pos[x]-1)*blo+1;i<=y;i++) cnt[a[i]]=0;
    return now;
}

void work()
{
    for(int i=1;i<=m;i++)
    {
        int l,r;scanf("%d%d",&l,&r);
        l=(l+ans-1)%n+1;r=(r+ans-1)%n+1;
        if(l>r) swap(l,r);
        ans=z[query(l,r)];
        printf("%d\n",ans);
    }
}

void init()
{
    scanf("%d%d",&n,&m);blo=sqrt(n);
    for(int i=1;i<=n;i++) pos[i]=(i-1)/blo+1;
    for(int i=1;i<=n;i++) scanf("%d",&a[i]),z[i]=a[i];
    sort(z+1,z+1+n);int cnt=unique(z+1,z+1+n)-z-1;
    for(int i=1;i<=n;i++) a[i]=lower_bound(z+1,z+1+cnt,a[i])-z;
    for(int i=1;i<=pos[n];i++) pre(i);
    for(int i=1;i<=n;i++) sum[a[i]][pos[i]]++;
    for(int i=1;i<=cnt;i++)
        for(int j=1;j<=pos[n];j++)
            sum[i][j]+=sum[i][j-1];
}

int main()
{
    init();
    work();
    return 0;
}