非旋treap套线段树

BZOJ3065。

去年用pascal 块链过了。。

今年来试了试非旋treap大法   注定被块链完爆

代码留这。

第一份 ：辣鸡的  垃圾回收做法  跑得极慢

 #include <bits/stdc++.h>
 #define M 70000
 using namespace std;
 struct W{int l,r,s,v,e;}a[];
 struct O{int l,r,s,k;}b[];
 int nex[],A[],n,m,rt,ans,l,r,x,k,t,c,T,d[]; char op[];
 void CL(int u){
     if (!u) return; --b[u].k;
     if (!b[u].k) nex[u]=t,t=u;
 }
 void NEW(int &u){
     int U=t; t=nex[t]; CL(u);
     CL(b[U].l); CL(b[U].r);
     b[U]=(O){b[u].l,b[u].r,b[u].s,};
     u=U; ++b[b[u].l].k; ++b[b[u].r].k;
 }
 void add(int &u,int x,int p,int q,int X){
     NEW(u); b[u].s+=X;
     if (p<q)
     if (x<=p+q>>) add(b[u].l,x,p,p+q>>,X);
         else add(b[u].r,x,(p+q>>)+,q,X);
     if (!b[u].s) {nex[u]=t; t=u; u=;}
 }
 void he(int &u,int v,int p,int q){
     if (!v) return;
     if (!u) {++b[u=v].k; return;}
     NEW(u); b[u].s+=b[v].s;
     he(b[u].l,b[v].l,p,p+q>>);
     he(b[u].r,b[v].r,(p+q>>)+,q);
 }
 void up(int u){
     a[u].s=a[a[u].l].s+a[a[u].r].s+;
     CL(a[u].e); a[u].e=;
     add(a[u].e,a[u].v,,M,);
     he(a[u].e,a[a[u].l].e,,M);
     he(a[u].e,a[a[u].r].e,,M);
 }
 void build(int &u,int p,int q){
     if (p>q) return;
     u=p+q>>; a[u].v=A[u];
     build(a[u].l,p,u-);
     build(a[u].r,u+,q);
     up(u);
 }
 void split(int u,int k,int &l,int &r){
     if (!k) {l=; r=u; return;}
     if (k==a[u].s) {l=u; r=; return;}
     if (k<=a[a[u].l].s) split(a[u].l,k,l,a[u].l),up(r=u);
      else split(a[u].r,k-a[a[u].l].s-,a[u].r,r),up(l=u);
 }
 void join(int &u,int k){
     if ((1ll+a[u].s)*rand()<=RAND_MAX){
         int A,B; split(u,k,A,B);
         a[T].l=A; a[T].r=B; up(u=T);
     }else{
         if (k<=a[a[u].l].s) join(a[u].l,k);
         else join(a[u].r,k-a[a[u].l].s-);
         add(a[u].e,a[T].v,,M,); ++a[u].s;
     }
 }
 void mody(int u,int k){
     if (a[a[u].l].s+==k)
         r=a[u].v,a[u].v=l; else
     if (k<=a[a[u].l].s) mody(a[u].l,k);
         else mody(a[u].r,k-a[a[u].l].s-);
     add(a[u].e,r,,M,-);
     add(a[u].e,l,,M,);
 }
 void find(int u,int l,int r){
     if (r<||l>a[u].s) return;
     if (l<=&&a[u].s<=r) {d[++c]=a[u].e; return;}
     if (l<=a[a[u].l].s+&&a[a[u].l].s<r)
         add(d[],a[u].v,,M,);
     find(a[u].l,l,r);
     find(a[u].r,l-a[a[u].l].s-,r-a[a[u].l].s-);
 }
 int main(){
     srand();
     scanf("%d",&n);
     for (int i=;i<=n;++i) scanf("%d",&A[i]);
     for (int i=;i<;++i) nex[i]=i+;
     t=; build(rt,,n); T=n;
     scanf("%d",&m);
     while (m--){
         scanf("%s",op);
         if (op[]=='Q'){
             scanf("%d%d%d",&l,&r,&k);
             l^=ans; r^=ans; k^=ans;
             c=; find(rt,l,r);
             l=; r=M; int tt=d[];
             while (l<r){
                 x=;
                 for (int i=;i<=c;++i) x+=b[b[d[i]].l].s;
                 if (x<k){
                     l=(l+r>>)+; k-=x;
                     for (int i=;i<=c;++i) d[i]=b[d[i]].r;
                 }else{
                     r=l+r>>;
                     for (int i=;i<=c;++i) d[i]=b[d[i]].l;
                 }
             }
             CL(tt); d[]=;
             printf("%d\n",ans=l);
         }else
         if (op[]=='M'){
             scanf("%d%d",&k,&l);
             l^=ans; k^=ans;
             mody(rt,k);
         }else{
             scanf("%d%d",&l,&k);
             l^=ans; k^=ans; --l;
             ++T; a[T].v=k;
             join(rt,l);
         }
     }
     return ;
 }

