待修改莫队与普通莫队优化

待修改莫队与普通莫队优化

普通莫队比较简单,先不讲了。

奇偶华排序

观察数据:

复制代码
1 1 
2 100
3 1
4 100

显然这么排会浪费许多时间,构造排列:

复制代码
1 1
3 1
2 100
4 100

这样可以极大地减少指针移动次数。

就是将 l l l 指针的所在块进行分类,如果是奇数块就从小到大排,否则从大到小排。

可以玄学优化 30 % 30\% 30% 的时间(逃

待修改莫队

普通莫队是不可以修改的,但是我们可以强行让其可以单点修改。

我们在普通莫队的基础上维护一维时间维。

我们定义时间是修改的顺序编号,即第几个修改操作。

具体地:

  • 单点修改,把某一位的数字修改掉。对于现在这个询问,如果当前区间的时间比它小,则暴力增加修改操作。否则撤销操作。
  • 增加操作:如果操作的位置在区间内,则删除原数字,增加改动的数字。如果不在,直接修改数组。
  • 撤销操作:即增加的逆操作。

例题

P1903 【模板】带修莫队 / [国家集训队] 数颜色 / 维护队列

代码:

cpp 复制代码
#include<bits/stdc++.h>
using namespace std;
typedef long long ljl;
#define FUP(i,x,y) for(auto i=(x);i<=(y);++i)
#define FDW(i,x,y) for(auto i=(x);i>=(y);--i)
inline void Rd(auto &num);
const int N=133335,K=1e6+5;
int n,m,cnt[K],cntc,cntq,pos[N],curl,curr,curc,ans[N],sum,a[N];
void blk()
{
	int blksz=(int)pow(n,2.0/3.0);
	FUP(i,1,n)
		pos[i]=(i-1)/blksz+1;
	return;
}
struct C{
	int p,c;
}c[N];
struct Q{
	int l,r,uid,tim;
	bool operator < (const Q &x)const{
		if(pos[l]!=pos[x.l])return pos[l]<pos[x.l];
		if(pos[r]!=pos[x.r])
		{
			if(pos[l]&1)
				return pos[r]<pos[x.r];
			return pos[r]>pos[x.r];
		}
		return tim<x.tim;
	}
}q[N];
void add(int x)
{
	if(!cnt[x])++sum;
	++cnt[x];
	return;
}
void del(int x)
{
	if(cnt[x]==1)--sum;
	--cnt[x];
	return;
}
void getch(char &ch)
{
	ch=getchar();
	while(ch!='Q'&&ch!='R')ch=getchar();
	return;
}
void upd(int i,int curc)
{
	if(q[i].l<=c[curc].p&&c[curc].p<=q[i].r)
	{
		del(a[c[curc].p]);
		add(c[curc].c);
	}
	swap(a[c[curc].p],c[curc].c);//小技巧:下次操作一定是反着的
	return;
}
int main(){
	Rd(n);Rd(m);
	blk();
	FUP(i,1,n)Rd(a[i]);
	for(int i=1,l,r;i<=m;++i)
	{
		char ch;getch(ch);
		Rd(l);Rd(r);
		if(ch=='Q')
		{
			q[++cntq].l=l;q[cntq].r=r;
			q[cntq].uid=cntq;q[cntq].tim=cntc;
		}
		else{
			c[++cntc].p=l;c[cntc].c=r;}
	}
	sort(q+1,q+cntq+1);curl=1;curr=0;curc=0;
	FUP(i,1,cntq)
	{
		while(q[i].l<curl)add(a[--curl]);
		while(curr<q[i].r)add(a[++curr]);
		while(curl<q[i].l)del(a[curl++]);
		while(q[i].r<curr)del(a[curr--]);
		while(curc<q[i].tim){
			++curc;upd(i,curc);}
		while(curc>q[i].tim){
			upd(i,curc);--curc;}
		ans[q[i].uid]=sum;
	}
	FUP(i,1,cntq)
		printf("%d\n",ans[i]);
	return 0;
}
inline void Rd(auto &num)
{
	num=0;char ch=getchar();bool f=0;
	while(ch<'0'||ch>'9')
	{
		if(ch=='-')f=1;
		ch=getchar();
	}
	while(ch>='0'&&ch<='9')
	{
		num=(num<<1)+(num<<3)+(ch-'0');
		ch=getchar();
	}
	if(f)num=-num;
	return;
}
相关推荐
renhongxia11 分钟前
如何基于知识图谱进行故障原因、事故原因推理,需要用到哪些算法
人工智能·深度学习·算法·机器学习·自然语言处理·transformer·知识图谱
坚持就完事了2 分钟前
数据结构之树(Java实现)
java·算法
WangYaolove13142 分钟前
基于python的在线水果销售系统(源码+文档)
python·mysql·django·毕业设计·源码
算法备案代理4 分钟前
大模型备案与算法备案,企业该如何选择?
人工智能·算法·大模型·算法备案
Monly215 分钟前
Java:修改打包配置文件
java·开发语言
roman_日积跬步-终至千里5 分钟前
【架构设计与实现】动态数据源切换:核心代码实现手册
java
AALoveTouch9 分钟前
大麦网协议分析
javascript·python
XiaoFan01214 分钟前
免密批量抓取日志并集中输出
java·linux·服务器
●VON18 分钟前
React Native for OpenHarmony:2048 小游戏的开发与跨平台适配实践
javascript·学习·react native·react.js·von
顾北1220 分钟前
MCP服务端开发:图片搜索助力旅游计划
java·spring boot·dubbo