《算法导论》读书笔记之第13章 红黑树

摘要:

  红黑树是一种二叉查找树,但在每个结点上增加了一个存储位表示结点的颜色,可以是RED或者BLACK。通过对任何一条从根到叶子的路径上各个着色方式的限制,红黑树确保没有一条路径会比其他路径长出两倍,因而是接近平衡的。本章主要介绍了红黑树的性质、左右旋转、插入和删除。重点分析了在红黑树中插入和删除元素的过程,分情况进行详细讨论。一棵高度为h的二叉查找树可以实现任何一种基本的动态集合操作,如SEARCH、PREDECESSOR、SUCCESSOR、MIMMUM、MAXMUM、INSERT、DELETE等。当二叉查找树的高度较低时,这些操作执行的比较快,但是当树的高度较高时,这些操作的性能可能不比用链表好。红黑树(red-black tree)是一种平衡的二叉查找树,它能保证在最坏情况下,基本的动态操作集合运行时间为O(lgn)。本章内容有些复杂,看了两天,才大概清楚其插入和删除过程,日后需要经常回顾,争取完全消化掉。红黑树的用途非常广泛,例如STL中的map就是采用红黑树实现的,效率非常之高,有机会可以研究一下STL的源代码。

1、红黑树的性质

  红黑树中的每个结点包含五个域:color、key、left、right和parent。如果某结点没有一个子结点或父结点,则该结点相应的指针parent域包含值为NIL(NIL并是是空指针,此处有些迷惑,一会解释)。把NIL视为指向红黑树的外结点(叶子)的指针,而把带关键字的结点视为红黑树的内结点。红黑树结点结构如下所示:

 1 #define RED  0
 2 #define BLACK 1
 3 struct RedBlackTreeNode
 4 { 
 5     T key;
 6     struct RedBlackTreeNode * parent;
 7     struct RedBlackTreeNode * left;
 8     struct RedBlackTreeNode * right;
 9     int color;
10 };

红黑树的性质如下:

(1)每个结点或是红色,或是黑色。

(2)根结点是黑色。

(3)每个叶子结点(NIL)是黑色。

(4)如果有一个结点是红色,则它的两个儿子都是黑色。

(5)对每个结点,从该结点到其孙子结点的所有路径上包含相同数目的黑色结点。

如下图是一棵红黑树:

从图可以看出NIL不是空指针,而是一个叶子结点。实际操作的时候可以将NIL视为哨兵,这样便于对黑红色进行操作。红黑树的操作主要是对内部结点操作,因为内部结点存储了关键字的值。书中为了便于讨论,忽略了叶子结点的,如是上图红黑树变成如下图所示:

  书中给出了黑高度的概念:从某个结点x出发(不包含该结点)到达一个叶子结点的任意一条路径上,黑色结点的个数称为该结点的黑高度。由红黑树的性质(5)可知,从该结点出发的所有下降路径都有相同的黑色结点个数。红黑树的黑高度定义为其根结点的黑高度。

  书中给出了一个引理来说明为什么红黑树是一种好的查找树,并对引理进行了证明(采用归纳法进行证明的,需要很强的归纳推理知识,正是我的不足之处,看书的痛苦在于此)。

引理:一棵有n个内结点的红黑树的高度之多为2lg(n+1)。

