jzoj5924 Queue 分块+链表

本文介绍了一个使用Splay树解决序列轮转及区间查询的问题,通过分块和链表维护,实现高效的区间轮转和特定元素计数查询。针对给定长度的数列和操作次数,提出了优化的解决方案。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

Description


Hack 国的居民人人都是 OI 大师,Hometown 得知便赶紧来到 Hack 国学习。可想要进入 Hack 国并不是件容易的事情,首先就必须通过 Hack 国海关小 B 的考验。小 B 觉得 Hometown 比较菜,于是就扔了一道小水题给 Hometown。
给定一个长度为 n 的数列 a i ,接下来会对这个序列进行 m 次操作。操作类型分为以下两种:

  1. l r,表示将区间 [l,r] 轮转一次,具体来说,a l ,a l+1 ,a l+2 ,··· ,a r−1 ,a r 经过一次轮转之后,会变为 a r ,a l ,a l+1 ,··· ,a r−1 ;
  2. l r k,询问区间 [l,r] 内 ai = k 的个数。

可惜 Hometown 还是不会做,他只能期待你能解决这个问题了。
n,m<=1e5

Solution


这很noip

注意到这个轮换操作非常的splay,我们可以用splay来维护这个操作
一个暴力的做法就是分块。我们预处理每个块内的桶,每次轮换操作等价于从每个块中取出最后一个位置,加入下一个块中作为下一块的第一个位置。这个在splay上面求第k个即可
分析一波复杂度,无论是修改还是查询都是nnlog⁡nn\sqrt n\log nnnlogn的,实测跑不过暴力QUQ,说多了都是泪

冷静分析一波可以发现我们这个splay没有什么用,注意到我们只需要资瓷求块内的开始元素和结束元素,因此我们只需要维护n\sqrt nn个链表即可,对于散块直接暴力把链表拉出来再重构
开了O2,各种stl直接上。

Code


#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <deque>
#include <math.h>
#define rep(i,st,ed) for (register int i=st;i<=ed;++i)
#define drp(i,st,ed) for (register int i=st;i>=ed;--i)

const int N=100005;

int bel[N],st[N],ed[N];
int b[350][N],c[N];

struct Block {
	std:: deque <int> q;

	void push_front(int x) {
		q.push_front(x);
	}

	void push_back(int x) {
		q.push_back(x);
	}

	int pop_front() {
		int x=q.front();
		q.pop_front();
		return x;
	}

	int pop_back() {
		int x=q.back();
		q.pop_back();
		return x;
	}

	void move(int l,int r) {
		for (c[0]=0;q.size();) c[++c[0]]=q.front(),q.pop_front();
		int tmp=c[r];
		drp(i,r,l+1) c[i]=c[i-1];
		c[l]=tmp;
		rep(i,1,c[0]) q.push_back(c[i]);
	}
} q[350];

inline int read() {
	int x=0,v=1; char ch=getchar();
	for (;ch<'0'||ch>'9';v=(ch=='-')?(-1):(v),ch=getchar());
	for (;ch<='9'&&ch>='0';x=x*10+ch-'0',ch=getchar());
	return x*v;
}

int main(void) {
	freopen("queue.in","r",stdin);
	freopen("queue.out","w",stdout);
	int n=read(),m=read(),B=sqrt(n)+1;
	rep(i,1,n) {
		int x=read();
		bel[i]=(i-1)/B+1;
		q[bel[i]].push_back(x);
		if (!st[bel[i]]) st[bel[i]]=i;
		ed[bel[i]]=std:: max(ed[bel[i]],i);
		++b[bel[i]][x];
	}
	for (;m--;) {
		int opt=read(),l=read(),r=read();
		if (opt==1) {
			if (bel[l]!=bel[r]) {
				int tmp=q[bel[l]].pop_back();
				b[bel[l]][tmp]--;
				rep(i,bel[l]+1,bel[r]-1) {
					b[i][tmp]++;
					q[i].push_front(tmp);
					tmp=q[i].pop_back();
					b[i][tmp]--;
				}
				b[bel[r]][tmp]++;
				c[0]=0;
				rep(i,st[bel[r]],r) c[++c[0]]=q[bel[r]].pop_front();
				b[bel[r]][c[c[0]]]--;
				drp(i,c[0]-1,1) q[bel[r]].push_front(c[i]);
				q[bel[r]].push_front(tmp);
				tmp=c[c[0]]; c[0]=0;
				rep(i,st[bel[l]],l-1) c[++c[0]]=q[bel[l]].pop_front();
				b[bel[l]][tmp]++;
				q[bel[l]].push_front(tmp);
				drp(i,c[0],1) q[bel[l]].push_front(c[i]);
			} else q[bel[l]].move(l-(bel[l]-1)*B,r-(bel[l]-1)*B);
		} else { int x=read(),ans=0,p;
			if (bel[l]==bel[r]) {
				c[0]=0;
				rep(i,st[bel[l]],ed[bel[l]]) c[++c[0]]=q[bel[l]].pop_front();
				rep(i,l-(bel[l]-1)*B,r-(bel[l]-1)*B) (c[i]==x)?(++ans):0;
				drp(i,c[0],1) q[bel[l]].push_front(c[i]);
			} else {
				c[0]=0;
				rep(i,l,ed[bel[l]]) {
					c[++c[0]]=q[bel[l]].pop_back();
					(c[c[0]]==x)?(++ans):0;
				}
				drp(i,c[0],1) q[bel[l]].push_back(c[i]);
				c[0]=0;
				rep(i,st[bel[r]],r) {
					c[++c[0]]=q[bel[r]].pop_front();
					(c[c[0]]==x)?(++ans):0;
				}
				drp(i,c[0],1) q[bel[r]].push_front(c[i]);
				rep(i,bel[l]+1,bel[r]-1) ans+=b[i][x];
			}
			printf("%d\n", ans);
		}
	}
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值