莫队算法入门

最新推荐文章于 2022-09-19 20:35:35 发布
最新推荐文章于 2022-09-19 20:35:35 发布
阅读量255 收藏 1
点赞数 1
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Emma_goddess/article/details/108677796
版权

昨天重温了一下CaptainMo的职业生涯 莫队的模板，看了下别人的博客，把三个板子打了，做练习前先小小总结了一下吧。
一. 基础莫队算法
莫队算法 = 离线 + 暴力 + 分块，它通常用于不修改只查询的一类区间问题，复杂度为在这里插入图片描述
主要就是通过排序过后再处理询问能优化暴力，排序则是利用分块，至于为什么更优，附张别人博客看到的图。
在这里插入图片描述
这就很显然图二的走法更短，然后 编码时，还可以对排序做一个小优化：奇偶性排序，让奇数块和偶数块的排序相反。例如左端点L都在奇数块，则对R从大到小排序；若L在偶数块，则对R从小到大排序（反过来也可以：奇数块从小到大，偶数块从大到小）。
luogu P1972 HH的项链代码，这代码不能拿全分(不会卡常的屑 )，但正确性无误

#include<bits/stdc++.h>
const int N=1e6+10;
using namespace std;
int n,m,a[N],cnt[N],ans[N],block,id[N],t,L,R,sum;
struct ask
{
	int l;
	int r;
	int Id;
}q[N];
inline int Read()
{
	int x=0,f=1;
	char ch;
	while(ch>'9'||ch<'0')
	{
		if(ch=='-')	f=-1;
		ch=getchar();
	}
	while(ch>='0'&&ch<='9')
	{
		x=x*10+ch-'0';
		ch=getchar();
	}
	return f*x;
}
bool comp(ask a,ask b)
{
	if(id[a.l]!=id[b.l]) return id[a.l]<id[b.l];
	if((id[a.l]%2)==1) return a.r>b.r;
	else return a.r<b.r; 
}
void add(int x)
{
	x=a[x];
	if(!cnt[x]) sum++;
	cnt[x]++;
}
void del(int x)
{
	x=a[x];
	cnt[x]--;
	if(!cnt[x]) sum--;
}
int main()
{
	n=Read();
	block=sqrt(n)+1;
	t=n/block;
	if(n%block) t++;
	for(int i=1;i<=n;i++) a[i]=Read(),id[i]=(i-1)/block+1;
	m=Read();
	for(int i=1;i<=m;i++)
	{
		q[i].l=Read(),q[i].r=Read();
		q[i].Id=i;
	}
	sort(q+1,q+1+m,comp);
	L=1,R=0;
	for(int i=1;i<=m;i++)
	{
		while(L<q[i].l) del(L),L++;
		while(L>q[i].l) L--,add(L);
		while(R<q[i].r) R++,add(R);
		while(R>q[i].r) del(R),R--;
		ans[q[i].Id]=sum;
	}
	for(int i=1;i<=m;i++) printf("%d\n",ans[i]);
	return 0;
}

二. 带修改的莫队
如果是比较简单的“单点修改”，也能应用莫队算法，复杂度在这里插入图片描述
还是有道例题luogu P1903数颜色
如果用莫队算法求解，必须离线，先把查询操作和修改操作分别记录下来。记录查询操作的时候，增加一个变量，记录本次查询前做少次修改。
.

如果没有修改，就是基础莫队，一个查询的左右端点是[L, R]。加上修改之后，一个查询表示为(L, R, t)，t表示在查询[L, R]前进行次修改操作。可以把t理解为“时间”，t的范围是1 ≤ t ≤ m，m是操作次数。
.
反正就是要多考虑修改操作对查询的影响，又因为是单点修改，可以直接暴力，查询中就多了这么一截

void CaptainMo()
{
	L=1,R=now=0;
	for(int i=1;i<=qnum;i++)
	{
		while(L<q[i].l) del(a[L]),L++;
		while(L>q[i].l) L--,add(a[L]);
		while(R<q[i].r) R++,add(a[R]);
		while(R>q[i].r) del(a[R]),R--;
		while(now>q[i].pre) update(now,i),now--;//需要还原修改
		while(now<q[i].pre) now++,update(now,i);//需要继续修改
		ans[q[i].Id]=sum;
	}
}

