AcWing 2521：数颜色 ← 带修改的莫队算法

【题目来源】
https://www.acwing.com/problem/content/2523/

【题目描述】
墨墨购买了一套 N 支彩色画笔（其中有些颜色可能相同），摆成一排，你需要回答墨墨的提问。
墨墨会像你发布如下指令：
    Q L R：代表询问你从第 L 支画笔到第 R 支画笔中共有几种不同颜色的画笔。
    R P Col：把第 P 支画笔替换为颜色 Col。
为了满足墨墨的要求，你知道你需要干什么了吗？

【输入格式】
第 1 行两个整数 N,M，分别代表初始画笔的数量以及墨墨会做的事情的个数。
第 2 行 N 个整数，分别代表初始画笔排中第 i 支画笔的颜色。
第 3 行到第 2+M 行，每行分别代表墨墨会做的一件事情，格式见题干部分。

【输出格式】
对于每一个 Query 的询问，你需要在对应的行中给出一个数字，代表第 L 支画笔到第 R 支画笔中共有几种不同颜色的画笔。

【数据范围】
1≤N,M≤10000,
修改操作不多于 1000 次，
所有的输入数据中出现的所有整数均大于等于 1 且不超过 10^6。

【输入样例】
6 5
1 2 3 4 5 5
Q 1 4
Q 2 6
R 1 2
Q 1 4
Q 2 6

【输出样例】
4
4
3
4

【算法分析】
● 基础莫队算法详见：https://blog.csdn.net/hnjzsyjyj/article/details/138976338  

● 基础莫队算法是“离线+暴力+分块”算法，通常用于“不修改、只查询”的区间问题。
“在线”是交互式的，一问一答。特别的，如果前面的答案用于后面的提问，称为“强制在线”。
“离线”是非交互的，一次性读取所有问题，一起回答。
如果只是涉及简单的“单点修改”的区间问题，也可以使用莫队算法。这便是带修改的莫队算法。

● 在带修改的莫队算法中，一个查询通常表示为 (L,R,t)，含义为在查询区间 [L,R] 前进行了 t 次修改操作。一般地，可以把 t 理解为“时间”。t 的范围为 1≤t≤m，m为操作（查询、修改）次数。两个查询的 t 相差越小，它们对应的数据差别越小，计算量也越小。所以，对 t 排序能减少计算量。

● 在带修改的莫队算法中，块的大小为 ，而不是基础莫队算法的 sqrt(n) 个块。其中，n 为元素个数。

莫队分块思想详见：https://blog.csdn.net/hnjzsyjyj/article/details/138955263

● 莫队算法的精髓在于奇偶性排序。因此，带修改的莫队算法依然会进行奇偶性排序。也就是依据左端点 L 位于奇数块还是偶数块，决定右端点 R 从小到大排序还是从大到小排序。
（1）首先，利用分块算法将给定的 n 个数分成  块，而不是基础莫队算法的 sqrt(n) 个块；
（2）然后，将多个询问的左端点 L 按块从小到大排序。
若 L 位于奇数块，则对右端点 R 从小到大排序。若 L 位于偶数块，则对右端点 R 从大到小排序。反之亦可。
示意图如下所示：

​​​​​
【算法代码】

#include<bits/stdc++.h>
using namespace std;

const int maxn=1e4+5;
const int maxs=1e6+5;

int a[maxn],id[maxn],ans[maxn];
int cnt[maxs];
int qt,mt; //query times, modify times
int num;
int n,m;

struct Query {
    int sn,x,y,t; //The order is the same as line 59
} q[maxn];

struct Change {
    int p,v; //position,val
} md[maxn];

bool cmp(Query a,Query b) {
    if(id[a.x]!=id[b.x]) return id[a.x]<id[b.x];
    if(id[a.y]!=id[b.y]) return id[a.y]<id[b.y];
    return a.t<b.t;
}

void add(int x) {
    cnt[x]++;
    if(cnt[x]==1) num++;
}

void del(int x) {
    cnt[x]--;
    if(!cnt[x]) num--;
}

void mod(int t,int i) {
    if(q[i].x<=md[t].p && md[t].p<=q[i].y) {
        add(md[t].v);
        del(a[md[t].p]);
    }
    swap(a[md[t].p],md[t].v);
}

int main() {
    scanf("%d%d",&n,&m);
    int block=pow(n,2.0/3.0);
    for(int i=1; i<=n; i++) {
        scanf("%d",&a[i]);
        id[i]=(i-1)/block+1;
    }

    char ch;
    int x,y;
    for(int i=1; i<=m; i++) {
        getchar();
        scanf("%c",&ch);
        scanf("%d%d",&x,&y);
        if(ch=='Q') q[++qt]= {qt,x,y,mt};
        else md[++mt]= {x,y};
    }

    sort(q+1,q+qt+1,cmp);

    int L=1,R=0,t=0;
    for(int i=1; i<=qt; i++) {
        while(L<q[i].x) del(a[L++]);
        while(L>q[i].x) add(a[--L]);
        while(R<q[i].y) add(a[++R]);
        while(R>q[i].y) del(a[R--]);
        while(t<q[i].t) mod(++t,i);
        while(t>q[i].t) mod(t--,i);
        ans[q[i].sn]=num;
    }

    for(int i=1; i<=qt; i++) printf("%d\n",ans[i]);

    return 0;
}


/*
in:
6 5
1 2 3 4 5 5
Q 1 4
Q 2 6
R 1 2
Q 1 4
Q 2 6

out:
4
4
3
4
*/

【参考文献】
https://blog.csdn.net/weixin_75161465/article/details/137195888
https://blog.csdn.net/GROZAX/article/details/130069889
http://fangkaipeng.com/?p=1504
https://www.acwing.com/solution/content/157290/