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

题目大意:

n个数，m次询问，每次询问区间众数，强制在线，1<=n<=40000,1<=m<=50000。

解题思路：

分块裸题。存一下模板。
用$O(n\sqrt n)$的时间预处理出sum[i][j]表示i在前j块的出现次数，ans[i][j]表示第i块到第j块的答案(枚举左端点)，询问时只用考虑边界上至多$2\sqrt n$个数的影响即可。

#include<bits/stdc++.h>
#define ll long long
using namespace std;

int getint()
{
    int i=0,f=1;char c;
    for(c=getchar();(c!='-')&&(c<'0'||c>'9');c=getchar());
    if(c=='-')f=-1,c=getchar();
    for(;c>='0'&&c<='9';c=getchar())i=(i<<3)+(i<<1)+c-'0';
    return i*f;
}

const int N=40005;
int n,m,c,x,mx,S,num,a[N],b[N],sum[N][205],ans[205][205],cnt[N];
bool vis[N];

int main()
{
    //freopen("lx.in","r",stdin);
    //freopen("lx.out","w",stdout);
    n=getint(),m=getint(),S=sqrt(n);
    for(int i=1;i<=n;i++)a[i]=b[i]=getint();
    sort(b+1,b+n+1),c=unique(b+1,b+n+1)-b-1;
    for(int i=1;i<=n;i++)a[i]=lower_bound(b+1,b+c+1,a[i])-b;
    for(int i=1;i<=n;i++)
    {
        cnt[a[i]]++;
        if(i%S==0||i==n)
        {
            ++num;
            for(int j=1;j<=c;j++)sum[j][num]=cnt[j];
        }
    }
    for(int i=1;i<=num;i++)
    {
        for(int j=1;j<=c;j++)cnt[j]=0;
        int id=0;
        for(int j=(i-1)*S+1;j<=n;j++)
        {
            cnt[a[j]]++;
            if(cnt[a[j]]>cnt[id]||cnt[a[j]]==cnt[id]&&a[j]<id)id=a[j];
            if(j%S==0||j==n)ans[i][(j-1)/S+1]=id;
        }
    }
    while(m--)
    {
        int l=(getint()+x-1)%n+1,r=(getint()+x-1)%n+1;
        if(l>r)swap(l,r);
        x=0;
        int lS=(l-1)/S+2,rS=(r-1)/S;
        if(lS>rS)
        {
            for(int i=l;i<=r;i++)
            {
                vis[a[i]]?++cnt[a[i]]:vis[a[i]]=cnt[a[i]]=1;
                if(cnt[a[i]]>cnt[x]||cnt[a[i]]==cnt[x]&&a[i]<x)x=a[i];
            }
            for(int i=l;i<=r;i++)vis[a[i]]=0;
        }
        else
        {
            x=ans[lS][rS],mx=sum[x][rS]-sum[x][lS-1];
            int led=(lS-1)*S,rst=rS*S+1;
            for(int i=l;i<=led;i++)
            {
                vis[a[i]]?++cnt[a[i]]:vis[a[i]]=cnt[a[i]]=1;
                int tmp=cnt[a[i]]+sum[a[i]][rS]-sum[a[i]][lS-1];
                if(tmp>mx||tmp==mx&&a[i]<x)x=a[i],mx=tmp;
            }
            for(int i=rst;i<=r;i++)
            {
                vis[a[i]]?++cnt[a[i]]:vis[a[i]]=cnt[a[i]]=1;
                int tmp=cnt[a[i]]+sum[a[i]][rS]-sum[a[i]][lS-1];
                if(tmp>mx||tmp==mx&&a[i]<x)x=a[i],mx=tmp;
            }
            for(int i=l;i<=led;i++)vis[a[i]]=0;
            for(int i=rst;i<=r;i++)vis[a[i]]=0;
        }
        printf("%d\n",x=b[x]);
    }
    return 0;
}