【BZOJ4241】【回滚莫队】历史研究 （非题解的学习分析）

回滚莫队：用于处理难以删除但是易于添加（其实易于删除难以添加也可以，但是没见过这样题-_-）的莫队，排序照常，如果左右端点在同一块直接暴力，这部分最多n sqrt n，否则把左端点在一块的一起处理，清空莫队，然后直接令莫队左端点在块尾，这部分n sqrtn，右端点照常走，这部分n sqrtn ，左端点每次走的时候记录更改了哪些量，走到地方记录完答案把修改回滚回去，这部分也是n sqrtn，所以复杂度还是根号的，但是把删除干掉了

neither_nor dalao 上课讲到回滚莫队，好像网上只有他写过。Orz%%%

例题：BZOJ4241
neither_nor 大佬的Blog
题目：
Description
IOI国历史研究的第一人——JOI教授，最近获得了一份被认为是古代IOI国的住民写下的日记。JOI教授为了通过这份日记来研究古代IOI国的生活，开始着手调查日记中记载的事件。
日记中记录了连续N天发生的时间，大约每天发生一件。
事件有种类之分。第i天(1<=i<=N)发生的事件的种类用一个整数Xi表示，Xi越大，事件的规模就越大。
JOI教授决定用如下的方法分析这些日记：
1. 选择日记中连续的一些天作为分析的时间段
2. 事件种类t的重要度为t*(这段时间内重要度为t的事件数)
3. 计算出所有事件种类的重要度，输出其中的最大值
现在你被要求制作一个帮助教授分析的程序，每次给出分析的区间，你需要输出重要度的最大值。
Input
第一行两个空格分隔的整数N和Q，表示日记一共记录了N天，询问有Q次。
接下来一行N个空格分隔的整数X1…XN，Xi表示第i天发生的事件的种类
接下来Q行，第i行(1<=i<=Q)有两个空格分隔整数Ai和Bi，表示第i次询问的区间为[Ai,Bi]。
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<=Xi<=10^9 (1<=i<=N)
Source
JOI 2013~2014 春季training合宿 竞技1 By PoPoQQQ

思路见上
分块和莫队套些奇怪的东西也可以过这题

//By neither_nor
#include<iostream>  
#include<cstdio>  
#include<cstdlib>  
#include<algorithm>  
#include<iomanip>  
#include<cstring>  
#include<cmath>  
#include<ctime>  
#include<vector>  
#include<stack>  
#include<queue>  
#include<set>  
#include<bitset>  
#include<map>  
using namespace std;  
#define MAXN 100010  
#define MAXM 100010  
#define INF 1000000000  
#define MOD 1000000007  
#define ll long long  
#define eps 1e-8  
char xB[1<<15],*xS=xB,*xT=xB;  
#define getc() (xS==xT&&(xT=(xS=xB)+fread(xB,1,1<<15,stdin),xS==xT)?0:*xS++)

inline int read() {
    int x=0,f=1;char xch=getc();
    while(xch<'0'|xch>'9'){if(xch=='-')f=-1;xch=getc();}
    while(xch>='0'&&xch<='9'){x=x*10+xch-'0';xch=getc();}
    return x*f;
}
struct que {  
    int l;
    int r;
    int num;
    int k;
    friend bool operator < (que x, que y){
        return x.k != y.k ? x.k < y.k : x.r < y.r;
    }
};
int siz;
que q[MAXM];
int a[MAXN];
int cnt[MAXN];
ll ans[MAXN];
int num[MAXN];
int visc[MAXN];
ll ANS;
map<int,int>h;
int tls[MAXN],tln,mx;
int g[MAXN];
int n,m;
int L,R;
int T;
int *st[MAXN];
int stv[MAXN],tp;
ll *sta[MAXN];
ll stva[MAXN];
int tpa;

void roolback() {
    while(tp) (*st[tp]) = stv[tp--];
    while(tpa) (*sta[tpa])=stva[tpa--];
}

int main() {  
    int i,j;  
    n = read();
    m = read();
    for(i = 1; i <= n; i++) a[i] = read(), tls[++tln] = a[i];
    sort(tls+1, tls+tln+1);
    tls[0] = -INF;
    for(i = 1; i <= tln; i++) if(tls[i] != tls[i-1]) h[tls[i]] = ++mx, g[mx] = tls[i];
    for(i = 1; i <= n; i++) a[i] = h[a[i]];
    siz = sqrt(n);
    for(i = 1; i <= m; i++) {
        q[i].l = read();
        q[i].r = read();
        q[i].k = (q[i].l-1)/siz+1;
        q[i].num = i;
    }
    sort(q+1, q+m+1);
    for(i = 1; i <= m; i++) {
        if(q[i].k != q[i-1].k) {
            T++;
            L = siz*q[i].k+1;
            R = L-1;
            ANS = 0;
        }
        if(q[i].r < L) {
            ANS = 0;
            for(j = q[i].l; j <= q[i].r; j++) {
                if(visc[a[j]] != T) cnt[a[j]] = 0;
                visc[a[j]] = T;
                cnt[a[j]]++;
                if((ll)g[a[j]]*cnt[a[j]] > ANS) ANS = (ll)g[a[j]]*cnt[a[j]];
            }
            ans[q[i].num] = ANS;
            T++;
            ANS = 0;
            continue ;
        }
        while(R < q[i].r) {
            R++;
            if(visc[a[R]] != T) cnt[a[R]] = 0;
            visc[a[R]] = T;
            cnt[a[R]]++;
            if((ll)g[a[R]]*cnt[a[R]] > ANS) ANS = (ll)g[a[R]]*cnt[a[R]];
        }
        while(L > q[i].l) {
            L--;
            if(visc[a[L]] != T) cnt[a[L]] = 0;
            st[++tp] = &visc[a[L]];
            stv[tp] = visc[a[L]];
            visc[a[L]] = T;
            st[++tp] = &cnt[a[L]];
            stv[tp] = cnt[a[L]];
            cnt[a[L]]++;
            if((ll)g[a[L]]*cnt[a[L]] > ANS) {
                sta[++tpa] = &ANS;
                stva[tpa] = ANS;
                ANS = (ll)g[a[L]]*cnt[a[L]];
            }
        }
        ans[q[i].num] = ANS;
        roolback();
        L = siz*q[i].k+1;
    }
    for(i = 1; i <= m; i++) printf("%lld\n",ans[i]);
    return 0;
}