题意
您需要写一种数据结构(可参考题目标题),来维护一个有序数列,其中需要提供以下操作:
1.查询k在区间内的排名
2.查询区间内排名为k的值
3.修改某一位值上的数值
4.查询k在区间内的前驱(前驱定义为小于x,且最大的数)
5.查询k在区间内的后继(后继定义为大于x,且最小的数)
题解
求区间第k大,树套树模板题。
外层线段树定序列,内层平衡树询问信息。
询问给定数的区间排名直接做就行了。复杂度O(log22n)
求区间排名为k的值还需要加一层二分答案,用上一问题验证。复杂度O(log32n)。
不知道有没有复杂度更低的方法……
还要注意一下常数问题,其他的不必说了。
#include<cstdio>
#include<cstdlib>
const int maxn=50005;
int max(int x,int y){ return x>y?x:y; }
int min(int x,int y){ return x<y?x:y; }
struct node{
node* ch[2];
int key,fix,size,cnt;
node(){}
node(int w,node* son){ key=w; ch[0]=ch[1]=son; fix=rand(); size=cnt=1; }
void maintain(){ size=ch[0]->size + ch[1]->size + cnt; }
};
typedef node* P_node;
struct Treap{
node nil;
P_node root,null,len;
Treap(){
root=null=&nil;
null->size=null->cnt=0;
null->ch[0]=null->ch[1]=null;
}
void rot(P_node &p,int d){
P_node k=p->ch[d^1]; p->ch[d^1]=k->ch[d]; k->ch[d]=p;
p->maintain(); k->maintain(); p=k;
}
void _Insert(P_node &p,int tkey){
if(p==null) p=new node(tkey,null); else
if(p->key==tkey) p->cnt++; else{
int d=tkey>p->key;
_Insert(p->ch[d],tkey); if(p->ch[d]->fix > p->fix) rot(p,d^1);
}
p->maintain();
}
void _Erase(P_node &p,int tkey){
if(p==null) return;
if(p->key==tkey){
if(p->cnt>1) p->cnt--; else
if(p->ch[0]==null||p->ch[1]==null){
P_node tem=p; p->ch[0]==null?p=p->ch[1]:p=p->ch[0];
delete tem;
}else{
int d=p->ch[0]->fix > p->ch[1]->fix;
rot(p,d); _Erase(p->ch[d],tkey);
}
} else _Erase(p->ch[tkey>p->key],tkey);
p->maintain();
}
int _Rank(P_node p,int tkey,int res){
if(p==null) return res+1;
if(p->key==tkey) return p->ch[0]->size+res+1;
if(tkey<p->key) return _Rank(p->ch[0],tkey,res);
return _Rank(p->ch[1],tkey,res+ p->ch[0]->size + p->cnt);
}
int _Pred(P_node p,int tkey){
if(p==null) return -2147483647;
if(tkey<=p->key) return _Pred(p->ch[0],tkey);
return max(p->key,_Pred(p->ch[1],tkey));
}
int _Succ(P_node p,int tkey){
if(p==null) return 2147483647;
if(tkey>=p->key) return _Succ(p->ch[1],tkey);
return min(p->key,_Succ(p->ch[0],tkey));
}
void Insert(int tkey){ _Insert(root,tkey); }
void Erase(int tkey){ _Erase(root,tkey); }
int Rank(int tkey){ return _Rank(root,tkey,0); }
int Pred(int tkey){ return _Pred(root,tkey); }
int Succ(int tkey){ return _Succ(root,tkey); }
};
struct Segment{
int L,R;
Segment(){ L=R=-1; }
Treap T;
} seg[maxn*4];
void Seg_Build(int p,int L,int R){
seg[p].L=L; seg[p].R=R;
if(L==R) return;
int mid=(L+R)>>1;
Seg_Build(p<<1,L,mid); Seg_Build(p<<1|1,mid+1,R);
}
void Seg_Insert(int p,int pos,int tkey){
if(p>=4*maxn||seg[p].R<pos||pos<seg[p].L) return;
if(seg[p].L<=pos&&pos<=seg[p].R) seg[p].T.Insert(tkey);
Seg_Insert(p<<1,pos,tkey); Seg_Insert(p<<1|1,pos,tkey);
}
void Seg_Erase(int p,int pos,int tkey){
if(p>=4*maxn||seg[p].R<pos||pos<seg[p].L) return;
if(seg[p].L<=pos&&pos<=seg[p].R) seg[p].T.Erase(tkey);
Seg_Erase(p<<1,pos,tkey); Seg_Erase(p<<1|1,pos,tkey);
}
int Seg_Rank(int p,int L,int R,int tkey){
if(p>=4*maxn||seg[p].R<L||R<seg[p].L) return 0;
if(L<=seg[p].L&&seg[p].R<=R) return seg[p].T.Rank(tkey)-1;
return Seg_Rank(p<<1,L,R,tkey)+Seg_Rank(p<<1|1,L,R,tkey);
}
int Seg_Kth(int q_L,int q_R,int k){
int L=0,R=1e+8,res;
while(L<=R){
int mid=(L+R)>>1;
int rk=Seg_Rank(1,q_L,q_R,mid)+1;
if(rk<=k) L=mid+1, res=mid;
else R=mid-1;
}
return res;
}
int Seg_Pred(int p,int L,int R,int tkey){
if(p>=4*maxn||seg[p].R<L||R<seg[p].L) return -2147483647;
if(L<=seg[p].L&&seg[p].R<=R) return seg[p].T.Pred(tkey);
return max(Seg_Pred(p<<1,L,R,tkey),Seg_Pred(p<<1|1,L,R,tkey));
}
int Seg_Succ(int p,int L,int R,int tkey){
if(p>=4*maxn||seg[p].R<L||R<seg[p].L) return 2147483647;
if(L<=seg[p].L&&seg[p].R<=R) return seg[p].T.Succ(tkey);
return min(Seg_Succ(p<<1,L,R,tkey),Seg_Succ(p<<1|1,L,R,tkey));
}
int getint(){
char ch=getchar(); int res=0,ff=1;
while(!('0'<=ch&&ch<='9')){ if(ch=='-') ff=-1; ch=getchar(); }
while('0'<=ch&&ch<='9') res=res*10+ch-'0', ch=getchar();
return res*ff;
}
int n,m,a[maxn];
int main(){
scanf("%d%d",&n,&m);
Seg_Build(1,1,n);
for(int i=1;i<=n;i++) Seg_Insert(1,i,a[i]=getint());
while(m--){
int pd=getint();
if(pd==1){
int x=getint(),y=getint(),z=getint();
printf("%d\n",Seg_Rank(1,x,y,z)+1);
} else
if(pd==2){
int x=getint(),y=getint(),z=getint();
printf("%d\n",Seg_Kth(x,y,z));
} else
if(pd==3){
int x=getint(),y=getint();
Seg_Erase(1,x,a[x]); Seg_Insert(1,x,y);
a[x]=y;
} else
if(pd==4){
int x=getint(),y=getint(),z=getint();
printf("%d\n",Seg_Pred(1,x,y,z));
} else
if(pd==5){
int x=getint(),y=getint(),z=getint();
printf("%d\n",Seg_Succ(1,x,y,z));
}
}
return 0;
}