inplace_merge.h 9.3 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231
  1. //===----------------------------------------------------------------------===//
  2. //
  3. // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
  4. // See https://llvm.org/LICENSE.txt for license information.
  5. // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
  6. //
  7. //===----------------------------------------------------------------------===//
  8. #ifndef _LIBCPP___ALGORITHM_INPLACE_MERGE_H
  9. #define _LIBCPP___ALGORITHM_INPLACE_MERGE_H
  10. #include <__algorithm/comp.h>
  11. #include <__algorithm/comp_ref_type.h>
  12. #include <__algorithm/lower_bound.h>
  13. #include <__algorithm/min.h>
  14. #include <__algorithm/move.h>
  15. #include <__algorithm/rotate.h>
  16. #include <__algorithm/upper_bound.h>
  17. #include <__config>
  18. #include <__iterator/iterator_traits.h>
  19. #include <__utility/swap.h>
  20. #include <memory>
  21. #if !defined(_LIBCPP_HAS_NO_PRAGMA_SYSTEM_HEADER)
  22. # pragma GCC system_header
  23. #endif
  24. _LIBCPP_PUSH_MACROS
  25. #include <__undef_macros>
  26. _LIBCPP_BEGIN_NAMESPACE_STD
  27. template <class _Predicate>
  28. class __invert // invert the sense of a comparison
  29. {
  30. private:
  31. _Predicate __p_;
  32. public:
  33. _LIBCPP_INLINE_VISIBILITY __invert() {}
  34. _LIBCPP_INLINE_VISIBILITY
  35. explicit __invert(_Predicate __p) : __p_(__p) {}
  36. template <class _T1>
  37. _LIBCPP_INLINE_VISIBILITY
  38. bool operator()(const _T1& __x) {return !__p_(__x);}
  39. template <class _T1, class _T2>
  40. _LIBCPP_INLINE_VISIBILITY
  41. bool operator()(const _T1& __x, const _T2& __y) {return __p_(__y, __x);}
  42. };
  43. template <class _Compare, class _InputIterator1, class _InputIterator2,
  44. class _OutputIterator>
  45. void __half_inplace_merge(_InputIterator1 __first1, _InputIterator1 __last1,
  46. _InputIterator2 __first2, _InputIterator2 __last2,
  47. _OutputIterator __result, _Compare __comp)
  48. {
  49. for (; __first1 != __last1; ++__result)
  50. {
  51. if (__first2 == __last2)
  52. {
  53. _VSTD::move(__first1, __last1, __result);
  54. return;
  55. }
  56. if (__comp(*__first2, *__first1))
  57. {
  58. *__result = _VSTD::move(*__first2);
  59. ++__first2;
  60. }
  61. else
  62. {
  63. *__result = _VSTD::move(*__first1);
  64. ++__first1;
  65. }
  66. }
  67. // __first2 through __last2 are already in the right spot.
  68. }
  69. template <class _Compare, class _BidirectionalIterator>
  70. void
  71. __buffered_inplace_merge(_BidirectionalIterator __first, _BidirectionalIterator __middle, _BidirectionalIterator __last,
  72. _Compare __comp, typename iterator_traits<_BidirectionalIterator>::difference_type __len1,
  73. typename iterator_traits<_BidirectionalIterator>::difference_type __len2,
  74. typename iterator_traits<_BidirectionalIterator>::value_type* __buff)
  75. {
  76. typedef typename iterator_traits<_BidirectionalIterator>::value_type value_type;
  77. __destruct_n __d(0);
  78. unique_ptr<value_type, __destruct_n&> __h2(__buff, __d);
  79. if (__len1 <= __len2)
  80. {
  81. value_type* __p = __buff;
  82. for (_BidirectionalIterator __i = __first; __i != __middle; __d.template __incr<value_type>(), (void) ++__i, (void) ++__p)
  83. ::new ((void*)__p) value_type(_VSTD::move(*__i));
  84. _VSTD::__half_inplace_merge<_Compare>(__buff, __p, __middle, __last, __first, __comp);
  85. }
  86. else
  87. {
  88. value_type* __p = __buff;
  89. for (_BidirectionalIterator __i = __middle; __i != __last; __d.template __incr<value_type>(), (void) ++__i, (void) ++__p)
  90. ::new ((void*)__p) value_type(_VSTD::move(*__i));
  91. typedef reverse_iterator<_BidirectionalIterator> _RBi;
  92. typedef reverse_iterator<value_type*> _Rv;
  93. typedef __invert<_Compare> _Inverted;
  94. _VSTD::__half_inplace_merge<_Inverted>(_Rv(__p), _Rv(__buff),
  95. _RBi(__middle), _RBi(__first),
  96. _RBi(__last), _Inverted(__comp));
  97. }
  98. }
  99. template <class _Compare, class _BidirectionalIterator>
  100. void
  101. __inplace_merge(_BidirectionalIterator __first, _BidirectionalIterator __middle, _BidirectionalIterator __last,
  102. _Compare __comp, typename iterator_traits<_BidirectionalIterator>::difference_type __len1,
  103. typename iterator_traits<_BidirectionalIterator>::difference_type __len2,
  104. typename iterator_traits<_BidirectionalIterator>::value_type* __buff, ptrdiff_t __buff_size)
  105. {
  106. typedef typename iterator_traits<_BidirectionalIterator>::difference_type difference_type;
  107. while (true)
  108. {
  109. // if __middle == __last, we're done
  110. if (__len2 == 0)
  111. return;
  112. if (__len1 <= __buff_size || __len2 <= __buff_size)
  113. return _VSTD::__buffered_inplace_merge<_Compare>
  114. (__first, __middle, __last, __comp, __len1, __len2, __buff);
  115. // shrink [__first, __middle) as much as possible (with no moves), returning if it shrinks to 0
  116. for (; true; ++__first, (void) --__len1)
  117. {
  118. if (__len1 == 0)
  119. return;
  120. if (__comp(*__middle, *__first))
  121. break;
  122. }
  123. // __first < __middle < __last
  124. // *__first > *__middle
  125. // partition [__first, __m1) [__m1, __middle) [__middle, __m2) [__m2, __last) such that
  126. // all elements in:
  127. // [__first, __m1) <= [__middle, __m2)
  128. // [__middle, __m2) < [__m1, __middle)
  129. // [__m1, __middle) <= [__m2, __last)
  130. // and __m1 or __m2 is in the middle of its range
  131. _BidirectionalIterator __m1; // "median" of [__first, __middle)
  132. _BidirectionalIterator __m2; // "median" of [__middle, __last)
  133. difference_type __len11; // distance(__first, __m1)
  134. difference_type __len21; // distance(__middle, __m2)
  135. // binary search smaller range
  136. if (__len1 < __len2)
  137. { // __len >= 1, __len2 >= 2
  138. __len21 = __len2 / 2;
  139. __m2 = __middle;
  140. _VSTD::advance(__m2, __len21);
  141. __m1 = _VSTD::__upper_bound<_Compare>(__first, __middle, *__m2, __comp);
  142. __len11 = _VSTD::distance(__first, __m1);
  143. }
  144. else
  145. {
  146. if (__len1 == 1)
  147. { // __len1 >= __len2 && __len2 > 0, therefore __len2 == 1
  148. // It is known *__first > *__middle
  149. swap(*__first, *__middle);
  150. return;
  151. }
  152. // __len1 >= 2, __len2 >= 1
  153. __len11 = __len1 / 2;
  154. __m1 = __first;
  155. _VSTD::advance(__m1, __len11);
  156. __m2 = _VSTD::__lower_bound<_Compare>(__middle, __last, *__m1, __comp);
  157. __len21 = _VSTD::distance(__middle, __m2);
  158. }
  159. difference_type __len12 = __len1 - __len11; // distance(__m1, __middle)
  160. difference_type __len22 = __len2 - __len21; // distance(__m2, __last)
  161. // [__first, __m1) [__m1, __middle) [__middle, __m2) [__m2, __last)
  162. // swap middle two partitions
  163. __middle = _VSTD::rotate(__m1, __middle, __m2);
  164. // __len12 and __len21 now have swapped meanings
  165. // merge smaller range with recursive call and larger with tail recursion elimination
  166. if (__len11 + __len21 < __len12 + __len22)
  167. {
  168. _VSTD::__inplace_merge<_Compare>(__first, __m1, __middle, __comp, __len11, __len21, __buff, __buff_size);
  169. // _VSTD::__inplace_merge<_Compare>(__middle, __m2, __last, __comp, __len12, __len22, __buff, __buff_size);
  170. __first = __middle;
  171. __middle = __m2;
  172. __len1 = __len12;
  173. __len2 = __len22;
  174. }
  175. else
  176. {
  177. _VSTD::__inplace_merge<_Compare>(__middle, __m2, __last, __comp, __len12, __len22, __buff, __buff_size);
  178. // _VSTD::__inplace_merge<_Compare>(__first, __m1, __middle, __comp, __len11, __len21, __buff, __buff_size);
  179. __last = __middle;
  180. __middle = __m1;
  181. __len1 = __len11;
  182. __len2 = __len21;
  183. }
  184. }
  185. }
  186. template <class _BidirectionalIterator, class _Compare>
  187. inline _LIBCPP_INLINE_VISIBILITY
  188. void
  189. inplace_merge(_BidirectionalIterator __first, _BidirectionalIterator __middle, _BidirectionalIterator __last,
  190. _Compare __comp)
  191. {
  192. typedef typename iterator_traits<_BidirectionalIterator>::value_type value_type;
  193. typedef typename iterator_traits<_BidirectionalIterator>::difference_type difference_type;
  194. difference_type __len1 = _VSTD::distance(__first, __middle);
  195. difference_type __len2 = _VSTD::distance(__middle, __last);
  196. difference_type __buf_size = _VSTD::min(__len1, __len2);
  197. pair<value_type*, ptrdiff_t> __buf = _VSTD::get_temporary_buffer<value_type>(__buf_size);
  198. unique_ptr<value_type, __return_temporary_buffer> __h(__buf.first);
  199. typedef typename __comp_ref_type<_Compare>::type _Comp_ref;
  200. return _VSTD::__inplace_merge<_Comp_ref>(__first, __middle, __last, __comp, __len1, __len2,
  201. __buf.first, __buf.second);
  202. }
  203. template <class _BidirectionalIterator>
  204. inline _LIBCPP_INLINE_VISIBILITY
  205. void
  206. inplace_merge(_BidirectionalIterator __first, _BidirectionalIterator __middle, _BidirectionalIterator __last)
  207. {
  208. _VSTD::inplace_merge(__first, __middle, __last,
  209. __less<typename iterator_traits<_BidirectionalIterator>::value_type>());
  210. }
  211. _LIBCPP_END_NAMESPACE_STD
  212. _LIBCPP_POP_MACROS
  213. #endif // _LIBCPP___ALGORITHM_INPLACE_MERGE_H