莫队、带修莫队、树上莫队详解

这几天学习了莫队算法，试着写一篇比较详细的莫队教程吧...

普通莫队

简介

莫队是一种基于分块思想的离线算法，用于解决区间问题，适用范围如下：

1. 只有询问没有修改。

2. 允许离线。

3. 在已知询问 ([l,r]) 答案的情况下可以 (O(1)) 得到 ([l,r-1],[l,r+1],[l-1,r],[l+1,r]) 的答案。

满足以上三个条件就可以在 (O(nsqrt{m}+mlogm)) 的时间复杂度下得到每个询问的解。

算法思想

莫队的精髓就在于通过对询问进行排序，并把询问的结果作为下一个询问求解的基础，使得暴力求解的复杂度得到保证。

上文中“适用范围”的第三点“在已知询问 ([l,r]) 答案的情况下可以 (O(1)) 得到 ([l,r-1],[l,r+1],[l-1,r],[l+1,r]) 的答案”即是“把询问的结果作为下一个询问求解的基础”的方法。

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

在这题中，用 (cnt_i) 表示当前处理的区间内颜色为i的袜子出现的次数，用 (mathrm{len}) 表示当前处理的区间的长度，用 (x) 表示新增的那只袜子的颜色。

以已知区间 ([l,r]) 的答案求解区间 ([l,r+1]) 为例。分别处理分子和分母：

• 分母为任选两只袜子的组合总数，原先是 (frac{mathrm{len}*(mathrm{len}-1)}{2})，现在是 (frac{mathrm{len}*(mathrm{len}+1)}{2})，增加了 (mathrm{len}) 。
• 分子为两只袜子颜色相同的组合总数，比原来增加了 (cnt_x)，即新增的这只袜子和原本就在当前区间内的相同颜色的袜子的组合。

因此，将一只颜色为x的袜子计入答案的函数就可以写出来了：

//fz代表分子，fm代表分母
void add(int x)
{
	fz+=cnt[x];
	++cnt[x];
	fm+=len;
	++len;
}

同理可以写出将一只颜色为x的袜子移出答案的函数：

void del(int x)
{
	--cnt[x];
	fz-=cnt[x];
	--len;
	fm-=len;
}

于是，我们就可以得到一个暴力的算法：用 (l) 和 (r) 分别记录当前区间的两个端点，然后用下面这段代码来更新答案（q[i].l,q[i].r代表正在处理的询问的两个端点，col[p]代表第 (p) 只袜子的颜色）：

while (l>q[i].l)
{
	add(col[--l]);
}
while (r<q[i].r)
{
	add(col[++r]);
}
while (l<q[i].l) 
{
	del(col[l++]);
}
while (r>q[i].r)
{
	del(col[r--]);
}

然而，这个算法的时间复杂度是 (O(nm)) 的（因为最坏情况下每次 (l) 和 (r) 两个指针都要走 (O(n)) 的距离，而一共有 (m) 次询问），和暴力完全一样甚至跑的更慢。

别忘了，之前我说过，莫队的精髓就在于通过对询问进行排序，使得暴力求解的复杂度得到保证。

我们的目的是使 (l) 和 (r) 两个指针走过的总距离尽量的小，这时候就要用到分块的思想了。

把整个区间 ([1,n]) 分成若干块，以询问的左端点所在块为第一关键字，以询问的右端点大小为第二关键字，对询问进行排序，那么：

• 对于同一块的询问，(l) 指针每次最多移动块的大小，(r) 指针的移动则是单调的，总共移动最多 (n) 。
• 对于不同块的询问，(l) 每次换块时最多移动两倍块的大小， (r) 每次换块时最多移动 (n) 。

总结：（用 (B) 表示块的大小）(l) 指针每次移动 (O(B))，(r) 指针每块移动 (O(n)) 。

所以：

• (l) 的移动次数最多为询问数×块的大小，即 (O(mB)) 。
• (r) 的移动次数最多为块的个数×总区间大小，即 (O(n^2/B)) 。

因此，总移动次数为 (O(mB+n^2/B)) 。

没错，这就是个双勾函数，所以当 (B=sqrt{frac{n^2}{m}}) 即 (frac{n}{sqrt{m}}) 时复杂度最小，为 (O(nsqrt{m})) 。

剩下的最后一个问题：初始的当前区间是什么？

只要任意指定一个空区间就好了，如 (l=1,r=0) 。

所以，整个莫队算法就可以概括为：

1. 将询问记录下来。

