平衡树之treap

treap是平衡树中的一种，为什么要介绍他是因为treap的代码量少，逻辑清晰简单，容易调试，并且性能也很不错.考场里相比起红黑树虽然稳定快速但庞大难调的代码，treap已经蛮够用了.
treap用的是一个随机的思想，给每个节点加入一个键值，这个键值是随机给的，可以rand但很慢，最好自己手写rand（代码如下）.在treap左旋右旋的操作里，就是通过保证子树的键值一定大于父亲来进行旋转操作，使整棵树保持平衡（旋转操作可以百度一下）.代码如下（bzoj3224）；

#include<stdio.h>int n,size,root,ans;inline const int read(){   register int x=0,f=1;   register char ch=getchar();   while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}   while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}   return f*x;}inline int randad(){//手写rand   static int seed=447;//seed可以随便取   return seed=int(seed*48271LL%2147483647);}struct data{   int l,r,v,rnd,size,w;}tr[100005];void update(int k){//更新操作   tr[k].size=tr[tr[k].l].size+tr[tr[k].r].size+tr[k].w;}void rturn(int &k){//右旋   int t=tr[k].l;tr[k].l=tr[t].r;tr[t].r=k;   tr[t].size=tr[k].size;update(k);k=t;}void lturn(int &k){//左旋   int t=tr[k].r;tr[k].r=tr[t].l;tr[t].l=k;   tr[t].size=tr[k].size;update(k);k=t;}void insert(int &k,int x){//插入操作   if(k==0){      size++;k=size;      tr[k].size=tr[k].w=1;tr[k].v=x;tr[k].rnd=randad();return;   }   tr[k].size++;   if(tr[k].v==x) tr[k].w++;   else if(x>tr[k].v) {insert(tr[k].r,x); if(tr[tr[k].r].rnd<tr[k].rnd) lturn(k);}   else {insert(tr[k].l,x); if(tr[tr[k].l].rnd<tr[k].rnd) rturn(k);}}//最后几行是使子树键值大于父亲键值void del(int &k,int x){   if(k==0) return;   if(tr[k].v==x){//w是这个数的个数      if(tr[k].w>1) {tr[k].w--;tr[k].size--;return;}      if(tr[k].l*tr[k].r==0) k=tr[k].l+tr[k].r;      else if(tr[tr[k].l].rnd<tr[tr[k].r].rnd) rturn(k),del(k,x);      else lturn(k),del(k,x);   }     else if(x>tr[k].v) tr[k].size--,del(tr[k].r,x);   else tr[k].size--,del(tr[k].l,x); }int query_rank(int k,int x){//求排名   if(k==0) return 0;   if(tr[k].v==x) return tr[tr[k].l].size+1;   else if(x>tr[k].v) return tr[tr[k].l].size+tr[k].w+query_rank(tr[k].r,x);    else return query_rank(tr[k].l,x);}int query_num(int k,int x){//求排名第几是谁   if(k==0) return 0;   if(x<=tr[tr[k].l].size) return query_num(tr[k].l,x);   else if(x>tr[tr[k].l].size+tr[k].w) return query_num(tr[k].r,x-tr[tr[k].l].size-tr[k].w);   else return tr[k].v;}void query_pro(int k,int x){   if(k==0) return;   if(tr[k].v<x) ans=k,query_pro(tr[k].r,x);   else query_pro(tr[k].l,x);}void query_sub(int k,int x){   if(k==0) return;   if(tr[k].v>x) ans=k,query_sub(tr[k].l,x);   else query_sub(tr[k].r,x);}//以上是求lower——bound和upper——boundint main(){   int opt,x;   n=read();   for(register int i=1;i<=n;i++){      opt=read(),x=read();      switch(opt){         case 1:insert(root,x);break;         case 2:del(root,x);break;         case 3:printf("%d\n",query_rank(root,x));break;         case 4:printf("%d\n",query_num(root,x));break;         case 5:ans=0,query_pro(root,x);printf("%d\n",tr[ans].v);break;         case 6:ans=0,query_sub(root,x);printf("%d\n",tr[ans].v);break;      }   } }