Fork me on GitHub

普通Splay详解

预备知识:

二叉搜索树(BST)

至于BST,随便看一下就可以,

我们知道二叉搜索树是O(logN)的,那我们为什么要用平衡树呢?

之前我们了解到,BST的插入是小的往左子树走,大的往右子树走,如果凉心出题人给出的序列是有序的呢

1

这样我们就只能O(N)的操作,GG

旋转(rotate):

Splay的经典操作就是旋转

在Splay中,我们用旋转来保持平衡,也就是保持是log(N)的量级

旋转就是将节点向上旋转到父亲节点的位置,同时保持平衡

有zig,和zag两种情况(其实都一个样

2                                 3

具体要怎么旋转呢

4

如图,X,Y,Z 是三个节点,A,B,C 是三颗子树,我们要把 Z 转到 Y 的位置

其实就只有3步

 

根据平衡树的性质

Z 是 Y 的左儿子,所以Z < Y

Y 是 X 的左儿子,所以Y < X

我们要把 Z 旋转上去的话,就把 Z 放到 Y 的位置

整完了长这样

5

因为我们还没操作 Y,所以Y还有连向其父亲和儿子的边

总结一下

Step1:把要旋转的节点放到父亲的位置

 

而 Y > Z 且  Y < X,所以这时Y就成了Z的右儿子

6

总结一下

Step2:把要旋转节点的父亲设为其儿子

 

这时会有三个节点(子树)连向 Z,而Y只有一个儿子,显然,Z,子树 B 和子树 C 都是小于 Y 的,子树 B 大于 Z,所以 B 成为 Y 的左儿子

7

这样就完成了

总结一下
Step3:把 父节点所占旋转节点的儿子 设为父节点的对应儿子

代码:

定义一波:

struct tree {
    int fa, cnt, sum, val;
    //父亲
    //计数(几个值为x的点)
    //以当前点为根节点的子树节点个数
    //当前点的值
    int ch[2];
    //左右儿子,0为左儿子,1为右儿子
} t[N];    

关于获得这个节点是左儿子还是右儿子:

inline int get(int x) {
    return t[t[x].fa].ch[0] == x ? 0 : 1;
}

 

更新:

inline void pushup(int x) {
    t[x].sum = t[t[x].ch[0]].sum + t[t[x].ch[1]].sum + t[x].cnt;
}

 

旋转:

inline void rotate(int x) {
    int fa = t[x].fa, gfa = t[fa].fa;//父亲和爷爷
    int k = get(x);//x是其父节点的那个儿子
    //step1  
    t[x].fa = gfa;
    t[gfa].ch[get(fa)] = x;
    //step2
    t[t[x].ch[k ^ 1]].fa = fa;
    t[fa].ch[k] = t[x].ch[k ^ 1];
    //step3
    t[fa].fa = x;
    t[x].ch[k ^ 1] = fa;
    pushup(fa), pushup(x);
    //因为旋转后父节点成了当前点的子节点,所以先更新父亲
}

关于为什么是 k ^ 1,假设我们要旋转的点是左儿子,那他的父亲一定会成为他的右儿子,同理,如果要旋转的点是左儿子,他的父节点一定会成为他的右儿子

伸展(splay)

 splay操作就是把一个点旋转到指定的点

最容易想到的,就是一直旋转到指定的节点,然而这样是错的

8

这时我们就要用到双旋,双旋有两大种四小种情况

1、zig-zig或zag-zag

当节点是父亲的左儿子且父节点是祖父节点的左儿子

或节点是父亲的右儿子且父节点是祖父节点的右儿子

先旋转父亲,再旋转自己

借用一下GeeksofrGeeks的图:

Zig-Zig (Left Left Case):
       G                        P                           X       
      / \                     /   \                        / \      
     P  T4   rightRotate(G)  X     G     rightRotate(P)  T1   P     
    / \      ============>  / \   / \    ============>       / \    
   X  T3                   T1 T2 T3 T4                      T2  G
  / \                                                          / \ 
 T1 T2                                                        T3  T4 

Zag-Zag (Right Right Case):
  G                          P                           X       
 /  \                      /   \                        / \      
T1   P     leftRotate(G)  G     X     leftRotate(P)    P   T4
    / \    ============> / \   / \    ============>   / \   
   T2   X               T1 T2 T3 T4                  G   T3
       / \                                          / \ 
      T3 T4                                        T1  T2
2.zig-zag或zag-zig

当节点是父亲的左儿子且父节点是祖父节点的右儿子

或节点是父亲的右儿子且父节点是祖父节点的左儿子

旋转两次自己

再次借用GeeksforGeeks的图:

Zag-Zig (Left Right Case):
       G                        G                            X       
      / \                     /   \                        /   \      
     P   T4  leftRotate(P)   X     T4    rightRotate(G)   P     G     
   /  \      ============>  / \          ============>   / \   /  \    
  T1   X                   P  T3                       T1  T2 T3  T4 
      / \                 / \                                       
    T2  T3              T1   T2                                     

Zig-Zag (Right Left Case):
  G                          G                           X       
 /  \                      /  \                        /   \      
T1   P    rightRotate(P)  T1   X     leftRotate(P)    G     P
    / \   =============>      / \    ============>   / \   / \   
   X  T4                    T2   P                 T1  T2 T3  T4
  / \                           / \                
 T2  T3                        T3  T4  

 

代码:

inline void splay(int x, int pos) {
    while (t[x].fa != pos) {//一直旋转成为目标位置的儿子
        int fa = t[x].fa, gfa = t[fa].fa;
        if (gfa != pos) (t[gfa].ch[0] == fa) ^ (t[fa].ch[0] == x) ? rotate(x) : rotate(fa);//判断是哪个儿子并旋转
        rotate(x);//无论哪种情况都要旋转x
    }
    if (pos == 0) root = x;
}

 

插入(insert)

对于一个新的值x

如果x等于根的值,从根节点开始比较节点的val值

如果x==val的话,这个点的计数器++,

x小于val的话向左搜,x大于val的话向右搜

如果不存在某个点的val是x,这时我们即使搜到最底端也没有找到,就直接新建这个节点

因为在插入时可能会形成一条链,在最后的时候还要splay一下把新插入的节点转为根

代码:

inline void insert(int x) {
    int u = root, fa = 0;    //当前位置u,父节点fa
    while (u && t[u].val != x) {    //当u不存在且u的值不等于x。······①
        fa = u;    //向下找u的儿子,父亲为u
        u = t[u].ch[x > t[u].val];    //大于当前位置u向右找,小于向左找
    }
    if (u) t[u].cnt++;    //如果有一个节点的值等于x,计数器++
    else {
        u = ++tot;    //新节点的位置
        if (fa) t[fa].ch[x > t[fa].val] = u;     //如果父节点非根
        t[u].ch[1] = t[u].ch[0] = 0;    //没有儿子
        t[u].fa = fa, t[u].val = x, t[u].cnt = 1, t[u].sum = 1;
    }
    splay(u, 0);    //旋转保持树的平衡
}   

查找(find)

与操作插入操作相似

只需要向左右子树找所查找的数

如果当前点的值等于所查找的数,把当前节点splay到根

inline void find(int x) {   //查找x的位置并旋转到根
    int u = root;
    if (!u) return ;    //空树
    while (t[u].ch[x > t[u].val] && x != t[u].val)    //存在儿子且当前节点的值不等于x。······②
        u = t[u].ch[x > t[u].val];//跳转到儿子
    splay(u, 0);    //旋转到根
}

在初学的时候在这里糊了一下,在这里稍微说明

在insert和find中,一个是当u存在(①),一个是当u的儿子存在(②),当时还试着改了一下代码,结果

bz

其实也很简单

在插入的时候,如果没有一个节点的值等于x,我们在找的时候u会找到树外(u为0,就表示了这个节点不存在),这时我们就新建节点

在查找的时候,不能找出树外,所以要判断u对应的子节点是否存在,不能让u跑到树外面

前驱/后继(nx)

先find一下,把要找的数先转到根

以后继为例,确定后继比x大,所以在右子树里找

有因为后继是右子树里最小的,就在右子树一直向左找,找到叶节点

前驱相反

inline int nx(int x, int f) {    //0 next;1 pre
    find(x);
    int u = root;
    if (t[u].val > x && f) return u;//如果当前节点的值大于x并且要查找的是后继
    if (t[u].val < x && !f) return u;//如果当前节点的值小于x并且要查找的是前驱
    u = t[u].ch[f];    //前驱在左子树里找,后继在右子树里找
    while (t[u].ch[f ^ 1]) u = t[u].ch[f ^ 1];//在另一个方向上找
    return u;
}

 

第k小的数(rank) 

先判断一下是不是有这么多数

看一下左子树的大小,如果k小于左子树大小的话就在左子树里找第k小

如果k大于(左子树大小+当前节点的个数),在右子树上找第(k-左子树大小-当前节点的个数)小

否则,就是根节点的值

 

inline int rank(int x) {
    int u = root;
    if (t[u].sum < x) return 0;    //没有这么多节点
    while (1) {
        int v = t[u].ch[0];    //左子树
        if (x > t[v].sum + t[u].cnt) {    //如果排名大于左子树的大小+当前节点的数量
            x -= t[v].sum + t[u].cnt;
            u = t[u].ch[1];    //当前排名的数一定在右儿子上
        } else if (t[v].sum >= x) u = v;    //在左子树上
        else return t[u].val;    //根节点
    }
}

 

删除(Del)

删除一个点的话

把前驱转到根,把后继转到前驱的下面

后继比前驱大,在前驱的右子树,当前数比前驱大,在前驱的右子树

而在右子树内比后继小的只有当前数,在后继的左子树,所以直接删去后继的左子树

 

inline void Del(int x) {
    int last = nx(x, 0), nxt = nx(x, 1);    //前驱,后继
    splay(last, 0), splay(nxt, last);
    int del = t[nxt].ch[0];    //后继的左子树
    if (t[del].cnt > 1) {    //超过一个
        t[del].cnt--;    //计数--
        splay(del, 0);
    } else t[nxt].ch[0] = 0;    //删除
}

 

模板:

P3369 【模板】普通平衡树

#include <bits/stdc++.h>
using namespace std;
const int N = 1e5 + 10;
int n, m, tot, root;
struct tree {
    int fa, cnt, sum, val;
    int ch[2];
} t[N];

template<class T>inline void read(T &x) {
    x = 0; int f = 0; char ch = getchar();
    while (!isdigit(ch)) f |= (ch == '-'), ch = getchar();
    while (isdigit(ch)) x = x * 10 + ch - '0', ch = getchar();
    x = f ? -x : x;
    return ;
}

inline int get(int x) {
    return t[t[x].fa].ch[0] == x ? 0 : 1;
}

inline void pushup(int x) {
    t[x].sum = t[t[x].ch[0]].sum + t[t[x].ch[1]].sum + t[x].cnt;
}

inline void rotate(int x) {
    int fa = t[x].fa, gfa = t[fa].fa;
    int k = get(x);

    t[x].fa = gfa;
    t[gfa].ch[get(fa)] = x;

    t[t[x].ch[k ^ 1]].fa = fa;
    t[fa].ch[k] = t[x].ch[k ^ 1];

    t[fa].fa = x;
    t[x].ch[k ^ 1] = fa;
    pushup(fa), pushup(x);
}

inline void splay(int x, int pos) {
    while (t[x].fa != pos) {
        int fa = t[x].fa, gfa = t[fa].fa;
        if (gfa != pos) (t[fa].ch[0] == x) ^ (t[gfa].ch[0] == fa) ? rotate(x) : rotate(fa);
        rotate(x);
    }
    if (pos == 0) root = x;
}

inline void find(int x) {                             //查找x的位置并旋转到根
    int u = root;
    if (!u) return ;
    while (t[u].ch[x > t[u].val] && x != t[u].val) u = t[u].ch[x > t[u].val];
    splay(u, 0);
}

inline void insert(int x) {
    int u = root, fa = 0;                            //当前位置u,u的父节点ff
    while (u && t[u].val != x) {
        fa = u;
        u = t[u].ch[x > t[u].val];
    }
    if (u) t[u].cnt++;
    else {
        u = ++tot;
        if (fa) t[fa].ch[x > t[fa].val] = u;     //如果父节点非根
        t[u].ch[1] = t[u].ch[0] = 0;                  //没有儿子
        t[u].fa = fa, t[u].val = x, t[u].cnt = 1, t[u].sum = 1;
    }
    splay(u, 0);
}

inline int nx(int x, int f) {                        //0 next;1 pre
    find(x);
    int u = root;
    if (t[u].val > x && f) return u;
    if (t[u].val < x && !f) return u;
    u = t[u].ch[f];                                   //后继往左找,前驱往右找
    while (t[u].ch[f ^ 1]) u = t[u].ch[f ^ 1];
    return u;
}


inline int rank(int x) {
    int u = root;
    if (t[u].sum < x) return 0;
    while (1) {
        int v = t[u].ch[0];
        if (x > t[v].sum + t[u].cnt) {                //如果排名比左儿子的大小和当前节点的数量要大
            x -= t[v].sum + t[u].cnt;                 //那么当前排名的数一定在右儿子上找
            u = t[u].ch[1];
        } else if (t[v].sum >= x) u = v;
        else return t[u].val;
    }
}

inline void Del(int x) {
    int last = nx(x, 0), nxt = nx(x, 1);
    splay(last, 0), splay(nxt, last);
    int del = t[nxt].ch[0];
    if (t[del].cnt > 1) {
        t[del].cnt--;
        splay(del, 0);
    } else t[nxt].ch[0] = 0;
}

int main(int argc, char const *argv[]) {
    insert(2147483647), insert(-2147483647);
    read(n);
    while (n --) {
        int opt, k;
        read(opt);
        if (opt == 1) read(k), insert(k);
        else if (opt == 2) read(k), Del(k);
        else if (opt == 3) {
            read(k);
            find(k);
            printf("%d\n", t[t[root].ch[0]].sum);
        }
        else if (opt == 4) {
            read(k);
            printf("%d\n", rank(k + 1));
        }
        else if (opt == 5) {
            read(k);
            printf("%d\n", t[nx(k, 0)].val);
        }
        else if (opt == 6) {
            read(k);
            printf("%d\n", t[nx(k, 1)].val);
        }
    }
    return 0;
}
Splay模板

 

posted @ 2019-02-13 22:23  Chrety  阅读(449)  评论(2编辑  收藏  举报