可持久化线段树/主席树:
顾名思义,该数据结构是可以访问历史版本的线段树。用于解决需要查询历史信息的区间问题。
在功能与时间复杂度上与开n棵线段树无异,然而空间复杂度从$O(n imes nlogn)$降到了$O(nlogn)$。
实现方法:
每次只更新有关的节点(每层一个,共$logn$个),其余节点不动。
用一个数组$rt[i]$记录第$i$个版本线段树的根节点(显然每次必更新根节点)。
查询时一路走下去,将两个需要查询的历史版本的节点信息差分。
没了。
模板题目:洛谷P3834
#include<algorithm> #include<iostream> #include<cstring> #include<cstdio> using namespace std; #define MAXN 200005 #define MAXM 500005 #define INF 0x7fffffff #define ll long long inline int read(){ int x=0,f=1; char c=getchar(); for(;!isdigit(c);c=getchar()) if(c=='-') f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f; } int N,M,Q,A[MAXN],B[MAXN],tot,rt[MAXN]; int ls[MAXN*20],rs[MAXN*20],sum[MAXN*20]; inline void build(int l,int r,int &k){ k=++tot; if(l==r) return; int mid=l+r>>1; build(l,mid,ls[k]); build(mid+1,r,rs[k]); } inline void update(int l,int r,int p,int las,int &k){ k=++tot; sum[k]=sum[las]+1; ls[k]=ls[las],rs[k]=rs[las]; if(l==r) return; int mid=l+r>>1; if(p<=mid) update(l,mid,p,ls[las],ls[k]); else update(mid+1,r,p,rs[las],rs[k]); } inline int query(int l,int r,int s,int u,int v){ if(l==r) return l; int mid=l+r>>1,x=sum[ls[v]]-sum[ls[u]]; if(s<=x) return query(l,mid,s,ls[u],ls[v]); else return query(mid+1,r,s-x,rs[u],rs[v]); } int main(){ N=read(),Q=read(); for(int i=1;i<=N;i++) A[i]=read(); memcpy(B,A,sizeof(A)); sort(B+1,B+1+N); M=unique(B+1,B+1+N)-B-1; build(1,M,rt[0]); for(int i=1;i<=N;i++){ int p=lower_bound(B+1,B+1+M,A[i])-B; update(1,M,p,rt[i-1],rt[i]); } while(Q--){ int l=read(),r=read(),k=read(); printf("%d ",B[query(1,M,k,rt[l-1],rt[r])]); } return 0; }