洛谷2617 动态可持久线段树

最新推荐文章于 2024-11-14 19:30:39 发布

mozeli

最新推荐文章于 2024-11-14 19:30:39 发布

阅读量81

点赞数

分类专栏：高级课程文章标签： c语言数据结构算法

本文链接：https://blog.csdn.net/mozelyr/article/details/120026627

版权

高级课程专栏收录该内容

13 篇文章 0 订阅

订阅专栏

本文介绍了一种使用C++ STL实现的数据结构，通过`tre`结构存储整数对，并实现了快速的更新和查询功能。博客详细讲解了`build`、`update`、`add`、`sums`和`query`等关键函数，重点在于如何高效地处理大规模数据的增删和查找操作。

摘要由CSDN通过智能技术生成

添加链接描述

#include<bits/stdc++.h>
using namespace std;
const int maxn=2e5+5;
struct st{
	int l,r,sum;
}tre[maxn<<9];
struct op{
	char s[2];
	int x,y,z;
}q[maxn];
int a[maxn],b[maxn];
int tot,ver[maxn];
int n,m,c;
int ndl[maxn],ndr[maxn],cntl,cntr;
int build(int l,int r){
	int root=tot;
	st &rt=tre[root];
	rt.l =root;
	rt.r =root;
	rt.sum =0;
	return root;
}
int update(int node,int l,int r,int p,int v){
	int root=++tot;
	st &rt=tre[root];
	rt=tre[node];
	if(l==r){
		rt.sum +=v;
		return root;
	}
	int mid=(l+r)>>1;
	if(p<=mid) rt.l =update(rt.l ,l,mid,p,v);
	else rt.r =update(rt.r ,mid+1,r,p,v);
	rt.sum =tre[rt.l ].sum +tre[rt.r ].sum ;
	return root;
}
void add(int x,int p,int f){
	x=lower_bound(b+1,b+c+1,x)-b;
	for(;p<=n;p+=p&-p){
		ver[p]=update(ver[p],1,c,x,f);
	}
}
int sums(int node[],int cnt){
	int res=0;
	for(int i=0;i<cnt;i++) res+=tre[tre[node[i]].l ].sum ;
	return res;
}
int query(int l,int r,int k){
	if(l==r) return l;
	int x=sums(ndr,cntr)-sums(ndl,cntl);
	int mid=(l+r)>>1;
	if(x>=k){
		for(int i=0;i<cntl;i++) ndl[i]=tre[ndl[i]].l ;
		for(int i=0;i<cntr;i++) ndr[i]=tre[ndr[i]].l ;
		return query(l,mid,k);
	}
	else{
		for(int i=0;i<cntl;i++) ndl[i]=tre[ndl[i]].r ;
		for(int i=0;i<cntr;i++) ndr[i]=tre[ndr[i]].r ;
		return query(mid+1,r,k-x);
	}
}
int ask(int l,int r,int k){
	cntr=cntl=0;
	for(int i=r;i>0;i-=i&-i) ndr[cntr++]=ver[i];
	for(int i=l;i>0;i-=i&-i) ndl[cntl++]=ver[i];
	return query(1,c,k);
}
int main()
{
	c=0;
	tot=0;
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++){
		scanf("%d",&a[i]);
		b[i]=a[i];
	}
	for(int i=1;i<=m;i++){
		scanf("%s%d%d",q[i].s ,&q[i].x ,&q[i].y );
		if(q[i].s [0]=='C') b[++c+n]=q[i].y ;
		else scanf("%d",&q[i].z );
	}
	sort(b+1,b+1+c+n);
	c=unique(b+1,b+n+c+1)-b-1;
	memset(ver,0,sizeof(ver));
	ver[0]=build(1,c);
	for(int i=1;i<=n;i++) add(a[i],i,1);
	for(int i=1;i<=m;i++){
		if(q[i].s[0] =='Q') printf("%d\n",b[ask(q[i].x -1,q[i].y ,q[i].z )]);
		else if(a[q[i].x ]!=q[i].y ){
			add(a[q[i].x ],q[i].x ,-1);
			add(q[i].y ,q[i].x ,1);
			a[q[i].x]=q[i].y;//
		}
	}
	return 0;
}
}