BZOJ4241 历史研究（莫队）

　　如果分块的话与区间众数没有本质区别。这里考虑莫队。

　　显然莫队时的删除可以用堆维护，但多了一个log不太跑得过。

　　有一种叫回滚莫队的trick，可以将问题变为只有加入操作。按莫队时分的块依次处理，一块中左端点的差不超过√n，右端点单调递增。首先将右端点也在该块中的询问暴力处理。然后令指针l在下一块开头，指针r在这一块结尾。暴力扩展右端点移动指针r，到达询问点时，移动指针l以回答询问，但不让指针l的移动对之后的询问产生影响，即回滚。这样就可以处理删除了。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
int read()
{
int x=0,f=1;char c=getchar();
while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}
while (c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
return x*f;
}
#define N 100010
#define ll long long
int n,m,a
,b
,cnt
;
ll ans
;
struct data
{
int x,y,k,i;
bool operator <(const data&a) const
{
return k<a.k||k==a.k&&y<a.y;
}
}q
;
int main()
{
#ifndef ONLINE_JUDGE
freopen("bzoj4241.in","r",stdin);
freopen("bzoj4241.out","w",stdout);
const char LL[]="%I64d\n";
#else
const char LL[]="%lld\n";
#endif
n=read(),m=read();
for (int i=1;i<=n;i++) b[i]=a[i]=read();
sort(b+1,b+n+1);
int t=unique(b+1,b+n+1)-b-1;
for (int i=1;i<=n;i++) a[i]=lower_bound(b+1,b+t+1,a[i])-b;
int block=sqrt(n);
for (int i=1;i<=m;i++) q[i].x=read(),q[i].y=read(),q[i].k=q[i].x/block,q[i].i=i;
sort(q+1,q+m+1);
for (int i=1;i<=m;i++)
{
int t=i;while (t<m&&q[t+1].k==q[i].k) t++;
while (i<=t&&q[i].y<(q[i].k+1)*block)
{
for (int j=q[i].x;j<=q[i].y;j++)
{
cnt[a[j]]++;
ans[q[i].i]=max(ans[q[i].i],1ll*cnt[a[j]]*b[a[j]]);
}
for (int j=q[i].x;j<=q[i].y;j++) cnt[a[j]]--;
i++;
}
int r=(q[i].k+1)*block-1;ll v=0;
for (int j=i;j<=t;j++)
{
while (r<q[j].y)
{
cnt[a[++r]]++;
v=max(v,1ll*cnt[a[r]]*b[a[r]]);
}
ans[q[j].i]=v;
for (int k=(q[i].k+1)*block-1;k>=q[j].x;k--)
{
cnt[a[k]]++;
ans[q[j].i]=max(ans[q[j].i],1ll*cnt[a[k]]*b[a[k]]);
}
for (int k=(q[i].k+1)*block-1;k>=q[j].x;k--) cnt[a[k]]--;
}
memset(cnt,0,sizeof(cnt));
i=t;
}
for (int i=1;i<=m;i++) printf(LL,ans[i]);
return 0;
}