奈奈

第二份：暴力回收，暴力合并，快了许多

 #include <bits/stdc++.h>
 #define M 70000
 #define S 15000000
 using namespace std;
 int D[S],rt,t,d[],m,x,l,r,k,n,ans,Q; char OP;
 struct O{int l,r,s,v,e;}a[];
 struct U{int l,r,s;}b[S];
 void cl(int &u){
     if (!u) return;
     cl(b[u].l); cl(b[u].r);
     D[++t]=u; b[u].s=; u=;
 }
 void add(int &u,int p,int q,int x){
     if (!u) u=D[t--]; ++b[u].s;
     if (p!=q)
         if (x<=p+q>>) add(b[u].l,p,p+q>>,x);
         else add(b[u].r,(p+q>>)+,q,x);
 }
 void del(int &u,int p,int q,int x){
     --b[u].s;
     if (p!=q)
         if (x<=p+q>>) del(b[u].l,p,p+q>>,x);
         else del(b[u].r,(p+q>>)+,q,x);
     if (!b[u].s) D[++t]=u,u=;
 }
 int he(int l,int r){
     if (!l&&!r) return ;
     int u=D[t--];
     b[u].s=b[l].s+b[r].s;
     b[u].l=he(b[l].l,b[r].l);
     b[u].r=he(b[l].r,b[r].r);
     return u;
 }
 void up(int u){
     a[u].s=a[a[u].l].s+a[a[u].r].s+;
     cl(a[u].e);
     a[u].e=he(a[a[u].l].e,a[a[u].r].e);
     add(a[u].e,,M,a[u].v);
 }
 void build(int &u,int p,int q){
     if (p>q) return;
     u=p+q>>;
     build(a[u].l,p,u-);
     build(a[u].r,u+,q);
     up(u);
 }
 void find(int u,int l,int r){
     if (r<||l>a[u].s) return;
     if (l<=&&a[u].s<=r) {d[++m]=a[u].e; return;}
     int i=a[a[u].l].s+;
     if (l<=i&&i<=r)
         add(d[],,M,a[u].v);
     find(a[u].l,l,r);
     find(a[u].r,l-i,r-i);
 }
 void mody(int u){
     int i=a[a[u].l].s+;
     if (i==k) l=a[u].v,a[u].v=r; else
     if (k<i) mody(a[u].l);
     else k-=i,mody(a[u].r);
     del(a[u].e,,M,l);
     add(a[u].e,,M,r);
 }
 void split(int u,int k,int &l,int &r){
     if (!k) {l=,r=u; return;}
     if (k==a[u].s) {l=u,r=; return;}
     if (k<=a[a[u].l].s){
         split(a[u].l,k,l,a[u].l); up(r=u);
     }else{
         split(a[u].r,k-a[a[u].l].s-,a[u].r,r); up(l=u);
     }
 }
 void join(int &u){
     if ((1ll+a[u].s)*rand()<=RAND_MAX){
         int A,B; split(u,k,A,B);
         a[n].l=A; a[n].r=B; up(u=n);
     }else{
         if (k<=a[a[u].l].s) join(a[u].l);
         else k-=a[a[u].l].s+,join(a[u].r);
         add(a[u].e,,M,a[n].v); ++a[u].s;
     }
 }
 int main(){
     srand();
     scanf("%d",&n);
     for (int i=;i<=n;++i) scanf("%d",&a[i].v);
     for (int i=S-;i;--i) D[i]=++t;
     build(rt,,n);
     scanf("%d",&Q);
     while (Q--){
         scanf("\n%c",&OP);
         if (OP=='Q'){
             scanf("%d%d%d",&l,&r,&k);
             l^=ans; r^=ans; k^=ans;
             m=; find(rt,l,r);
             l=; r=M; int tt=d[];
             while (l<r){
                 x=;
                 for (int i=;i<=m;++i) x+=b[b[d[i]].l].s;
                 if (x<k){
                     l=(l+r>>)+; k-=x;
                     for (int i=;i<=m;++i) d[i]=b[d[i]].r;
                 }else{
                     r=l+r>>;
                     for (int i=;i<=m;++i) d[i]=b[d[i]].l;
                 }
             }
             cl(d[]=tt); printf("%d\n",ans=l);
         }else
         if (OP=='M'){
             scanf("%d%d",&k,&r);
             k^=ans; r^=ans;
             mody(rt);
         }else{
             scanf("%d%d",&k,&x);
             k^=ans; a[++n].v=x^ans;
             --k; join(rt);
         }
     }
     return ;
 }

雲珠桜