然后就是一些小细节的不同，比如分块大小应该为在这里插入图片描述
因为这样更快。排序时右端点r也按所在块的序号大小排，不是r大小。还有就是在修改时要注意写这句swap(a[x],p[pos].val);为什么的话模拟一下就知道了，当时理解了好久才看懂。
代码

#include<bits/stdc++.h>
using namespace std;
const int N=1400000;
int n,m,L,R,x,y,now,id[N],block,a[N],ans[N],cnt[N],qnum,pnum,sum;
struct ask
{
	int l;
	int r;
	int Id;
	int pre;
}q[N];
struct upd
{
	int st;
	int val;
}p[N];
inline int Read()
{
	int x=0,f=1;
	char ch;
	while(ch>'9'||ch<'0')
	{
		if(ch=='-')	f=-1;
		ch=getchar();
	}
	while(ch>='0'&&ch<='9')
	{
		x=x*10+ch-'0';
		ch=getchar();
	}
	return f*x;
}
bool comp(ask A,ask b)
{
	if(id[A.l]!=id[b.l]) return id[A.l]<id[b.l];
	if(id[A.r]!=id[b.r]) return id[A.r]<id[b.r];
	return A.pre<b.pre;
}
void add(int x)
{
	if(!cnt[x]) sum++;
	cnt[x]++;
}
void del(int x)
{
	cnt[x]--;
	if(!cnt[x]) sum--;
}
void update(int pos,int i)
{
	int x=p[pos].st;
	if(x>=q[i].l&&x<=q[i].r)
	{
    	cnt[a[x]]--;
     	if(!cnt[a[x]]) sum--;
		if(!cnt[p[pos].val]) sum++;
		cnt[p[pos].val]++;		
	}
	swap(a[x],p[pos].val);
}
void CaptainMo()
{
	L=1,R=now=0;
	for(int i=1;i<=qnum;i++)
	{
		while(L<q[i].l) del(a[L]),L++;
		while(L>q[i].l) L--,add(a[L]);
		while(R<q[i].r) R++,add(a[R]);
		while(R>q[i].r) del(a[R]),R--;
		while(now>q[i].pre) update(now,i),now--;
		while(now<q[i].pre) now++,update(now,i);
		ans[q[i].Id]=sum;
	}
}
int main()
{
	n=Read(),m=Read();
	double X=pow(n,2.0/3.0);
	block=(int)X;
	for(int i=1;i<=n;i++) a[i]=Read(),id[i]=(i-1)/block+1;
	for(int i=1;i<=m;i++)
	{
		char s;
		cin>>s,x=Read(),y=Read();
		if(s=='Q')
		q[++qnum].l=x,q[qnum].r=y,q[qnum].Id=qnum,q[qnum].pre=pnum;
		else
		p[++pnum].st=x,p[pnum].val=y;
	}
	sort(q+1,q+1+qnum,comp);
	CaptainMo();
	for(int i=1;i<=qnum;i++) printf("%d\n",ans[i]);
	return 0;
}

三. 树上莫队
基础莫队和带修改的莫队操作的都是一维数组。基于其他的数据结构的问题，如果能转换成一维数组而且是区间问题，那么也能应用莫队算法。
典型的例子是树形结构上的路径问题，可以利用“欧拉序”把整棵树的结点顺序转化为一个一维数组，路径问题也变成了区间问题，就能利用莫队算法求解。下面还是有道例题。
SP10707 COT2 - Count on a tree II

相信学过树剖后对树上的节点存到队列已经熟悉，然后就考虑(u, v)上的路径有哪些结点？首先计算出u、v的lca(u, v)（最近公共祖先），然后讨论两种情况：
（1）lca(u, v) = u或lca(u, v) = v，这种情况最简单，可以直接做，不多说了。
（2）lca(u, v) ≠ u且lca(u, v) ≠ v，这种情况其实也不复杂，就在查询时多查询一个lca再计算答案，然后再查询一次消除lca对目前的sum值影响。
至于如何忽略掉区间内出现了两次的点，这个很简单，我们多记录一个vis[x]，表示x这个点有没有被加入，每次处理的时候如果vis[x]=0则需要添加节点；如果vis[x]=1则需要删除节点。
代码

