洛谷P4168 [Violet]蒲公英【分块】

时空限制 2000ms / 512MB

题目描述

在乡下的小路旁种着许多蒲公英，而我们的问题正是与这些蒲公英有关。

为了简化起见，我们把所有的蒲公英看成一个长度为n的序列 (a_1,a_2…a_n)，其中ai为一个正整数，表示第i棵蒲公英的种类编号。

而每次询问一个区间 [l,r]，你需要回答区间里出现次数最多的是哪种蒲公英，如果有若干种蒲公英出现次数相同，则输出种类编号最小的那个。

注意，你的算法必须是在线的

输入格式：

第一行两个整数 n,m ，表示有n株蒲公英，m 次询问。

接下来一行n个空格分隔的整数 ai，表示蒲公英的种类

再接下来m 行每行两个整数l0,r0 ，我们令上次询问的结果为 x（如果这是第一次询问， 则 x=0）。

令$l=(l_0+x-1)\mathrm{m}\mathrm{o}\mathrm{d}n+1,r=(r_0+x-1)\mathrm{m}\mathrm{o}\mathrm{d}n+1$，如果 l>r，则交换 l,r 。

最终的询问区间为[l,r]。

输出格式：

输出m 行。每行一个整数，表示每次询问的结果。

说明

对于 100% 的数据，保证 $1\le n\le 40000,1\le m\le 50000,1\le a_i\le 10^9$

---

题目分析

套路分块，黑题评分十分鬼畜
$sum[i][x]$记录从第i块到序列末尾每个数出现次数，mx[i][j]记录块i到块j的区间众数
一个[l,r]的询问，假如其中整块的区间为[L,R]
那么众数要么出现在[L,R]中，要么出现在[l,L)或(R,r]中
遍历这两部分不是整块的部分，记录每个数出现次数，结合sum数组更新即可
注意数值离散化，出现次数相同更新编号小的

#include<iostream>
#include<cstdio>
#include<cmath>
#include<queue>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long lt;

int read()
{
    int f=1,x=0;
    char ss=getchar();
    while(ss<'0'||ss>'9'){if(ss=='-')f=-1;ss=getchar();}
    while(ss>='0'&&ss<='9'){x=x*10+ss-'0';ss=getchar();}
    return f*x;
}

const int maxn=50010;
int n,m,t;
int a[maxn],b[maxn],pos[maxn],cnt;
int lft[310],rht[310],bl[maxn];
int sum[310][maxn],mx[310][310];
int rem[maxn],ans;

void build()
{
    for(int i=1;i<=t;++i)
    lft[i]=(i-1)*t+1,rht[i]=i*t;
    if(rht[t]<n) t++,lft[t]=rht[t-1]+1,rht[t]=n;
    
    for(int i=1;i<=t;++i)
    for(int j=lft[i];j<=rht[i];++j)
    bl[j]=i;
    
    for(int i=1;i<=t;++i)
    {
        int tt=0;
        for(int j=lft[i];j<=n;++j)
        {
            sum[i][a[j]]++;
            if(sum[i][a[j]]>sum[i][tt]) tt=a[j];
            else if(sum[i][a[j]]==sum[i][tt]) tt=pos[a[j]]<pos[tt]?a[j]:tt;
            if(bl[j]!=bl[j+1]) mx[i][bl[j]]=tt;
        }
    }
}

void cls(int ll,int rr)
{
    for(int i=ll;i<=rr;++i) rem[a[i]]=0;
}

int solve(int ll,int rr)
{
    int p=bl[ll],q=bl[rr],res=0;
    if(q-p<=1)
    {
        for(int i=ll;i<=rr;++i) rem[a[i]]++;
        for(int i=ll;i<=rr;++i)
        {
            if(rem[a[i]]>rem[res]) res=a[i];
            else if(rem[a[i]]==rem[res]) res=pos[a[i]]<pos[res]?a[i]:res;
        }
        cls(ll,rr);
    }
    else
    {
        res=mx[p+1][q-1];
        for(int i=ll;i<=rht[p];++i) rem[a[i]]++;
        for(int i=lft[q];i<=rr;++i) rem[a[i]]++;
        
        for(int i=ll;i<=rht[p];++i)
        {
            int tt1=rem[a[i]]+sum[p+1][a[i]]-sum[q][a[i]];
            int tt2=rem[res]+sum[p+1][res]-sum[q][res];
            if(tt1>tt2) res=a[i];
            else if(tt1==tt2) res=pos[a[i]]<pos[res]?a[i]:res;
        }
        for(int i=lft[q];i<=rr;++i)
        {
            int tt1=rem[a[i]]+sum[p+1][a[i]]-sum[q][a[i]];
            int tt2=rem[res]+sum[p+1][res]-sum[q][res];
            if(tt1>tt2) res=a[i];
            else if(tt1==tt2) res=pos[a[i]]<pos[res]?a[i]:res;
        }
        cls(ll,rht[p]); cls(lft[q],rr);
    }
    return res;
}

int main()
{
    n=read();m=read();
    for(int i=1;i<=n;++i)
    a[i]=b[i]=read();
    
    sort(b+1,b+1+n);
    for(int i=1;i<=n;++i)
    if(i==1||b[i]!=b[i-1])
    pos[++cnt]=b[i];
    
    for(int i=1;i<=n;++i)
    a[i]=lower_bound(pos+1,pos+1+cnt,a[i])-pos;
    
    t=sqrt(n); build();
    while(m--)
    {
    	int ll=(read()+ans-1)%n+1,rr=(read()+ans-1)%n+1;
    	if(ll>rr) swap(ll,rr);
    	ans=pos[solve(ll,rr)];
    	printf("%d\n",ans);
    }
    return 0;
}