非旋treap套线段树

BZOJ3065。

去年用pascal 块链过了。。

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

代码留这。

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

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

 

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

#include <bits/stdc++.h>
#define M 70000
#define S 15000000
using namespace std;
int D[S],rt,t,d[70005],m,x,l,r,k,n,ans,Q; char OP;
struct O{int l,r,s,v,e;}a[70005];
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=0; u=0;
}
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>>1) add(b[u].l,p,p+q>>1,x);
        else add(b[u].r,(p+q>>1)+1,q,x);
}
void del(int &u,int p,int q,int x){
    --b[u].s;
    if (p!=q)
        if (x<=p+q>>1) del(b[u].l,p,p+q>>1,x);
        else del(b[u].r,(p+q>>1)+1,q,x);
    if (!b[u].s) D[++t]=u,u=0;
}
int he(int l,int r){
    if (!l&&!r) return 0;
    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+1;
    cl(a[u].e);
    a[u].e=he(a[a[u].l].e,a[a[u].r].e);
    add(a[u].e,0,M,a[u].v);
}
void build(int &u,int p,int q){
    if (p>q) return;
    u=p+q>>1;
    build(a[u].l,p,u-1);
    build(a[u].r,u+1,q);
    up(u);
}
void find(int u,int l,int r){
    if (r<1||l>a[u].s) return;
    if (l<=1&&a[u].s<=r) {d[++m]=a[u].e; return;}
    int i=a[a[u].l].s+1;
    if (l<=i&&i<=r)
        add(d[1],0,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+1;
    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,0,M,l);
    add(a[u].e,0,M,r);
}
void split(int u,int k,int &l,int &r){
    if (!k) {l=0,r=u; return;}
    if (k==a[u].s) {l=u,r=0; 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-1,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+1,join(a[u].r);
        add(a[u].e,0,M,a[n].v); ++a[u].s;
    }
}
int main(){
    srand(23);
    scanf("%d",&n);
    for (int i=1;i<=n;++i) scanf("%d",&a[i].v);
    for (int i=S-1;i;--i) D[i]=++t;
    build(rt,1,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=1; find(rt,l,r);
            l=0; r=M; int tt=d[1];
            while (l<r){
                x=0;
                for (int i=1;i<=m;++i) x+=b[b[d[i]].l].s;
                if (x<k){
                    l=(l+r>>1)+1; k-=x;
                    for (int i=1;i<=m;++i) d[i]=b[d[i]].r;
                }else{
                    r=l+r>>1;
                    for (int i=1;i<=m;++i) d[i]=b[d[i]].l;
                }
            }
            cl(d[1]=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 0;
}