2、旋转

  在红黑树上进行结点插入和删除操作时,会改变树的结构形状,导致结果可能不满足了红黑树的某些性质,为了保证每次插入和删除操作后,仍然能报维持红黑树的性质,需要改变树中某些结点的颜色和指针结构。其中的指针结构的改变通过旋转完成的。书中给出了两种旋转:左旋转和右旋转。如下图是旋转过程:

  从图可以得出左右旋转的过程,假设对某个结点x进行左旋转,y是x的右孩子,则左旋转过程为:以x和y之间的链为“支轴”进行的,使得x的右孩子为y的左孩子,y的父节点为x的父节点,y的左孩子为x。书中给出了左旋转的伪代码如下:

 1 LEFT_ROTATE(T,x)
 2    y = right[x]   //获取右孩子
 3    rihgt[x] = left[y]  //设置x的右孩子为y的左孩子
 4    if left[y] != NIL
 5        then parent[left[x]] = x
 6     parent[y] = parent[x]  //设置y的父节点为x的父节点
 7     if parent[x] == NIL
 8        then root[T] = y
 9        else if x==left[parent[x]
10               then left[parent[x]] = y
11               else  right[[parent[x]] = y
12     left[y] = x  //设置y的左孩子为x
13     parent[x] =y
14 
15    

假设对某个结点y进行右旋转,x是y的左孩子,则左旋转过程为:y的左孩子设置为x的右孩子,将x的父节点设置为y的父节点,x的右孩子设置为y。书中并没有给出右旋转的伪代码,参照左旋转的伪代码很容易实现:

 1  RIGHT_ROTATE(T,y)
 2      x = left[y]    //获取左孩子
 3      left[y] = right[x] //设置y的左孩子为x的右孩子
 4      if right[x] != NIL
 5         then parent[right[x]] = y
 6      parent[x] = parent[y]  //设为x的父节点为y的父结点
 7      if parent[y] == NIL
 8          then root = x
 9          else if y== left[parent[y]]
10                then left[parent[y]] = x
11                else  right[parent[y]] = x
12      right[x] = y //设置x的右孩子为y
13      parent[y] = x

 为了更好的理解旋转操作,书中给出了一个左旋转的例如,如下图所示:

3、插入

  红黑树插入一个新结点的过程RB_INSERT是在二叉查找树插入过程的基础上改进的,先按照二叉排序的插入过程插入到红黑树中,然后将新插入的结点标记为红色(疑问:为什么是红色,而不是黑色呢?),然后调用一个辅助的过程RB_INSERT_FIXUP来调整结点并重新着色,使得满足红黑树的性质。关于二叉查找树的插入过程可以参考上一篇日志:http://www.cnblogs.com/Anker/archive/2013/01/28/2880581.html。书中给出了RB_INSERT的伪代码:

 1 RB_INSERT(T,z)
 2   y = NIL
 3   x =root(T)
 4   while x != NIL
 5        do y=x
 6            if key[z]<key[x]
 7              then x=left[x]
 8              else  x=right[x]
 9   parent[z] = y
10   if y =NIL
11      then root =z
12      else if key[z] < key[y]
13             then left[y] =z
14             else  right[y] =z
15    left[z] = NIL
16    right[z] =NIL
17    color[z] = RED  //新插入结点标记为红色
18    RB_INSERT_FIXUP(T,z)  //进行调整,使得满足红黑树性质

  红黑树的插入过程最主要的是RB_INSERT_FIXUP过程,书中发了很大的篇幅进行介绍。首先分析了当插入一个新的结点后,会破坏红黑树的哪些性质,然后针对可能的破坏性质进行分类讨论并给出了给出了解决办法。因为每次插入的新元素标记为RED,这样可能性质2(根节点为黑色)和性质4(一个红结点的左右孩子都是黑色的)被破坏。例如下图插入一个新结点,破坏了性质4。

      如果每次插入新的结点z导致红黑树性质被破坏,则之多只有一个性质被破坏,并且不是性质2就是性质4。违反性质2是因为z是根且为红色,违反性质4是因为z和其父节点parent[z]都是红色的。

  如果性质2被违反了,则红色的根必定是新增的结点z,它是树中唯一的内结点,由于z的父接点和两个子女都是NIL(黑色),不违反性质4。违反性质2在整个插入过程中只有这一次。所以对于违反性质2的结点,直接将其结点变成黑色即可。

  剩下的问题是对于违反性质4的处理,在插入新结点z之前,红黑树的性质没有被破坏。插入结点z后违反性质4,必定是因为z和其父亲结点parent[z]都是红色的,此时只违反性质4,而没有违反其他性质。假设新插入结点z,导致红黑树性质4被破坏,此时z和其父节点parent[z]都是红色,由于在插入结点z之前红黑树的性质没有被破坏,parent[z]是红色,很容易推出z的祖父结点parent[parent[z]]必定是黑色。此时根据parent[z]是parent[parent[z]]的左孩子还是右孩子进行讨论。因为左右之间是对称的,书中只给出了parent[z]作为parent[parent[z]]的左孩子进行讨论的,然后给出了三种可能的情况。

情况1):z的叔叔结点y是红色的

  此时parent[z]和y都是红色的,解决办法是将z的父节点parent[z]和叔叔结点y都着为黑色,而将z的祖父结点parent[parent[z]]着为红色,然后从祖父结点parent[parent[z]]继续向上判断是否破坏红黑树的性质。处理过程如下图所示:

情况2):z的叔叔y是黑色的,而且z是右孩子