#include<bits/stdc++.h>
using namespace std;
const int N=1e5+100;
int vis[N],L,R,cnt[N],f[40010][20],n,m,sum,b[N],a[N],len,ans[N],idx[N],Idx[N],ref[N],index_,id[N],block,dep[N],first[N],Next[N],to[N],tot;
struct ask
{
	int l;
	int r;
	int id;
	int zx;
}q[N];
int Read()
{
	int x=0,f=1;
	char ch=getchar();
	while(!isdigit(ch))
	{
		if(ch=='-') f=-1;
		ch=getchar();
	} 
	while(isdigit(ch))
	{
		x=(x<<1)+(x<<3)+ch-'0';
		ch=getchar();
	}
	return f*x;
}
void add(int a,int b)
{
	tot++;
	Next[tot]=first[a];
	to[tot]=b;
	first[a]=tot;
}
void dfs(int u,int fa)
{
	index_++;
	idx[u]=index_;
	ref[index_]=u;
	dep[u]=dep[fa]+1;
	for(int i=1;i<20;i++) f[u][i]=f[f[u][i-1]][i-1];
	for(int i=first[u];i;i=Next[i])
	{
		int v=to[i];
		if(v==fa) continue;
		f[v][0]=u;
		dfs(v,u);
	}
	Idx[u]=++index_,ref[index_]=u;
}
int lca(int u,int v)
{
	if(dep[u]<dep[v]) swap(u,v);
	for(int i=19;i>=0;i--) if(dep[f[u][i]]>=dep[v]) u=f[u][i];
	if(u==v) return u;
	for(int i=19;i>=0;i--) if(f[u][i]!=f[v][i]) u=f[u][i],v=f[v][i];
	return f[u][0];
}
bool comp(ask a,ask b)
{
	if(id[a.l]!=id[b.l]) return id[a.l]<id[b.l];
	if(id[a.l]%2==1) return a.r>b.r;
	else return a.r<b.r;
}
void update(int x)
{
	vis[x]^=1;
	if(vis[x]) 
	{
		cnt[a[x]]++;
		if(cnt[a[x]]==1) sum++;
	}
	else
	{
		cnt[a[x]]--;
		if(!cnt[a[x]]) sum--;
	}
}
void CaptainMo()
{
	L=1,R=0;
	for(int i=1;i<=m;i++)
	{
		while(L>q[i].l) L--,update(ref[L]);
		while(L<q[i].l) update(ref[L]),L++;
		while(R>q[i].r) update(ref[R]),R--;
		while(R<q[i].r) R++,update(ref[R]);
		if(q[i].zx) update(q[i].zx);
		ans[q[i].id]=sum;
		if(q[i].zx) update(q[i].zx);		
	}
}
int main()
{
	n=Read(),m=Read();
	block=(int)pow(n,2.0/3.0);
	for(int i=1;i<=n;i++) a[i]=b[i]=Read();
	for(int i=1;i<=2*n;i++ ) id[i]=(i-1)/block+1;
	sort(b+1,b+1+n);
	len=unique(b+1,b+1+n)-b-1;
	for(int i=1;i<=n;i++) a[i]=lower_bound(b+1,b+1+len,a[i])-b;
	for(int i=1;i<n;i++) 
	{
		int x=Read(),y=Read();
		add(x,y),add(y,x);
	}
	dfs(1,0);
//	for(int i=1;i<=index_;i++) cout<<ref[i]<<" ";
	for(int i=1;i<=m;i++)
	{
		int x=Read(),y=Read();
		int z=lca(x,y);
		if(z==x||z==y)
		{
			if(idx[x]>idx[y]) swap(x,y);
			q[i].l=idx[x],q[i].r=idx[y];
		}
		else
		{
			if(idx[x]>Idx[y]) swap(x,y);
			q[i].l=Idx[x],q[i].r=idx[y];
			q[i].zx=z;
		}
		q[i].id=i;
	}
	sort(q+1,q+1+m,comp);
	CaptainMo();
	for(int i=1;i<=m;i++) printf("%d\n",ans[i]);
	return 0;
}

该去做练习了啊…