2. 以 (frac{n}{sqrt{m}}) 为块的大小，以询问的左端点所在块为第一关键字，以询问的右端点大小为第二关键字，对询问进行排序。

3. 暴力处理每个询问。

4. 输出答案。

总的复杂度为 (O(nsqrt{m}+mlogm)) 。

P.S. 网上很多教程说分块大小取 (sqrt{n}) 最优，复杂度为 (O(nsqrt{n}))，这是不严谨的，当n、m差别较大时使用 (sqrt{n}) 作为分块大小效率会明显偏低。

例题代码

[国家集训队]小Z的袜子 AC代码：

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>

using namespace std;

const int N=50010;

void add(int x);
void del(int x);
int gcd(int a,int b);

int n,m,B,fz,fm,len,col[N],cnt[N],ans[N][2];

struct Query
{
    int l,r,id;
    bool operator<(Query& b)
    {
        return l/B==b.l/B?r<b.r:l<b.l;
    }
} q[N];

int main()
{
    int i,l=1,r=0,g;
    
    cin>>n>>m;
    
    B=n/sqrt(m);
    
    for (i=1;i<=n;++i)
    {
        cin>>col[i];
    }
    
    for (i=0;i<m;++i)
    {
        cin>>q[i].l>>q[i].r;
        q[i].id=i;
    }
    
    sort(q,q+m);
    
    for (i=0;i<m;++i)
    {
        if (q[i].l==q[i].r)
        {
            ans[q[i].id][0]=0;
            ans[q[i].id][1]=1;
            continue;
        }
        while (l>q[i].l)
        {
            add(col[--l]);
        }
        while (r<q[i].r)
        {
            add(col[++r]);
        }
        while (l<q[i].l) 
        {
            del(col[l++]);
        }
        while (r>q[i].r)
        {
            del(col[r--]);
        }
        g=gcd(fz,fm);
        ans[q[i].id][0]=fz/g;
        ans[q[i].id][1]=fm/g;
    }
    
    for (i=0;i<m;++i)
    {
        printf("%d/%d
",ans[i][0],ans[i][1]);
    }
    
    return 0;
}

void add(int x)
{
    fz+=cnt[x];
    ++cnt[x];
    fm+=len;
    ++len;
}

void del(int x)
{
    --cnt[x];
    fz-=cnt[x];
    --len;
    fm-=len;
}

int gcd(int a,int b)
{
    return b==0?a:gcd(b,a%b);
}

其它例题

小B的询问

带修莫队

前面说过，普通的莫队只能解决没有修改的问题，那么带修改的问题怎么解决呢？带修莫队就是一种支持单点修改的莫队算法。

算法简介

还是对询问进行排序，每个询问除了左端点和右端点还要记录这次询问是在第几次修改之后（时间），以左端点所在块为第一关键字，以右端点所在块为第二关键字，以时间为第三关键字进行排序。

暴力查询时，如果当前修改数比询问的修改数少就把没修改的进行修改，反之回退。

需要注意的是，修改分为两部分：

1. 若修改的位置在当前区间内，需要更新答案（del原颜色，add修改后的颜色）。

2. 无论修改的位置是否在当前区间内，都要进行修改（以供add和del函数在以后更新答案）。

分块大小的选择以及复杂度证明

（用 (B) 表示分块大小，(c) 表示修改个数，(q) 表示询问个数，l块表示以 (l/B) 分的块，r块表示以 (r/B) 分的块，每个l块包含 (n/B) 个r块）

1. 对于时间指针 (now)：对于每个r块，最坏情况下会移动 (c)，共有 (left(frac{n}{B} ight)^2) 个r块，所以总移动次数为 (frac{cn^2}{B^2}) 。

2. 对于左端点指针 (l) ：l块内移动每次最多 (B)，换l块每次最多 (2B)，所以总移动次数为 (O(qB)) 。

3. 对于右端点指针 (r)：r块内移动每次最多 (B)，换r块每次最多 (2B)，所有l块内移动次数之和为 (O(qB))；换l块时最多移动 (n)，总的换l块时移动次数为 (Oleft(frac{n^2}{B} ight))；所以总的移动次数为 (Oleft(qB+frac{n^2}{B} ight)) 。

所以：总移动次数为 (Oleft(frac{cn^2}{B^2}+qB+frac{n^2}{B} ight)) 。

由于一般的题目都不会告诉你修改和询问分别的个数，所以统一用 (m) 表示，即 (Oleft(frac{mn^2}{B^2}+mB+frac{n^2}{B} ight)) 。

那么 (B) 取多少呢...Mathematica告诉我大约是这个 （图炸不补，就是一个很复杂的式子）。

所以还是不要纠结带修莫队的最佳分块大小好了...视作 (n=m) 的话，就可以得到总移动次数为 (Oleft(frac{n^3}{B^2}+nB+frac{n^2}{B} ight))，那么 (B=n^{frac{2}{3}}) 时取最小值 (Oleft(n^{frac{5}{3}} ight)) 。

所以：带修莫队的渐进时间复杂度为 (Oleft(nlogn+n^{frac{5}{3}} ight)) （视作 (n=m)）。

例题代码

这次就不详细分析例题了，直接上代码。

[国家集训队]数颜色 AC代码：

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>

using namespace std;

void add(int x);
void del(int x);
void modify(int x,int ti);  //这个函数会执行或回退修改ti（执行还是回退取决于是否执行过，具体通过swap实现），x表明当前的询问是x，即若修改了区间[q[x].l,q[x].r]便要更新答案

int n,m,B,cnt[1000010],a[50010],ans,ccnt,qcnt,now,out[50010];

struct Change
{
    int p,col;
} c[50010];

struct Query
{
    int l,r,t,id;
    bool operator<(Query& b)
    {
        return l/B==b.l/B?(r/B==b.r/B?t<b.t:r<b.r):l<b.l;
    }
} q[50010];

int main()
{
    int i,l=2,r=1;
    char type[10];
    
    cin>>n>>m;
    B=pow(n,0.66666);
        
    for (i=1;i<=n;++i)
    {
        cin>>a[i];
    }
    
    for (i=1;i<=m;++i)
    {
        scanf("%s",type);
        if (type[0]=='Q')
        {
            ++qcnt;
            cin>>q[qcnt].l>>q[qcnt].r;
            q[qcnt].t=ccnt;
            q[qcnt].id=qcnt;
        }
        else
        {
        	++ccnt;
        	cin>>c[ccnt].p>>c[ccnt].col;
        }
    }
    
    sort(q+1,q+qcnt+1);
    
    for (i=1;i<=qcnt;++i)
    {
        while (l>q[i].l)
        {
            add(a[--l]);
        }
        while (r<q[i].r)
        {
            add(a[++r]);
        }
        while (l<q[i].l)
        {
            del(a[l++]);
        }
        while (r>q[i].r)
        {
            del(a[r--]);
        }
        while (now<q[i].t)
        {
            modify(i,++now);
        }
        while (now>q[i].t)
        {
            modify(i,now--);
        }
        out[q[i].id]=ans;
    }
    
    for (i=1;i<=qcnt;++i)
    {
    	cout<<out[i]<<endl;
    }
    
    return 0;
}

void add(int x)
{
    if (cnt[x]++==0)
    {
        ++ans;
    }
}

void del(int x)
{
    if (--cnt[x]==0)
    {
        --ans;
    }
}

void modify(int x,int ti)
{
    if (c[ti].p>=q[x].l&&c[ti].p<=q[x].r)
    {
        del(a[c[ti].p]);
    	add(c[ti].col);
    }
    swap(a[c[ti].p],c[ti].col); //下次执行时必定是回退这次操作，直接互换就可以了 
}

其它例题

CF940F Machine Learning

树上莫队

其实，莫队算法除了序列还可以用于树。复杂度同序列上的莫队（不带修 (O(nsqrt{m}+mlogm))，带修 (Oleft(nlogn+n^{frac{5}{3}} ight))）。

例题：[WC2013]糖果公园

分块方式

这里需要看一道专门为树上莫队设计的题目 [SCOI2005]王室联邦。

用这道题所要求的方式进行分块，并用后文的方式更新答案，就能保证复杂度（复杂度分析见后文）。

那么如何满足每块大小在 ([B,3B])，块内每个点到核心点路径上的所有点都在块内呢？

这里先提供一种构造方式，再予以证明：

dfs，并创建一个栈，dfs一个点时先记录初始栈顶高度，每dfs完当前节点的一棵子树就判断栈内（相对于刚开始dfs时）新增节点的数量是否>=B，是则将栈内所有新增点分为同一块，核心点为当前dfs的点，当前节点结束dfs时将当前节点入栈，整个dfs结束后将栈内所有剩余节点归入已经分好的最后一个块。

参考代码：

void dfs(int u,int fa)
{
    int t=top;
    for (int i=head[u];i;i=nxt[i])
    {
        int v=to[i];
        if (v!=fa)
        {
            dfs(v,u);
            if (top-t>=B)
            {
                ++tot;
                while (top>t)
                {
                    bl[sta[top--]]=tot;
                }
            }
        }
    }
    sta[++top]=u;
}

dfs(1,0);

while (top)
{
    bl[sta[top--]]=tot;
}

如果你看懂了这个方法的话，每块大小>=B是显然的，下面证明为何每块大小<=3B：

对于当前节点的每一棵子树：

• 若未被分块的节点数>B，那么在dfs这棵子树的根节点时就一定会把这棵子树的一部分分为一块直至这棵子树的剩余节点数<=B，所以这种情况不存在。

• 若未被分块的节点数=B，这些节点一定会和栈中所有节点分为一块，栈中之前还剩 ([0,B-1]) 个节点，那么这一块的大小为 ([B,2B-1]) 。

• 若未被分块的节点数<B，当未被分块的节点数+栈中剩余节点数>=B时，这一块的大小为 ([B,2B-1))，否则继续进行下一棵子树。

对于dfs结束后栈内剩余节点，数量一定在 ([1,B]) 内，而已经分好的每一块的大小为 ([B,2B-1])，所以每块的大小都在 ([B,3B)) 内（我看有的博客写的剩余节点数量在 ([1,B+1]) 内，所以最后一块可能达到 (3B)...然而我觉得最多 (3B-1)啊QAQ）。

修改方式

所谓“修改”，就是由询问 ((cu,cv)) 更新至询问 ((tu,tv)) 。

如果把两条路径上的点全部修改..显然是和暴力一样的嘛！

这里直接给出结论好了...

（下文中 (T(u,v)) 表示 (u) 到 (v) 的路径上除 (lca(u,v)) 外的所有点构成的集合，(S(u,v)) 代表u到v的路径，(xor) 表示集合对称差（就跟异或差不多））

• 两个指针 (cu,cv) （相当于序列莫队的 (l,r) 两个指针）， (ans)记录(T(cu,cv)) 的答案，(vis) 数组记录每个节点是否在 (T(cu,cv)) 内；

• 由 (T(cu,cv)) 更新至 (T(tu,tv)) 时，将 (T(cu,tu)) 和 (T(cv,tv)) 的 (vis) 分别取反，并相应地更新答案；

• 将答案记录到 (out) 数组（离线后用于输出那个）时对 (lca(cu,cv)) （此时的 (cu,cv) 已更新为上一步中的 (tu,tv)） 的 (vis) 取反并更新答案，记录完再改回来（因为lca比较烦，所以就这样做了QAQ）。

第二步证明如下：

(quad\,T(cu,cv) xor T(tu,tv))

(=[S(cu,root) xor S(cv,root)] xor [S(tu,root) xor S(tv,root)]) （lca及以上相消）

(=[S(cu,root) xor S(tu,root)] xor [S(cv,root) xor S(tv,root)]) （交换律、结合律）

(=T(cu,tu) xor T(cv,tv))

之所以要把 (T(cu,cv) xor T(tu,tv)) 转化成 (T(cu,tu) xor T(cv,tv))，是因为这样的话就能通过对询问排序来保证复杂度。

关于单点修改

树上莫队的单点修改和序列莫队类似，唯一不同就是，修改后是否更新答案通过vis数组判断。

复杂度分析

每块大小在 ([B,3B))，所以两点间路径长度也在 ([B,3B))，块内移动就是 (O(B)) 的；编号相邻的块位置必然是相邻的，所以两块间路径长度也是 (O(B))；然后就和序列莫队的复杂度分析类似了...

例题代码

[WC2013]糖果公园

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>

using namespace std;

const int N=100010;

void pathmodify(int u,int v); //将T(u,v)取反并更新答案
void opp(int x); //将节点x取反并更新答案
void modify(int ti); //进行或回退修改ti
int lca(int u,int v);
void dfs(int u); //进行分块并记录dep数组、f数组（用于求lca、两点间路径）
void add(int u,int v);

int head[N],nxt[N<<1],to[N<<1],cnt;
int n,m,Q,B,bl[N],tot,V[N],W[N],a[N],sta[N],top,qcnt,ccnt,dep[N],f[20][N],num[N],now;
long long ans,out[N];
bool vis[N];

struct Query
{
	int u,v,t,id;
	bool operator<(Query& y)
	{
		return bl[u]==bl[y.u]?(bl[v]==bl[y.v]?t<y.t:bl[v]<bl[y.v]):bl[u]<bl[y.u];
	}
} q[N];

struct Change
{
	int p,x;
} c[N];

int main()
{
	int i,j,u,v,lc,type;
	
	cin>>n>>m>>Q;
	B=pow(n,0.666);
	
	for (i=1;i<=m;++i)
	{
		cin>>V[i];
	}
	
	for (i=1;i<=n;++i)
	{
		cin>>W[i];
	}
	
	for (i=1;i<n;++i)
	{
		cin>>u>>v;
		add(u,v);
		add(v,u);
	}
	
	dfs(1);
	
	for (i=1;i<=16;++i)
	{
		for (j=1;j<=n;++j)
		{
			f[i][j]=f[i-1][f[i-1][j]];
		}
	}
	
	while (top)
	{
		bl[sta[top--]]=tot;
	}
	
	for (i=1;i<=n;++i)
	{
		cin>>a[i];
	}
	
	for (i=0;i<Q;++i)
	{
		cin>>type;
		if (type==0)
		{
			++ccnt;
			cin>>c[ccnt].p>>c[ccnt].x;
		}
		else
		{
			cin>>q[qcnt].u>>q[qcnt].v;
			q[qcnt].t=ccnt;
			q[qcnt].id=qcnt;
			++qcnt;
		}
	}
	
	sort(q,q+qcnt);
	
	u=v=1;
	
	for (i=0;i<qcnt;++i)
	{
		pathmodify(u,q[i].u);
		pathmodify(v,q[i].v);
		u=q[i].u;
		v=q[i].v;
		while (now<q[i].t)
		{
			modify(++now);
		}
		while (now>q[i].t)
		{
			modify(now--);
		}
		lc=lca(u,v);
		opp(lc);
		out[q[i].id]=ans;
		opp(lc);
	}
	
	for (i=0;i<qcnt;++i)
	{
		cout<<out[i]<<endl;
	}
	
	return 0;
}

void pathmodify(int u,int v)
{
	if (dep[u]<dep[v])
	{
		swap(u,v);
	}
	while (dep[u]>dep[v])
	{
		opp(u);
		u=f[0][u];
	}
	while (u!=v)
	{
		opp(u);
		opp(v);
		u=f[0][u];
		v=f[0][v];
	}
}

void opp(int x)
{
	if (vis[x])
	{
		ans-=1ll*V[a[x]]*W[num[a[x]]--];
	}
	else
	{
		ans+=1ll*V[a[x]]*W[++num[a[x]]];
	}
	vis[x]^=1;
}

void modify(int ti)
{
	if (vis[c[ti].p])
	{
		opp(c[ti].p);
		swap(a[c[ti].p],c[ti].x);
		opp(c[ti].p);
	}
	else
	{
		swap(a[c[ti].p],c[ti].x);
	}
}

int lca(int u,int v)
{
	if (dep[u]<dep[v])
	{
		swap(u,v);
	}
	int i;
	for (i=0;i<=16;++i)
	{
		if ((dep[u]-dep[v])&(1<<i))
		{
			u=f[i][u];
		}
	}
	if (u==v)
	{
		return u;
	}
	for (i=16;i>=0;--i)
	{
		if (f[i][u]!=f[i][v])
		{
			u=f[i][u];
			v=f[i][v];
		}
	}
	return f[0][u];
}

void dfs(int u)
{
	int t=top;
	for (int i=head[u];i;i=nxt[i])
	{
		int v=to[i];
		if (v!=f[0][u])
		{
			f[0][v]=u;
			dep[v]=dep[u]+1;
			dfs(v);
			if (top-t>=B)
			{
				++tot;
				while (top>t)
				{
					bl[sta[top--]]=tot;
				}
			}
		}
	}
	sta[++top]=u;
}

void add(int u,int v)
{
	nxt[++cnt]=head[u];
	head[u]=cnt;
	to[cnt]=v;
}

莫队的扩展

其实莫队可以扩展到高维，参见二维莫队解题报告。

更一般地，若 (Q(x_1,x_2,cdots,x_k)) 为一个询问，(forall iin[1,k])，(x_i) 的规模都为 (n)，可以在时间 (mathrm{T}) 内求解 (Q(x_1,x_2,cdots,x_ipm 1,cdots,x_n))，共有 (m) 个询问，那么就可以在 (Oleft(kmlogm+nTm^frac{k-1}{k} ight)) 的时间复杂度下离线求解。

（蒟蒻的大胆猜想而已..并没有严格证明）

To be finished:回滚莫队（只增莫队）..（有时间再填坑吧Orz）