ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

带修莫队例题详解

2022-08-24 01:03:23  阅读:131  来源: 互联网

标签:dui int ir read tm cntq 例题 详解 修莫队


带修莫队

[P1903 国家集训队] 数颜色 / 维护队列

版本更新内容:

  • 在普通莫队基础上增加时间坐标,提高游戏难度;
  • 排序时以时间坐标为第三关键字,奇偶排序玄学值上调 \(20\%\);
  • 代码常数加大,请玩家将分块大小调至 \(n^{\frac{2}{3}}\) 以抵消常数因子;
  • 莫队函数主体内容增加:双指针操作完,请将时间轴调至正确位置,并关注区间内时间扰动;
  • 在时间跳跃时,若目标是未来,请先抵达未来再处理时间扰动,若目标是过去,请先处理时间扰动再回到过去。\(swap\) 是个神奇且玄妙的操作,请学会使用;详见代码中 \(upd()\) 函数的注释;
  • 请分别储存时间节点询问节点,提问者会感谢你;
  • 请勿忘记 \(I/O\) 优化 , 祝您游戏愉快。
#include<bits/stdc++.h>
#define re register
#define rep(i,a,n) for(re int i=a;i<=n;++i)
using namespace std;
inline int read() {
    re int x=0,f=1;re char c=getchar();//getchar()yyds!
    while(c<'0'||c>'9') {if(c=='-') f=-1;c=getchar();}
    while (c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
    return x*f;
}
const int maxn = 1e6;
struct area{int l,r,id,tm;}dui[maxn]; 
int a[maxn],ans,n,m,bel[maxn],ens[maxn];
int cntr,cntq;
int cnt[maxn];
bool cmp(area x,area y) {
	if(bel[x.l]==bel[y.l]) {
		if(bel[x.r]==bel[y.r]) return x.tm < y.tm;
		return x.r < y.r;
	}
	return x.l < y.l;
}
bool cmp1(area x,area y)
{
	if(bel[x.l] != bel[y.l]) return x.l<y.l;
	if(bel[x.r] != bel[y.r]) if(bel[x.l]^1) return x.r < y.r;
	return x.tm > y.tm; 
} 
struct Memory {int pos,val;}mem[maxn];
void upd(int l,int r,int tm)
{
	re int pos = mem[tm].pos ,val = mem[tm].val;
	if(pos >= l && pos <= r)
	{
		ans -= !(--cnt[a[pos]]);
		ans += !(cnt[val]); cnt[val]++;
	}
	swap(a[pos],mem[tm].val); 
	//这里有个很巧妙的操作
    //对于一个操作,下一次需要为的颜色是本次被改变的颜色
    //比如,我把颜色3改为了7,那么再进行这次修改的时候就是把7改为3
    //所以直接交换两种颜色就好 
    // 假如经过两次,则swap操作就奇妙地把原数组还原了 
}
void solv() {
	ans=0;
	re int il=1,ir=0,nt=0;
	rep(i,1,cntq)
	{
		int l=dui[i].l,r=dui[i].r,tm=dui[i].tm;
		while(il < l){--cnt[a[il]];ans -= !(cnt[a[il]]);++il;}
		while(il > l){il--; ans += !cnt[a[il]]; ++cnt[a[il]]; } 
		while(ir < r){++ir; ans += !cnt[a[ir]]; ++cnt[a[ir]];}
		while(ir > r){--cnt[a[ir]]; ans -= !cnt[a[ir]]; --ir;}
		while(nt > tm){upd(il,ir,nt);--nt;} // 先更新再回溯
		while(nt < tm){++nt;upd(il,ir,nt);} // 先回溯再更新,原因是upd里的 swap 
		ens[dui[i].id] = ans;
	}	
}
int main()
{
	n=read(),m=read();
	int blk = (int)pow(n,(double)2*1.0/(double)3); // 2/3次方块长优化
	//int blk = (int)sqrt(n); // 根号已落后于版本
	rep(i,1,n) a[i] = read(),bel[i]=(i-1)/blk;
	char ch;
	int l,r,p,co;
	rep(i,1,m) {
		scanf("%c",&ch);
		if(ch == 'R')
		{
			++cntr; p=read(),co=read();
			mem[cntr].pos = p, mem[cntr].val = co;
		}
		else
		{
			dui[++cntq].l=read(),dui[cntq].r=read();
			dui[cntq].id = cntq;
			dui[cntq].tm = cntr;
		}
	}
    sort(dui+1,dui+cntq+1,cmp);
    solv();
    rep(i,1,cntq) printf("%d\n",ens[i]);
	return 0;
} 

作者:@魔幻世界魔幻人生
本文为作者原创,转载请注明出处:https://www.cnblogs.com/subtlemaple/p/16618356.html

标签:dui,int,ir,read,tm,cntq,例题,详解,修莫队
来源: https://www.cnblogs.com/subtlemaple/p/16618387.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有