BZOJ3224: Tyvj 1728 普通平衡树[treap]

平衡树

sdsc时gty讲的指针版真把我弄烦了，反人类啊根本打不快并且->一大片更容易混乱

于是学习了一下数组版的

1.用l和r，操作需要分类讨论左还是右

2.旋转传引用，操作后x还是根结点

3.insert和del传引用因为要修改树的样子

左旋是根向左转，右边的上来

右旋相反

注意：前驱和后继时=v不能用这个更新ans

//
//  main.cpp
//  treap数组
//
//  Created by Candy on 26/11/2016.
//  Copyright © 2016 Candy. All rights reserved.
//

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int N=1e5+5;
int read(){
    char c=getchar();int x=0,f=1;
    while(c<‘0‘||c>‘9‘){if(c==‘-‘)f=-1; c=getchar();}
    while(c>=‘0‘&&c<=‘9‘){x=x*10+c-‘0‘; c=getchar();}
    return x*f;
}
struct node{
    int l,r,size,v,w,rnd;
}t[N];
int cnt,root;
inline void update(int x){t[x].size=t[t[x].l].size+t[t[x].r].size+t[x].w;}
inline void rturn(int &x){
    int c=t[x].l;t[x].l=t[c].r;t[c].r=x;
    t[c].size=t[x].size;update(x);x=c;
}
inline void lturn(int &x){
    int c=t[x].r;t[x].r=t[c].l;t[c].l=x;
    t[c].size=t[x].size;update(x);x=c;
}
void insert(int &x,int v){
    if(x==0){
        cnt++;x=cnt;
        t[cnt].l=t[cnt].r=0;
        t[cnt].v=v;t[cnt].w=t[cnt].size=1;t[cnt].rnd=rand();
        return;
    }
    t[x].size++;
    if(t[x].v==v) t[x].w++;
    else if(v<t[x].v){
        insert(t[x].l,v);
        if(t[t[x].l].rnd<t[x].rnd) rturn(x);
    }else{
        insert(t[x].r,v);
        if(t[t[x].r].rnd<t[x].rnd) lturn(x);
    }
}
void del(int &x,int v){
    if(x==0) return;
    if(t[x].v==v){
        if(t[x].w>1){t[x].w--;t[x].size--;return;}
        if(t[x].l*t[x].r==0) x=t[x].l+t[x].r;//一个孩子为空，直接用另一个代替
        else if(t[t[x].l].rnd<t[t[x].r].rnd)
            rturn(x),del(x,v);//del(t[x].r,v);需要额外处理size，不如等下一步del的else
        else lturn(x),del(x,v);//del(t[x].l,v);
    }else {
        t[x].size--;
        if(v<t[x].v) del(t[x].l,v);
        else del(t[x].r,v);
    }
}
int rnk(int x,int v){
    if(x==0) return 0;
    if(t[x].v==v) return t[t[x].l].size+1;
    else if(v<t[x].v) return rnk(t[x].l,v);
    else return t[t[x].l].size+t[x].w+rnk(t[x].r,v);
}
int kth(int x,int k){
    if(x==0) return 0;
    if(k<=t[t[x].l].size) return kth(t[x].l,k);
    else if(k>t[t[x].l].size+t[x].w)
        return kth(t[x].r,k-t[t[x].l].size-t[x].w);
    else return t[x].v;
}
int ans;
void pre(int x,int v){
    if(x==0) return;
    if(v>t[x].v) ans=x,pre(t[x].r,v);
    else pre(t[x].l,v);
}
void suf(int x,int v){
    if(x==0) return;
    if(v<t[x].v) ans=x,suf(t[x].l,v);
    else suf(t[x].r,v);
}
int n,op,x;
int main(int argc, const char * argv[]){
    srand(222);
    n=read();
    while(n--){
        op=read();x=read();
        switch(op){
            case 1:insert(root,x);break;
            case 2:del(root,x);break;
            case 3:printf("%d\n",rnk(root,x));break;
            case 4:printf("%d\n",kth(root,x));break;
            case 5:ans=0;pre(root,x);printf("%d\n",t[ans].v);break;
            case 6:ans=0;suf(root,x);printf("%d\n",t[ans].v);break;
        }
    }
    return 0;
}

以及使用一个诡异的宏来简化代码版

//
//  main.cpp
//  treap_prac1
//
//  Created by Candy on 26/11/2016.
//  Copyright © 2016 Candy. All rights reserved.
//

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
#define lc t[x].l
#define rc t[x].r
const int N=1e5+5;
int read(){
    char c=getchar();int x=0,f=1;
    while(c<‘0‘||c>‘9‘){if(c==‘-‘)f=-1; c=getchar();}
    while(c>=‘0‘&&c<=‘9‘){x=x*10+c-‘0‘; c=getchar();}
    return x*f;
}
struct node{
    int l,r,v,w,size,rnd;
}t[N];
int cnt=0,root;
inline void update(int x){t[x].size=t[lc].size+t[rc].size+t[x].w;}
inline void rturn(int &x){
    int c=lc;lc=t[c].r;t[c].r=x;
    t[c].size=t[x].size;update(x);x=c;
}
inline void lturn(int &x){
    int c=rc;rc=t[c].l;t[c].l=x;
    t[c].size=t[x].size;update(x);x=c;
}
void ins(int &x,int v){
    if(x==0){
        cnt++;
        t[cnt].l=t[cnt].r=0;t[cnt].size=t[cnt].w=1;
        t[cnt].v=v;t[cnt].rnd=rand();
        return;
    }
    t[x].size++;
    if(t[x].v==v) {t[x].w++;return;}
    if(v<t[x].v){
        ins(lc,v);
        if(t[lc].rnd<t[x].rnd) rturn(x);
    }else{
        ins(rc,v);
        if(t[rc].rnd<t[x].rnd) lturn(x);
    }
}
void del(int &x,int v){
    if(x==0) return;
    if(t[x].v==v){
        if(t[x].w>1){t[x].w--,t[x].size--;return;}
        if(lc*rc==0) x=lc+rc;
        else if(t[lc].rnd<t[rc].rnd) rturn(x),del(x,v);
        else lturn(x),del(x,v);
    }else{
        t[x].size--;
        if(v<t[x].v) del(lc,v);
        else del(rc,v);
    }
}
int rnk(int x,int v){
    if(t[x].v==v) return t[lc].size+1;
    if(v<t[x].v) return rnk(lc,v);
    else return t[lc].size+t[x].w+rnk(rc,v);
}
int kth(int x,int k){
    if(x==0) return 0;
    if(k<=t[lc].size) return kth(lc,k);
    else if(k>t[lc].size+t[x].w) return kth(rc,k-t[lc].size-t[x].w);
    else return t[x].v;
}
int ans;
void pre(int x,int v){
    if(x==0) return;
    if(v>t[x].v) ans=x,pre(rc,v);
    else pre(lc,v);
}
void suf(int x,int v){
    if(x==0) return;
    if(v<t[x].v) ans=x,suf(lc,v);
    else suf(rc,v);
}
int main(int argc, const char * argv[]) {
    srand(222);
    int n=read(),op,x;
    while(n--){
        op=read();x=read();
        switch(op){
            case 1:ins(root,x);break;
            case 2:del(root,x);break;
            case 3:printf("%d\n",rnk(root,x));break;
            case 4:printf("%d\n",kth(root,x));break;
            case 5:ans=0;pre(root,x);printf("%d\n",t[ans].v);break;
            case 6:ans=0;suf(root,x);printf("%d\n",t[ans].v);break;
        }
    }
    return 0;
}

BZOJ3224: Tyvj 1728 普通平衡树[treap]

3224: Tyvj 1728 普通平衡树

Description

Input

Output

Sample Input

Sample Output

HINT

Source