浅谈莫队算法--洛谷P1972 HH的项链

2020.5.18
莫队算法之前有说过，区域赛那道题，但是当时我仍然不太理解莫队算法这个算法，只能寄希望于带板子上场到时候撕了这道题。然而，这并不太现实，要是想在区域赛干掉其它队伍去一趟nac光荣退役，靠板子显然不太现实。

区间众数出现次数请参见我上一篇farming mars的题解

最近终于能集中精力搞一下算法了（住嘴，你mysql还没开始学呢），cf现在能做到d了，待我打回蓝名。莫队这边，之前一直听人强调说我们能在O(1)时间内得知区间l和r到l-1和r+1的情况，一直在揣摩这句话，到底是个什么意思，没搞懂。后来有一天跑步的时候突然想到了，那不就是只用搞定子状态就行了？比如查众数，修整区间的时候就是如果增大区间，那么看看这个区间里相关的数字出现次数+1就行了，反之亦然。然后回去看了看farming mars那道题，哦哦，我明白了，然后照着上面那个，之前用的是树状数组和主席树，重写了一遍莫队，然而不幸被毒瘤出题人卡掉了嘤嘤嘤。

不过算是了解了吧，如果当时我在场上，现在的我也会用莫队去思考问题，这大概是我能带给队友们的一些帮助，奥利给。题解如下

#include <bits/stdc++.h>
using namespace std;
#define limit (1000000 + 5)//防止溢出
#define INF 0x3f3f3f3f
#define inf 0x3f3f3f3f3f
#define lowbit(i) i&(-i)//一步两步
#define EPS 1e-6
#define FASTIO  ios::sync_with_stdio(false);cin.tie(0);
#define ff(a) printf("%lld\n",a );
#define pi(a,b) pair<a,b>
#define rep(i, a, b) for(int i = a ; i <= b ; ++i)
#define per(i, a, b) for(int i = b ; i >= a ; --i)
#define mint(a,b,c) min(min(a,b), c)
#define MOD 998244353
#define FOPEN freopen("C:\\Users\\administrator01\\CLionProjects\\untitled24\\data.txt", "rt", stdin)
typedef long long ll;
typedef unsigned long long ull;
ll read(){
    ll sign = 1, x = 0;char s = getchar();
    while(s > '9' || s < '0' ){if(s == '-')sign = -1;s = getchar();}
    while(s >= '0' && s <= '9'){x = x * 10 + s - '0';s = getchar();}
    return x * sign;
}//快读
void write(ll x){
    if(x / 10) write(x / 10);
    putchar(x % 10 + '0');
}
int n, k;
struct node{
    int l, r, qid, blo;

    bool operator<(const node &rhs)const{
        if(blo != rhs.blo)return l < rhs.l;
        else return blo & 1 ? r < rhs.r : rhs.r < r;
    }
}query[limit];
int  cnt[limit],a[limit],block,answer[limit],ans;
void add(int x){
    if(!cnt[a[x]]){
        ++ans;
    }
    ++cnt[a[x]];
}
void del(int x){
    --cnt[a[x]];
    if(!cnt[a[x]])--ans;
}
int main(){
#ifdef LOCAL
    FOPEN;
    //freopen("C:\\Users\\administrator01\\CLionProjects\\untitled24\\out.txt", "w", stdout);
#endif
    n = read();

    rep(i ,1, n)a[i] = read();
    k = read();
    block = int(sqrt(n >= 3 ? n * (2.0 / 3) : n));//分块
    block = (k / block) ? k /block : block;
    rep(i ,1, k){
        query[i].l = read(), query[i].r = read(), query[i].qid = i;
        query[i].blo = query[i].l / block;//分块
    }
    sort(query + 1 , query + 1 + k);
    int l = 0 , r = 0;
    ans = 0;
    rep(i ,1, k){
        if(query[i].l == query[i].r){
            answer[query[i].qid] = 1;
        }
        while(l < query[i].l)del(l++);//缩进
        while(l > query[i].l)add(--l);
        while(r < query[i].r)add(++r);
        while(r > query[i].r)del(r--);
        answer[query[i].qid] = ans;
    }
    rep( i,1,k)write(answer[i]), puts("");
    return 0;
}