情况3):z的叔叔y是黑色的,而且z是左孩子

  情况2和情况3中y都是黑色的,通过z是左孩子还是右孩子进行区分的。可以将情况2通过旋转为情况3。情况2中z是右孩子,旋转后成为情况3,使得z变为左孩子,可以在parent[z]结点出使用一次左旋转来完成。无论是间接还是直接的通过情况2进入到情况3,z的叔叔y总是黑色的。在情况3中,将parent[z]着为黑色,parent[parent[z]]着为红色,然后从parent[parent[z]]处进行一次右旋转。情况2、3修正了对性质4的违反,修正过程不会导致其他的红黑性质被破坏。修正过程如下图所示:

  给一个完整的例子来说明插入过程,如下图所示:

  书中给出了RB_INSERT_FIXUP的伪代码,伪代码中只给出了z的父节点为左孩子的情况,为右孩子的情况与左孩子的情况是对称的,只需将左孩子中的right换成left即可。

 1 RB_INSERT_FIXUP(T,z)
 2          while color[parent[z]] = RED
 3            do if parent[z] == left[parent[parent[z]]]
 4                   then y = right[parent[parent[z]]]
 5                        if color[y] == RED    //情况1,z的叔叔为红色
 6                             then color[parent[z]] = BLACK
 7                                  color[y] = BLACK
 8                                  color[parent[parent[z]]=RED 
9
z= parent[parent[z]] 10 else if z == right[parent[z]] //情况2,z的叔叔为黑色,z为右孩子 11 then z = parent[z] 12 LEFT_ROTATE(T,z) 13 color[parent[z]]=BLACK //情况3,z的叔叔为黑色,z为左孩子 14 color[parent[parent[z]] = RED 15 RIGHT_ROTATE(T, parent[parent[z]]) 16 else (same as then clause with “right” and “left” exchanged) 17 color(root(T)) = BLACK; //将根结点设置为黑色

4、删除

   删除过程最复杂,看了好多遍才明白个大概,需要反复看,多想删除过程中会破坏哪些性质,然后又针对性的去调整。

  红黑树删除结点过程是在二叉查找树删除结点过程的基础改进的。与二叉查找树类似,删除的结点分为三种情况:<1>无左右孩子、<2>有左孩子或者右孩子、<3>既有树=左孩子又有右孩子。删除过程可以参考前一篇日志:http://www.cnblogs.com/Anker/archive/2013/01/28/2880581.html。红黑树在删除结点后需要检查是否破坏了红黑树的性质。如果删除的结点y是红色的,则删除后的树仍然是保持红黑树的性质,因为树中各个结点的黑高度没有改变,不存在两个相邻(父结点和子结点)的红色结点,y是红色不可能是根,所有根仍然是黑色。如果删除的结点z是黑色的,则这个是破坏了红黑树的性质,需要调用RB_DELETE_FIXUP进行调整。从删除结点y的唯一孩子结点x或者是NIL处开始调整。书中给出了RB_DELETE的伪代码:

 1 RB_DELETE(T,z)
 2      if left[z] ==NIL or right[z] == NIL
 3         then y=z
 4         else  y=TREE_SUCCESSOR(z)
 5     if left[y] != NIL
 6         then x=left[y]
 7         else  x=right[y]
 8     parent[x] = parent[y]
 9     if p[y] ==NIL
10        then root[T] =x
11        else if y = left[[prarnt[y]]
12                    then left[parent[y]] = x
13                    else  right[parent[y]] =x
14      if y!=z
15          then key[z] = key[y]
16                copy y's data into z
17      if color[y] == BLACK    //当被删除结点为黑色时候进行调整
18          then RB_DELETE_FIXUP(T,x)
19       return y

  书中分析了被删除结点y是黑色会产生的问题:首先,如果y是根,而y的一个红色孩子变成了新根,则违反了性质2。其次,如果x和parent[y](此时parent[x] = parent[y])都是红色,就违反了性质4。第三,删除y将会导致先前包含y的任何路径上黑结点个数减少1,违反了性质5。书中给出了解决第三个问题的办法:将结点x设为还有额外的一重黑色(此处看的不是很明白,我的理解是是不管是x是什么颜色,将x增加了额外一重黑色,这样可以保证黑结点数目增加1个),即将任意包含结点x的路径上黑结点个数加1,这样可以保证性质5成立。当将黑色结点y被删除时,将其黑色“下推”至其子结点,导致问题变成为结点x可能即不是红,又不是黑,从而违反性质1。因为给x增加了一种颜色,即结点x是双重黑色或者是红黑色。这样就分别给包含x的路径上黑结点个数贡献2个或1个。但是x的color属性仍然是RED(如果x是红黑的)或BLACK(如果x是双重黑色)。换而言之,一个结点额外的黑色反映在x指向它,而不是它的color属性。

  过程RB_DELETE_FIXUP恢复性质1,2,4。对于恢复性质2、4很简单,因为x是红色,所有直接将x结点着为黑色即可。书中着重介绍了如何恢复性质1。此时x是黑色,需要根据x是左孩子还是右孩子两种情况进行恢复,因为左右是对称的,书中只给出了x是左孩子的恢复过程。将x作为第一个额外的黑色结点,从x结点开始循环,将额外的黑色结点沿着树向上移,直到:

(1)x指向一个红黑结点,此时将x单独着为黑色。

(2)x指向根,这时可以简单地消除那个额外的黑色,或者

(3)做必要的旋转和颜色改变

在循环过程中,x总是指向具有双重黑色的那个非根结点。设w是x的兄弟结点,因为x是双重黑色的,故w不可能是NIL。书中分四种情况讨论:

情况1:x的兄弟w是红色的

          此时因为x是双重黑色,贡献两个黑色结点,所有w必有黑色孩子。此时将w着为黑色,parent[x]为红色,在对parent[x]做一次左旋转。此时x的新兄弟w是黑色,这样将情况1转换为情况2、3或4。情况1的处理过程下图所示:

情况2:x的兄弟w是黑色的,而且w的两个孩子都是黑色的。

     处理过程是从x和w上去掉一重黑色,即x只有一重黑色而w着为红色,给x的父节点parent[x]添加额外黑色。处理过程如下图所示:

 

情况3:x的兄弟w是黑色的,w的左孩子是红色的,右孩子是黑色的

       交换w和其左孩子left[w]的颜色,并对w进行右旋转。旋转后x的新兄弟w是一个有红色右孩子的黑结点,转换成了情况4。处理过程如下图所示:

情况4:x的兄弟w是黑色的,而且w的右孩子是红色的。

  执行过程是将w的颜色设置为parent[x]的颜色,将parent[x]的颜色设置为黑色,将w的右孩子着为黑色,然后在parent[x]做一次右旋,最后将x设置为根root。处理过程如下图所示:

书中给出了RB_DELETE_FIXUP的伪代码:

 1 RB_DELETE_FIXUP(T,x)
 2    while x!= root[T] and color[x] ==BLACK
 3          do if x == left[parent[x]]
 4                then w = right[parent[x]]
 5                     if color[w] == RED  //case 1 x的兄弟w是红色的
 6                        then color[w] = BLACK
 7                             color[parent[x]] = RED
 8                             LEFT_ROTATE(T,PARENT[x])
 9                             w = right[parent[x]]
10                        if color[left[w]] == BLACK and color[right[w]] = BLACK
11                           then color[w] = RED  //case 2
12                             x = parent[x]
13                            else if color[right[w]] =BLACK
14                                   then color[left[w]] = BLACK //case 3
15                                        color[w] = RED
16                                        RIGHT_ROTATE(T,w)
17                                        w = right[parent[x]]
18                                color[w] = color[parent[x]] //case 4
19                                color[parent[x]] = BLACK
20                                color[right[w]] = BLACK
21                                LEFT_ROTATE(T,parent[x])
22                                x=root(T)
23              else(same as then clasue with “right” and “left” exchanged)
24      color[x]=BLACK

5、编程实现

  这一章看了两天,宏观上把握了红黑树的插入和删除操作,中间还有细节问题需要思考。看完后要实现才能消化,于是我采用C++语言设计了简单的红黑树结点和红黑树类,设计的类如下所示:

 1 static const int RED = 0;
 2 static const int BLACK = 1;
 3 
 4 template <class T>
 5 class RedBlackTreeNode
 6 {
 7 public:
 8     RedBlackTreeNode():key(T()),parent(NULL),left(NULL),right(NULL),color(BLACK){}
 9     T key;
10     RedBlackTreeNode<T>* parent;
11     RedBlackTreeNode<T>* left;
12     RedBlackTreeNode<T>* right;
13     int color;
14 };
15 
16 template <class T>
17 class RedBlackTree
18 {
19 public:
20     RedBlackTree();
21     int search_element(const T& k) const;
22     int get_minmum(T& retmin)const;
23     int get_maxmum(T& retmax)const;
24     int get_successor(const T& k,T& ret) const;
25     int get_predecessor(const T& k,T& ret) const;
26     int insert_key(const T& k);
27     int delete_key(const T& k);
28     void inorder_tree_walk()const;
29     RedBlackTreeNode<T>* get_root() const;
30     ~RedBlackTree();
31 private:
32     RedBlackTreeNode<T>* root;
33     static  RedBlackTreeNode<T> *NIL;
34     RedBlackTreeNode<T>* get_parent(RedBlackTreeNode<T>* pnode) const;
35     RedBlackTreeNode<T>* get_left(RedBlackTreeNode<T>* pnode) const;
36     RedBlackTreeNode<T>* get_right(RedBlackTreeNode<T>* pnode) const;
37     T get_key(RedBlackTreeNode<T>* pnode) const;
38     int get_color(RedBlackTreeNode<T>* pnode) const;
39     void set_color(RedBlackTreeNode<T>* pnode,int color);
40     void left_rotate(RedBlackTreeNode<T> *pnode);
41     void right_rotate(RedBlackTreeNode<T> *pnode);
42     void rb_insert_fixup(RedBlackTreeNode<T> *pnode);
43     void rb_delete_fixup(RedBlackTreeNode<T> *pnode);
44     RedBlackTreeNode<T>* get_maxmum(RedBlackTreeNode<T> *root) const;
45     RedBlackTreeNode<T>* get_minmum(RedBlackTreeNode<T> *root) const;
46     RedBlackTreeNode<T>* get_successor(RedBlackTreeNode<T> *pnode) const;
47     RedBlackTreeNode<T>* get_predecessor(RedBlackTreeNode<T> *pnode) const;
48     RedBlackTreeNode<T>* search_tree_node(const T& k)const;
49     void make_empty(RedBlackTreeNode<T>* root);
50 };

设计过程中采用了C++的模板类型,这样可以支持多种数据类型,使得程序具备扩展性,完整的程序实现如下所示:

View Code
  1 #include <iostream>
  2 #include <stack>
  3 using namespace std;
  4 
  5 static const int RED = 0;
  6 static const int BLACK = 1;
  7 
  8 template <class T>
  9 class RedBlackTreeNode
 10 {
 11 public:
 12     RedBlackTreeNode():key(T()),parent(NULL),left(NULL),right(NULL),color(BLACK){}
 13     T key;
 14     RedBlackTreeNode<T>* parent;
 15     RedBlackTreeNode<T>* left;
 16     RedBlackTreeNode<T>* right;
 17     int color;
 18 };
 19 
 20 template <class T>
 21 class RedBlackTree
 22 {
 23 public:
 24     RedBlackTree();
 25     int search_element(const T& k) const;
 26     int get_minmum(T& retmin)const;
 27     int get_maxmum(T& retmax)const;
 28     int get_successor(const T& k,T& ret) const;
 29     int get_predecessor(const T& k,T& ret) const;
 30     int insert_key(const T& k);
 31     int delete_key(const T& k);
 32     void inorder_tree_walk()const;
 33     RedBlackTreeNode<T>* get_root() const;
 34     ~RedBlackTree();
 35 private:
 36     RedBlackTreeNode<T>* root;
 37     static  RedBlackTreeNode<T> *NIL;
 38     RedBlackTreeNode<T>* get_parent(RedBlackTreeNode<T>* pnode) const;
 39     RedBlackTreeNode<T>* get_left(RedBlackTreeNode<T>* pnode) const;
 40     RedBlackTreeNode<T>* get_right(RedBlackTreeNode<T>* pnode) const;
 41     T get_key(RedBlackTreeNode<T>* pnode) const;
 42     int get_color(RedBlackTreeNode<T>* pnode) const;
 43     void set_color(RedBlackTreeNode<T>* pnode,int color);
 44     void left_rotate(RedBlackTreeNode<T> *pnode);
 45     void right_rotate(RedBlackTreeNode<T> *pnode);
 46     void rb_insert_fixup(RedBlackTreeNode<T> *pnode);
 47     void rb_delete_fixup(RedBlackTreeNode<T> *pnode);
 48     RedBlackTreeNode<T>* get_maxmum(RedBlackTreeNode<T> *root) const;
 49     RedBlackTreeNode<T>* get_minmum(RedBlackTreeNode<T> *root) const;
 50     RedBlackTreeNode<T>* get_successor(RedBlackTreeNode<T> *pnode) const;
 51     RedBlackTreeNode<T>* get_predecessor(RedBlackTreeNode<T> *pnode) const;
 52     RedBlackTreeNode<T>* search_tree_node(const T& k)const;
 53     void make_empty(RedBlackTreeNode<T>* root);
 54 };
 55 
 56 template <class T>
 57 RedBlackTreeNode<T>* RedBlackTree<T>::NIL = new RedBlackTreeNode<T>;
 58 
 59 template <class T>
 60 RedBlackTree<T>::RedBlackTree()
 61 {
 62     root = NULL;
 63 }
 64 
 65 template <class T>
 66 int RedBlackTree<T>::search_element(const T& k) const
 67 {
 68     return (NIL != search_tree_node(k));
 69 }
 70 
 71 template <class T>
 72 int RedBlackTree<T>::get_minmum(T& retmin)const
 73 {
 74     if(root)
 75     {
 76         retmin = get_minmum(root)->key;
 77         return 0;
 78     }
 79     return -1;
 80 }
 81 
 82 template <class T>
 83 int RedBlackTree<T>::get_maxmum(T& retmax)const
 84 {
 85     if(root)
 86     {
 87         retmax = get_maxmum(root)->key;
 88         return 0;
 89     }
 90     return -1;
 91 }
 92 
 93 template <class T>
 94 int RedBlackTree<T>::get_successor(const T& k,T& ret) const
 95 {
 96     RedBlackTreeNode<T>* pnode = search_tree_node(k);
 97 
 98     if(pnode != NIL)
 99     {
100         pnode = get_successor(pnode);
101         if(pnode != NIL)
102         {
103             ret = pnode->key;
104             return 0;
105         }
106         return -1;
107     }
108     return -1;
109 }
110 template <class T>
111 int RedBlackTree<T>::get_predecessor(const T& k,T& ret) const
112 {
113     RedBlackTreeNode<T>* pnode = search_tree_node(k);
114     if(pnode != NIL)
115     {
116         pnode = get_predecessor(pnode);
117         if(pnode != NIL)
118         {
119             ret = pnode->key;
120             return 0;
121         }
122         return -1;
123     }
124     return -1;
125 }
126 
127 template <class T>
128 int RedBlackTree<T>::insert_key(const T& k)
129 {
130     RedBlackTreeNode<T> *newnode = new RedBlackTreeNode<T>;
131     newnode->key = k;
132     newnode->color = RED;
133     newnode->left = NIL;
134     newnode->right = NIL;
135     newnode->parent = NIL;
136 
137     if(NULL == root)
138         root = newnode;
139     else
140     {
141         RedBlackTreeNode<T>* pnode = root;
142         RedBlackTreeNode<T>* qnode;
143         while(pnode != NIL)
144         {
145             qnode = pnode;
146             if(pnode->key > newnode->key)
147                 pnode = pnode->left;
148             else
149                 pnode = pnode->right;
150         }
151         newnode->parent = qnode;
152         if(qnode->key > newnode->key)
153             qnode->left = newnode;
154         else
155             qnode->right = newnode;
156     }
157     rb_insert_fixup(newnode);
158     return 0;
159 }
160 
161 template <class T>
162 int RedBlackTree<T>::delete_key(const T& k)
163 {
164     RedBlackTreeNode<T>* pnode = search_tree_node(k);
165     if(NIL != pnode)
166     {
167         RedBlackTreeNode<T>* qnode,*tnode;
168         if(get_left(pnode) == NIL || get_right(pnode) == NIL)
169             qnode = pnode;
170         else
171             qnode = get_successor(pnode);
172         if(get_left(qnode) != NIL)
173             tnode = get_left(qnode);
174         else
175             tnode = get_right(qnode);
176         tnode->parent = get_parent(qnode);
177         if(get_parent(qnode) == NIL)
178             root = tnode;
179         else if(qnode == get_left(get_parent(qnode)))
180             qnode->parent->left = tnode;
181         else
182             qnode->parent->right = tnode;
183         if(qnode != pnode)
184              pnode->key = get_key(qnode);
185         if(get_color(qnode) == BLACK)
186             rb_delete_fixup(tnode);
187         delete qnode;
188         return 0;
189     }
190     return -1;
191 }
192 
193 template <class T>
194 RedBlackTree<T>::~RedBlackTree()
195 {
196     make_empty(root);
197 }
198 template <class T>
199 RedBlackTreeNode<T>* RedBlackTree<T>:: get_root() const
200 {
201     return root;
202 }
203 template <class T>
204 RedBlackTreeNode<T>* RedBlackTree<T>::get_parent(RedBlackTreeNode<T>* pnode) const
205 {
206     return pnode->parent;
207 }
208 template <class T>
209 RedBlackTreeNode<T>* RedBlackTree<T>::get_left(RedBlackTreeNode<T>* pnode) const
210 {
211     return pnode->left;
212 }
213 template <class T>
214 RedBlackTreeNode<T>* RedBlackTree<T>::get_right(RedBlackTreeNode<T>* pnode) const
215 {
216     return pnode->right;
217 }
218 template <class T>
219 T RedBlackTree<T>::get_key(RedBlackTreeNode<T>* pnode) const
220 {
221     return pnode->key;
222 }
223 template <class T>
224 int RedBlackTree<T>::get_color(RedBlackTreeNode<T>* pnode) const
225 {
226     return pnode->color;
227 }
228 
229 template <class T>
230 void RedBlackTree<T>::set_color(RedBlackTreeNode<T>* pnode,int color)
231 {
232     pnode->color = color;
233 }
234 
235 template <class T>
236 void RedBlackTree<T>::left_rotate(RedBlackTreeNode<T> *pnode)
237 {
238     RedBlackTreeNode<T>* rightnode = pnode->right;
239     pnode->right = rightnode->left;
240     if(rightnode->left != NIL)
241         rightnode->left->parent = pnode;
242     rightnode->parent = pnode->parent;
243     if(pnode->parent == NIL)
244         root = rightnode;
245     else if(pnode == pnode->parent->left)
246         pnode->parent->left = rightnode;
247     else
248         pnode->parent->right = rightnode;
249     rightnode->left = pnode;
250     pnode->parent = rightnode;
251 }
252 
253 template <class T>
254 void RedBlackTree<T>::right_rotate(RedBlackTreeNode<T> *pnode)
255 {
256     RedBlackTreeNode<T>* leftnode = pnode->left;
257     pnode->left = leftnode->right;
258     if(leftnode->right != NIL)
259         leftnode->right->parent = pnode;
260     leftnode->parent = pnode->parent;
261     if(pnode->parent == NIL)
262         root = leftnode;
263     else if(pnode == pnode->parent->left)
264         pnode->parent->left = leftnode;
265     else
266         pnode->parent->right = leftnode;
267     leftnode->right = pnode;
268     pnode->parent = leftnode;
269 }
270 template <class T>
271 void RedBlackTree<T>::rb_insert_fixup(RedBlackTreeNode<T>*pnode)
272 {
273     RedBlackTreeNode<T> *qnode,*tnode;
274     //当pnode的父节点为红色时,破坏性质4
275     while(get_color(get_parent(pnode))== RED)
276     {
277         qnode = get_parent(get_parent(pnode));//祖父结点
278         if(get_parent(pnode) == get_left(qnode))
279         {
280             tnode = get_right(qnode);//pnode的叔叔结点
281             if(get_color(tnode) == RED) //case1 叔叔结点为红色
282             {
283                 set_color(get_parent(pnode),BLACK);
284                 set_color(tnode,BLACK);
285                 set_color(qnode,RED);
286                 pnode = qnode;
287             }
288             else  //case 2 or case 3
289             {
290                 if(pnode == get_right(get_parent(pnode)))  //case2 pnode为右孩子
291                 {
292                     pnode = get_parent(pnode);
293                     left_rotate(pnode);
294                 }
295                 //case3 pnode为左孩子
296                 set_color(get_parent(pnode),BLACK);
297                 qnode = get_parent(get_parent(pnode));
298                 set_color(qnode,RED);
299                 right_rotate(qnode);
300             }
301         }
302         else
303         {
304             tnode = get_left(qnode);
305             if(get_color(tnode) == RED)
306             {
307                 set_color(get_parent(pnode),BLACK);
308                 set_color(tnode,BLACK);
309                 set_color(qnode,RED);
310                 pnode = qnode;
311             }
312             else
313             {
314                 if(pnode == get_left(get_parent(pnode)))
315                 {
316                     pnode = get_parent(pnode);
317                     right_rotate(pnode);
318                 }
319                 set_color(get_parent(pnode),BLACK);
320                 qnode = get_parent(get_parent(pnode));
321                 set_color(qnode,RED);
322                 left_rotate(qnode);
323             }
324         }
325     }
326     set_color(root,BLACK);
327 }
328 
329 template <class T>
330 void RedBlackTree<T>::rb_delete_fixup(RedBlackTreeNode<T> *pnode)
331 {
332     while(pnode != root && get_color(pnode) == BLACK)
333     {
334         RedBlackTreeNode<T> *qnode,*tnode;
335         if(pnode == get_left(get_parent(pnode)))
336         {
337             qnode = get_right(get_parent(pnode));
338             if(get_color(qnode) == RED)
339             {
340                 set_color(qnode,BLACK);
341                 set_color(get_parent(pnode),RED);
342                 left_rotate(get_parent(pnode));
343                 qnode = get_right(get_parent(pnode));
344             }
345             if(get_color(get_left(qnode)) == BLACK && get_color(get_right(qnode)) == BLACK)
346             {
347                 set_color(qnode,RED);
348                 pnode = get_parent(pnode);
349             }
350             else
351             {
352                 if(get_color(get_right(qnode)) == BLACK)
353                 {
354                     set_color(get_left(qnode),BLACK);
355                     set_color(qnode,RED);
356                     right_rotate(qnode);
357                     qnode = get_right(get_parent(pnode));
358                 }
359                 set_color(qnode,get_color(get_parent(pnode)));
360                 set_color(get_parent(pnode),BLACK);
361                 set_color(get_right(qnode),BLACK);
362                 left_rotate(get_parent(pnode));
363                 pnode = root;
364             }
365         }
366         else
367         {
368             qnode = get_left(get_parent(pnode));
369             if(get_color(qnode) == RED)
370             {
371                 set_color(qnode,BLACK);
372                 set_color(get_parent(pnode),RED);
373                 right_rotate(get_parent(pnode));
374                 qnode = get_left(get_parent(pnode));
375             }
376             if(get_color(get_right(qnode)) == BLACK && get_color(get_left(qnode)) == BLACK)
377             {
378                 set_color(qnode,RED);
379                 pnode = get_parent(pnode);
380             }
381             else
382             {
383                 if(get_color(get_left(qnode)) == BLACK)
384                 {
385                     set_color(get_right(qnode),BLACK);
386                     set_color(qnode,RED);
387                     left_rotate(qnode);
388                     qnode = get_left(get_parent(pnode));
389                 }
390                 set_color(qnode,get_color(get_parent(pnode)));
391                 set_color(get_parent(pnode),BLACK);
392                 set_color(get_left(qnode),BLACK);
393                 right_rotate(get_parent(pnode));
394                 pnode = root;
395             }
396         }
397     }
398     set_color(pnode,BLACK);
399 }
400 
401 template <class T>
402 RedBlackTreeNode<T>* RedBlackTree<T>::get_maxmum(RedBlackTreeNode<T> *root) const
403 {
404     RedBlackTreeNode<T> *pnode = root;
405     while(pnode->right != NIL)
406         pnode = pnode->right;
407     return pnode;
408 }
409 
410 template <class T>
411 RedBlackTreeNode<T>* RedBlackTree<T>::get_minmum(RedBlackTreeNode<T> *root) const
412 {
413     RedBlackTreeNode<T> *pnode = root;
414     while(pnode->left != NIL)
415         pnode = pnode->left;
416     return pnode;
417 }
418 
419 template <class T>
420 RedBlackTreeNode<T>* RedBlackTree<T>:: get_successor(RedBlackTreeNode<T> *pnode) const
421 {
422     if(pnode->right != NIL)
423         return get_minmum(pnode->right);
424     RedBlackTreeNode<T>* parentnode = get_parent(pnode);
425     while(parentnode != NIL && get_right(parentnode) == pnode)
426     {
427             pnode = parentnode;
428             parentnode = get_parent(pnode);
429     }
430     return parentnode;
431 }
432 
433 template <class T>
434 RedBlackTreeNode<T>* RedBlackTree<T>::get_predecessor(RedBlackTreeNode<T> *pnode) const
435 {
436     if(pnode->left != NIL)
437         return get_maxmum(pnode->left);
438     RedBlackTreeNode<T>* parentnode = get_parent(pnode);
439     while(parentnode != NIL && get_left(parentnode) == pnode)
440     {
441         pnode = parentnode;
442         parentnode = get_parent(pnode);
443     }
444     return parentnode;
445 }
446 
447 template <class T>
448 RedBlackTreeNode<T>* RedBlackTree<T>:: search_tree_node(const T& k)const
449 {
450     RedBlackTreeNode<T>* pnode = root;
451     while(pnode != NIL)
452     {
453         if(pnode->key == k)
454             break;
455         else if(pnode->key > k)
456             pnode = pnode->left;
457         else
458             pnode = pnode->right;
459     }
460     return pnode;
461 }
462 
463 template <class T>
464 void RedBlackTree<T>::make_empty(RedBlackTreeNode<T>* root)
465 {
466     if(root)
467     {
468         RedBlackTreeNode<T> *pleft = root->left;
469         RedBlackTreeNode<T>* pright = root->right;
470         delete root;
471         if(pleft != NIL)
472             make_empty(pleft);
473         if(pright != NIL)
474             make_empty(pright);
475     }
476 }
477 template <class T>
478 void RedBlackTree<T>::inorder_tree_walk()const
479 {
480     if(NULL != root)
481      {
482          stack<RedBlackTreeNode<T>* > s;
483          RedBlackTreeNode<T> *ptmpnode;
484          ptmpnode = root;
485          while(NIL != ptmpnode || !s.empty())
486          {
487              if(NIL != ptmpnode)
488              {
489                  s.push(ptmpnode);
490                  ptmpnode = ptmpnode->left;
491              }
492              else
493              {
494                  ptmpnode = s.top();
495                  s.pop();
496                  cout<<ptmpnode->key<<":";
497                  if(ptmpnode->color == BLACK)
498                     cout<<"Black"<<endl;
499                  else
500                     cout<<"Red"<<endl;
501                  ptmpnode = ptmpnode->right;
502              }
503          }
504      }
505 }
506 int main()
507 {
508     RedBlackTree<int> rbtree;
509     int value;
510     rbtree.insert_key(41);
511     rbtree.insert_key(38);
512     rbtree.insert_key(31);
513     rbtree.insert_key(12);
514     rbtree.insert_key(19);
515     rbtree.insert_key(8);
516     cout<<"root is: "<<rbtree.get_root()->key<<endl;
517     cout<<"Inorder walk red black tree:"<<endl;
518     rbtree.inorder_tree_walk();
519     if(rbtree.get_minmum(value) == 0)
520         cout<<"minmum is: "<<value<<endl;
521     if(rbtree.get_maxmum(value) == 0)
522         cout<<"maxmum is: "<<value<<endl;
523     if(rbtree.get_successor(19,value) == 0)
524         cout<<"19 successor is: "<<value<<endl;
525     if(rbtree.get_predecessor(19,value) == 0)
526         cout<<"19 predecessor is: "<<value<<endl;
527     if(rbtree.delete_key(38)==0)
528         cout<<"delete 38 successfully"<<endl;
529      cout<<"root is: "<<rbtree.get_root()->key<<endl;
530     cout<<"Inorder walk red black tree:"<<endl;
531     rbtree.inorder_tree_walk();
532     return 0;
533 }

程序测试结果如下所示:


  实现过程中感触非常多,需要很大的耐心去调试程序,关键的是insert和delete操作。

posted @ 2013-01-30 12:25  Rabbit_Dale  阅读(15541)  评论(11编辑  收藏  举报