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

3224: Tyvj 1728 普通平衡树

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 9046  Solved: 3840
[Submit][Status][Discuss]

Description

您需要写一种数据结构（可参考题目标题），来维护一些数，其中需要提供以下操作：
1. 插入x数
2. 删除x数(若有多个相同的数，因只删除一个)
3. 查询x数的排名(若有多个相同的数，因输出最小的排名)
4. 查询排名为x的数
5. 求x的前驱(前驱定义为小于x，且最大的数)
6. 求x的后继(后继定义为大于x，且最小的数)

Input

第一行为n，表示操作的个数,下面n行每行有两个数opt和x，opt表示操作的序号(1<=opt<=6)

Output

对于操作3,4,5,6每行输出一个数，表示对应答案

Sample Input

10
1 106465
4 1
1 317721
1 460929
1 644985
1 84185
1 89851
6 81968
1 492737
5 493598

Sample Output

106465
84185
492737

HINT

1.n的数据范围：n<=100000

2.每个数的数据范围：[-1e7,1e7]

数据如下http://pan.baidu.com/s/1jHMJwO2

Source

平衡树

于是学习了一下数组版的

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;
}