bzoj 4241: 历史研究 回滚莫队

Description
IOI国历史研究的第一人——JOI教授，最近获得了一份被认为是古代IOI国的住民写下的日记。JOI教授为了通过这份日记来研究古代IOI国的生活，开始着手调查日记中记载的事件。
日记中记录了连续N天发生的时间，大约每天发生一件。
事件有种类之分。第$i$天$(1<=i<=N)$发生的事件的种类用一个整数$X_i$表示，$X_i$越大，事件的规模就越大。
JOI教授决定用如下的方法分析这些日记：
1. 选择日记中连续的一些天作为分析的时间段
2. 事件种类t的重要度为t*(这段时间内重要度为t的事件数)
3. 计算出所有事件种类的重要度，输出其中的最大值
现在你被要求制作一个帮助教授分析的程序，每次给出分析的区间，你需要输出重要度的最大值。

Input
第一行两个空格分隔的整数$N$和$Q$，表示日记一共记录了$N$天，询问有$Q$次。
接下来一行N个空格分隔的整数$X_1...X_N$，$X_i$表示第i天发生的事件的种类。
接下来$Q$行，第$i$行$(1<=i<=Q)$有两个空格分隔整数$A_i$和$B_i$，表示第i次询问的区间为$[A_i,B_i]$。

Output
输出$Q$行，第$i$行$(1<=i<=Q)$一个整数，表示第i次询问的最大重要度

Sample Input
5 5
9 8 7 8 9
1 2
3 4
4 4
1 4
2 4

Sample Output
9
8
8
16
16

HINT
$1<=N<=10^5$
$1<=Q<=10^5$
$1<=X_i<=10^9 (1<=i<=N)$

分析：
听dalao讲课时，知道有这个东西，今天终于可以填坑了。
莫队很容易维护一些东西，但是也有缺陷，就是加入和删除的速度取决于题目所求。
比如说容易插入删除较难的有$max,min$等；删除容易插入较难的有$mex$等。
所以我们就可以使用回滚莫队来搞。
对于每一个询问，假如左右端点都在同一区间，直接暴力计算；
如果不是，我们可以先把左端点固定在当前块的最右端。先移动右端点，到达询问的右端点后，记录一个$cur$为当前答案，然后再移动左端点到询问左端点。搞出询问答案后，左端点移回当前块最右端，但是并不用修改，直接取回备份的$cur$值。

代码：

/**************************************************************
    Problem: 4241
    User: liangzihao
    Language: C++
    Result: Accepted
    Time:12908 ms
    Memory:5604 kb
****************************************************************/

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <map>
#define LL long long

const int maxn=1e5+7;

using namespace std;

int n,m,num,t;
int l[maxn],r[maxn],belong[maxn],h[maxn],b[maxn],a[maxn];
LL ans,cur;
int block;

struct query{
    int l,r,num;
    LL ans;
}q[maxn];

bool cmp1(query a,query b)
{
    if (belong[a.l]==belong[b.l]) return a.r<b.r;
    return a.l<b.l;
}

bool cmp2(query a,query b)
{
    return a.num<b.num;
}

void updata(int x,int d)
{
    if (d==1)
    {
        int k=++h[a[x]];
        ans=max(ans,(LL)k*(LL)b[a[x]]);
    }
    else
    {
        h[a[x]]--;
    }
}

int main()
{
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++)
    {
        scanf("%d",&a[i]);
        b[i]=a[i];
    }
    sort(b+1,b+n+1);
    int size=unique(b+1,b+n+1)-b-1;
    for (int i=1;i<=n;i++) a[i]=lower_bound(b+1,b+size+1,a[i])-b;
    block=trunc(sqrt(n));
    num=n/block+(n%block!=0);
    for (int i=1;i<=n;i++) belong[i]=(i-1)/block+1;  
    for (int i=1;i<=num;i++)
    {
        l[i]=(i-1)*block+1;
        r[i]=i*block;
    }
    r[num]=n;
    for (int i=1;i<=m;i++)
    {
        scanf("%d%d",&q[i].l,&q[i].r);
        q[i].num=i;
    }   
    sort(q+1,q+m+1,cmp1);
    int j=1;
    for (int i=1,last;i<=m;i=last+1,j++)
    {       
        last=i-1;
        while (belong[q[last+1].l]==j) last++;
        cur=ans=0;
        t=r[j];     
        for (int k=i;k<=last;k++)
        {
            if (belong[q[k].r]==j)
            {
                for (int d=q[k].l;d<=q[k].r;d++) updata(d,1);
                q[k].ans=ans;
                ans=0;
                for (int d=q[k].l;d<=q[k].r;d++) updata(d,-1);
            }
            else
            {
                for (;t<q[k].r;t++) updata(t+1,1);
                cur=ans;
                for (int d=q[k].l;d<=r[j];d++) updata(d,1);
                q[k].ans=ans;
                ans=cur;
                for (int d=q[k].l;d<=r[j];d++) updata(d,-1);
            }
        }
        for (int d=r[j]+1;d<=t;d++) updata(d,-1);
    }
    sort(q+1,q+m+1,cmp2);
    for (int i=1;i<=m;i++) printf("%lld\n",q[i].ans);
}