kmouyito
关注
莫队算法——从入门到黑题
weixin_30512043的博客
12-14 1918
众所周知，莫队是由莫涛大神提出的，一种玄学毒瘤暴力骗分区间操作算法，它以简短的框架、简单易记的板子和优秀的复杂度闻名于世。然而由于莫队算法应用的毒瘤，很多可做的莫队模板题都有着较高的难度评级，令很多初学者望而却步。然而，如果你真正理解了莫队的算法原理，那么它用起来还是很简单的。当然某些套左套右的毒瘤除外 0、前置芝士： 莫队算法还是比较独立的。不过你还是得了解了解以下的一些知识： \(1\)、分...
莫队算法入门 + 模板 Codeforces 617E
Rain Stopped的博客
05-13 2378
题意已知一个长度为n的数列 (0 ≤ ai ≤ 1 000 000) ，给m个区间，问每个区间有多少个子区间xor和为k  (1 ≤ n, m ≤ 100 000, 0 ≤ k ≤ 1 000 000) 莫队算法  如果你知道了[L,R]的答案。你可以在O(1)的时间下得到[L,R-1]和[L,R+1]和[L-1,R]和[L+1,R]的答案的话。就可以使用莫队算法。 先对序列分块。然后
莫队算法入门 BZOJ 2038
fufck
04-01 1154
目前的题型概括为三种：普通莫队，树形莫队以及带修莫队。下面介绍普通莫队算法：普通莫队算法是一种离线算法，不带修改 操作。其通过对询问操作的执行顺序进行更改，然后使用最暴力的方法，可以达到很好的复杂度。首先,如果要用莫队算法，则必须满足已知ans[l,r]可以得到ans[l,r+1]，ans[l,r-1],ans[l+1,r],ans[l-1,r]莫队算法的实现步骤为： 1、先对原序列进行分块。(s...
莫队算法 入门理解
立志变成java巨头
08-15 396
目前的题型概括为三种：普通莫队，树形莫队以及带修莫队。 今天主要讲的是普通莫队算法，莫队算法本来就是一个莽夫算法，但是他就是对查询区间排个序后就优化很多了。         莫队算法巧妙地将询问离线排序，使得其复杂度无比美妙……”在一般做题时我们时常遇到使用排序来优化枚举时间消耗的例子。莫队的优化基于分块思想：对于两个询问，若在其l在同块，那么将其r作为排序关键字，若l不在同块，就将l作为关键...
莫队算法入门 Codeforces617E
LinzhiQQQ的博客
09-03 540
莫队算法： 莫队算法的用处是，对于一个区间内的查询，当我们已经知道了 [ L , R ] 的答案的时候，有莫队算法可以在很短的时间内得到 [ L - 1 , R ] 或者是 [ L ,R + 1 ]的答案，前提是可以离线处理。 步骤： 1.输入 N 个点的时候，对每个点的进行分块 ，一般 N 个点分成 sqrt(n) 块 ，用一个pos[]数组记录分块 int s = sqrt(n); f
普通、带修(可持久化)莫队算法入门例题详解
菜鸡成长史
05-24 724
目录 【莫队算法】 【普通莫队】 【代码】 【题面】 【带修莫队】 【代码】 【题面】 【总结】 【莫队算法】 参考大米饼的莫队算法，目前的题型概括为三种：普通莫队，带修莫队以及树形莫队。 【普通莫队】 例题：2038: [2009国家集训队]小Z的袜子(hose) 题意：给定编号1-n的n只袜子的颜色，输出从询问的区间[L,R]中随机抽出两只袜子颜色相同的概率。...
莫队算法总结&专题训练1
Plozia 的小窝
12-08 530
1.概述 莫队算法，是由莫涛队长提出的一种，能够以玄学的复杂度来处理区间查询类的问题。 甲：区间查询类的问题不是可以用线段树等数据结构解决的吗？ 乙：那如果要求某个区间的区间众数要怎么办呢？ 甲：啊这。。。。。。 所以，莫队算法就是用来解决这种线段树等数据结构不好维护的区间查询问题。 这里先安利几个博客，写的非常好，建议各位读者可以看一看，写的比我好多了： 莫队算法——从入门到黑题（强烈推荐） 洛谷日报 #48 期：莫队算法初探 洛谷日报 #183 期：你以为莫队只能离线？莫队的在线化改造（难度较高，建议
莫队 从零基础到入门 超详细
glorious_dream的博客
09-19 800
莫队算法 带修莫队 树上莫队 回滚莫队 树上带修莫队 知识点总结
莫队算法初
一个蒟蒻的博客
05-15 263
emmmm，上周的比赛遇到莫队了，然而只会线段树。 特地学习一下。 参考博文： 洛谷日报 普通莫队(参考洛谷日报) 莫队算法一般分为两类，一是莫队维护区间答案，二是维护区间内的数据结构。 先看一个比较简单的问题: 给一个序列，m次询问，每次询问你区间[l,r][l,r][l,r]有多少种不同的颜色。 $ n,m \leq 100000 $ 最简单的暴力。最极致的超时? 对于每一次询问query...
莫队算法
ke_yi_的博客
07-21 937
莫队算法是一种解决区间问题的算法，是一种很优的暴力算法。 莫队算法的基本思想是，知道一个区间[l,r]的一些信息，就可以在很短的时间内求出它的扩展区间的一些信息。 算法讲解挺详细的链接点击跳转到百度文库pdf文档 题目描述 Given a sequence of integers a1, a2, ..., an and q pairs of integers (l1, r1), (l2, ...
莫队算法良心讲解
weixin_30483013的博客
09-24 155
问题：有n个数组成一个序列，有m个形如询问L, R的询问，每次询问需要回答区间内至少出现2次的数有哪些。 　　朴素的解法需要读取O(nm)次数。如果数据范围小，可以用数组，时间复杂度为O(nm)。如果使用STL的Map来保存出现的次数，则需要O(nmlogn)的复杂度。有没有更快的方法呢？ 　　注意到询问并没有强制在线，因此我们可以使用离线方法。注意到一点，如果我们有计算完[L, R]时的“中...
【莫队算法】
weixin_34358365的博客
06-02 111
·排序巧妙优化复杂度，带来NOIP前的最后一丝宁静。几个活蹦乱跳的指针的跳跃次数，决定着莫队算法的优劣…… ·目前的题型概括为三种：普通莫队，树形莫队以及带修莫队。 若谈及入门，那么BZOJ2038的美妙袜子一题堪称顶尖。 【例题一】袜子 ·述大意： 进行区间询问[l,r]，输出该区间内随机抽两次抽到相同颜色袜子的概率。 ·分析： 首先考虑对于一个长度为n区间内...
莫队入门
weixin_30457551的博客
12-28 80
模板代码HH的项链 学习笔记 时间复杂度 时间复杂度2 带修改莫队复杂度 1 #include <iostream> 2 #include <cstdio> 3 #include <cmath> 4 #include <algorithm> 5 #define reg register 6 using nam...
莫队分块算法模板[BZOJ2038]
snowy_smile的博客
10-28 1657
#include<stdio.h> #include<string.h> #include<ctype.h> #include<math.h> #include<iostream> #include<string> #include<set> #include<map> #include<vector> #include<queue> #include<bitset> #include<algori
莫队算法及应用
剑锋OI SharpSwordOI
11-08 532
该文章非原创，原创地址为：https://www.cnblogs.com/Paul-Guderian/p/6933799.html 莫队算法 ·目前的题型概括为三种：普通莫队，树形莫队以及带修莫队。 若谈及入门，那么BZOJ2038的美妙袜子一题堪称顶尖。 【例题一】袜子 传输门 ·述大意： 进行区间询问[l,r]，输出该区间内随机抽两次抽到相同颜色袜子的概率...
小Z的袜子【莫队算法】
mfcheer
01-28 647
[2009国家集训队]小Z的袜子(hose)Time Limit: 20 Sec Memory Limit: 259 MB Submit: 5259 Solved: 2426 [Submit][Status][Discuss] Description作为一个生活散漫的人，小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿。终于有一天，小Z再也无法忍受这恼人的找袜子过程，于是他决定听
莫队 入门（离线基础板子）
HappyCode
02-03 164
解决类型题：给定一个数列，进行区间的查询. 对于给定的数列，如果我们已经知道了一段区间要求的结果，设左右边界为L和R，然后让我们求qL到qR区间的结果ans，此时，我们可以让R和L移动，到qL和qR上，在移动的过程中，对ans进行修改即可，有点像尺取法对吧。我们要优化时间，就要尽量减少移动的次数，所以此时对整个数列进行分块，进行一个sort排序，如果左端点在一个块内，就以右端点排序，否则以左端点...
理解循环：从‘直到型’到‘当型’——编程基础与算法入门
主讲人莫影在潭州教育分享了结构化程序设计、面向对象程序设计以及流程图等编程启蒙思维。 在结构化程序设计中，有三种基本结构：顺序、选择和循环。循环结构包括当型循环和直到型循环。直到型循环的特点是先执行...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值