专栏文章

莫队学习笔记

GGYHL
算法·理论参与者 1已保存评论 0

文章操作

快速查看文章及其快照的属性，并进行相关操作。

当前评论
0 条
当前快照
1 份
快照标识符
@mipb7c3t
此快照首次捕获于
2025/12/03 09:09
3 个月前
此快照最后确认于
2025/12/03 09:09
3 个月前
查看原文

莫队

莫队是个好东西

莫队简介

莫队算法是由莫队提出的。虽然他提出之前该算法已经在Codeforces中小范围流传，但是莫涛是第一个详细总结莫队算法的人。

莫队作用

莫队是一种离线处理静态区间查询的一类方法。时间复杂度为 O(nq+qlog⁡q)O(n\sqrt{q}+q\log q)，当 n,qn,q 同阶时复杂度为 O(nn)O(n\sqrt{n})。

普通莫队

形式

钦定 n,mn,m 同阶，那么对于区间询问问题，从 [l,r][l,r] 的答案 O(1)O(1) 扩展到 [l±1,r±1][l\pm1,r\pm1]（即与 [l,r][l,r] 相邻的区间）的答案，那么可以 O(nn)O(n\sqrt{n}) 求出所有询问的答案。详细看解释。

解释

离线后按照特定的方式排序，顺序处理每个询问，暴力从上一个区间的答案转移到下一个区间的答案（一步一步的将 l,rl,r 指针移动即可）。
那么有人就要问了：这样子的时间复杂度在 n,qn,q 同阶时不就是 O(n2)O(n^2) 吗？
实际上不是。注意到上面描述作用时，有一个 O(qlog⁡q)O(q\log q)。容易看出是排序的时间复杂度，也就是说，我们在指针移动之前要排序。

排序方法

我们把问题变成一个简单的理解就是在平面坐标系上游走，有 qq 个目标点 (li,ri)(l_i,r_i)，需要找到总长度最短的路径经过每个目标点一次。
很容易想到，如果 li,ril_i,r_i 分别单调递增，则用双指针即可 O(n)O(n) 维护这个问题。但是，如果我们放宽一下，允许 ll 在小范围内移动，rr 单增。其中，小范围看成分块。那么，排序方法就为按 ll 所在的块为第一关键字，rr 为第二关键字，将所有询问排序。
但是我们的时间复杂度得到优化了吗？如果块长为 TT，那么每次询问，ll 指针的移动量为 O(T)O(T) 的，总移动量 O(qT)O(qT)，rr 指针在 ll 不换块的时候是单调的，总移动量是 O(nT×n)O(\frac{n}{T}\times n)。根据均值不等式求出当 T=nqT=\frac{n}{\sqrt{q}} 时复杂度为 O(nq+qlog⁡q)O(n\sqrt{q}+q\log q)。
tips：关键字的意思是指从小到大排序时的依据，也就是按照谁的值来排序的。

代码

处理每个询问（先扩张在缩小，避免区间为空）

CPP
// 回答 [l,r]，当前在 [L,R]
while(L > l) add(-- L);
while(R < r) add(++ R);
while(L < l) del(L ++);
while(R > r) del(R --);

排序方法（重载）

CPP
struct node {
	int l, r, id;
	//b[x] 表示 x 所在的块
	bool operator < (const node &x) const {
		if (b[l] ^ b[x.l]) return l < x.l;
		else return r < x.r;
	}
};

优化

过程

这里我们要用到奇偶化排序。
什么是奇偶化排序？奇偶化排序即对于属于奇数块的询问，rr 按从小到大排序，对于属于偶数块的排序，rr 从大到小排序，这样我们的 rr 指针在处理完这个奇数块的问题后，将在返回的途中处理偶数块的问题，再向 nn 移动处理下一个奇数块的问题，优化了 rr 指针的移动次数。
建议自己模拟一下，加深理解。

代码

CPP
struct node {
  int l, r, id;
  bool operator<(const node &x) const {
    if (l / unit ^ x.l / unit) return l < x.l;
    if ((l / unit) & 1) return r < x.r;
    return r > x.r;
  }
};

习题

P1494 [国家集训队] 小 Z 的袜子

带修改莫队

普通莫队处理的都是静态问题，如果加上单点修改还能做吗？
充分发挥人类智慧。 我们可以认为是多了“时间”这一个维度，也就是询问形如 (l,r,t)(l,r,t) 表示询问进行了 tt 次单点修改之后，区间 [l,r][l,r] 的答案。
也就是变成了在一个三维的坐标系中游走，但是我们需要对这样的询问再后再一个排序的方式使其复杂度优于暴力。
其实我们只需要把时间也看成指针。特别地，如果加入的修改在当前区间里面，应当立刻对区间答案进行修改处理影响。
先按 lil_i 所在块的编号排序，第二关键字是 rir_i 所在块的编号，第三关键字才是 tit_i。
那么复杂度又是多少呢？
钦定 n,mn,m 同阶时，令块长 T=n23T=n^{\frac 23}，则有 n13n^{\frac 13} 个块。当 l,rl,r 所在的块固定时，时间指针只会移动 O(n)O(n) 次。而 l,rl,r 所在的块的组合仅有 O((n13)2)O((n^{\frac13})^2) 个，于是时间指针的中移动量为 O(n53)O(n^{\frac53})，l,rl,r 的总移动量都是 n13×(n23)2=n53n^{\frac13}\times(n^{\frac23})^2=n^{\frac53}，因此总时间复杂度是 O(n53)O(n^{\frac53})。
P1903 [国家集训队] 数颜色 / 维护队列

回滚莫队

我们发现无论是普通莫队还是带修改莫队，都可以快速地对区间进行扩展或者收缩。
如果强行想用莫队做区间最大值等增加或删除无法实现的问题怎么办？
当只增加和删除中的一个不可实现时，就可以使用回滚莫队 O(nq)O(n\sqrt{q}) 时间解决。

核心

只能使用一种操作，就用一种操作，剩下的用回滚莫队来解决。

形式

假设现在只能区间扩展，不能删除（只删除不扩展是类似的）。
那么，我们怎么样让 ll 这一坨只加入不删除呢？我们始终让 ll 保持在当前它所在的区间的右端点，然后 rr 已知移动到询问的 rir_i，然后 ll 再移动到 lil_i。
走了之后我们所做的操作不叫删除，叫回退。回退就是我们知道我们走的时候的状态，然后一步一步回到这个状态。
举个栗子：比如你修改了 [l,r][l,r] 的 max⁡\max 值，你在移动到 [l,r+1][l,r+1] 前备份你的最大值。

参考文献及鸣谢

本文部分内容参考oi-wiki上的内容。
特别感谢摸鱼酱Aurelia_Veilchenyixuan18提供的宝贵意见。
欢迎在下方留言哦。

评论

共 0 条评论，欢迎与作者交流。

正在加载评论...