【莫队算法】总结&CF940F Machine Learning

题意：

给出一个长度为N序列（为什么D、E、F都是序列）
需要支持两种操作：
1、定义一个区间的值为：这段区间任意元素出现次数的集合的mex，给出l,r求原串中[l,r]这段区间的值
2、修改某个点的值
对mex的定义与SG函数中是相同的，表示一个自然数集中未出现的最小的整数。
例如：$1、3、2、1、2、2、2$这个序列的值为3：
3出现了1次，1出现了2次，2出现了4次，集合中的数为$0,1,2,4$（所有除了这三个数以外的$[0，10^9]$的数均未出现，所以为0）
$N，M（操作次数）\le 10^5$
时限：4s

---

分析：

很直观的带修改的莫队算法题（我居然没学过就猜中了）

因为只学过普通的莫队算法，而且还是16年准备省选的时候学的，早就忘的差不多了，于是又回去补了补，现在重新写一个总结，以免以后再忘。

莫队算法：

莫队算法被称为区间问题的神器，但这神器的限制其实也颇多：必须离线操作，不支持区间修改，不支持无法在$O(1)$或近似的复杂度内，将区间$[l,r]$的答案通过加入一个值/删去一个值以得到$[l-1,r]$,$[l,r-1]$,$[l+1,r]$,$[l,r+1]$的答案的算法。并且，其复杂度为$(N+M)\ast \sqrt{N}$，所以当$N\ge 10^6$时，莫队算法就有TLE的可能。

说了那么多缺点，那么接下来就该大吹特吹了：
其优点是代码简单，考场上实用性强，是许多区间问题中bug级的存在，令出题人焦头烂额（谁也不想自己琢磨了半天搞出来的区间题，被人用莫队水过去）

先说说不带修改的莫队算法：
首先，莫队算法是基于分块进行的，将原串拆分为$\sqrt{N}$个块，每个块的大小自然也是$\sqrt{N}$级的，现在我们将所有询问，按照其左端点所在的块排序，左端在相同块中的，按照右端点排序。

假设当前已经得到了$[L_i,R_i]$的答案，我们通过依次插入和删除的方式，来得到$[L_{i+1},R_{i+1}]$
这看似十分暴力的算法，实则是有严格的复杂度证明的。

我们从$L$与$R$的偏移量来计算其复杂度：
首先是$L$的偏移量：因为我们总共需要处理M个询问，每两个相邻的询问之间，若$L$在同一个块中，那么其偏移量为$\sqrt{N}$，若不在同一个块中，其偏移量为$2\sqrt{N}$，所以$L$的复杂度为$M\sqrt{N}$

然后是$R$的偏移量：其实只要稍加分析就会发现，对于同一个块中，$R$是升序排列的，所以每个块的$R$的偏移量均为$N$，总计$\sqrt{N}$个块，所以$R$的偏移量为$N\sqrt{N}$

所以，莫队算法的复杂度就可以得到了：$(M+N)\sqrt{N}$

模板题：BZOJ2038小Z的袜子

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define SF scanf
#define PF printf
#define MAXN 50010
using namespace std;
int s[MAXN],n,m;
long long ans;
int b[MAXN];
int pos[MAXN];
struct node{
   
    int id;
    long long a,b,l,r;
}a[MAXN];
long long gcd(long long x,long long y){
   
    if(x==0)
        return y;
    return gcd(y%x,x);
}
bool cmp(node x,node y){
   
    return pos[x.l]<pos[y.l]||(pos[x.l]==pos[y.l]&&x.r<y.r);