inplace_merge.h 10.0 KB

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