BZOJ2724 [Violet]蒲公英(分块)

区间众数。分块,预处理任意两块间所有数的众数,和每块中所有数的出现次数的前缀和。查询时对不是整块的部分暴力,显然只有这里出现的数可能更新答案。于是可以优美地做到O(n√n)。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
int read()
{
int x=,f=;char c=getchar();
while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
return x*f;
}
#define N 50010
#define BLOCK 250
int n,m,a[N],b[N],lastans=;
int block,tot,L[N],R[N],pos[N];
int cnt[N],f[BLOCK][BLOCK],sum[BLOCK][N];
int main()
{
freopen("bzoj2724.in","r",stdin);
freopen("bzoj2724.out","w",stdout);
n=read(),m=read();
for (int i=;i<=n;i++) b[i]=a[i]=read();
sort(b+,b+n+);
int t=unique(b+,b+n+)-b;
for (int i=;i<=n;i++) a[i]=lower_bound(b+,b+t,a[i])-b;
block=sqrt(n);tot=n/block+(n%block>);
for (int i=;i<=n/block;i++)
L[i]=(i-)*block+,R[i]=(i-)*block+block;
if (n/block<tot) L[tot]=n/block*block+,R[tot]=n;
for (int i=;i<=tot;i++)
{
memset(cnt,,sizeof(cnt));
int num=;
for (int j=i;j<=tot;j++)
{
for (int k=L[j];k<=R[j];k++)
{
cnt[a[k]]++;
if (cnt[a[k]]>cnt[num]||cnt[a[k]]==cnt[num]&&a[k]<num) num=a[k];
}
f[i][j]=num;
}
memcpy(sum[i],sum[i-],sizeof(sum[i]));
for (int j=L[i];j<=R[i];j++)
pos[j]=i,sum[i][a[j]]++;
}
memset(cnt,,sizeof(cnt));
while (m--)
{
int x=read(),y=read();
x=(x+lastans-)%n+,y=(y+lastans-)%n+;
if (x>y) swap(x,y);
int num=;
if (pos[x]==pos[y])
{
for (int i=x;i<=y;i++)
{
cnt[a[i]]++;
if (cnt[a[i]]>cnt[num]||cnt[a[i]]==cnt[num]&&a[i]<num) num=a[i];
}
for (int i=x;i<=y;i++) cnt[a[i]]--;
}
else
{
num=f[pos[x]+][pos[y]-];
for (int i=x;i<=R[pos[x]];i++)
{
cnt[a[i]]++;
if (cnt[a[i]]+sum[pos[y]-][a[i]]-sum[pos[x]][a[i]]>cnt[num]+sum[pos[y]-][num]-sum[pos[x]][num]
||cnt[a[i]]+sum[pos[y]-][a[i]]-sum[pos[x]][a[i]]==cnt[num]+sum[pos[y]-][num]-sum[pos[x]][num]&&a[i]<num)
num=a[i];
}
for (int i=L[pos[y]];i<=y;i++)
{
cnt[a[i]]++;
if (cnt[a[i]]+sum[pos[y]-][a[i]]-sum[pos[x]][a[i]]>cnt[num]+sum[pos[y]-][num]-sum[pos[x]][num]
||cnt[a[i]]+sum[pos[y]-][a[i]]-sum[pos[x]][a[i]]==cnt[num]+sum[pos[y]-][num]-sum[pos[x]][num]&&a[i]<num)
num=a[i];
}
for (int i=x;i<=R[pos[x]];i++) cnt[a[i]]--;
for (int i=L[pos[y]];i<=y;i++) cnt[a[i]]--;
}
lastans=b[num];
printf("%d\n",b[num]);
}
fclose(stdin);fclose(stdout);
return ;
}
上一篇:(转)IC验证概述


下一篇:NumberSort选择法、冒泡法排序