stl_tree.h
上传用户:kellyonhid
上传日期:2013-10-12
资源大小:932k
文件大小:43k
源码类别:

3D图形编程

开发平台:

Visual C++

  1. /*
  2.  *
  3.  * Copyright (c) 1996,1997
  4.  * Silicon Graphics Computer Systems, Inc.
  5.  *
  6.  * Permission to use, copy, modify, distribute and sell this software
  7.  * and its documentation for any purpose is hereby granted without fee,
  8.  * provided that the above copyright notice appear in all copies and
  9.  * that both that copyright notice and this permission notice appear
  10.  * in supporting documentation.  Silicon Graphics makes no
  11.  * representations about the suitability of this software for any
  12.  * purpose.  It is provided "as is" without express or implied warranty.
  13.  *
  14.  *
  15.  * Copyright (c) 1994
  16.  * Hewlett-Packard Company
  17.  *
  18.  * Permission to use, copy, modify, distribute and sell this software
  19.  * and its documentation for any purpose is hereby granted without fee,
  20.  * provided that the above copyright notice appear in all copies and
  21.  * that both that copyright notice and this permission notice appear
  22.  * in supporting documentation.  Hewlett-Packard Company makes no
  23.  * representations about the suitability of this software for any
  24.  * purpose.  It is provided "as is" without express or implied warranty.
  25.  *
  26.  *
  27.  */
  28. /* NOTE: This is an internal header file, included by other STL headers.
  29.  *   You should not attempt to use it directly.
  30.  */
  31. #ifndef __SGI_STL_INTERNAL_TREE_H
  32. #define __SGI_STL_INTERNAL_TREE_H
  33. /*
  34. Red-black tree class, designed for use in implementing STL
  35. associative containers (set, multiset, map, and multimap). The
  36. insertion and deletion algorithms are based on those in Cormen,
  37. Leiserson, and Rivest, Introduction to Algorithms (MIT Press, 1990),
  38. except that
  39. (1) the header cell is maintained with links not only to the root
  40. but also to the leftmost node of the tree, to enable constant time
  41. begin(), and to the rightmost node of the tree, to enable linear time
  42. performance when used with the generic set algorithms (set_union,
  43. etc.);
  44. (2) when a node being deleted has two children its successor node is
  45. relinked into its place, rather than copied, so that the only
  46. iterators invalidated are those referring to the deleted node.
  47. */
  48. #include <stl_algobase.h>
  49. #include <stl_alloc.h>
  50. #include <stl_construct.h>
  51. #include <stl_function.h>
  52. __STL_BEGIN_NAMESPACE 
  53. #if defined(__sgi) && !defined(__GNUC__) && (_MIPS_SIM != _MIPS_SIM_ABI32)
  54. #pragma set woff 1375
  55. #endif
  56. typedef bool _Rb_tree_Color_type;
  57. const _Rb_tree_Color_type _S_rb_tree_red = false;
  58. const _Rb_tree_Color_type _S_rb_tree_black = true;
  59. struct _Rb_tree_node_base
  60. {
  61.   typedef _Rb_tree_Color_type _Color_type;
  62.   typedef _Rb_tree_node_base* _Base_ptr;
  63.   _Color_type _M_color; 
  64.   _Base_ptr _M_parent;
  65.   _Base_ptr _M_left;
  66.   _Base_ptr _M_right;
  67.   static _Base_ptr _S_minimum(_Base_ptr __x)
  68.   {
  69.     while (__x->_M_left != 0) __x = __x->_M_left;
  70.     return __x;
  71.   }
  72.   static _Base_ptr _S_maximum(_Base_ptr __x)
  73.   {
  74.     while (__x->_M_right != 0) __x = __x->_M_right;
  75.     return __x;
  76.   }
  77. };
  78. template <class _Value>
  79. struct _Rb_tree_node : public _Rb_tree_node_base
  80. {
  81.   typedef _Rb_tree_node<_Value>* _Link_type;
  82.   _Value _M_value_field;
  83. };
  84. struct _Rb_tree_base_iterator
  85. {
  86.   typedef _Rb_tree_node_base::_Base_ptr _Base_ptr;
  87.   typedef bidirectional_iterator_tag iterator_category;
  88.   typedef ptrdiff_t difference_type;
  89.   _Base_ptr _M_node;
  90.   void _M_increment()
  91.   {
  92.     if (_M_node->_M_right != 0) {
  93.       _M_node = _M_node->_M_right;
  94.       while (_M_node->_M_left != 0)
  95.         _M_node = _M_node->_M_left;
  96.     }
  97.     else {
  98.       _Base_ptr __y = _M_node->_M_parent;
  99.       while (_M_node == __y->_M_right) {
  100.         _M_node = __y;
  101.         __y = __y->_M_parent;
  102.       }
  103.       if (_M_node->_M_right != __y)
  104.         _M_node = __y;
  105.     }
  106.   }
  107.   void _M_decrement()
  108.   {
  109.     if (_M_node->_M_color == _S_rb_tree_red &&
  110.         _M_node->_M_parent->_M_parent == _M_node)
  111.       _M_node = _M_node->_M_right;
  112.     else if (_M_node->_M_left != 0) {
  113.       _Base_ptr __y = _M_node->_M_left;
  114.       while (__y->_M_right != 0)
  115.         __y = __y->_M_right;
  116.       _M_node = __y;
  117.     }
  118.     else {
  119.       _Base_ptr __y = _M_node->_M_parent;
  120.       while (_M_node == __y->_M_left) {
  121.         _M_node = __y;
  122.         __y = __y->_M_parent;
  123.       }
  124.       _M_node = __y;
  125.     }
  126.   }
  127. };
  128. template <class _Value, class _Ref, class _Ptr>
  129. struct _Rb_tree_iterator : public _Rb_tree_base_iterator
  130. {
  131.   typedef _Value value_type;
  132.   typedef _Ref reference;
  133.   typedef _Ptr pointer;
  134.   typedef _Rb_tree_iterator<_Value, _Value&, _Value*>             
  135.     iterator;
  136.   typedef _Rb_tree_iterator<_Value, const _Value&, const _Value*> 
  137.     const_iterator;
  138.   typedef _Rb_tree_iterator<_Value, _Ref, _Ptr>                   
  139.     _Self;
  140.   typedef _Rb_tree_node<_Value>* _Link_type;
  141.   _Rb_tree_iterator() {}
  142.   _Rb_tree_iterator(_Link_type __x) { _M_node = __x; }
  143.   _Rb_tree_iterator(const iterator& __it) { _M_node = __it._M_node; }
  144.   reference operator*() const { return _Link_type(_M_node)->_M_value_field; }
  145. #ifndef __SGI_STL_NO_ARROW_OPERATOR
  146.   pointer operator->() const { return &(operator*()); }
  147. #endif /* __SGI_STL_NO_ARROW_OPERATOR */
  148.   _Self& operator++() { _M_increment(); return *this; }
  149.   _Self operator++(int) {
  150.     _Self __tmp = *this;
  151.     _M_increment();
  152.     return __tmp;
  153.   }
  154.     
  155.   _Self& operator--() { _M_decrement(); return *this; }
  156.   _Self operator--(int) {
  157.     _Self __tmp = *this;
  158.     _M_decrement();
  159.     return __tmp;
  160.   }
  161. };
  162. inline bool operator==(const _Rb_tree_base_iterator& __x,
  163.                        const _Rb_tree_base_iterator& __y) {
  164.   return __x._M_node == __y._M_node;
  165. }
  166. inline bool operator!=(const _Rb_tree_base_iterator& __x,
  167.                        const _Rb_tree_base_iterator& __y) {
  168.   return __x._M_node != __y._M_node;
  169. }
  170. #ifndef __STL_CLASS_PARTIAL_SPECIALIZATION
  171. inline bidirectional_iterator_tag
  172. iterator_category(const _Rb_tree_base_iterator&) {
  173.   return bidirectional_iterator_tag();
  174. }
  175. inline _Rb_tree_base_iterator::difference_type*
  176. distance_type(const _Rb_tree_base_iterator&) {
  177.   return (_Rb_tree_base_iterator::difference_type*) 0;
  178. }
  179. template <class _Value, class _Ref, class _Ptr>
  180. inline _Value* value_type(const _Rb_tree_iterator<_Value, _Ref, _Ptr>&) {
  181.   return (_Value*) 0;
  182. }
  183. #endif /* __STL_CLASS_PARTIAL_SPECIALIZATION */
  184. inline void 
  185. _Rb_tree_rotate_left(_Rb_tree_node_base* __x, _Rb_tree_node_base*& __root)
  186. {
  187.   _Rb_tree_node_base* __y = __x->_M_right;
  188.   __x->_M_right = __y->_M_left;
  189.   if (__y->_M_left !=0)
  190.     __y->_M_left->_M_parent = __x;
  191.   __y->_M_parent = __x->_M_parent;
  192.   if (__x == __root)
  193.     __root = __y;
  194.   else if (__x == __x->_M_parent->_M_left)
  195.     __x->_M_parent->_M_left = __y;
  196.   else
  197.     __x->_M_parent->_M_right = __y;
  198.   __y->_M_left = __x;
  199.   __x->_M_parent = __y;
  200. }
  201. inline void 
  202. _Rb_tree_rotate_right(_Rb_tree_node_base* __x, _Rb_tree_node_base*& __root)
  203. {
  204.   _Rb_tree_node_base* __y = __x->_M_left;
  205.   __x->_M_left = __y->_M_right;
  206.   if (__y->_M_right != 0)
  207.     __y->_M_right->_M_parent = __x;
  208.   __y->_M_parent = __x->_M_parent;
  209.   if (__x == __root)
  210.     __root = __y;
  211.   else if (__x == __x->_M_parent->_M_right)
  212.     __x->_M_parent->_M_right = __y;
  213.   else
  214.     __x->_M_parent->_M_left = __y;
  215.   __y->_M_right = __x;
  216.   __x->_M_parent = __y;
  217. }
  218. inline void 
  219. _Rb_tree_rebalance(_Rb_tree_node_base* __x, _Rb_tree_node_base*& __root)
  220. {
  221.   __x->_M_color = _S_rb_tree_red;
  222.   while (__x != __root && __x->_M_parent->_M_color == _S_rb_tree_red) {
  223.     if (__x->_M_parent == __x->_M_parent->_M_parent->_M_left) {
  224.       _Rb_tree_node_base* __y = __x->_M_parent->_M_parent->_M_right;
  225.       if (__y && __y->_M_color == _S_rb_tree_red) {
  226.         __x->_M_parent->_M_color = _S_rb_tree_black;
  227.         __y->_M_color = _S_rb_tree_black;
  228.         __x->_M_parent->_M_parent->_M_color = _S_rb_tree_red;
  229.         __x = __x->_M_parent->_M_parent;
  230.       }
  231.       else {
  232.         if (__x == __x->_M_parent->_M_right) {
  233.           __x = __x->_M_parent;
  234.           _Rb_tree_rotate_left(__x, __root);
  235.         }
  236.         __x->_M_parent->_M_color = _S_rb_tree_black;
  237.         __x->_M_parent->_M_parent->_M_color = _S_rb_tree_red;
  238.         _Rb_tree_rotate_right(__x->_M_parent->_M_parent, __root);
  239.       }
  240.     }
  241.     else {
  242.       _Rb_tree_node_base* __y = __x->_M_parent->_M_parent->_M_left;
  243.       if (__y && __y->_M_color == _S_rb_tree_red) {
  244.         __x->_M_parent->_M_color = _S_rb_tree_black;
  245.         __y->_M_color = _S_rb_tree_black;
  246.         __x->_M_parent->_M_parent->_M_color = _S_rb_tree_red;
  247.         __x = __x->_M_parent->_M_parent;
  248.       }
  249.       else {
  250.         if (__x == __x->_M_parent->_M_left) {
  251.           __x = __x->_M_parent;
  252.           _Rb_tree_rotate_right(__x, __root);
  253.         }
  254.         __x->_M_parent->_M_color = _S_rb_tree_black;
  255.         __x->_M_parent->_M_parent->_M_color = _S_rb_tree_red;
  256.         _Rb_tree_rotate_left(__x->_M_parent->_M_parent, __root);
  257.       }
  258.     }
  259.   }
  260.   __root->_M_color = _S_rb_tree_black;
  261. }
  262. inline _Rb_tree_node_base*
  263. _Rb_tree_rebalance_for_erase(_Rb_tree_node_base* __z,
  264.                              _Rb_tree_node_base*& __root,
  265.                              _Rb_tree_node_base*& __leftmost,
  266.                              _Rb_tree_node_base*& __rightmost)
  267. {
  268.   _Rb_tree_node_base* __y = __z;
  269.   _Rb_tree_node_base* __x = 0;
  270.   _Rb_tree_node_base* __x_parent = 0;
  271.   if (__y->_M_left == 0)     // __z has at most one non-null child. y == z.
  272.     __x = __y->_M_right;     // __x might be null.
  273.   else
  274.     if (__y->_M_right == 0)  // __z has exactly one non-null child. y == z.
  275.       __x = __y->_M_left;    // __x is not null.
  276.     else {                   // __z has two non-null children.  Set __y to
  277.       __y = __y->_M_right;   //   __z's successor.  __x might be null.
  278.       while (__y->_M_left != 0)
  279.         __y = __y->_M_left;
  280.       __x = __y->_M_right;
  281.     }
  282.   if (__y != __z) {          // relink y in place of z.  y is z's successor
  283.     __z->_M_left->_M_parent = __y; 
  284.     __y->_M_left = __z->_M_left;
  285.     if (__y != __z->_M_right) {
  286.       __x_parent = __y->_M_parent;
  287.       if (__x) __x->_M_parent = __y->_M_parent;
  288.       __y->_M_parent->_M_left = __x;      // __y must be a child of _M_left
  289.       __y->_M_right = __z->_M_right;
  290.       __z->_M_right->_M_parent = __y;
  291.     }
  292.     else
  293.       __x_parent = __y;  
  294.     if (__root == __z)
  295.       __root = __y;
  296.     else if (__z->_M_parent->_M_left == __z)
  297.       __z->_M_parent->_M_left = __y;
  298.     else 
  299.       __z->_M_parent->_M_right = __y;
  300.     __y->_M_parent = __z->_M_parent;
  301.     __STD::swap(__y->_M_color, __z->_M_color);
  302.     __y = __z;
  303.     // __y now points to node to be actually deleted
  304.   }
  305.   else {                        // __y == __z
  306.     __x_parent = __y->_M_parent;
  307.     if (__x) __x->_M_parent = __y->_M_parent;   
  308.     if (__root == __z)
  309.       __root = __x;
  310.     else 
  311.       if (__z->_M_parent->_M_left == __z)
  312.         __z->_M_parent->_M_left = __x;
  313.       else
  314.         __z->_M_parent->_M_right = __x;
  315.     if (__leftmost == __z) 
  316.       if (__z->_M_right == 0)        // __z->_M_left must be null also
  317.         __leftmost = __z->_M_parent;
  318.     // makes __leftmost == _M_header if __z == __root
  319.       else
  320.         __leftmost = _Rb_tree_node_base::_S_minimum(__x);
  321.     if (__rightmost == __z)  
  322.       if (__z->_M_left == 0)         // __z->_M_right must be null also
  323.         __rightmost = __z->_M_parent;  
  324.     // makes __rightmost == _M_header if __z == __root
  325.       else                      // __x == __z->_M_left
  326.         __rightmost = _Rb_tree_node_base::_S_maximum(__x);
  327.   }
  328.   if (__y->_M_color != _S_rb_tree_red) { 
  329.     while (__x != __root && (__x == 0 || __x->_M_color == _S_rb_tree_black))
  330.       if (__x == __x_parent->_M_left) {
  331.         _Rb_tree_node_base* __w = __x_parent->_M_right;
  332.         if (__w->_M_color == _S_rb_tree_red) {
  333.           __w->_M_color = _S_rb_tree_black;
  334.           __x_parent->_M_color = _S_rb_tree_red;
  335.           _Rb_tree_rotate_left(__x_parent, __root);
  336.           __w = __x_parent->_M_right;
  337.         }
  338.         if ((__w->_M_left == 0 || 
  339.              __w->_M_left->_M_color == _S_rb_tree_black) &&
  340.             (__w->_M_right == 0 || 
  341.              __w->_M_right->_M_color == _S_rb_tree_black)) {
  342.           __w->_M_color = _S_rb_tree_red;
  343.           __x = __x_parent;
  344.           __x_parent = __x_parent->_M_parent;
  345.         } else {
  346.           if (__w->_M_right == 0 || 
  347.               __w->_M_right->_M_color == _S_rb_tree_black) {
  348.             if (__w->_M_left) __w->_M_left->_M_color = _S_rb_tree_black;
  349.             __w->_M_color = _S_rb_tree_red;
  350.             _Rb_tree_rotate_right(__w, __root);
  351.             __w = __x_parent->_M_right;
  352.           }
  353.           __w->_M_color = __x_parent->_M_color;
  354.           __x_parent->_M_color = _S_rb_tree_black;
  355.           if (__w->_M_right) __w->_M_right->_M_color = _S_rb_tree_black;
  356.           _Rb_tree_rotate_left(__x_parent, __root);
  357.           break;
  358.         }
  359.       } else {                  // same as above, with _M_right <-> _M_left.
  360.         _Rb_tree_node_base* __w = __x_parent->_M_left;
  361.         if (__w->_M_color == _S_rb_tree_red) {
  362.           __w->_M_color = _S_rb_tree_black;
  363.           __x_parent->_M_color = _S_rb_tree_red;
  364.           _Rb_tree_rotate_right(__x_parent, __root);
  365.           __w = __x_parent->_M_left;
  366.         }
  367.         if ((__w->_M_right == 0 || 
  368.              __w->_M_right->_M_color == _S_rb_tree_black) &&
  369.             (__w->_M_left == 0 || 
  370.              __w->_M_left->_M_color == _S_rb_tree_black)) {
  371.           __w->_M_color = _S_rb_tree_red;
  372.           __x = __x_parent;
  373.           __x_parent = __x_parent->_M_parent;
  374.         } else {
  375.           if (__w->_M_left == 0 || 
  376.               __w->_M_left->_M_color == _S_rb_tree_black) {
  377.             if (__w->_M_right) __w->_M_right->_M_color = _S_rb_tree_black;
  378.             __w->_M_color = _S_rb_tree_red;
  379.             _Rb_tree_rotate_left(__w, __root);
  380.             __w = __x_parent->_M_left;
  381.           }
  382.           __w->_M_color = __x_parent->_M_color;
  383.           __x_parent->_M_color = _S_rb_tree_black;
  384.           if (__w->_M_left) __w->_M_left->_M_color = _S_rb_tree_black;
  385.           _Rb_tree_rotate_right(__x_parent, __root);
  386.           break;
  387.         }
  388.       }
  389.     if (__x) __x->_M_color = _S_rb_tree_black;
  390.   }
  391.   return __y;
  392. }
  393. // Base class to encapsulate the differences between old SGI-style
  394. // allocators and standard-conforming allocators.  In order to avoid
  395. // having an empty base class, we arbitrarily move one of rb_tree's
  396. // data members into the base class.
  397. #ifdef __STL_USE_STD_ALLOCATORS
  398. // _Base for general standard-conforming allocators.
  399. template <class _T, class _Alloc, bool _S_instanceless>
  400. class _Rb_tree_alloc_base {
  401. public:
  402.   typedef typename _Alloc_traits<_T, _Alloc>::allocator_type allocator_type;
  403.   allocator_type get_allocator() const { return _M_node_allocator; }
  404.   _Rb_tree_alloc_base(const allocator_type& __a)
  405.     : _M_node_allocator(__a), _M_header(0) {}
  406. protected:
  407.   typename _Alloc_traits<_Rb_tree_node<_T>, _Alloc>::allocator_type
  408.            _M_node_allocator;
  409.   _Rb_tree_node<_T>* _M_header;
  410.   _Rb_tree_node<_T>* _M_get_node() 
  411.     { return _M_node_allocator.allocate(1); }
  412.   void _M_put_node(_Rb_tree_node<_T>* __p) 
  413.     { _M_node_allocator.deallocate(__p, 1); }
  414. };
  415. // Specialization for instanceless allocators.
  416. template <class _T, class _Alloc>
  417. class _Rb_tree_alloc_base<_T, _Alloc, true> {
  418. public:
  419.   typedef typename _Alloc_traits<_T, _Alloc>::allocator_type allocator_type;
  420.   allocator_type get_allocator() const { return allocator_type(); }
  421.   _Rb_tree_alloc_base(const allocator_type&) : _M_header(0) {}
  422. protected:
  423.   _Rb_tree_node<_T>* _M_header;
  424.   typedef typename _Alloc_traits<_Rb_tree_node<_T>, _Alloc>::_Alloc_type
  425.           _Alloc_type;
  426.   _Rb_tree_node<_T>* _M_get_node()
  427.     { return _Alloc_type::allocate(1); }
  428.   void _M_put_node(_Rb_tree_node<_T>* __p)
  429.     { _Alloc_type::deallocate(__p, 1); }
  430. };
  431. template <class _T, class _Alloc>
  432. struct _Rb_tree_base
  433.   : public _Rb_tree_alloc_base<_T, _Alloc,
  434.                                _Alloc_traits<_T, _Alloc>::_S_instanceless>
  435. {
  436.   typedef _Rb_tree_alloc_base<_T, _Alloc,
  437.                               _Alloc_traits<_T, _Alloc>::_S_instanceless>
  438.           _Base;
  439.   typedef typename _Base::allocator_type allocator_type;
  440.   _Rb_tree_base(const allocator_type& __a) 
  441.     : _Base(__a) { _M_header = _M_get_node(); }
  442.   ~_Rb_tree_base() { _M_put_node(_M_header); }
  443. };
  444. #else /* __STL_USE_STD_ALLOCATORS */
  445. template <class _T, class _Alloc>
  446. struct _Rb_tree_base
  447. {
  448.   typedef _Alloc allocator_type;
  449.   allocator_type get_allocator() const { return allocator_type(); }
  450.   _Rb_tree_base(const allocator_type&) 
  451.     : _M_header(0) { _M_header = _M_get_node(); }
  452.   ~_Rb_tree_base() { _M_put_node(_M_header); }
  453. protected:
  454.   _Rb_tree_node<_T>* _M_header;
  455.   typedef simple_alloc<_Rb_tree_node<_T>, _Alloc> _Alloc_type;
  456.   _Rb_tree_node<_T>* _M_get_node()
  457.     { return _Alloc_type::allocate(1); }
  458.   void _M_put_node(_Rb_tree_node<_T>* __p)
  459.     { _Alloc_type::deallocate(__p, 1); }
  460. };
  461. #endif /* __STL_USE_STD_ALLOCATORS */
  462. template <class _Key, class _Value, class _KeyOfValue, class _Compare,
  463.           class _Alloc = __STL_DEFAULT_ALLOCATOR(_Value) >
  464. class _Rb_tree : protected _Rb_tree_base<_Value, _Alloc> {
  465.   typedef _Rb_tree_base<_Value, _Alloc> _Base;
  466. protected:
  467.   typedef _Rb_tree_node_base* _Base_ptr;
  468.   typedef _Rb_tree_node<_Value> _Rb_tree_node;
  469.   typedef _Rb_tree_Color_type _Color_type;
  470. public:
  471.   typedef _Key key_type;
  472.   typedef _Value value_type;
  473.   typedef value_type* pointer;
  474.   typedef const value_type* const_pointer;
  475.   typedef value_type& reference;
  476.   typedef const value_type& const_reference;
  477.   typedef _Rb_tree_node* _Link_type;
  478.   typedef size_t size_type;
  479.   typedef ptrdiff_t difference_type;
  480.   typedef typename _Base::allocator_type allocator_type;
  481.   allocator_type get_allocator() const { return _Base::get_allocator(); }
  482. protected:
  483. #ifdef __STL_USE_NAMESPACES
  484.   using _Base::_M_get_node;
  485.   using _Base::_M_put_node;
  486.   using _Base::_M_header;
  487. #endif /* __STL_USE_NAMESPACES */
  488. protected:
  489.   _Link_type _M_create_node(const value_type& __x)
  490.   {
  491.     _Link_type __tmp = _M_get_node();
  492.     __STL_TRY {
  493.       construct(&__tmp->_M_value_field, __x);
  494.     }
  495.     __STL_UNWIND(_M_put_node(__tmp));
  496.     return __tmp;
  497.   }
  498.   _Link_type _M_clone_node(_Link_type __x)
  499.   {
  500.     _Link_type __tmp = _M_create_node(__x->_M_value_field);
  501.     __tmp->_M_color = __x->_M_color;
  502.     __tmp->_M_left = 0;
  503.     __tmp->_M_right = 0;
  504.     return __tmp;
  505.   }
  506.   void destroy_node(_Link_type __p)
  507.   {
  508.     destroy(&__p->_M_value_field);
  509.     _M_put_node(__p);
  510.   }
  511. protected:
  512.   size_type _M_node_count; // keeps track of size of tree
  513.   _Compare _M_key_compare;
  514.   _Link_type& _M_root() const 
  515.     { return (_Link_type&) _M_header->_M_parent; }
  516.   _Link_type& _M_leftmost() const 
  517.     { return (_Link_type&) _M_header->_M_left; }
  518.   _Link_type& _M_rightmost() const 
  519.     { return (_Link_type&) _M_header->_M_right; }
  520.   static _Link_type& _S_left(_Link_type __x)
  521.     { return (_Link_type&)(__x->_M_left); }
  522.   static _Link_type& _S_right(_Link_type __x)
  523.     { return (_Link_type&)(__x->_M_right); }
  524.   static _Link_type& _S_parent(_Link_type __x)
  525.     { return (_Link_type&)(__x->_M_parent); }
  526.   static reference _S_value(_Link_type __x)
  527.     { return __x->_M_value_field; }
  528.   static const _Key& _S_key(_Link_type __x)
  529.     { return _KeyOfValue()(_S_value(__x)); }
  530.   static _Color_type& _S_color(_Link_type __x)
  531.     { return (_Color_type&)(__x->_M_color); }
  532.   static _Link_type& _S_left(_Base_ptr __x)
  533.     { return (_Link_type&)(__x->_M_left); }
  534.   static _Link_type& _S_right(_Base_ptr __x)
  535.     { return (_Link_type&)(__x->_M_right); }
  536.   static _Link_type& _S_parent(_Base_ptr __x)
  537.     { return (_Link_type&)(__x->_M_parent); }
  538.   static reference _S_value(_Base_ptr __x)
  539.     { return ((_Link_type)__x)->_M_value_field; }
  540.   static const _Key& _S_key(_Base_ptr __x)
  541.     { return _KeyOfValue()(_S_value(_Link_type(__x)));} 
  542.   static _Color_type& _S_color(_Base_ptr __x)
  543.     { return (_Color_type&)(_Link_type(__x)->_M_color); }
  544.   static _Link_type _S_minimum(_Link_type __x) 
  545.     { return (_Link_type)  _Rb_tree_node_base::_S_minimum(__x); }
  546.   static _Link_type _S_maximum(_Link_type __x)
  547.     { return (_Link_type) _Rb_tree_node_base::_S_maximum(__x); }
  548. public:
  549.   typedef _Rb_tree_iterator<value_type, reference, pointer> iterator;
  550.   typedef _Rb_tree_iterator<value_type, const_reference, const_pointer> 
  551.           const_iterator;
  552. #ifdef __STL_CLASS_PARTIAL_SPECIALIZATION
  553.   typedef reverse_iterator<const_iterator> const_reverse_iterator;
  554.   typedef reverse_iterator<iterator> reverse_iterator;
  555. #else /* __STL_CLASS_PARTIAL_SPECIALIZATION */
  556.   typedef reverse_bidirectional_iterator<iterator, value_type, reference,
  557.                                          difference_type>
  558.           reverse_iterator; 
  559.   typedef reverse_bidirectional_iterator<const_iterator, value_type,
  560.                                          const_reference, difference_type>
  561.           const_reverse_iterator;
  562. #endif /* __STL_CLASS_PARTIAL_SPECIALIZATION */ 
  563. private:
  564.   iterator _M_insert(_Base_ptr __x, _Base_ptr __y, const value_type& __v);
  565.   _Link_type _M_copy(_Link_type __x, _Link_type __p);
  566.   void _M_erase(_Link_type __x);
  567. public:
  568.                                 // allocation/deallocation
  569.   _Rb_tree()
  570.     : _Base(allocator_type()), _M_node_count(0), _M_key_compare()
  571.     { _M_empty_initialize(); }
  572.   _Rb_tree(const _Compare& __comp)
  573.     : _Base(allocator_type()), _M_node_count(0), _M_key_compare(__comp) 
  574.     { _M_empty_initialize(); }
  575.   _Rb_tree(const _Compare& __comp, const allocator_type& __a)
  576.     : _Base(__a), _M_node_count(0), _M_key_compare(__comp) 
  577.     { _M_empty_initialize(); }
  578.   _Rb_tree(const _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>& __x) 
  579.     : _Base(__x.get_allocator()),
  580.       _M_node_count(0), _M_key_compare(__x._M_key_compare)
  581.   { 
  582.     if (__x._M_root() == 0)
  583.       _M_empty_initialize();
  584.     else {
  585.       _S_color(_M_header) = _S_rb_tree_red;
  586.       _M_root() = _M_copy(__x._M_root(), _M_header);
  587.       _M_leftmost() = _S_minimum(_M_root());
  588.       _M_rightmost() = _S_maximum(_M_root());
  589.     }
  590.     _M_node_count = __x._M_node_count;
  591.   }
  592.   ~_Rb_tree() { clear(); }
  593.   _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>& 
  594.   operator=(const _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>& __x);
  595. private:
  596.   void _M_empty_initialize() {
  597.     _S_color(_M_header) = _S_rb_tree_red; // used to distinguish header from 
  598.                                           // __root, in iterator.operator++
  599.     _M_root() = 0;
  600.     _M_leftmost() = _M_header;
  601.     _M_rightmost() = _M_header;
  602.   }
  603. public:    
  604.                                 // accessors:
  605.   _Compare key_comp() const { return _M_key_compare; }
  606.   iterator begin() { return _M_leftmost(); }
  607.   const_iterator begin() const { return _M_leftmost(); }
  608.   iterator end() { return _M_header; }
  609.   const_iterator end() const { return _M_header; }
  610.   reverse_iterator rbegin() { return reverse_iterator(end()); }
  611.   const_reverse_iterator rbegin() const { 
  612.     return const_reverse_iterator(end()); 
  613.   }
  614.   reverse_iterator rend() { return reverse_iterator(begin()); }
  615.   const_reverse_iterator rend() const { 
  616.     return const_reverse_iterator(begin());
  617.   } 
  618.   bool empty() const { return _M_node_count == 0; }
  619.   size_type size() const { return _M_node_count; }
  620.   size_type max_size() const { return size_type(-1); }
  621.   void swap(_Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>& __t) {
  622.     __STD::swap(_M_header, __t._M_header);
  623.     __STD::swap(_M_node_count, __t._M_node_count);
  624.     __STD::swap(_M_key_compare, __t._M_key_compare);
  625.   }
  626.     
  627. public:
  628.                                 // insert/erase
  629.   pair<iterator,bool> insert_unique(const value_type& __x);
  630.   iterator insert_equal(const value_type& __x);
  631.   iterator insert_unique(iterator __position, const value_type& __x);
  632.   iterator insert_equal(iterator __position, const value_type& __x);
  633. #ifdef __STL_MEMBER_TEMPLATES  
  634.   template <class _InputIterator>
  635.   void insert_unique(_InputIterator __first, _InputIterator __last);
  636.   template <class _InputIterator>
  637.   void insert_equal(_InputIterator __first, _InputIterator __last);
  638. #else /* __STL_MEMBER_TEMPLATES */
  639.   void insert_unique(const_iterator __first, const_iterator __last);
  640.   void insert_unique(const value_type* __first, const value_type* __last);
  641.   void insert_equal(const_iterator __first, const_iterator __last);
  642.   void insert_equal(const value_type* __first, const value_type* __last);
  643. #endif /* __STL_MEMBER_TEMPLATES */
  644.   void erase(iterator __position);
  645.   size_type erase(const key_type& __x);
  646.   void erase(iterator __first, iterator __last);
  647.   void erase(const key_type* __first, const key_type* __last);
  648.   void clear() {
  649.     if (_M_node_count != 0) {
  650.       _M_erase(_M_root());
  651.       _M_leftmost() = _M_header;
  652.       _M_root() = 0;
  653.       _M_rightmost() = _M_header;
  654.       _M_node_count = 0;
  655.     }
  656.   }      
  657. public:
  658.                                 // set operations:
  659.   iterator find(const key_type& __x);
  660.   const_iterator find(const key_type& __x) const;
  661.   size_type count(const key_type& __x) const;
  662.   iterator lower_bound(const key_type& __x);
  663.   const_iterator lower_bound(const key_type& __x) const;
  664.   iterator upper_bound(const key_type& __x);
  665.   const_iterator upper_bound(const key_type& __x) const;
  666.   pair<iterator,iterator> equal_range(const key_type& __x);
  667.   pair<const_iterator, const_iterator> equal_range(const key_type& __x) const;
  668. public:
  669.                                 // Debugging.
  670.   bool __rb_verify() const;
  671. };
  672. template <class _Key, class _Value, class _KeyOfValue, 
  673.           class _Compare, class _Alloc>
  674. inline bool 
  675. operator==(const _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>& __x, 
  676.            const _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>& __y)
  677. {
  678.   return __x.size() == __y.size() &&
  679.          equal(__x.begin(), __x.end(), __y.begin());
  680. }
  681. template <class _Key, class _Value, class _KeyOfValue, 
  682.           class _Compare, class _Alloc>
  683. inline bool 
  684. operator<(const _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>& __x, 
  685.           const _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>& __y)
  686. {
  687.   return lexicographical_compare(__x.begin(), __x.end(), 
  688.                                  __y.begin(), __y.end());
  689. }
  690. #ifdef __STL_FUNCTION_TMPL_PARTIAL_ORDER
  691. template <class _Key, class _Value, class _KeyOfValue, 
  692.           class _Compare, class _Alloc>
  693. inline void 
  694. swap(_Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>& __x, 
  695.      _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>& __y)
  696. {
  697.   __x.swap(__y);
  698. }
  699. #endif /* __STL_FUNCTION_TMPL_PARTIAL_ORDER */
  700. template <class _Key, class _Value, class _KeyOfValue, 
  701.           class _Compare, class _Alloc>
  702. _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>& 
  703. _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>
  704.   ::operator=(const _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>& __x)
  705. {
  706.   if (this != &__x) {
  707.                                 // Note that _Key may be a constant type.
  708.     clear();
  709.     _M_node_count = 0;
  710.     _M_key_compare = __x._M_key_compare;        
  711.     if (__x._M_root() == 0) {
  712.       _M_root() = 0;
  713.       _M_leftmost() = _M_header;
  714.       _M_rightmost() = _M_header;
  715.     }
  716.     else {
  717.       _M_root() = _M_copy(__x._M_root(), _M_header);
  718.       _M_leftmost() = _S_minimum(_M_root());
  719.       _M_rightmost() = _S_maximum(_M_root());
  720.       _M_node_count = __x._M_node_count;
  721.     }
  722.   }
  723.   return *this;
  724. }
  725. template <class _Key, class _Value, class _KeyOfValue, 
  726.           class _Compare, class _Alloc>
  727. typename _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>::iterator
  728. _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>
  729.   ::_M_insert(_Base_ptr __x_, _Base_ptr __y_, const _Value& __v)
  730. {
  731.   _Link_type __x = (_Link_type) __x_;
  732.   _Link_type __y = (_Link_type) __y_;
  733.   _Link_type __z;
  734.   if (__y == _M_header || __x != 0 || 
  735.       _M_key_compare(_KeyOfValue()(__v), _S_key(__y))) {
  736.     __z = _M_create_node(__v);
  737.     _S_left(__y) = __z;               // also makes _M_leftmost() = __z 
  738.                                       //    when __y == _M_header
  739.     if (__y == _M_header) {
  740.       _M_root() = __z;
  741.       _M_rightmost() = __z;
  742.     }
  743.     else if (__y == _M_leftmost())
  744.       _M_leftmost() = __z;   // maintain _M_leftmost() pointing to min node
  745.   }
  746.   else {
  747.     __z = _M_create_node(__v);
  748.     _S_right(__y) = __z;
  749.     if (__y == _M_rightmost())
  750.       _M_rightmost() = __z;  // maintain _M_rightmost() pointing to max node
  751.   }
  752.   _S_parent(__z) = __y;
  753.   _S_left(__z) = 0;
  754.   _S_right(__z) = 0;
  755.   _Rb_tree_rebalance(__z, _M_header->_M_parent);
  756.   ++_M_node_count;
  757.   return iterator(__z);
  758. }
  759. template <class _Key, class _Value, class _KeyOfValue, 
  760.           class _Compare, class _Alloc>
  761. typename _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>::iterator
  762. _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>
  763.   ::insert_equal(const _Value& __v)
  764. {
  765.   _Link_type __y = _M_header;
  766.   _Link_type __x = _M_root();
  767.   while (__x != 0) {
  768.     __y = __x;
  769.     __x = _M_key_compare(_KeyOfValue()(__v), _S_key(__x)) ? 
  770.             _S_left(__x) : _S_right(__x);
  771.   }
  772.   return _M_insert(__x, __y, __v);
  773. }
  774. template <class _Key, class _Value, class _KeyOfValue, 
  775.           class _Compare, class _Alloc>
  776. pair<typename _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>::iterator, 
  777.      bool>
  778. _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>
  779.   ::insert_unique(const _Value& __v)
  780. {
  781.   _Link_type __y = _M_header;
  782.   _Link_type __x = _M_root();
  783.   bool __comp = true;
  784.   while (__x != 0) {
  785.     __y = __x;
  786.     __comp = _M_key_compare(_KeyOfValue()(__v), _S_key(__x));
  787.     __x = __comp ? _S_left(__x) : _S_right(__x);
  788.   }
  789.   iterator __j = iterator(__y);   
  790.   if (__comp)
  791.     if (__j == begin())     
  792.       return pair<iterator,bool>(_M_insert(__x, __y, __v), true);
  793.     else
  794.       --__j;
  795.   if (_M_key_compare(_S_key(__j._M_node), _KeyOfValue()(__v)))
  796.     return pair<iterator,bool>(_M_insert(__x, __y, __v), true);
  797.   return pair<iterator,bool>(__j, false);
  798. }
  799. template <class _Key, class _Val, class _KeyOfValue, 
  800.           class _Compare, class _Alloc>
  801. typename _Rb_tree<_Key, _Val, _KeyOfValue, _Compare, _Alloc>::iterator 
  802. _Rb_tree<_Key, _Val, _KeyOfValue, _Compare, _Alloc>
  803.   ::insert_unique(iterator __position, const _Val& __v)
  804. {
  805.   if (__position._M_node == _M_header->_M_left) { // begin()
  806.     if (size() > 0 && 
  807.         _M_key_compare(_KeyOfValue()(__v), _S_key(__position._M_node)))
  808.       return _M_insert(__position._M_node, __position._M_node, __v);
  809.     // first argument just needs to be non-null 
  810.     else
  811.       return insert_unique(__v).first;
  812.   } else if (__position._M_node == _M_header) { // end()
  813.     if (_M_key_compare(_S_key(_M_rightmost()), _KeyOfValue()(__v)))
  814.       return _M_insert(0, _M_rightmost(), __v);
  815.     else
  816.       return insert_unique(__v).first;
  817.   } else {
  818.     iterator __before = __position;
  819.     --__before;
  820.     if (_M_key_compare(_S_key(__before._M_node), _KeyOfValue()(__v)) 
  821.         && _M_key_compare(_KeyOfValue()(__v), _S_key(__position._M_node))) {
  822.       if (_S_right(__before._M_node) == 0)
  823.         return _M_insert(0, __before._M_node, __v); 
  824.       else
  825.         return _M_insert(__position._M_node, __position._M_node, __v);
  826.     // first argument just needs to be non-null 
  827.     } else
  828.       return insert_unique(__v).first;
  829.   }
  830. }
  831. template <class _Key, class _Val, class _KeyOfValue, 
  832.           class _Compare, class _Alloc>
  833. typename _Rb_tree<_Key,_Val,_KeyOfValue,_Compare,_Alloc>::iterator 
  834. _Rb_tree<_Key,_Val,_KeyOfValue,_Compare,_Alloc>
  835.   ::insert_equal(iterator __position, const _Val& __v)
  836. {
  837.   if (__position._M_node == _M_header->_M_left) { // begin()
  838.     if (size() > 0 && 
  839.         _M_key_compare(_KeyOfValue()(__v), _S_key(__position._M_node)))
  840.       return _M_insert(__position._M_node, __position._M_node, __v);
  841.     // first argument just needs to be non-null 
  842.     else
  843.       return insert_equal(__v);
  844.   } else if (__position._M_node == _M_header) {// end()
  845.     if (!_M_key_compare(_KeyOfValue()(__v), _S_key(_M_rightmost())))
  846.       return _M_insert(0, _M_rightmost(), __v);
  847.     else
  848.       return insert_equal(__v);
  849.   } else {
  850.     iterator __before = __position;
  851.     --__before;
  852.     if (!_M_key_compare(_KeyOfValue()(__v), _S_key(__before._M_node))
  853.         && !_M_key_compare(_S_key(__position._M_node), _KeyOfValue()(__v))) {
  854.       if (_S_right(__before._M_node) == 0)
  855.         return _M_insert(0, __before._M_node, __v); 
  856.       else
  857.         return _M_insert(__position._M_node, __position._M_node, __v);
  858.     // first argument just needs to be non-null 
  859.     } else
  860.       return insert_equal(__v);
  861.   }
  862. }
  863. #ifdef __STL_MEMBER_TEMPLATES  
  864. template <class _Key, class _Val, class _KoV, class _Cmp, class _Alloc>
  865.   template<class _II>
  866. void _Rb_tree<_Key,_Val,_KoV,_Cmp,_Alloc>
  867.   ::insert_equal(_II __first, _II __last)
  868. {
  869.   for ( ; __first != __last; ++__first)
  870.     insert_equal(*__first);
  871. }
  872. template <class _Key, class _Val, class _KoV, class _Cmp, class _Alloc> 
  873.   template<class _II>
  874. void _Rb_tree<_Key,_Val,_KoV,_Cmp,_Alloc>
  875.   ::insert_unique(_II __first, _II __last) {
  876.   for ( ; __first != __last; ++__first)
  877.     insert_unique(*__first);
  878. }
  879. #else /* __STL_MEMBER_TEMPLATES */
  880. template <class _Key, class _Val, class _KoV, class _Cmp, class _Alloc>
  881. void
  882. _Rb_tree<_Key,_Val,_KoV,_Cmp,_Alloc>
  883.   ::insert_equal(const _Val* __first, const _Val* __last)
  884. {
  885.   for ( ; __first != __last; ++__first)
  886.     insert_equal(*__first);
  887. }
  888. template <class _Key, class _Val, class _KoV, class _Cmp, class _Alloc>
  889. void
  890. _Rb_tree<_Key,_Val,_KoV,_Cmp,_Alloc>
  891.   ::insert_equal(const_iterator __first, const_iterator __last)
  892. {
  893.   for ( ; __first != __last; ++__first)
  894.     insert_equal(*__first);
  895. }
  896. template <class _Key, class _Val, class _KoV, class _Cmp, class _Alloc>
  897. void 
  898. _Rb_tree<_Key,_Val,_KoV,_Cmp,_Alloc>
  899.   ::insert_unique(const _Val* __first, const _Val* __last)
  900. {
  901.   for ( ; __first != __last; ++__first)
  902.     insert_unique(*__first);
  903. }
  904. template <class _Key, class _Val, class _KoV, class _Cmp, class _Alloc>
  905. void _Rb_tree<_Key,_Val,_KoV,_Cmp,_Alloc>
  906.   ::insert_unique(const_iterator __first, const_iterator __last)
  907. {
  908.   for ( ; __first != __last; ++__first)
  909.     insert_unique(*__first);
  910. }
  911. #endif /* __STL_MEMBER_TEMPLATES */
  912.          
  913. template <class _Key, class _Value, class _KeyOfValue, 
  914.           class _Compare, class _Alloc>
  915. inline void _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>
  916.   ::erase(iterator __position)
  917. {
  918.   _Link_type __y = 
  919.     (_Link_type) _Rb_tree_rebalance_for_erase(__position._M_node,
  920.                                               _M_header->_M_parent,
  921.                                               _M_header->_M_left,
  922.                                               _M_header->_M_right);
  923.   destroy_node(__y);
  924.   --_M_node_count;
  925. }
  926. template <class _Key, class _Value, class _KeyOfValue, 
  927.           class _Compare, class _Alloc>
  928. typename _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>::size_type 
  929. _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>::erase(const _Key& __x)
  930. {
  931.   pair<iterator,iterator> __p = equal_range(__x);
  932.   size_type __n = 0;
  933.   distance(__p.first, __p.second, __n);
  934.   erase(__p.first, __p.second);
  935.   return __n;
  936. }
  937. template <class _Key, class _Val, class _KoV, class _Compare, class _Alloc>
  938. typename _Rb_tree<_Key, _Val, _KoV, _Compare, _Alloc>::_Link_type 
  939. _Rb_tree<_Key,_Val,_KoV,_Compare,_Alloc>
  940.   ::_M_copy(_Link_type __x, _Link_type __p)
  941. {
  942.                         // structural copy.  __x and __p must be non-null.
  943.   _Link_type __top = _M_clone_node(__x);
  944.   __top->_M_parent = __p;
  945.  
  946.   __STL_TRY {
  947.     if (__x->_M_right)
  948.       __top->_M_right = _M_copy(_S_right(__x), __top);
  949.     __p = __top;
  950.     __x = _S_left(__x);
  951.     while (__x != 0) {
  952.       _Link_type __y = _M_clone_node(__x);
  953.       __p->_M_left = __y;
  954.       __y->_M_parent = __p;
  955.       if (__x->_M_right)
  956.         __y->_M_right = _M_copy(_S_right(__x), __y);
  957.       __p = __y;
  958.       __x = _S_left(__x);
  959.     }
  960.   }
  961.   __STL_UNWIND(_M_erase(__top));
  962.   return __top;
  963. }
  964. template <class _Key, class _Value, class _KeyOfValue, 
  965.           class _Compare, class _Alloc>
  966. void _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>
  967.   ::_M_erase(_Link_type __x)
  968. {
  969.                                 // erase without rebalancing
  970.   while (__x != 0) {
  971.     _M_erase(_S_right(__x));
  972.     _Link_type __y = _S_left(__x);
  973.     destroy_node(__x);
  974.     __x = __y;
  975.   }
  976. }
  977. template <class _Key, class _Value, class _KeyOfValue, 
  978.           class _Compare, class _Alloc>
  979. void _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>
  980.   ::erase(iterator __first, iterator __last)
  981. {
  982.   if (__first == begin() && __last == end())
  983.     clear();
  984.   else
  985.     while (__first != __last) erase(__first++);
  986. }
  987. template <class _Key, class _Value, class _KeyOfValue, 
  988.           class _Compare, class _Alloc>
  989. void _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>
  990.   ::erase(const _Key* __first, const _Key* __last) 
  991. {
  992.   while (__first != __last) erase(*__first++);
  993. }
  994. template <class _Key, class _Value, class _KeyOfValue, 
  995.           class _Compare, class _Alloc>
  996. typename _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>::iterator 
  997. _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>::find(const _Key& __k)
  998. {
  999.   _Link_type __y = _M_header;      // Last node which is not less than __k. 
  1000.   _Link_type __x = _M_root();      // Current node. 
  1001.   while (__x != 0) 
  1002.     if (!_M_key_compare(_S_key(__x), __k))
  1003.       __y = __x, __x = _S_left(__x);
  1004.     else
  1005.       __x = _S_right(__x);
  1006.   iterator __j = iterator(__y);   
  1007.   return (__j == end() || _M_key_compare(__k, _S_key(__j._M_node))) ? 
  1008.      end() : __j;
  1009. }
  1010. template <class _Key, class _Value, class _KeyOfValue, 
  1011.           class _Compare, class _Alloc>
  1012. typename _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>::const_iterator 
  1013. _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>::find(const _Key& __k) const
  1014. {
  1015.   _Link_type __y = _M_header; /* Last node which is not less than __k. */
  1016.   _Link_type __x = _M_root(); /* Current node. */
  1017.   while (__x != 0) {
  1018.     if (!_M_key_compare(_S_key(__x), __k))
  1019.       __y = __x, __x = _S_left(__x);
  1020.     else
  1021.       __x = _S_right(__x);
  1022.   }
  1023.   const_iterator __j = const_iterator(__y);   
  1024.   return (__j == end() || _M_key_compare(__k, _S_key(__j._M_node))) ?
  1025.     end() : __j;
  1026. }
  1027. template <class _Key, class _Value, class _KeyOfValue, 
  1028.           class _Compare, class _Alloc>
  1029. typename _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>::size_type 
  1030. _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>
  1031.   ::count(const _Key& __k) const
  1032. {
  1033.   pair<const_iterator, const_iterator> __p = equal_range(__k);
  1034.   size_type __n = 0;
  1035.   distance(__p.first, __p.second, __n);
  1036.   return __n;
  1037. }
  1038. template <class _Key, class _Value, class _KeyOfValue, 
  1039.           class _Compare, class _Alloc>
  1040. typename _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>::iterator 
  1041. _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>
  1042.   ::lower_bound(const _Key& __k)
  1043. {
  1044.   _Link_type __y = _M_header; /* Last node which is not less than __k. */
  1045.   _Link_type __x = _M_root(); /* Current node. */
  1046.   while (__x != 0) 
  1047.     if (!_M_key_compare(_S_key(__x), __k))
  1048.       __y = __x, __x = _S_left(__x);
  1049.     else
  1050.       __x = _S_right(__x);
  1051.   return iterator(__y);
  1052. }
  1053. template <class _Key, class _Value, class _KeyOfValue, 
  1054.           class _Compare, class _Alloc>
  1055. typename _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>::const_iterator 
  1056. _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>
  1057.   ::lower_bound(const _Key& __k) const
  1058. {
  1059.   _Link_type __y = _M_header; /* Last node which is not less than __k. */
  1060.   _Link_type __x = _M_root(); /* Current node. */
  1061.   while (__x != 0) 
  1062.     if (!_M_key_compare(_S_key(__x), __k))
  1063.       __y = __x, __x = _S_left(__x);
  1064.     else
  1065.       __x = _S_right(__x);
  1066.   return const_iterator(__y);
  1067. }
  1068. template <class _Key, class _Value, class _KeyOfValue, 
  1069.           class _Compare, class _Alloc>
  1070. typename _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>::iterator 
  1071. _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>
  1072.   ::upper_bound(const _Key& __k)
  1073. {
  1074.   _Link_type __y = _M_header; /* Last node which is greater than __k. */
  1075.   _Link_type __x = _M_root(); /* Current node. */
  1076.    while (__x != 0) 
  1077.      if (_M_key_compare(__k, _S_key(__x)))
  1078.        __y = __x, __x = _S_left(__x);
  1079.      else
  1080.        __x = _S_right(__x);
  1081.    return iterator(__y);
  1082. }
  1083. template <class _Key, class _Value, class _KeyOfValue, 
  1084.           class _Compare, class _Alloc>
  1085. typename _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>::const_iterator 
  1086. _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>
  1087.   ::upper_bound(const _Key& __k) const
  1088. {
  1089.   _Link_type __y = _M_header; /* Last node which is greater than __k. */
  1090.   _Link_type __x = _M_root(); /* Current node. */
  1091.    while (__x != 0) 
  1092.      if (_M_key_compare(__k, _S_key(__x)))
  1093.        __y = __x, __x = _S_left(__x);
  1094.      else
  1095.        __x = _S_right(__x);
  1096.    return const_iterator(__y);
  1097. }
  1098. template <class _Key, class _Value, class _KeyOfValue, 
  1099.           class _Compare, class _Alloc>
  1100. inline 
  1101. pair<typename _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>::iterator,
  1102.      typename _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>::iterator>
  1103. _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>
  1104.   ::equal_range(const _Key& __k)
  1105. {
  1106.   return pair<iterator, iterator>(lower_bound(__k), upper_bound(__k));
  1107. }
  1108. template <class _Key, class _Value, class _KoV, class _Compare, class _Alloc>
  1109. inline 
  1110. pair<typename _Rb_tree<_Key, _Value, _KoV, _Compare, _Alloc>::const_iterator,
  1111.      typename _Rb_tree<_Key, _Value, _KoV, _Compare, _Alloc>::const_iterator>
  1112. _Rb_tree<_Key, _Value, _KoV, _Compare, _Alloc>
  1113.   ::equal_range(const _Key& __k) const
  1114. {
  1115.   return pair<const_iterator,const_iterator>(lower_bound(__k),
  1116.                                              upper_bound(__k));
  1117. }
  1118. inline int 
  1119. __black_count(_Rb_tree_node_base* __node, _Rb_tree_node_base* __root)
  1120. {
  1121.   if (__node == 0)
  1122.     return 0;
  1123.   else {
  1124.     int __bc = __node->_M_color == _S_rb_tree_black ? 1 : 0;
  1125.     if (__node == __root)
  1126.       return __bc;
  1127.     else
  1128.       return __bc + __black_count(__node->_M_parent, __root);
  1129.   }
  1130. }
  1131. template <class _Key, class _Value, class _KeyOfValue, 
  1132.           class _Compare, class _Alloc>
  1133. bool _Rb_tree<_Key,_Value,_KeyOfValue,_Compare,_Alloc>::__rb_verify() const
  1134. {
  1135.   if (_M_node_count == 0 || begin() == end())
  1136.     return _M_node_count == 0 && begin() == end() &&
  1137.       _M_header->_M_left == _M_header && _M_header->_M_right == _M_header;
  1138.   
  1139.   int __len = __black_count(_M_leftmost(), _M_root());
  1140.   for (const_iterator __it = begin(); __it != end(); ++__it) {
  1141.     _Link_type __x = (_Link_type) __it._M_node;
  1142.     _Link_type __L = _S_left(__x);
  1143.     _Link_type __R = _S_right(__x);
  1144.     if (__x->_M_color == _S_rb_tree_red)
  1145.       if ((__L && __L->_M_color == _S_rb_tree_red) ||
  1146.           (__R && __R->_M_color == _S_rb_tree_red))
  1147.         return false;
  1148.     if (__L && _M_key_compare(_S_key(__x), _S_key(__L)))
  1149.       return false;
  1150.     if (__R && _M_key_compare(_S_key(__R), _S_key(__x)))
  1151.       return false;
  1152.     if (!__L && !__R && __black_count(__x, _M_root()) != __len)
  1153.       return false;
  1154.   }
  1155.   if (_M_leftmost() != _Rb_tree_node_base::_S_minimum(_M_root()))
  1156.     return false;
  1157.   if (_M_rightmost() != _Rb_tree_node_base::_S_maximum(_M_root()))
  1158.     return false;
  1159.   return true;
  1160. }
  1161. // Class rb_tree is not part of the C++ standard.  It is provided for
  1162. // compatibility with the HP STL.
  1163. template <class _Key, class _Value, class _KeyOfValue, class _Compare,
  1164.           class _Alloc = __STL_DEFAULT_ALLOCATOR(_Value) >
  1165. struct rb_tree : public _Rb_tree<_Key, _Value, _KeyOfValue, _Compare, _Alloc>
  1166. {
  1167.   typedef _Rb_tree<_Key, _Value, _KeyOfValue, _Compare, _Alloc> _Base;
  1168.   typedef _Base::allocator_type allocator_type;
  1169.   rb_tree(const _Compare& __comp = _Compare(),
  1170.           const allocator_type& __a = allocator_type())
  1171.     : _Base(__comp, __a) {}
  1172.   
  1173.   ~rb_tree() {}
  1174. };
  1175. #if defined(__sgi) && !defined(__GNUC__) && (_MIPS_SIM != _MIPS_SIM_ABI32)
  1176. #pragma reset woff 1375
  1177. #endif
  1178. __STL_END_NAMESPACE 
  1179. #endif /* __SGI_STL_INTERNAL_TREE_H */
  1180. // Local Variables:
  1181. // mode:C++
  1182. // End: