git
/
gcc
mirror of git://gcc.gnu.org/git/gcc.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity
gcc/libstdc++-v3/include/bits/ranges_algo.h

5453 lines
172 KiB
C++
Raw Blame History

// Core algorithmic facilities -*- C++ -*-
// Copyright (C) 2020-2025 Free Software Foundation, Inc.
//
// This file is part of the GNU ISO C++ Library. This library is free
// software; you can redistribute it and/or modify it under the
// terms of the GNU General Public License as published by the
// Free Software Foundation; either version 3, or (at your option)
// any later version.
// This library is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU General Public License for more details.
// Under Section 7 of GPL version 3, you are granted additional
// permissions described in the GCC Runtime Library Exception, version
// 3.1, as published by the Free Software Foundation.
// You should have received a copy of the GNU General Public License and
// a copy of the GCC Runtime Library Exception along with this program;
// see the files COPYING3 and COPYING.RUNTIME respectively. If not, see
// <http://www.gnu.org/licenses/>.
/** @file bits/ranges_algo.h
* This is an internal header file, included by other library headers.
* Do not attempt to use it directly. @headername{algorithm}
*/
#ifndef _RANGES_ALGO_H
#define _RANGES_ALGO_H 1
#if __cplusplus > 201703L
#include <bit> // __bit_width
#if __cplusplus > 202002L
#include <optional>
#endif
#include <bits/ranges_algobase.h>
#include <bits/ranges_util.h>
#include <bits/uniform_int_dist.h> // concept uniform_random_bit_generator
#if __glibcxx_concepts
namespace std _GLIBCXX_VISIBILITY(default)
{
_GLIBCXX_BEGIN_NAMESPACE_VERSION
namespace ranges
{
namespace __detail
{
template<typename _Fp>
using __by_ref_or_value_fn
= __conditional_t<is_scalar_v<_Fp> || is_empty_v<_Fp>, _Fp, _Fp&>;
template<typename _Comp, typename _Proj>
struct _Comp_proj
{
[[no_unique_address]] __by_ref_or_value_fn<_Comp> _M_comp;
[[no_unique_address]] __by_ref_or_value_fn<_Proj> _M_proj;
constexpr
_Comp_proj(_Comp& __comp, _Proj& __proj)
: _M_comp(__comp), _M_proj(__proj)
{ }
template<typename _Tp, typename _Up>
constexpr bool
operator()(_Tp&& __x, _Up&& __y)
{
return std::__invoke(_M_comp,
std::__invoke(_M_proj, std::forward<_Tp>(__x)),
std::__invoke(_M_proj, std::forward<_Up>(__y)));
}
};
template<typename _Comp, typename _Proj>
constexpr _Comp_proj<_Comp, _Proj>
__make_comp_proj(_Comp& __comp, _Proj& __proj)
{ return {__comp, __proj}; }
template<typename _Pred, typename _Proj>
struct _Pred_proj
{
[[no_unique_address]] __by_ref_or_value_fn<_Pred> _M_pred;
[[no_unique_address]] __by_ref_or_value_fn<_Proj> _M_proj;
constexpr
_Pred_proj(_Pred& __pred, _Proj& __proj)
: _M_pred(__pred), _M_proj(__proj)
{ }
template<typename _Tp>
constexpr bool
operator()(_Tp&& __x)
{
return std::__invoke(_M_pred,
std::__invoke(_M_proj, std::forward<_Tp>(__x)));
}
};
template<typename _Pred, typename _Proj>
constexpr _Pred_proj<_Pred, _Proj>
__make_pred_proj(_Pred& __pred, _Proj& __proj)
{ return {__pred, __proj}; }
} // namespace __detail
struct __all_of_fn
{
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
indirect_unary_predicate<projected<_Iter, _Proj>> _Pred>
[[nodiscard]] constexpr bool
operator()(_Iter __first, _Sent __last,
_Pred __pred, _Proj __proj = {}) const
{
for (; __first != __last; ++__first)
if (!(bool)std::__invoke(__pred, std::__invoke(__proj, *__first)))
return false;
return true;
}
template<input_range _Range, typename _Proj = identity,
indirect_unary_predicate<projected<iterator_t<_Range>, _Proj>>
_Pred>
[[nodiscard]] constexpr bool
operator()(_Range&& __r, _Pred __pred, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__pred), std::move(__proj));
}
};
inline constexpr __all_of_fn all_of{};
struct __any_of_fn
{
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
indirect_unary_predicate<projected<_Iter, _Proj>> _Pred>
[[nodiscard]] constexpr bool
operator()(_Iter __first, _Sent __last,
_Pred __pred, _Proj __proj = {}) const
{
for (; __first != __last; ++__first)
if (std::__invoke(__pred, std::__invoke(__proj, *__first)))
return true;
return false;
}
template<input_range _Range, typename _Proj = identity,
indirect_unary_predicate<projected<iterator_t<_Range>, _Proj>>
_Pred>
[[nodiscard]] constexpr bool
operator()(_Range&& __r, _Pred __pred, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__pred), std::move(__proj));
}
};
inline constexpr __any_of_fn any_of{};
struct __none_of_fn
{
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
indirect_unary_predicate<projected<_Iter, _Proj>> _Pred>
[[nodiscard]] constexpr bool
operator()(_Iter __first, _Sent __last,
_Pred __pred, _Proj __proj = {}) const
{
for (; __first != __last; ++__first)
if (std::__invoke(__pred, std::__invoke(__proj, *__first)))
return false;
return true;
}
template<input_range _Range, typename _Proj = identity,
indirect_unary_predicate<projected<iterator_t<_Range>, _Proj>>
_Pred>
[[nodiscard]] constexpr bool
operator()(_Range&& __r, _Pred __pred, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__pred), std::move(__proj));
}
};
inline constexpr __none_of_fn none_of{};
template<typename _Iter, typename _Fp>
struct in_fun_result
{
[[no_unique_address]] _Iter in;
[[no_unique_address]] _Fp fun;
template<typename _Iter2, typename _F2p>
requires convertible_to<const _Iter&, _Iter2>
&& convertible_to<const _Fp&, _F2p>
constexpr
operator in_fun_result<_Iter2, _F2p>() const &
{ return {in, fun}; }
template<typename _Iter2, typename _F2p>
requires convertible_to<_Iter, _Iter2> && convertible_to<_Fp, _F2p>
constexpr
operator in_fun_result<_Iter2, _F2p>() &&
{ return {std::move(in), std::move(fun)}; }
};
template<typename _Iter, typename _Fp>
using for_each_result = in_fun_result<_Iter, _Fp>;
struct __for_each_fn
{
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
indirectly_unary_invocable<projected<_Iter, _Proj>> _Fun>
constexpr for_each_result<_Iter, _Fun>
operator()(_Iter __first, _Sent __last, _Fun __f, _Proj __proj = {}) const
{
for (; __first != __last; ++__first)
std::__invoke(__f, std::__invoke(__proj, *__first));
return { std::move(__first), std::move(__f) };
}
template<input_range _Range, typename _Proj = identity,
indirectly_unary_invocable<projected<iterator_t<_Range>, _Proj>>
_Fun>
constexpr for_each_result<borrowed_iterator_t<_Range>, _Fun>
operator()(_Range&& __r, _Fun __f, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__f), std::move(__proj));
}
};
inline constexpr __for_each_fn for_each{};
template<typename _Iter, typename _Fp>
using for_each_n_result = in_fun_result<_Iter, _Fp>;
struct __for_each_n_fn
{
template<input_iterator _Iter, typename _Proj = identity,
indirectly_unary_invocable<projected<_Iter, _Proj>> _Fun>
constexpr for_each_n_result<_Iter, _Fun>
operator()(_Iter __first, iter_difference_t<_Iter> __n,
_Fun __f, _Proj __proj = {}) const
{
if constexpr (random_access_iterator<_Iter>)
{
if (__n <= 0)
return {std::move(__first), std::move(__f)};
auto __last = __first + __n;
return ranges::for_each(std::move(__first), std::move(__last),
std::move(__f), std::move(__proj));
}
else
{
while (__n-- > 0)
{
std::__invoke(__f, std::__invoke(__proj, *__first));
++__first;
}
return {std::move(__first), std::move(__f)};
}
}
};
inline constexpr __for_each_n_fn for_each_n{};
// find, find_if and find_if_not are defined in <bits/ranges_util.h>.
struct __find_first_of_fn
{
template<input_iterator _Iter1, sentinel_for<_Iter1> _Sent1,
forward_iterator _Iter2, sentinel_for<_Iter2> _Sent2,
typename _Pred = ranges::equal_to,
typename _Proj1 = identity, typename _Proj2 = identity>
requires indirectly_comparable<_Iter1, _Iter2, _Pred, _Proj1, _Proj2>
[[nodiscard]] constexpr _Iter1
operator()(_Iter1 __first1, _Sent1 __last1,
_Iter2 __first2, _Sent2 __last2, _Pred __pred = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
for (; __first1 != __last1; ++__first1)
for (auto __iter = __first2; __iter != __last2; ++__iter)
if (std::__invoke(__pred,
std::__invoke(__proj1, *__first1),
std::__invoke(__proj2, *__iter)))
return __first1;
return __first1;
}
template<input_range _Range1, forward_range _Range2,
typename _Pred = ranges::equal_to,
typename _Proj1 = identity, typename _Proj2 = identity>
requires indirectly_comparable<iterator_t<_Range1>, iterator_t<_Range2>,
_Pred, _Proj1, _Proj2>
[[nodiscard]] constexpr borrowed_iterator_t<_Range1>
operator()(_Range1&& __r1, _Range2&& __r2, _Pred __pred = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
return (*this)(ranges::begin(__r1), ranges::end(__r1),
ranges::begin(__r2), ranges::end(__r2),
std::move(__pred),
std::move(__proj1), std::move(__proj2));
}
};
inline constexpr __find_first_of_fn find_first_of{};
struct __count_fn
{
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
typename _Tp _GLIBCXX26_RANGE_ALGO_DEF_VAL_T(_Iter, _Proj)>
requires indirect_binary_predicate<ranges::equal_to,
projected<_Iter, _Proj>,
const _Tp*>
[[nodiscard]] constexpr iter_difference_t<_Iter>
operator()(_Iter __first, _Sent __last,
const _Tp& __value, _Proj __proj = {}) const
{
iter_difference_t<_Iter> __n = 0;
for (; __first != __last; ++__first)
if (std::__invoke(__proj, *__first) == __value)
++__n;
return __n;
}
template<input_range _Range, typename _Proj = identity,
typename _Tp
_GLIBCXX26_RANGE_ALGO_DEF_VAL_T(iterator_t<_Range>, _Proj)>
requires indirect_binary_predicate<ranges::equal_to,
projected<iterator_t<_Range>, _Proj>,
const _Tp*>
[[nodiscard]] constexpr range_difference_t<_Range>
operator()(_Range&& __r, const _Tp& __value, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
__value, std::move(__proj));
}
};
inline constexpr __count_fn count{};
struct __count_if_fn
{
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
indirect_unary_predicate<projected<_Iter, _Proj>> _Pred>
constexpr iter_difference_t<_Iter>
operator()(_Iter __first, _Sent __last,
_Pred __pred, _Proj __proj = {}) const
{
iter_difference_t<_Iter> __n = 0;
for (; __first != __last; ++__first)
if (std::__invoke(__pred, std::__invoke(__proj, *__first)))
++__n;
return __n;
}
template<input_range _Range,
typename _Proj = identity,
indirect_unary_predicate<projected<iterator_t<_Range>, _Proj>>
_Pred>
constexpr range_difference_t<_Range>
operator()(_Range&& __r, _Pred __pred, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__pred), std::move(__proj));
}
};
inline constexpr __count_if_fn count_if{};
// in_in_result, mismatch and search are defined in <bits/ranges_util.h>.
struct __search_n_fn
{
template<forward_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Pred = ranges::equal_to, typename _Proj = identity,
typename _Tp _GLIBCXX26_RANGE_ALGO_DEF_VAL_T(_Iter, _Proj)>
requires indirectly_comparable<_Iter, const _Tp*, _Pred, _Proj>
constexpr subrange<_Iter>
operator()(_Iter __first, _Sent __last, iter_difference_t<_Iter> __count,
const _Tp& __value, _Pred __pred = {}, _Proj __proj = {}) const
{
if (__count <= 0)
return {__first, __first};
auto __value_comp = [&] <typename _Rp> (_Rp&& __arg) -> bool {
return std::__invoke(__pred, std::forward<_Rp>(__arg), __value);
};
if (__count == 1)
{
__first = ranges::find_if(std::move(__first), __last,
std::move(__value_comp),
std::move(__proj));
if (__first == __last)
return {__first, __first};
else
{
auto __end = __first;
return {__first, ++__end};
}
}
if constexpr (sized_sentinel_for<_Sent, _Iter>
&& random_access_iterator<_Iter>)
{
auto __tail_size = __last - __first;
auto __remainder = __count;
while (__remainder <= __tail_size)
{
__first += __remainder;
__tail_size -= __remainder;
auto __backtrack = __first;
while (__value_comp(std::__invoke(__proj, *--__backtrack)))
{
if (--__remainder == 0)
return {__first - __count, __first};
}
__remainder = __count + 1 - (__first - __backtrack);
}
auto __i = __first + __tail_size;
return {__i, __i};
}
else
{
__first = ranges::find_if(__first, __last, __value_comp, __proj);
while (__first != __last)
{
auto __n = __count;
auto __i = __first;
++__i;
while (__i != __last && __n != 1
&& __value_comp(std::__invoke(__proj, *__i)))
{
++__i;
--__n;
}
if (__n == 1)
return {__first, __i};
if (__i == __last)
return {__i, __i};
__first = ranges::find_if(++__i, __last, __value_comp, __proj);
}
return {__first, __first};
}
}
template<forward_range _Range,
typename _Pred = ranges::equal_to, typename _Proj = identity,
typename _Tp
_GLIBCXX26_RANGE_ALGO_DEF_VAL_T(iterator_t<_Range>, _Proj)>
requires indirectly_comparable<iterator_t<_Range>, const _Tp*,
_Pred, _Proj>
constexpr borrowed_subrange_t<_Range>
operator()(_Range&& __r, range_difference_t<_Range> __count,
const _Tp& __value, _Pred __pred = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__count), __value,
std::move(__pred), std::move(__proj));
}
};
inline constexpr __search_n_fn search_n{};
#if __glibcxx_ranges_starts_ends_with // C++ >= 23
struct __starts_with_fn
{
template<input_iterator _Iter1, sentinel_for<_Iter1> _Sent1,
input_iterator _Iter2, sentinel_for<_Iter2> _Sent2,
typename _Pred = ranges::equal_to,
typename _Proj1 = identity, typename _Proj2 = identity>
requires indirectly_comparable<_Iter1, _Iter2, _Pred, _Proj1, _Proj2>
constexpr bool
operator()(_Iter1 __first1, _Sent1 __last1,
_Iter2 __first2, _Sent2 __last2, _Pred __pred = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
iter_difference_t<_Iter1> __n1 = -1;
iter_difference_t<_Iter2> __n2 = -1;
if constexpr (sized_sentinel_for<_Sent1, _Iter1>)
__n1 = __last1 - __first1;
if constexpr (sized_sentinel_for<_Sent2, _Iter2>)
__n2 = __last2 - __first2;
return _S_impl(std::move(__first1), __last1, __n1,
std::move(__first2), __last2, __n2,
std::move(__pred),
std::move(__proj1), std::move(__proj2));
}
template<input_range _Range1, input_range _Range2,
typename _Pred = ranges::equal_to,
typename _Proj1 = identity, typename _Proj2 = identity>
requires indirectly_comparable<iterator_t<_Range1>, iterator_t<_Range2>,
_Pred, _Proj1, _Proj2>
constexpr bool
operator()(_Range1&& __r1, _Range2&& __r2, _Pred __pred = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
range_difference_t<_Range1> __n1 = -1;
range_difference_t<_Range2> __n2 = -1;
if constexpr (sized_range<_Range1>)
__n1 = ranges::size(__r1);
if constexpr (sized_range<_Range2>)
__n2 = ranges::size(__r2);
return _S_impl(ranges::begin(__r1), ranges::end(__r1), __n1,
ranges::begin(__r2), ranges::end(__r2), __n2,
std::move(__pred),
std::move(__proj1), std::move(__proj2));
}
private:
template<typename _Iter1, typename _Sent1, typename _Iter2, typename _Sent2,
typename _Pred,
typename _Proj1, typename _Proj2>
static constexpr bool
_S_impl(_Iter1 __first1, _Sent1 __last1, iter_difference_t<_Iter1> __n1,
_Iter2 __first2, _Sent2 __last2, iter_difference_t<_Iter2> __n2,
_Pred __pred, _Proj1 __proj1, _Proj2 __proj2)
{
if (__first2 == __last2) [[unlikely]]
return true;
else if (__n1 == -1 || __n2 == -1)
return ranges::mismatch(std::move(__first1), __last1,
std::move(__first2), __last2,
std::move(__pred),
std::move(__proj1), std::move(__proj2)).in2 == __last2;
else if (__n1 < __n2)
return false;
else if constexpr (random_access_iterator<_Iter1>)
return ranges::equal(__first1, __first1 + iter_difference_t<_Iter1>(__n2),
std::move(__first2), __last2,
std::move(__pred),
std::move(__proj1), std::move(__proj2));
else
return ranges::equal(counted_iterator(std::move(__first1),
iter_difference_t<_Iter1>(__n2)),
default_sentinel,
std::move(__first2), __last2,
std::move(__pred),
std::move(__proj1), std::move(__proj2));
}
friend struct __ends_with_fn;
};
inline constexpr __starts_with_fn starts_with{};
struct __ends_with_fn
{
template<input_iterator _Iter1, sentinel_for<_Iter1> _Sent1,
input_iterator _Iter2, sentinel_for<_Iter2> _Sent2,
typename _Pred = ranges::equal_to,
typename _Proj1 = identity, typename _Proj2 = identity>
requires (forward_iterator<_Iter1> || sized_sentinel_for<_Sent1, _Iter1>)
&& (forward_iterator<_Iter2> || sized_sentinel_for<_Sent2, _Iter2>)
&& indirectly_comparable<_Iter1, _Iter2, _Pred, _Proj1, _Proj2>
constexpr bool
operator()(_Iter1 __first1, _Sent1 __last1,
_Iter2 __first2, _Sent2 __last2, _Pred __pred = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
iter_difference_t<_Iter1> __n1 = -1;
iter_difference_t<_Iter2> __n2 = -1;
if constexpr (sized_sentinel_for<_Sent1, _Iter1>)
__n1 = __last1 - __first1;
if constexpr (sized_sentinel_for<_Sent2, _Iter2>)
__n2 = __last2 - __first2;
return _S_impl(std::move(__first1), __last1, __n1,
std::move(__first2), __last2, __n2,
std::move(__pred),
std::move(__proj1), std::move(__proj2));
}
template<input_range _Range1, input_range _Range2,
typename _Pred = ranges::equal_to,
typename _Proj1 = identity, typename _Proj2 = identity>
requires (forward_range<_Range1> || sized_range<_Range1>)
&& (forward_range<_Range2> || sized_range<_Range2>)
&& indirectly_comparable<iterator_t<_Range1>, iterator_t<_Range2>,
_Pred, _Proj1, _Proj2>
constexpr bool
operator()(_Range1&& __r1, _Range2&& __r2, _Pred __pred = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
range_difference_t<_Range1> __n1 = -1;
range_difference_t<_Range2> __n2 = -1;
if constexpr (sized_range<_Range1>)
__n1 = ranges::size(__r1);
if constexpr (sized_range<_Range2>)
__n2 = ranges::size(__r2);
return _S_impl(ranges::begin(__r1), ranges::end(__r1), __n1,
ranges::begin(__r2), ranges::end(__r2), __n2,
std::move(__pred),
std::move(__proj1), std::move(__proj2));
}
private:
template<typename _Iter1, typename _Sent1,
typename _Iter2, typename _Sent2,
typename _Pred,
typename _Proj1, typename _Proj2>
static constexpr bool
_S_impl(_Iter1 __first1, _Sent1 __last1, iter_difference_t<_Iter1> __n1,
_Iter2 __first2, _Sent2 __last2, iter_difference_t<_Iter2> __n2,
_Pred __pred, _Proj1 __proj1, _Proj2 __proj2)
{
if constexpr (!random_access_iterator<_Iter1>
&& bidirectional_iterator<_Iter1> && same_as<_Iter1, _Sent1>
&& bidirectional_iterator<_Iter2> && same_as<_Iter2, _Sent2>)
return starts_with._S_impl(std::make_reverse_iterator(__last1),
std::make_reverse_iterator(__first1),
__n1,
std::make_reverse_iterator(__last2),
std::make_reverse_iterator(__first2),
__n2,
std::move(__pred),
std::move(__proj1), std::move(__proj2));
if (__first2 == __last2) [[unlikely]]
return true;
if constexpr (forward_iterator<_Iter2>)
if (__n2 == -1)
__n2 = ranges::distance(__first2, __last2);
// __glibcxx_assert(__n2 != -1);
if (__n1 != -1)
{
if (__n1 < __n2)
return false;
auto __shift = __n1 - iter_difference_t<_Iter1>(__n2);
if (random_access_iterator<_Iter1>
|| !bidirectional_iterator<_Iter1>
|| !same_as<_Iter1, _Sent1>
|| __shift < __n2)
{
ranges::advance(__first1, __shift);
return ranges::equal(std::move(__first1), __last1,
std::move(__first2), __last2,
std::move(__pred),
std::move(__proj1), std::move(__proj2));
}
}
if constexpr (bidirectional_iterator<_Iter1> && same_as<_Iter1, _Sent1>)
{
_Iter1 __it1 = __last1;
if (__n1 != -1)
ranges::advance(__it1, -iter_difference_t<_Iter1>(__n2));
else
{
// We can't use ranges::advance if the haystack size is
// unknown, since we need to detect and return false if
// it's smaller than the needle.
iter_difference_t<_Iter2> __m = __n2;
while (__m != 0 && __it1 != __first1)
{
--__m;
--__it1;
}
if (__m != 0)
return false;
}
return ranges::equal(__it1, __last1,
std::move(__first2), __last2,
std::move(__pred),
std::move(__proj1), std::move(__proj2));
}
else if constexpr (forward_iterator<_Iter1>)
{
// __glibcxx_assert(__n1 == -1);
_Iter1 __prev_first1;
__n1 = 0;
while (true)
{
iter_difference_t<_Iter2> __m = __n2;
_Iter1 __it1 = __first1;
while (__m != 0 && __it1 != __last1)
{
++__n1;
--__m;
++__it1;
}
if (__m != 0)
{
// __glibcxx_assert(__it1 == __last1);
if (__n1 < __n2)
return false;
__first1 = ranges::next(__prev_first1,
iter_difference_t<_Iter1>(__n2 - __m));
break;
}
__prev_first1 = __first1;
__first1 = __it1;
}
return ranges::equal(__first1, __last1,
std::move(__first2), __last2,
std::move(__pred),
std::move(__proj1), std::move(__proj2));
}
else
// If the haystack is non-forward then it must be sized, in which case
// we already returned via the __n1 != 1 case.
__builtin_unreachable();
}
};
inline constexpr __ends_with_fn ends_with{};
#endif // __glibcxx_ranges_starts_ends_with
struct __find_end_fn
{
template<forward_iterator _Iter1, sentinel_for<_Iter1> _Sent1,
forward_iterator _Iter2, sentinel_for<_Iter2> _Sent2,
typename _Pred = ranges::equal_to,
typename _Proj1 = identity, typename _Proj2 = identity>
requires indirectly_comparable<_Iter1, _Iter2, _Pred, _Proj1, _Proj2>
[[nodiscard]] constexpr subrange<_Iter1>
operator()(_Iter1 __first1, _Sent1 __last1,
_Iter2 __first2, _Sent2 __last2, _Pred __pred = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
if constexpr (bidirectional_iterator<_Iter1>
&& bidirectional_iterator<_Iter2>)
{
auto __i1 = ranges::next(__first1, __last1);
auto __i2 = ranges::next(__first2, __last2);
auto __rresult
= ranges::search(reverse_iterator<_Iter1>{__i1},
reverse_iterator<_Iter1>{__first1},
reverse_iterator<_Iter2>{__i2},
reverse_iterator<_Iter2>{__first2},
std::move(__pred),
std::move(__proj1), std::move(__proj2));
auto __result_first = ranges::end(__rresult).base();
auto __result_last = ranges::begin(__rresult).base();
if (__result_last == __first1)
return {__i1, __i1};
else
return {__result_first, __result_last};
}
else
{
auto __i = ranges::next(__first1, __last1);
if (__first2 == __last2)
return {__i, __i};
auto __result_begin = __i;
auto __result_end = __i;
for (;;)
{
auto __new_range = ranges::search(__first1, __last1,
__first2, __last2,
__pred, __proj1, __proj2);
auto __new_result_begin = ranges::begin(__new_range);
auto __new_result_end = ranges::end(__new_range);
if (__new_result_begin == __last1)
return {__result_begin, __result_end};
else
{
__result_begin = __new_result_begin;
__result_end = __new_result_end;
__first1 = __result_begin;
++__first1;
}
}
}
}
template<forward_range _Range1, forward_range _Range2,
typename _Pred = ranges::equal_to,
typename _Proj1 = identity, typename _Proj2 = identity>
requires indirectly_comparable<iterator_t<_Range1>, iterator_t<_Range2>,
_Pred, _Proj1, _Proj2>
[[nodiscard]] constexpr borrowed_subrange_t<_Range1>
operator()(_Range1&& __r1, _Range2&& __r2, _Pred __pred = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
return (*this)(ranges::begin(__r1), ranges::end(__r1),
ranges::begin(__r2), ranges::end(__r2),
std::move(__pred),
std::move(__proj1), std::move(__proj2));
}
};
inline constexpr __find_end_fn find_end{};
// adjacent_find is defined in <bits/ranges_util.h>.
struct __is_permutation_fn
{
template<forward_iterator _Iter1, sentinel_for<_Iter1> _Sent1,
forward_iterator _Iter2, sentinel_for<_Iter2> _Sent2,
typename _Proj1 = identity, typename _Proj2 = identity,
indirect_equivalence_relation<projected<_Iter1, _Proj1>,
projected<_Iter2, _Proj2>> _Pred
= ranges::equal_to>
[[nodiscard]] constexpr bool
operator()(_Iter1 __first1, _Sent1 __last1,
_Iter2 __first2, _Sent2 __last2, _Pred __pred = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
constexpr bool __sized_iters
= (sized_sentinel_for<_Sent1, _Iter1>
&& sized_sentinel_for<_Sent2, _Iter2>);
if constexpr (__sized_iters)
{
auto __d1 = ranges::distance(__first1, __last1);
auto __d2 = ranges::distance(__first2, __last2);
if (__d1 != __d2)
return false;
}
// Efficiently compare identical prefixes: O(N) if sequences
// have the same elements in the same order.
for (; __first1 != __last1 && __first2 != __last2;
++__first1, (void)++__first2)
if (!(bool)std::__invoke(__pred,
std::__invoke(__proj1, *__first1),
std::__invoke(__proj2, *__first2)))
break;
if constexpr (__sized_iters)
{
if (__first1 == __last1)
return true;
}
else
{
auto __d1 = ranges::distance(__first1, __last1);
auto __d2 = ranges::distance(__first2, __last2);
if (__d1 == 0 && __d2 == 0)
return true;
if (__d1 != __d2)
return false;
}
for (auto __scan = __first1; __scan != __last1; ++__scan)
{
auto&& __scan_deref = *__scan;
auto&& __proj_scan =
std::__invoke(__proj1, std::forward<decltype(__scan_deref)>(__scan_deref));
auto __comp_scan = [&] <typename _Tp> (_Tp&& __arg) -> bool {
return std::__invoke(__pred,
std::forward<decltype(__proj_scan)>(__proj_scan),
std::forward<_Tp>(__arg));
};
if (__scan != ranges::find_if(__first1, __scan,
__comp_scan, __proj1))
continue; // We've seen this one before.
auto __matches = ranges::count_if(__first2, __last2,
__comp_scan, __proj2);
if (__matches == 0
|| ranges::count_if(__scan, __last1,
__comp_scan, __proj1) != __matches)
return false;
}
return true;
}
template<forward_range _Range1, forward_range _Range2,
typename _Proj1 = identity, typename _Proj2 = identity,
indirect_equivalence_relation<
projected<iterator_t<_Range1>, _Proj1>,
projected<iterator_t<_Range2>, _Proj2>> _Pred = ranges::equal_to>
[[nodiscard]] constexpr bool
operator()(_Range1&& __r1, _Range2&& __r2, _Pred __pred = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
// _GLIBCXX_RESOLVE_LIB_DEFECTS
// 3560. ranges::is_permutation should short-circuit for sized_ranges
if constexpr (sized_range<_Range1>)
if constexpr (sized_range<_Range2>)
if (ranges::distance(__r1) != ranges::distance(__r2))
return false;
return (*this)(ranges::begin(__r1), ranges::end(__r1),
ranges::begin(__r2), ranges::end(__r2),
std::move(__pred),
std::move(__proj1), std::move(__proj2));
}
};
inline constexpr __is_permutation_fn is_permutation{};
template<typename _Iter, typename _Out>
using copy_if_result = in_out_result<_Iter, _Out>;
struct __copy_if_fn
{
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
weakly_incrementable _Out, typename _Proj = identity,
indirect_unary_predicate<projected<_Iter, _Proj>> _Pred>
requires indirectly_copyable<_Iter, _Out>
constexpr copy_if_result<_Iter, _Out>
operator()(_Iter __first, _Sent __last, _Out __result,
_Pred __pred, _Proj __proj = {}) const
{
for (; __first != __last; ++__first)
if (std::__invoke(__pred, std::__invoke(__proj, *__first)))
{
*__result = *__first;
++__result;
}
return {std::move(__first), std::move(__result)};
}
template<input_range _Range, weakly_incrementable _Out,
typename _Proj = identity,
indirect_unary_predicate<projected<iterator_t<_Range>, _Proj>>
_Pred>
requires indirectly_copyable<iterator_t<_Range>, _Out>
constexpr copy_if_result<borrowed_iterator_t<_Range>, _Out>
operator()(_Range&& __r, _Out __result,
_Pred __pred, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__result),
std::move(__pred), std::move(__proj));
}
};
inline constexpr __copy_if_fn copy_if{};
template<typename _Iter1, typename _Iter2>
using swap_ranges_result = in_in_result<_Iter1, _Iter2>;
struct __swap_ranges_fn
{
template<input_iterator _Iter1, sentinel_for<_Iter1> _Sent1,
input_iterator _Iter2, sentinel_for<_Iter2> _Sent2>
requires indirectly_swappable<_Iter1, _Iter2>
constexpr swap_ranges_result<_Iter1, _Iter2>
operator()(_Iter1 __first1, _Sent1 __last1,
_Iter2 __first2, _Sent2 __last2) const
{
for (; __first1 != __last1 && __first2 != __last2;
++__first1, (void)++__first2)
ranges::iter_swap(__first1, __first2);
return {std::move(__first1), std::move(__first2)};
}
template<input_range _Range1, input_range _Range2>
requires indirectly_swappable<iterator_t<_Range1>, iterator_t<_Range2>>
constexpr swap_ranges_result<borrowed_iterator_t<_Range1>,
borrowed_iterator_t<_Range2>>
operator()(_Range1&& __r1, _Range2&& __r2) const
{
return (*this)(ranges::begin(__r1), ranges::end(__r1),
ranges::begin(__r2), ranges::end(__r2));
}
};
inline constexpr __swap_ranges_fn swap_ranges{};
template<typename _Iter, typename _Out>
using unary_transform_result = in_out_result<_Iter, _Out>;
template<typename _Iter1, typename _Iter2, typename _Out>
struct in_in_out_result
{
[[no_unique_address]] _Iter1 in1;
[[no_unique_address]] _Iter2 in2;
[[no_unique_address]] _Out out;
template<typename _IIter1, typename _IIter2, typename _OOut>
requires convertible_to<const _Iter1&, _IIter1>
&& convertible_to<const _Iter2&, _IIter2>
&& convertible_to<const _Out&, _OOut>
constexpr
operator in_in_out_result<_IIter1, _IIter2, _OOut>() const &
{ return {in1, in2, out}; }
template<typename _IIter1, typename _IIter2, typename _OOut>
requires convertible_to<_Iter1, _IIter1>
&& convertible_to<_Iter2, _IIter2>
&& convertible_to<_Out, _OOut>
constexpr
operator in_in_out_result<_IIter1, _IIter2, _OOut>() &&
{ return {std::move(in1), std::move(in2), std::move(out)}; }
};
template<typename _Iter1, typename _Iter2, typename _Out>
using binary_transform_result = in_in_out_result<_Iter1, _Iter2, _Out>;
struct __transform_fn
{
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
weakly_incrementable _Out,
copy_constructible _Fp, typename _Proj = identity>
requires indirectly_writable<_Out,
indirect_result_t<_Fp&,
projected<_Iter, _Proj>>>
constexpr unary_transform_result<_Iter, _Out>
operator()(_Iter __first1, _Sent __last1, _Out __result,
_Fp __op, _Proj __proj = {}) const
{
for (; __first1 != __last1; ++__first1, (void)++__result)
*__result = std::__invoke(__op, std::__invoke(__proj, *__first1));
return {std::move(__first1), std::move(__result)};
}
template<input_range _Range, weakly_incrementable _Out,
copy_constructible _Fp, typename _Proj = identity>
requires indirectly_writable<_Out,
indirect_result_t<_Fp&,
projected<iterator_t<_Range>, _Proj>>>
constexpr unary_transform_result<borrowed_iterator_t<_Range>, _Out>
operator()(_Range&& __r, _Out __result, _Fp __op, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__result),
std::move(__op), std::move(__proj));
}
template<input_iterator _Iter1, sentinel_for<_Iter1> _Sent1,
input_iterator _Iter2, sentinel_for<_Iter2> _Sent2,
weakly_incrementable _Out, copy_constructible _Fp,
typename _Proj1 = identity, typename _Proj2 = identity>
requires indirectly_writable<_Out,
indirect_result_t<_Fp&,
projected<_Iter1, _Proj1>,
projected<_Iter2, _Proj2>>>
constexpr binary_transform_result<_Iter1, _Iter2, _Out>
operator()(_Iter1 __first1, _Sent1 __last1,
_Iter2 __first2, _Sent2 __last2,
_Out __result, _Fp __binary_op,
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
for (; __first1 != __last1 && __first2 != __last2;
++__first1, (void)++__first2, ++__result)
*__result = std::__invoke(__binary_op,
std::__invoke(__proj1, *__first1),
std::__invoke(__proj2, *__first2));
return {std::move(__first1), std::move(__first2), std::move(__result)};
}
template<input_range _Range1, input_range _Range2,
weakly_incrementable _Out, copy_constructible _Fp,
typename _Proj1 = identity, typename _Proj2 = identity>
requires indirectly_writable<_Out,
indirect_result_t<_Fp&,
projected<iterator_t<_Range1>, _Proj1>,
projected<iterator_t<_Range2>, _Proj2>>>
constexpr binary_transform_result<borrowed_iterator_t<_Range1>,
borrowed_iterator_t<_Range2>, _Out>
operator()(_Range1&& __r1, _Range2&& __r2, _Out __result, _Fp __binary_op,
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
return (*this)(ranges::begin(__r1), ranges::end(__r1),
ranges::begin(__r2), ranges::end(__r2),
std::move(__result), std::move(__binary_op),
std::move(__proj1), std::move(__proj2));
}
};
inline constexpr __transform_fn transform{};
struct __replace_fn
{
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
typename _Tp1 _GLIBCXX26_RANGE_ALGO_DEF_VAL_T(_Iter, _Proj),
typename _Tp2 _GLIBCXX26_DEF_VAL_T(_Tp1)>
requires indirectly_writable<_Iter, const _Tp2&>
&& indirect_binary_predicate<ranges::equal_to, projected<_Iter, _Proj>,
const _Tp1*>
constexpr _Iter
operator()(_Iter __first, _Sent __last,
const _Tp1& __old_value, const _Tp2& __new_value,
_Proj __proj = {}) const
{
for (; __first != __last; ++__first)
if (std::__invoke(__proj, *__first) == __old_value)
*__first = __new_value;
return __first;
}
template<input_range _Range, typename _Proj = identity,
typename _Tp1
_GLIBCXX26_RANGE_ALGO_DEF_VAL_T(iterator_t<_Range>, _Proj),
typename _Tp2 _GLIBCXX26_DEF_VAL_T(_Tp1)>
requires indirectly_writable<iterator_t<_Range>, const _Tp2&>
&& indirect_binary_predicate<ranges::equal_to,
projected<iterator_t<_Range>, _Proj>,
const _Tp1*>
constexpr borrowed_iterator_t<_Range>
operator()(_Range&& __r,
const _Tp1& __old_value, const _Tp2& __new_value,
_Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
__old_value, __new_value, std::move(__proj));
}
};
inline constexpr __replace_fn replace{};
struct __replace_if_fn
{
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
typename _Tp _GLIBCXX26_RANGE_ALGO_DEF_VAL_T(_Iter, _Proj),
indirect_unary_predicate<projected<_Iter, _Proj>> _Pred>
requires indirectly_writable<_Iter, const _Tp&>
constexpr _Iter
operator()(_Iter __first, _Sent __last,
_Pred __pred, const _Tp& __new_value, _Proj __proj = {}) const
{
for (; __first != __last; ++__first)
if (std::__invoke(__pred, std::__invoke(__proj, *__first)))
*__first = __new_value;
return std::move(__first);
}
template<input_range _Range, typename _Proj = identity,
typename _Tp
_GLIBCXX26_RANGE_ALGO_DEF_VAL_T(iterator_t<_Range>, _Proj),
indirect_unary_predicate<projected<iterator_t<_Range>, _Proj>>
_Pred>
requires indirectly_writable<iterator_t<_Range>, const _Tp&>
constexpr borrowed_iterator_t<_Range>
operator()(_Range&& __r,
_Pred __pred, const _Tp& __new_value, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__pred), __new_value, std::move(__proj));
}
};
inline constexpr __replace_if_fn replace_if{};
template<typename _Iter, typename _Out>
using replace_copy_result = in_out_result<_Iter, _Out>;
struct __replace_copy_fn
{
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Out, typename _Proj = identity,
typename _Tp1 _GLIBCXX26_RANGE_ALGO_DEF_VAL_T(_Iter, _Proj),
typename _Tp2 _GLIBCXX26_DEF_VAL_T(iter_value_t<_Out>)>
requires indirectly_copyable<_Iter, _Out>
&& indirect_binary_predicate<ranges::equal_to,
projected<_Iter, _Proj>, const _Tp1*>
&& output_iterator<_Out, const _Tp2&>
constexpr replace_copy_result<_Iter, _Out>
operator()(_Iter __first, _Sent __last, _Out __result,
const _Tp1& __old_value, const _Tp2& __new_value,
_Proj __proj = {}) const
{
for (; __first != __last; ++__first, (void)++__result)
if (std::__invoke(__proj, *__first) == __old_value)
*__result = __new_value;
else
*__result = *__first;
return {std::move(__first), std::move(__result)};
}
template<input_range _Range, typename _Out,
typename _Proj = identity,
typename _Tp1
_GLIBCXX26_RANGE_ALGO_DEF_VAL_T(iterator_t<_Range>, _Proj),
typename _Tp2 _GLIBCXX26_DEF_VAL_T(iter_value_t<_Out>)>
requires indirectly_copyable<iterator_t<_Range>, _Out>
&& indirect_binary_predicate<ranges::equal_to,
projected<iterator_t<_Range>, _Proj>,
const _Tp1*>
&& output_iterator<_Out, const _Tp2&>
constexpr replace_copy_result<borrowed_iterator_t<_Range>, _Out>
operator()(_Range&& __r, _Out __result,
const _Tp1& __old_value, const _Tp2& __new_value,
_Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__result), __old_value,
__new_value, std::move(__proj));
}
};
inline constexpr __replace_copy_fn replace_copy{};
template<typename _Iter, typename _Out>
using replace_copy_if_result = in_out_result<_Iter, _Out>;
struct __replace_copy_if_fn
{
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Out,
typename _Tp _GLIBCXX26_DEF_VAL_T(iter_value_t<_Out>),
typename _Proj = identity,
indirect_unary_predicate<projected<_Iter, _Proj>> _Pred>
requires indirectly_copyable<_Iter, _Out>
&& output_iterator<_Out, const _Tp&>
constexpr replace_copy_if_result<_Iter, _Out>
operator()(_Iter __first, _Sent __last, _Out __result,
_Pred __pred, const _Tp& __new_value, _Proj __proj = {}) const
{
for (; __first != __last; ++__first, (void)++__result)
if (std::__invoke(__pred, std::__invoke(__proj, *__first)))
*__result = __new_value;
else
*__result = *__first;
return {std::move(__first), std::move(__result)};
}
template<input_range _Range,
typename _Out,
typename _Tp _GLIBCXX26_DEF_VAL_T(iter_value_t<_Out>),
typename _Proj = identity,
indirect_unary_predicate<projected<iterator_t<_Range>, _Proj>>
_Pred>
requires indirectly_copyable<iterator_t<_Range>, _Out>
&& output_iterator<_Out, const _Tp&>
constexpr replace_copy_if_result<borrowed_iterator_t<_Range>, _Out>
operator()(_Range&& __r, _Out __result,
_Pred __pred, const _Tp& __new_value, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__result), std::move(__pred),
__new_value, std::move(__proj));
}
};
inline constexpr __replace_copy_if_fn replace_copy_if{};
struct __generate_n_fn
{
template<input_or_output_iterator _Out, copy_constructible _Fp>
requires invocable<_Fp&>
&& indirectly_writable<_Out, invoke_result_t<_Fp&>>
constexpr _Out
operator()(_Out __first, iter_difference_t<_Out> __n, _Fp __gen) const
{
for (; __n > 0; --__n, (void)++__first)
*__first = std::__invoke(__gen);
return __first;
}
};
inline constexpr __generate_n_fn generate_n{};
struct __generate_fn
{
template<input_or_output_iterator _Out, sentinel_for<_Out> _Sent,
copy_constructible _Fp>
requires invocable<_Fp&>
&& indirectly_writable<_Out, invoke_result_t<_Fp&>>
constexpr _Out
operator()(_Out __first, _Sent __last, _Fp __gen) const
{
for (; __first != __last; ++__first)
*__first = std::__invoke(__gen);
return __first;
}
template<typename _Range, copy_constructible _Fp>
requires invocable<_Fp&> && output_range<_Range, invoke_result_t<_Fp&>>
constexpr borrowed_iterator_t<_Range>
operator()(_Range&& __r, _Fp __gen) const
{
return (*this)(ranges::begin(__r), ranges::end(__r), std::move(__gen));
}
};
inline constexpr __generate_fn generate{};
struct __remove_if_fn
{
template<permutable _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
indirect_unary_predicate<projected<_Iter, _Proj>> _Pred>
[[nodiscard]] constexpr subrange<_Iter>
operator()(_Iter __first, _Sent __last,
_Pred __pred, _Proj __proj = {}) const
{
__first = ranges::find_if(__first, __last, __pred, __proj);
if (__first == __last)
return {__first, __first};
auto __result = __first;
++__first;
for (; __first != __last; ++__first)
if (!std::__invoke(__pred, std::__invoke(__proj, *__first)))
{
*__result = ranges::iter_move(__first);
++__result;
}
return {__result, __first};
}
template<forward_range _Range, typename _Proj = identity,
indirect_unary_predicate<projected<iterator_t<_Range>, _Proj>>
_Pred>
requires permutable<iterator_t<_Range>>
[[nodiscard]] constexpr borrowed_subrange_t<_Range>
operator()(_Range&& __r, _Pred __pred, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__pred), std::move(__proj));
}
};
inline constexpr __remove_if_fn remove_if{};
struct __remove_fn
{
template<permutable _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
typename _Tp _GLIBCXX26_RANGE_ALGO_DEF_VAL_T(_Iter, _Proj)>
requires indirect_binary_predicate<ranges::equal_to,
projected<_Iter, _Proj>,
const _Tp*>
[[nodiscard]] constexpr subrange<_Iter>
operator()(_Iter __first, _Sent __last,
const _Tp& __value, _Proj __proj = {}) const
{
auto __pred = [&] (auto&& __arg) -> bool {
return std::forward<decltype(__arg)>(__arg) == __value;
};
return ranges::remove_if(__first, __last,
std::move(__pred), std::move(__proj));
}
template<forward_range _Range, typename _Proj = identity,
typename _Tp
_GLIBCXX26_RANGE_ALGO_DEF_VAL_T(iterator_t<_Range>, _Proj)>
requires permutable<iterator_t<_Range>>
&& indirect_binary_predicate<ranges::equal_to,
projected<iterator_t<_Range>, _Proj>,
const _Tp*>
[[nodiscard]] constexpr borrowed_subrange_t<_Range>
operator()(_Range&& __r, const _Tp& __value, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
__value, std::move(__proj));
}
};
inline constexpr __remove_fn remove{};
template<typename _Iter, typename _Out>
using remove_copy_if_result = in_out_result<_Iter, _Out>;
struct __remove_copy_if_fn
{
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
weakly_incrementable _Out, typename _Proj = identity,
indirect_unary_predicate<projected<_Iter, _Proj>> _Pred>
requires indirectly_copyable<_Iter, _Out>
constexpr remove_copy_if_result<_Iter, _Out>
operator()(_Iter __first, _Sent __last, _Out __result,
_Pred __pred, _Proj __proj = {}) const
{
for (; __first != __last; ++__first)
if (!std::__invoke(__pred, std::__invoke(__proj, *__first)))
{
*__result = *__first;
++__result;
}
return {std::move(__first), std::move(__result)};
}
template<input_range _Range, weakly_incrementable _Out,
typename _Proj = identity,
indirect_unary_predicate<projected<iterator_t<_Range>, _Proj>>
_Pred>
requires indirectly_copyable<iterator_t<_Range>, _Out>
constexpr remove_copy_if_result<borrowed_iterator_t<_Range>, _Out>
operator()(_Range&& __r, _Out __result,
_Pred __pred, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__result),
std::move(__pred), std::move(__proj));
}
};
inline constexpr __remove_copy_if_fn remove_copy_if{};
template<typename _Iter, typename _Out>
using remove_copy_result = in_out_result<_Iter, _Out>;
struct __remove_copy_fn
{
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
weakly_incrementable _Out, typename _Proj = identity,
typename _Tp _GLIBCXX26_RANGE_ALGO_DEF_VAL_T(_Iter, _Proj)>
requires indirectly_copyable<_Iter, _Out>
&& indirect_binary_predicate<ranges::equal_to,
projected<_Iter, _Proj>,
const _Tp*>
constexpr remove_copy_result<_Iter, _Out>
operator()(_Iter __first, _Sent __last, _Out __result,
const _Tp& __value, _Proj __proj = {}) const
{
for (; __first != __last; ++__first)
if (!(std::__invoke(__proj, *__first) == __value))
{
*__result = *__first;
++__result;
}
return {std::move(__first), std::move(__result)};
}
template<input_range _Range, weakly_incrementable _Out,
typename _Proj = identity,
typename _Tp
_GLIBCXX26_RANGE_ALGO_DEF_VAL_T(iterator_t<_Range>, _Proj)>
requires indirectly_copyable<iterator_t<_Range>, _Out>
&& indirect_binary_predicate<ranges::equal_to,
projected<iterator_t<_Range>, _Proj>,
const _Tp*>
constexpr remove_copy_result<borrowed_iterator_t<_Range>, _Out>
operator()(_Range&& __r, _Out __result,
const _Tp& __value, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__result), __value, std::move(__proj));
}
};
inline constexpr __remove_copy_fn remove_copy{};
struct __unique_fn
{
template<permutable _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
indirect_equivalence_relation<
projected<_Iter, _Proj>> _Comp = ranges::equal_to>
[[nodiscard]] constexpr subrange<_Iter>
operator()(_Iter __first, _Sent __last,
_Comp __comp = {}, _Proj __proj = {}) const
{
__first = ranges::adjacent_find(__first, __last, __comp, __proj);
if (__first == __last)
return {__first, __first};
auto __dest = __first;
++__first;
while (++__first != __last)
if (!std::__invoke(__comp,
std::__invoke(__proj, *__dest),
std::__invoke(__proj, *__first)))
*++__dest = ranges::iter_move(__first);
return {++__dest, __first};
}
template<forward_range _Range, typename _Proj = identity,
indirect_equivalence_relation<
projected<iterator_t<_Range>, _Proj>> _Comp = ranges::equal_to>
requires permutable<iterator_t<_Range>>
[[nodiscard]] constexpr borrowed_subrange_t<_Range>
operator()(_Range&& __r, _Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__comp), std::move(__proj));
}
};
inline constexpr __unique_fn unique{};
namespace __detail
{
template<typename _Out, typename _Tp>
concept __can_reread_output = input_iterator<_Out>
&& same_as<_Tp, iter_value_t<_Out>>;
}
template<typename _Iter, typename _Out>
using unique_copy_result = in_out_result<_Iter, _Out>;
struct __unique_copy_fn
{
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
weakly_incrementable _Out, typename _Proj = identity,
indirect_equivalence_relation<
projected<_Iter, _Proj>> _Comp = ranges::equal_to>
requires indirectly_copyable<_Iter, _Out>
&& (forward_iterator<_Iter>
|| __detail::__can_reread_output<_Out, iter_value_t<_Iter>>
|| indirectly_copyable_storable<_Iter, _Out>)
constexpr unique_copy_result<_Iter, _Out>
operator()(_Iter __first, _Sent __last, _Out __result,
_Comp __comp = {}, _Proj __proj = {}) const
{
if (__first == __last)
return {std::move(__first), std::move(__result)};
// _GLIBCXX_RESOLVE_LIB_DEFECTS
// 4269. unique_copy passes arguments to its predicate backwards
// TODO: perform a closer comparison with reference implementations
if constexpr (forward_iterator<_Iter>)
{
auto __next = __first;
*__result = *__next;
while (++__next != __last)
if (!std::__invoke(__comp,
std::__invoke(__proj, *__first),
std::__invoke(__proj, *__next)))
{
__first = __next;
*++__result = *__first;
}
return {__next, std::move(++__result)};
}
else if constexpr (__detail::__can_reread_output<_Out, iter_value_t<_Iter>>)
{
*__result = *__first;
while (++__first != __last)
if (!std::__invoke(__comp,
std::__invoke(__proj, *__result),
std::__invoke(__proj, *__first)))
*++__result = *__first;
return {std::move(__first), std::move(++__result)};
}
else // indirectly_copyable_storable<_Iter, _Out>
{
iter_value_t<_Iter> __value(*__first);
*__result = __value;
while (++__first != __last)
{
if (!(bool)std::__invoke(__comp,
std::__invoke(__proj, __value),
std::__invoke(__proj, *__first)))
{
__value = *__first;
*++__result = __value;
}
}
return {std::move(__first), std::move(++__result)};
}
}
template<input_range _Range,
weakly_incrementable _Out, typename _Proj = identity,
indirect_equivalence_relation<
projected<iterator_t<_Range>, _Proj>> _Comp = ranges::equal_to>
requires indirectly_copyable<iterator_t<_Range>, _Out>
&& (forward_iterator<iterator_t<_Range>>
|| __detail::__can_reread_output<_Out, range_value_t<_Range>>
|| indirectly_copyable_storable<iterator_t<_Range>, _Out>)
constexpr unique_copy_result<borrowed_iterator_t<_Range>, _Out>
operator()(_Range&& __r, _Out __result,
_Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__result),
std::move(__comp), std::move(__proj));
}
};
inline constexpr __unique_copy_fn unique_copy{};
struct __reverse_fn
{
template<bidirectional_iterator _Iter, sentinel_for<_Iter> _Sent>
requires permutable<_Iter>
constexpr _Iter
operator()(_Iter __first, _Sent __last) const
{
auto __i = ranges::next(__first, __last);
auto __tail = __i;
if constexpr (random_access_iterator<_Iter>)
{
if (__first != __last)
{
--__tail;
while (__first < __tail)
{
ranges::iter_swap(__first, __tail);
++__first;
--__tail;
}
}
return __i;
}
else
{
for (;;)
if (__first == __tail || __first == --__tail)
break;
else
{
ranges::iter_swap(__first, __tail);
++__first;
}
return __i;
}
}
template<bidirectional_range _Range>
requires permutable<iterator_t<_Range>>
constexpr borrowed_iterator_t<_Range>
operator()(_Range&& __r) const
{
return (*this)(ranges::begin(__r), ranges::end(__r));
}
};
inline constexpr __reverse_fn reverse{};
template<typename _Iter, typename _Out>
using reverse_copy_result = in_out_result<_Iter, _Out>;
struct __reverse_copy_fn
{
template<bidirectional_iterator _Iter, sentinel_for<_Iter> _Sent,
weakly_incrementable _Out>
requires indirectly_copyable<_Iter, _Out>
constexpr reverse_copy_result<_Iter, _Out>
operator()(_Iter __first, _Sent __last, _Out __result) const
{
auto __i = ranges::next(__first, __last);
auto __tail = __i;
while (__first != __tail)
{
--__tail;
*__result = *__tail;
++__result;
}
return {__i, std::move(__result)};
}
template<bidirectional_range _Range, weakly_incrementable _Out>
requires indirectly_copyable<iterator_t<_Range>, _Out>
constexpr reverse_copy_result<borrowed_iterator_t<_Range>, _Out>
operator()(_Range&& __r, _Out __result) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__result));
}
};
inline constexpr __reverse_copy_fn reverse_copy{};
struct __rotate_fn
{
template<permutable _Iter, sentinel_for<_Iter> _Sent>
constexpr subrange<_Iter>
operator()(_Iter __first, _Iter __middle, _Sent __last) const
{
auto __lasti = ranges::next(__first, __last);
if (__first == __middle)
return {__lasti, __lasti};
if (__last == __middle)
return {std::move(__first), std::move(__lasti)};
if constexpr (random_access_iterator<_Iter>)
{
auto __n = __lasti - __first;
auto __k = __middle - __first;
if (__k == __n - __k)
{
ranges::swap_ranges(__first, __middle, __middle, __middle + __k);
return {std::move(__middle), std::move(__lasti)};
}
auto __p = __first;
auto __ret = __first + (__lasti - __middle);
for (;;)
{
if (__k < __n - __k)
{
// TODO: is_pod is deprecated, but this condition is
// consistent with the STL implementation.
if constexpr (__is_pod(iter_value_t<_Iter>))
if (__k == 1)
{
auto __mid = ranges::next(__p, __n - 1);
auto __end = ranges::next(__mid);
iter_value_t<_Iter> __t(ranges::iter_move(__p));
ranges::move(ranges::next(__p), __end, __p);
*__mid = std::move(__t);
return {std::move(__ret), std::move(__lasti)};
}
auto __q = __p + __k;
for (decltype(__n) __i = 0; __i < __n - __k; ++ __i)
{
ranges::iter_swap(__p, __q);
++__p;
++__q;
}
__n %= __k;
if (__n == 0)
return {std::move(__ret), std::move(__lasti)};
ranges::swap(__n, __k);
__k = __n - __k;
}
else
{
__k = __n - __k;
// TODO: is_pod is deprecated, but this condition is
// consistent with the STL implementation.
if constexpr (__is_pod(iter_value_t<_Iter>))
if (__k == 1)
{
auto __mid = ranges::next(__p, __n - 1);
auto __end = ranges::next(__mid);
iter_value_t<_Iter> __t(ranges::iter_move(__mid));
ranges::move_backward(__p, __mid, __end);
*__p = std::move(__t);
return {std::move(__ret), std::move(__lasti)};
}
auto __q = __p + __n;
__p = __q - __k;
for (decltype(__n) __i = 0; __i < __n - __k; ++ __i)
{
--__p;
--__q;
ranges::iter_swap(__p, __q);
}
__n %= __k;
if (__n == 0)
return {std::move(__ret), std::move(__lasti)};
std::swap(__n, __k);
}
}
}
else if constexpr (bidirectional_iterator<_Iter>)
{
auto __tail = __lasti;
ranges::reverse(__first, __middle);
ranges::reverse(__middle, __tail);
while (__first != __middle && __middle != __tail)
{
ranges::iter_swap(__first, --__tail);
++__first;
}
if (__first == __middle)
{
ranges::reverse(__middle, __tail);
return {std::move(__tail), std::move(__lasti)};
}
else
{
ranges::reverse(__first, __middle);
return {std::move(__first), std::move(__lasti)};
}
}
else
{
auto __first2 = __middle;
do
{
ranges::iter_swap(__first, __first2);
++__first;
++__first2;
if (__first == __middle)
__middle = __first2;
} while (__first2 != __last);
auto __ret = __first;
__first2 = __middle;
while (__first2 != __last)
{
ranges::iter_swap(__first, __first2);
++__first;
++__first2;
if (__first == __middle)
__middle = __first2;
else if (__first2 == __last)
__first2 = __middle;
}
return {std::move(__ret), std::move(__lasti)};
}
}
template<forward_range _Range>
requires permutable<iterator_t<_Range>>
constexpr borrowed_subrange_t<_Range>
operator()(_Range&& __r, iterator_t<_Range> __middle) const
{
return (*this)(ranges::begin(__r), std::move(__middle),
ranges::end(__r));
}
};
inline constexpr __rotate_fn rotate{};
template<typename _Iter, typename _Out>
using rotate_copy_result = in_out_result<_Iter, _Out>;
struct __rotate_copy_fn
{
template<forward_iterator _Iter, sentinel_for<_Iter> _Sent,
weakly_incrementable _Out>
requires indirectly_copyable<_Iter, _Out>
constexpr rotate_copy_result<_Iter, _Out>
operator()(_Iter __first, _Iter __middle, _Sent __last,
_Out __result) const
{
auto __copy1 = ranges::copy(__middle,
std::move(__last),
std::move(__result));
auto __copy2 = ranges::copy(std::move(__first),
std::move(__middle),
std::move(__copy1.out));
return { std::move(__copy1.in), std::move(__copy2.out) };
}
template<forward_range _Range, weakly_incrementable _Out>
requires indirectly_copyable<iterator_t<_Range>, _Out>
constexpr rotate_copy_result<borrowed_iterator_t<_Range>, _Out>
operator()(_Range&& __r, iterator_t<_Range> __middle, _Out __result) const
{
return (*this)(ranges::begin(__r), std::move(__middle),
ranges::end(__r), std::move(__result));
}
};
inline constexpr __rotate_copy_fn rotate_copy{};
struct __sample_fn
{
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
weakly_incrementable _Out, typename _Gen>
requires (forward_iterator<_Iter> || random_access_iterator<_Out>)
&& indirectly_copyable<_Iter, _Out>
&& uniform_random_bit_generator<remove_reference_t<_Gen>>
_Out
operator()(_Iter __first, _Sent __last, _Out __out,
iter_difference_t<_Iter> __n, _Gen&& __g) const
{
// FIXME: Correctly handle integer-class difference types.
if constexpr (forward_iterator<_Iter>)
{
using _Size = iter_difference_t<_Iter>;
using __distrib_type = uniform_int_distribution<_Size>;
using __param_type = typename __distrib_type::param_type;
using _USize = __detail::__make_unsigned_like_t<_Size>;
using __uc_type
= common_type_t<typename remove_reference_t<_Gen>::result_type, _USize>;
if (__first == __last)
return __out;
__distrib_type __d{};
_Size __unsampled_sz = ranges::distance(__first, __last);
__n = std::min(__n, __unsampled_sz);
// If possible, we use __gen_two_uniform_ints to efficiently produce
// two random numbers using a single distribution invocation:
const __uc_type __urngrange = __g.max() - __g.min();
if (__urngrange / __uc_type(__unsampled_sz) >= __uc_type(__unsampled_sz))
// I.e. (__urngrange >= __unsampled_sz * __unsampled_sz) but without
// wrapping issues.
{
while (__n != 0 && __unsampled_sz >= 2)
{
const pair<_Size, _Size> __p =
__gen_two_uniform_ints(__unsampled_sz, __unsampled_sz - 1, __g);
--__unsampled_sz;
if (__p.first < __n)
{
*__out = *__first;
++__out;
--__n;
}
++__first;
if (__n == 0) break;
--__unsampled_sz;
if (__p.second < __n)
{
*__out = *__first;
++__out;
--__n;
}
++__first;
}
}
// The loop above is otherwise equivalent to this one-at-a-time version:
for (; __n != 0; ++__first)
if (__d(__g, __param_type{0, --__unsampled_sz}) < __n)
{
*__out = *__first;
++__out;
--__n;
}
return __out;
}
else
{
using __distrib_type
= uniform_int_distribution<iter_difference_t<_Iter>>;
using __param_type = typename __distrib_type::param_type;
__distrib_type __d{};
iter_difference_t<_Iter> __sample_sz = 0;
while (__first != __last && __sample_sz != __n)
{
__out[__sample_sz++] = *__first;
++__first;
}
for (auto __pop_sz = __sample_sz; __first != __last;
++__first, (void) ++__pop_sz)
{
const auto __k = __d(__g, __param_type{0, __pop_sz});
if (__k < __n)
__out[__k] = *__first;
}
return __out + iter_difference_t<_Out>(__sample_sz);
}
}
template<input_range _Range, weakly_incrementable _Out, typename _Gen>
requires (forward_range<_Range> || random_access_iterator<_Out>)
&& indirectly_copyable<iterator_t<_Range>, _Out>
&& uniform_random_bit_generator<remove_reference_t<_Gen>>
_Out
operator()(_Range&& __r, _Out __out,
range_difference_t<_Range> __n, _Gen&& __g) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__out), __n,
std::forward<_Gen>(__g));
}
};
inline constexpr __sample_fn sample{};
struct __shuffle_fn
{
template<random_access_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Gen>
requires permutable<_Iter>
&& uniform_random_bit_generator<remove_reference_t<_Gen>>
_Iter
operator()(_Iter __first, _Sent __last, _Gen&& __g) const
{
// FIXME: Correctly handle integer-class difference types.
if (__first == __last)
return __first;
using _DistanceType = iter_difference_t<_Iter>;
using __ud_type = __detail::__make_unsigned_like_t<_DistanceType>;
using __distr_type = std::uniform_int_distribution<__ud_type>;
using __p_type = typename __distr_type::param_type;
using __uc_type
= common_type_t<typename remove_reference_t<_Gen>::result_type, __ud_type>;
if constexpr (sized_sentinel_for<_Sent, _Iter>)
{
const __uc_type __urngrange = __g.max() - __g.min();
const __uc_type __urange = __uc_type(__last - __first);
if (__urngrange / __urange >= __urange)
// I.e. (__urngrange >= __urange * __urange) but without wrap issues.
{
_Iter __i = ranges::next(__first);
// Since we know the range isn't empty, an even number of elements
// means an uneven number of elements /to swap/, in which case we
// do the first one up front:
if ((__urange % 2) == 0)
{
__distr_type __d{0, 1};
ranges::iter_swap(__i++, ranges::next(__first, __d(__g)));
}
// Now we know that __last - __i is even, so we do the rest in pairs,
// using a single distribution invocation to produce swap positions
// for two successive elements at a time:
while (__i != __last)
{
const __uc_type __swap_range = __uc_type(__i - __first) + 1;
const pair<_DistanceType, _DistanceType> __pospos =
__gen_two_uniform_ints(__swap_range, __swap_range + 1, __g);
ranges::iter_swap(__i++, ranges::next(__first, __pospos.first));
ranges::iter_swap(__i++, ranges::next(__first, __pospos.second));
}
return __i;
}
}
__distr_type __d;
_Iter __i = ranges::next(__first);
for (; __i != __last; ++__i)
ranges::iter_swap(__i,
ranges::next(__first,
__d(__g, __p_type(0, __i - __first))));
return __i;
}
template<random_access_range _Range, typename _Gen>
requires permutable<iterator_t<_Range>>
&& uniform_random_bit_generator<remove_reference_t<_Gen>>
borrowed_iterator_t<_Range>
operator()(_Range&& __r, _Gen&& __g) const
{
if constexpr (sized_range<_Range>
&& !sized_sentinel_for<sentinel_t<_Range>,
iterator_t<_Range>>)
return (*this)(ranges::begin(__r),
ranges::begin(__r) + ranges::distance(__r),
std::forward<_Gen>(__g));
else
return (*this)(ranges::begin(__r), ranges::end(__r),
std::forward<_Gen>(__g));
}
};
inline constexpr __shuffle_fn shuffle{};
namespace __detail
{
template<typename _Iter, typename _Comp>
constexpr void
__push_heap(_Iter __first,
iter_difference_t<_Iter> __holeIndex,
iter_difference_t<_Iter> __topIndex,
iter_value_t<_Iter> __value,
_Comp __comp)
{
auto __parent = (__holeIndex - 1) / 2;
while (__holeIndex > __topIndex
&& __comp(*(__first + __parent), __value))
{
*(__first + __holeIndex) = ranges::iter_move(__first + __parent);
__holeIndex = __parent;
__parent = (__holeIndex - 1) / 2;
}
*(__first + __holeIndex) = std::move(__value);
}
} // namespace __detail
struct __push_heap_fn
{
template<random_access_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Comp = ranges::less, typename _Proj = identity>
requires sortable<_Iter, _Comp, _Proj>
constexpr _Iter
operator()(_Iter __first, _Sent __last,
_Comp __comp = {}, _Proj __proj = {}) const
{
if constexpr (!same_as<_Iter, _Sent>)
return (*this)(__first, ranges::next(__first, __last),
std::move(__comp), std::move(__proj));
else
{
auto __comp_proj = __detail::__make_comp_proj(__comp, __proj);
iter_value_t<_Iter> __value(ranges::iter_move(ranges::prev(__last)));
__detail::__push_heap(__first, (__last - __first) - 1,
0, std::move(__value), __comp_proj);
return __last;
}
}
template<random_access_range _Range,
typename _Comp = ranges::less, typename _Proj = identity>
requires sortable<iterator_t<_Range>, _Comp, _Proj>
constexpr borrowed_iterator_t<_Range>
operator()(_Range&& __r, _Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__comp), std::move(__proj));
}
};
inline constexpr __push_heap_fn push_heap{};
namespace __detail
{
template<typename _Iter, typename _Comp>
constexpr void
__adjust_heap(_Iter __first,
iter_difference_t<_Iter> __holeIndex,
iter_difference_t<_Iter> __len,
iter_value_t<_Iter> __value,
_Comp __comp)
{
auto __topIndex = __holeIndex;
auto __secondChild = __holeIndex;
while (__secondChild < (__len - 1) / 2)
{
__secondChild = 2 * (__secondChild + 1);
if (__comp(*(__first + __secondChild),
*(__first + (__secondChild - 1))))
__secondChild--;
*(__first + __holeIndex) = ranges::iter_move(__first + __secondChild);
__holeIndex = __secondChild;
}
if ((__len & 1) == 0 && __secondChild == (__len - 2) / 2)
{
__secondChild = 2 * (__secondChild + 1);
*(__first + __holeIndex) = ranges::iter_move(__first + (__secondChild - 1));
__holeIndex = __secondChild - 1;
}
__detail::__push_heap(__first, __holeIndex, __topIndex,
std::move(__value), __comp);
}
template<typename _Iter, typename _Comp>
constexpr void
__pop_heap(_Iter __first, _Iter __last, _Iter __result, _Comp __comp)
{
iter_value_t<_Iter> __value = ranges::iter_move(__result);
*__result = ranges::iter_move(__first);
__detail::__adjust_heap(__first, 0, __last - __first,
std::move(__value), __comp);
}
} // namespace __detail
struct __pop_heap_fn
{
template<random_access_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Comp = ranges::less, typename _Proj = identity>
requires sortable<_Iter, _Comp, _Proj>
constexpr _Iter
operator()(_Iter __first, _Sent __last,
_Comp __comp = {}, _Proj __proj = {}) const
{
if constexpr (!same_as<_Iter, _Sent>)
return (*this)(__first, ranges::next(__first, __last),
std::move(__comp), std::move(__proj));
else
{
if (__last - __first > 1)
{
auto __back = ranges::prev(__last);
auto __comp_proj = __detail::__make_comp_proj(__comp, __proj);
__detail::__pop_heap(__first, __back, __back, __comp_proj);
}
return __last;
}
}
template<random_access_range _Range,
typename _Comp = ranges::less, typename _Proj = identity>
requires sortable<iterator_t<_Range>, _Comp, _Proj>
constexpr borrowed_iterator_t<_Range>
operator()(_Range&& __r, _Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__comp), std::move(__proj));
}
};
inline constexpr __pop_heap_fn pop_heap{};
struct __make_heap_fn
{
template<random_access_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Comp = ranges::less, typename _Proj = identity>
requires sortable<_Iter, _Comp, _Proj>
constexpr _Iter
operator()(_Iter __first, _Sent __last,
_Comp __comp = {}, _Proj __proj = {}) const
{
if constexpr (!same_as<_Iter, _Sent>)
return (*this)(__first, ranges::next(__first, __last),
std::move(__comp), std::move(__proj));
else
{
const auto __len = __last - __first;
if (__len < 2)
return __last;
auto __comp_proj = __detail::__make_comp_proj(__comp, __proj);
auto __parent = (__len - 2) / 2;
while (true)
{
iter_value_t<_Iter> __value = ranges::iter_move(__first + __parent);
__detail::__adjust_heap(__first, __parent, __len,
std::move(__value),
__comp_proj);
if (__parent == 0)
break;
__parent--;
}
return __last;
}
}
template<random_access_range _Range,
typename _Comp = ranges::less, typename _Proj = identity>
requires sortable<iterator_t<_Range>, _Comp, _Proj>
constexpr borrowed_iterator_t<_Range>
operator()(_Range&& __r, _Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__comp), std::move(__proj));
}
};
inline constexpr __make_heap_fn make_heap{};
struct __sort_heap_fn
{
template<random_access_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Comp = ranges::less, typename _Proj = identity>
requires sortable<_Iter, _Comp, _Proj>
constexpr _Iter
operator()(_Iter __first, _Sent __last,
_Comp __comp = {}, _Proj __proj = {}) const
{
if constexpr (!same_as<_Iter, _Sent>)
return (*this)(__first, ranges::next(__first, __last),
std::move(__comp), std::move(__proj));
else
{
auto __comp_proj = __detail::__make_comp_proj(__comp, __proj);
_Iter __ret = __last;
while (__last - __first > 1)
{
--__last;
__detail::__pop_heap(__first, __last, __last, __comp_proj);
}
return __ret;
}
}
template<random_access_range _Range,
typename _Comp = ranges::less, typename _Proj = identity>
requires sortable<iterator_t<_Range>, _Comp, _Proj>
constexpr borrowed_iterator_t<_Range>
operator()(_Range&& __r, _Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__comp), std::move(__proj));
}
};
inline constexpr __sort_heap_fn sort_heap{};
struct __is_heap_until_fn
{
template<random_access_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
indirect_strict_weak_order<projected<_Iter, _Proj>>
_Comp = ranges::less>
constexpr _Iter
operator()(_Iter __first, _Sent __last,
_Comp __comp = {}, _Proj __proj = {}) const
{
iter_difference_t<_Iter> __n = ranges::distance(__first, __last);
iter_difference_t<_Iter> __parent = 0, __child = 1;
for (; __child < __n; ++__child)
if (std::__invoke(__comp,
std::__invoke(__proj, *(__first + __parent)),
std::__invoke(__proj, *(__first + __child))))
return __first + __child;
else if ((__child & 1) == 0)
++__parent;
return __first + __n;
}
template<random_access_range _Range,
typename _Proj = identity,
indirect_strict_weak_order<projected<iterator_t<_Range>, _Proj>>
_Comp = ranges::less>
constexpr borrowed_iterator_t<_Range>
operator()(_Range&& __r, _Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__comp), std::move(__proj));
}
};
inline constexpr __is_heap_until_fn is_heap_until{};
struct __is_heap_fn
{
template<random_access_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
indirect_strict_weak_order<projected<_Iter, _Proj>>
_Comp = ranges::less>
constexpr bool
operator()(_Iter __first, _Sent __last,
_Comp __comp = {}, _Proj __proj = {}) const
{
return (__last
== ranges::is_heap_until(__first, __last,
std::move(__comp),
std::move(__proj)));
}
template<random_access_range _Range,
typename _Proj = identity,
indirect_strict_weak_order<projected<iterator_t<_Range>, _Proj>>
_Comp = ranges::less>
constexpr bool
operator()(_Range&& __r, _Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__comp), std::move(__proj));
}
};
inline constexpr __is_heap_fn is_heap{};
namespace __detail
{
template<typename _Iter, typename _Comp>
constexpr void
__move_median_to_first(_Iter __result, _Iter __a, _Iter __b, _Iter __c,
_Comp __comp)
{
if (__comp(*__a, *__b))
{
if (__comp(*__b, *__c))
ranges::iter_swap(__result, __b);
else if (__comp(*__a, *__c))
ranges::iter_swap(__result, __c);
else
ranges::iter_swap(__result, __a);
}
else if (__comp(*__a, *__c))
ranges::iter_swap(__result, __a);
else if (__comp(*__b, *__c))
ranges::iter_swap(__result, __c);
else
ranges::iter_swap(__result, __b);
}
template<typename _Iter, typename _Comp>
constexpr void
__unguarded_linear_insert(_Iter __last, _Comp __comp)
{
iter_value_t<_Iter> __val = ranges::iter_move(__last);
_Iter __next = __last;
--__next;
while (__comp(__val, *__next))
{
*__last = ranges::iter_move(__next);
__last = __next;
--__next;
}
*__last = std::move(__val);
}
template<typename _Iter, typename _Comp>
constexpr void
__insertion_sort(_Iter __first, _Iter __last, _Comp __comp)
{
if (__first == __last)
return;
for (_Iter __i = ranges::next(__first); __i != __last; ++__i)
{
if (__comp(*__i, *__first))
{
iter_value_t<_Iter> __val = ranges::iter_move(__i);
ranges::move_backward(__first, __i, ranges::next(__i));
*__first = std::move(__val);
}
else
__detail::__unguarded_linear_insert(__i, __comp);
}
}
template<typename _Iter, typename _Comp>
constexpr void
__unguarded_insertion_sort(_Iter __first, _Iter __last, _Comp __comp)
{
for (_Iter __i = __first; __i != __last; ++__i)
__detail::__unguarded_linear_insert(__i, __comp);
}
inline constexpr int __sort_threshold = 16;
template<typename _Iter, typename _Comp>
constexpr void
__final_insertion_sort(_Iter __first, _Iter __last, _Comp __comp)
{
constexpr iter_difference_t<_Iter> __threshold = __sort_threshold;
if (__last - __first > __threshold)
{
__detail::__insertion_sort(__first, __first + __threshold, __comp);
__detail::__unguarded_insertion_sort(__first + __threshold, __last,
__comp);
}
else
__detail::__insertion_sort(__first, __last, __comp);
}
template<typename _Iter, typename _Comp>
constexpr _Iter
__unguarded_partition(_Iter __first, _Iter __last, _Iter __pivot, _Comp __comp)
{
while (true)
{
while (__comp(*__first, *__pivot))
++__first;
--__last;
while (__comp(*__pivot, *__last))
--__last;
if (!(__first < __last))
return __first;
ranges::iter_swap(__first, __last);
++__first;
}
}
template<typename _Iter, typename _Comp>
constexpr _Iter
__unguarded_partition_pivot(_Iter __first, _Iter __last, _Comp __comp)
{
_Iter __mid = __first + (__last - __first) / 2;
__detail::__move_median_to_first(__first, ranges::next(__first), __mid,
ranges::prev(__last), __comp);
return __detail::__unguarded_partition(ranges::next(__first), __last,
__first, __comp);
}
template<typename _Iter, typename _Comp>
constexpr void
__heap_select(_Iter __first, _Iter __middle, _Iter __last, _Comp __comp)
{
ranges::make_heap(__first, __middle, __comp);
for (_Iter __i = __middle; __i < __last; ++__i)
if (__comp(*__i, *__first))
__detail::__pop_heap(__first, __middle, __i, __comp);
}
template<typename _Iter, typename _Comp>
constexpr void
__partial_sort(_Iter __first, _Iter __middle, _Iter __last, _Comp __comp)
{
__detail::__heap_select(__first, __middle, __last, __comp);
ranges::sort_heap(__first, __middle, __comp);
}
template<typename _Iter, typename _Comp>
constexpr void
__introsort_loop(_Iter __first, _Iter __last, unsigned __depth_limit, _Comp __comp)
{
while (__last - __first > __sort_threshold)
{
if (__depth_limit == 0)
{
__detail::__partial_sort(__first, __last, __last, __comp);
return;
}
--__depth_limit;
_Iter __cut = __detail::__unguarded_partition_pivot(__first, __last, __comp);
__detail::__introsort_loop(__cut, __last, __depth_limit, __comp);
__last = __cut;
}
}
} // namespace __detail
struct __sort_fn
{
template<random_access_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Comp = ranges::less, typename _Proj = identity>
requires sortable<_Iter, _Comp, _Proj>
constexpr _Iter
operator()(_Iter __first, _Sent __last,
_Comp __comp = {}, _Proj __proj = {}) const
{
if constexpr (!same_as<_Iter, _Sent>)
return (*this)(__first, ranges::next(__first, __last),
std::move(__comp), std::move(__proj));
else
{
if (__first != __last)
{
auto __comp_proj = __detail::__make_comp_proj(__comp, __proj);
auto __n = __detail::__to_unsigned_like(__last - __first);
unsigned __depth_limit = (std::__bit_width(__n) - 1) * 2;
__detail::__introsort_loop(__first, __last, __depth_limit, __comp_proj);
__detail::__final_insertion_sort(__first, __last, __comp_proj);
}
return __last;
}
}
template<random_access_range _Range,
typename _Comp = ranges::less, typename _Proj = identity>
requires sortable<iterator_t<_Range>, _Comp, _Proj>
constexpr borrowed_iterator_t<_Range>
operator()(_Range&& __r, _Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__comp), std::move(__proj));
}
};
inline constexpr __sort_fn sort{};
namespace __detail
{
// This is a helper function for the __merge_sort_loop routines.
template<typename _Iter, typename _Out, typename _Comp>
_Out
__move_merge(_Iter __first1, _Iter __last1,
_Iter __first2, _Iter __last2,
_Out __result, _Comp __comp)
{
while (__first1 != __last1 && __first2 != __last2)
{
if (__comp(*__first2, *__first1))
{
*__result = ranges::iter_move(__first2);
++__first2;
}
else
{
*__result = ranges::iter_move(__first1);
++__first1;
}
++__result;
}
return ranges::move(__first2, __last2,
ranges::move(__first1, __last1, __result).out).out;
}
template<typename _Iter, typename _Out, typename _Distance, typename _Comp>
void
__merge_sort_loop(_Iter __first, _Iter __last, _Out __result,
_Distance __step_size, _Comp __comp)
{
const _Distance __two_step = 2 * __step_size;
while (__last - __first >= __two_step)
{
__result = __detail::__move_merge(__first, __first + __step_size,
__first + __step_size,
__first + __two_step,
__result, __comp);
__first += __two_step;
}
__step_size = ranges::min(_Distance(__last - __first), __step_size);
__detail::__move_merge(__first, __first + __step_size,
__first + __step_size, __last, __result, __comp);
}
template<typename _Iter, typename _Distance, typename _Compare>
constexpr void
__chunk_insertion_sort(_Iter __first, _Iter __last,
_Distance __chunk_size, _Compare __comp)
{
while (__last - __first >= __chunk_size)
{
__detail::__insertion_sort(__first, __first + __chunk_size, __comp);
__first += __chunk_size;
}
__detail::__insertion_sort(__first, __last, __comp);
}
template<typename _Iter, typename _Pointer, typename _Comp>
void
__merge_sort_with_buffer(_Iter __first, _Iter __last,
_Pointer __buffer, _Comp __comp)
{
using _Distance = iter_difference_t<_Iter>;
const _Distance __len = __last - __first;
const _Pointer __buffer_last = __buffer + ptrdiff_t(__len);
constexpr int __chunk_size = 7;
_Distance __step_size = __chunk_size;
__detail::__chunk_insertion_sort(__first, __last, __step_size, __comp);
while (__step_size < __len)
{
__detail::__merge_sort_loop(__first, __last, __buffer,
__step_size, __comp);
__step_size *= 2;
__detail::__merge_sort_loop(__buffer, __buffer_last, __first,
ptrdiff_t(__step_size), __comp);
__step_size *= 2;
}
}
template<typename _Iter, typename _Pointer, typename _Comp>
void
__merge_adaptive(_Iter __first, _Iter __middle, _Iter __last,
iter_difference_t<_Iter> __len1,
iter_difference_t<_Iter> __len2,
_Pointer __buffer, _Comp __comp); // defined near inplace_merge
template<typename _Iter, typename _Distance, typename _Pointer, typename _Comp>
void
__merge_adaptive_resize(_Iter __first, _Iter __middle, _Iter __last,
_Distance __len1, _Distance __len2,
_Pointer __buffer, _Distance __buffer_size,
_Comp __comp); // defined near inplace_merge
template<typename _Iter, typename _Distance, typename _Comp>
constexpr void
__merge_without_buffer(_Iter __first, _Iter __middle, _Iter __last,
_Distance __len1, _Distance __len2,
_Comp __comp); // defined near inplace_merge
template<typename _Iter, typename _Pointer, typename _Comp>
void
__stable_sort_adaptive(_Iter __first, _Iter __middle, _Iter __last,
_Pointer __buffer, _Comp __comp)
{
__detail::__merge_sort_with_buffer(__first, __middle, __buffer, __comp);
__detail::__merge_sort_with_buffer(__middle, __last, __buffer, __comp);
__detail::__merge_adaptive(__first, __middle, __last,
__middle - __first, __last - __middle,
__buffer, __comp);
}
template<typename _Iter, typename _Pointer, typename _Distance, typename _Comp>
void
__stable_sort_adaptive_resize(_Iter __first, _Iter __last,
_Pointer __buffer, _Distance __buffer_size,
_Comp __comp)
{
const _Distance __len = (__last - __first + 1) / 2;
const _Iter __middle = __first + __len;
if (__len > __buffer_size)
{
__detail::__stable_sort_adaptive_resize(__first, __middle, __buffer,
__buffer_size, __comp);
__detail::__stable_sort_adaptive_resize(__middle, __last, __buffer,
__buffer_size, __comp);
__detail::__merge_adaptive_resize(__first, __middle, __last,
_Distance(__middle - __first),
_Distance(__last - __middle),
__buffer, __buffer_size,
__comp);
}
else
__detail::__stable_sort_adaptive(__first, __middle, __last,
__buffer, __comp);
}
template<typename _Iter, typename _Comp>
constexpr void
__inplace_stable_sort(_Iter __first, _Iter __last, _Comp __comp)
{
if (__last - __first < 15)
{
__detail::__insertion_sort(__first, __last, __comp);
return;
}
_Iter __middle = __first + (__last - __first) / 2;
__detail::__inplace_stable_sort(__first, __middle, __comp);
__detail::__inplace_stable_sort(__middle, __last, __comp);
__detail::__merge_without_buffer(__first, __middle, __last,
__middle - __first,
__last - __middle,
__comp);
}
} // namespace __detail
struct __stable_sort_fn
{
template<random_access_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Comp = ranges::less, typename _Proj = identity>
requires sortable<_Iter, _Comp, _Proj>
_GLIBCXX26_CONSTEXPR
_Iter
operator()(_Iter __first, _Sent __last,
_Comp __comp = {}, _Proj __proj = {}) const
{
if constexpr (!same_as<_Iter, _Sent>)
return (*this)(__first, ranges::next(__first, __last),
std::move(__comp), std::move(__proj));
else
{
using _DistanceType = iter_difference_t<_Iter>;
if (__first == __last)
return __last;
auto __comp_proj = __detail::__make_comp_proj(__comp, __proj);
#if _GLIBCXX_HOSTED
# if __glibcxx_constexpr_algorithms >= 202306L // >= C++26
if consteval {
__detail::__inplace_stable_sort(__first, __last, __comp_proj);
return __last;
}
# endif
typedef _Temporary_buffer<_Iter, iter_value_t<_Iter>> _TmpBuf;
// __stable_sort_adaptive sorts the range in two halves,
// so the buffer only needs to fit half the range at once.
_TmpBuf __buf(__first, ptrdiff_t((__last - __first + 1) / 2));
if (__buf._M_requested_size() == __buf.size()) [[likely]]
__detail::__stable_sort_adaptive(__first,
__first + _DistanceType(__buf.size()),
__last, __buf.begin(), __comp_proj);
else if (__buf.begin()) [[unlikely]]
__detail::__inplace_stable_sort(__first, __last, __comp_proj);
else
__detail::__stable_sort_adaptive_resize(__first, __last, __buf.begin(),
_DistanceType(__buf.size()),
__comp_proj);
#else
__detail::__inplace_stable_sort(__first, __last, __comp_proj);
#endif
return __last;
}
}
template<random_access_range _Range,
typename _Comp = ranges::less, typename _Proj = identity>
requires sortable<iterator_t<_Range>, _Comp, _Proj>
_GLIBCXX26_CONSTEXPR
borrowed_iterator_t<_Range>
operator()(_Range&& __r, _Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__comp), std::move(__proj));
}
};
inline constexpr __stable_sort_fn stable_sort{};
struct __partial_sort_fn
{
template<random_access_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Comp = ranges::less, typename _Proj = identity>
requires sortable<_Iter, _Comp, _Proj>
constexpr _Iter
operator()(_Iter __first, _Iter __middle, _Sent __last,
_Comp __comp = {}, _Proj __proj = {}) const
{
if (__first == __middle)
return ranges::next(__first, __last);
ranges::make_heap(__first, __middle, __comp, __proj);
auto __i = __middle;
for (; __i != __last; ++__i)
if (std::__invoke(__comp,
std::__invoke(__proj, *__i),
std::__invoke(__proj, *__first)))
{
ranges::pop_heap(__first, __middle, __comp, __proj);
ranges::iter_swap(std::prev(__middle), __i);
ranges::push_heap(__first, __middle, __comp, __proj);
}
ranges::sort_heap(__first, __middle, __comp, __proj);
return __i;
}
template<random_access_range _Range,
typename _Comp = ranges::less, typename _Proj = identity>
requires sortable<iterator_t<_Range>, _Comp, _Proj>
constexpr borrowed_iterator_t<_Range>
operator()(_Range&& __r, iterator_t<_Range> __middle,
_Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), std::move(__middle),
ranges::end(__r),
std::move(__comp), std::move(__proj));
}
};
inline constexpr __partial_sort_fn partial_sort{};
template<typename _Iter, typename _Out>
using partial_sort_copy_result = in_out_result<_Iter, _Out>;
struct __partial_sort_copy_fn
{
template<input_iterator _Iter1, sentinel_for<_Iter1> _Sent1,
random_access_iterator _Iter2, sentinel_for<_Iter2> _Sent2,
typename _Comp = ranges::less,
typename _Proj1 = identity, typename _Proj2 = identity>
requires indirectly_copyable<_Iter1, _Iter2>
&& sortable<_Iter2, _Comp, _Proj2>
&& indirect_strict_weak_order<_Comp,
projected<_Iter1, _Proj1>,
projected<_Iter2, _Proj2>>
constexpr partial_sort_copy_result<_Iter1, _Iter2>
operator()(_Iter1 __first, _Sent1 __last,
_Iter2 __result_first, _Sent2 __result_last,
_Comp __comp = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
if (__result_first == __result_last)
{
// TODO: Eliminating the variable __lasti triggers an ICE.
auto __lasti = ranges::next(std::move(__first),
std::move(__last));
return {std::move(__lasti), std::move(__result_first)};
}
auto __result_real_last = __result_first;
while (__first != __last && __result_real_last != __result_last)
{
*__result_real_last = *__first;
++__result_real_last;
++__first;
}
ranges::make_heap(__result_first, __result_real_last, __comp, __proj2);
for (; __first != __last; ++__first)
if (std::__invoke(__comp,
std::__invoke(__proj1, *__first),
std::__invoke(__proj2, *__result_first)))
{
ranges::pop_heap(__result_first, __result_real_last,
__comp, __proj2);
*ranges::prev(__result_real_last) = *__first;
ranges::push_heap(__result_first, __result_real_last,
__comp, __proj2);
}
ranges::sort_heap(__result_first, __result_real_last, __comp, __proj2);
return {std::move(__first), std::move(__result_real_last)};
}
template<input_range _Range1, random_access_range _Range2,
typename _Comp = ranges::less,
typename _Proj1 = identity, typename _Proj2 = identity>
requires indirectly_copyable<iterator_t<_Range1>, iterator_t<_Range2>>
&& sortable<iterator_t<_Range2>, _Comp, _Proj2>
&& indirect_strict_weak_order<_Comp,
projected<iterator_t<_Range1>, _Proj1>,
projected<iterator_t<_Range2>, _Proj2>>
constexpr partial_sort_copy_result<borrowed_iterator_t<_Range1>,
borrowed_iterator_t<_Range2>>
operator()(_Range1&& __r, _Range2&& __out, _Comp __comp = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
ranges::begin(__out), ranges::end(__out),
std::move(__comp),
std::move(__proj1), std::move(__proj2));
}
};
inline constexpr __partial_sort_copy_fn partial_sort_copy{};
struct __is_sorted_until_fn
{
template<forward_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
indirect_strict_weak_order<projected<_Iter, _Proj>>
_Comp = ranges::less>
[[nodiscard]] constexpr _Iter
operator()(_Iter __first, _Sent __last,
_Comp __comp = {}, _Proj __proj = {}) const
{
if (__first == __last)
return __first;
auto __next = __first;
for (++__next; __next != __last; __first = __next, (void)++__next)
if (std::__invoke(__comp,
std::__invoke(__proj, *__next),
std::__invoke(__proj, *__first)))
return __next;
return __next;
}
template<forward_range _Range, typename _Proj = identity,
indirect_strict_weak_order<projected<iterator_t<_Range>, _Proj>>
_Comp = ranges::less>
[[nodiscard]] constexpr borrowed_iterator_t<_Range>
operator()(_Range&& __r, _Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__comp), std::move(__proj));
}
};
inline constexpr __is_sorted_until_fn is_sorted_until{};
struct __is_sorted_fn
{
template<forward_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
indirect_strict_weak_order<projected<_Iter, _Proj>>
_Comp = ranges::less>
[[nodiscard]] constexpr bool
operator()(_Iter __first, _Sent __last,
_Comp __comp = {}, _Proj __proj = {}) const
{
if (__first == __last)
return true;
auto __next = __first;
for (++__next; __next != __last; __first = __next, (void)++__next)
if (std::__invoke(__comp,
std::__invoke(__proj, *__next),
std::__invoke(__proj, *__first)))
return false;
return true;
}
template<forward_range _Range, typename _Proj = identity,
indirect_strict_weak_order<projected<iterator_t<_Range>, _Proj>>
_Comp = ranges::less>
[[nodiscard]] constexpr bool
operator()(_Range&& __r, _Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__comp), std::move(__proj));
}
};
inline constexpr __is_sorted_fn is_sorted{};
namespace __detail
{
template<typename _Iter, typename _Comp>
constexpr void
__introselect(_Iter __first, _Iter __nth, _Iter __last,
iter_difference_t<_Iter> __depth_limit, _Comp __comp)
{
while (__last - __first > 3)
{
if (__depth_limit == 0)
{
__detail::__heap_select(__first, ranges::next(__nth), __last,
__comp);
// Place the nth largest element in its final position.
ranges::iter_swap(__first, __nth);
return;
}
--__depth_limit;
_Iter __cut = __detail::__unguarded_partition_pivot(__first, __last, __comp);
if (__cut <= __nth)
__first = __cut;
else
__last = __cut;
}
__detail::__insertion_sort(__first, __last, __comp);
}
} // namespace __detail
struct __nth_element_fn
{
template<random_access_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Comp = ranges::less, typename _Proj = identity>
requires sortable<_Iter, _Comp, _Proj>
constexpr _Iter
operator()(_Iter __first, _Iter __nth, _Sent __last,
_Comp __comp = {}, _Proj __proj = {}) const
{
if constexpr (!same_as<_Iter, _Sent>)
return (*this)(__first, __nth, ranges::next(__first, __last),
std::move(__comp), std::move(__proj));
else
{
if (__first == __last || __nth == __last)
return __last;
auto __comp_proj = __detail::__make_comp_proj(__comp, __proj);
auto __n = __detail::__to_unsigned_like(__last - __first);
__detail::__introselect(__first, __nth, __last,
std::__bit_width(__n) * 2,
__comp_proj);
return __last;
}
}
template<random_access_range _Range,
typename _Comp = ranges::less, typename _Proj = identity>
requires sortable<iterator_t<_Range>, _Comp, _Proj>
constexpr borrowed_iterator_t<_Range>
operator()(_Range&& __r, iterator_t<_Range> __nth,
_Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), std::move(__nth),
ranges::end(__r), std::move(__comp), std::move(__proj));
}
};
inline constexpr __nth_element_fn nth_element{};
struct __lower_bound_fn
{
template<forward_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
typename _Tp _GLIBCXX26_RANGE_ALGO_DEF_VAL_T(_Iter, _Proj),
indirect_strict_weak_order<const _Tp*, projected<_Iter, _Proj>>
_Comp = ranges::less>
[[nodiscard]] constexpr _Iter
operator()(_Iter __first, _Sent __last,
const _Tp& __value, _Comp __comp = {}, _Proj __proj = {}) const
{
auto __len = ranges::distance(__first, __last);
while (__len > 0)
{
auto __half = __len / 2;
auto __middle = __first;
ranges::advance(__middle, __half);
if (std::__invoke(__comp, std::__invoke(__proj, *__middle), __value))
{
__first = __middle;
++__first;
__len = __len - __half - 1;
}
else
__len = __half;
}
return __first;
}
template<forward_range _Range,
typename _Proj = identity,
typename _Tp
_GLIBCXX26_RANGE_ALGO_DEF_VAL_T(iterator_t<_Range>, _Proj),
indirect_strict_weak_order<const _Tp*,
projected<iterator_t<_Range>, _Proj>>
_Comp = ranges::less>
[[nodiscard]] constexpr borrowed_iterator_t<_Range>
operator()(_Range&& __r,
const _Tp& __value, _Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
__value, std::move(__comp), std::move(__proj));
}
};
inline constexpr __lower_bound_fn lower_bound{};
struct __upper_bound_fn
{
template<forward_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
typename _Tp _GLIBCXX26_RANGE_ALGO_DEF_VAL_T(_Iter, _Proj),
indirect_strict_weak_order<const _Tp*, projected<_Iter, _Proj>>
_Comp = ranges::less>
[[nodiscard]] constexpr _Iter
operator()(_Iter __first, _Sent __last,
const _Tp& __value, _Comp __comp = {}, _Proj __proj = {}) const
{
auto __len = ranges::distance(__first, __last);
while (__len > 0)
{
auto __half = __len / 2;
auto __middle = __first;
ranges::advance(__middle, __half);
if (std::__invoke(__comp, __value, std::__invoke(__proj, *__middle)))
__len = __half;
else
{
__first = __middle;
++__first;
__len = __len - __half - 1;
}
}
return __first;
}
template<forward_range _Range,
typename _Proj = identity,
typename _Tp
_GLIBCXX26_RANGE_ALGO_DEF_VAL_T(iterator_t<_Range>, _Proj),
indirect_strict_weak_order<const _Tp*,
projected<iterator_t<_Range>, _Proj>>
_Comp = ranges::less>
[[nodiscard]] constexpr borrowed_iterator_t<_Range>
operator()(_Range&& __r,
const _Tp& __value, _Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
__value, std::move(__comp), std::move(__proj));
}
};
inline constexpr __upper_bound_fn upper_bound{};
struct __equal_range_fn
{
template<forward_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
typename _Tp _GLIBCXX26_RANGE_ALGO_DEF_VAL_T(_Iter, _Proj),
indirect_strict_weak_order<const _Tp*, projected<_Iter, _Proj>>
_Comp = ranges::less>
[[nodiscard]] constexpr subrange<_Iter>
operator()(_Iter __first, _Sent __last,
const _Tp& __value, _Comp __comp = {}, _Proj __proj = {}) const
{
auto __len = ranges::distance(__first, __last);
while (__len > 0)
{
auto __half = __len / 2;
auto __middle = __first;
ranges::advance(__middle, __half);
if (std::__invoke(__comp,
std::__invoke(__proj, *__middle),
__value))
{
__first = __middle;
++__first;
__len = __len - __half - 1;
}
else if (std::__invoke(__comp,
__value,
std::__invoke(__proj, *__middle)))
__len = __half;
else
{
auto __left
= ranges::lower_bound(__first, __middle,
__value, __comp, __proj);
ranges::advance(__first, __len);
auto __right
= ranges::upper_bound(++__middle, __first,
__value, __comp, __proj);
return {__left, __right};
}
}
return {__first, __first};
}
template<forward_range _Range,
typename _Proj = identity,
typename _Tp
_GLIBCXX26_RANGE_ALGO_DEF_VAL_T(iterator_t<_Range>, _Proj),
indirect_strict_weak_order<const _Tp*,
projected<iterator_t<_Range>, _Proj>>
_Comp = ranges::less>
[[nodiscard]] constexpr borrowed_subrange_t<_Range>
operator()(_Range&& __r, const _Tp& __value,
_Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
__value, std::move(__comp), std::move(__proj));
}
};
inline constexpr __equal_range_fn equal_range{};
struct __binary_search_fn
{
template<forward_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
typename _Tp _GLIBCXX26_RANGE_ALGO_DEF_VAL_T(_Iter, _Proj),
indirect_strict_weak_order<const _Tp*, projected<_Iter, _Proj>>
_Comp = ranges::less>
[[nodiscard]] constexpr bool
operator()(_Iter __first, _Sent __last,
const _Tp& __value, _Comp __comp = {}, _Proj __proj = {}) const
{
auto __i = ranges::lower_bound(__first, __last, __value, __comp, __proj);
if (__i == __last)
return false;
return !(bool)std::__invoke(__comp, __value,
std::__invoke(__proj, *__i));
}
template<forward_range _Range,
typename _Proj = identity,
typename _Tp
_GLIBCXX26_RANGE_ALGO_DEF_VAL_T(iterator_t<_Range>, _Proj),
indirect_strict_weak_order<const _Tp*,
projected<iterator_t<_Range>, _Proj>>
_Comp = ranges::less>
[[nodiscard]] constexpr bool
operator()(_Range&& __r, const _Tp& __value, _Comp __comp = {},
_Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
__value, std::move(__comp), std::move(__proj));
}
};
inline constexpr __binary_search_fn binary_search{};
struct __is_partitioned_fn
{
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
indirect_unary_predicate<projected<_Iter, _Proj>> _Pred>
[[nodiscard]] constexpr bool
operator()(_Iter __first, _Sent __last,
_Pred __pred, _Proj __proj = {}) const
{
__first = ranges::find_if_not(std::move(__first), __last,
__pred, __proj);
if (__first == __last)
return true;
++__first;
return ranges::none_of(std::move(__first), std::move(__last),
std::move(__pred), std::move(__proj));
}
template<input_range _Range, typename _Proj = identity,
indirect_unary_predicate<projected<iterator_t<_Range>, _Proj>>
_Pred>
[[nodiscard]] constexpr bool
operator()(_Range&& __r, _Pred __pred, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__pred), std::move(__proj));
}
};
inline constexpr __is_partitioned_fn is_partitioned{};
struct __partition_fn
{
template<permutable _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
indirect_unary_predicate<projected<_Iter, _Proj>> _Pred>
constexpr subrange<_Iter>
operator()(_Iter __first, _Sent __last,
_Pred __pred, _Proj __proj = {}) const
{
if constexpr (bidirectional_iterator<_Iter>)
{
auto __lasti = ranges::next(__first, __last);
auto __tail = __lasti;
for (;;)
{
for (;;)
if (__first == __tail)
return {std::move(__first), std::move(__lasti)};
else if (std::__invoke(__pred,
std::__invoke(__proj, *__first)))
++__first;
else
break;
--__tail;
for (;;)
if (__first == __tail)
return {std::move(__first), std::move(__lasti)};
else if (!(bool)std::__invoke(__pred,
std::__invoke(__proj, *__tail)))
--__tail;
else
break;
ranges::iter_swap(__first, __tail);
++__first;
}
}
else
{
if (__first == __last)
return {__first, __first};
while (std::__invoke(__pred, std::__invoke(__proj, *__first)))
if (++__first == __last)
return {__first, __first};
auto __next = __first;
while (++__next != __last)
if (std::__invoke(__pred, std::__invoke(__proj, *__next)))
{
ranges::iter_swap(__first, __next);
++__first;
}
return {std::move(__first), std::move(__next)};
}
}
template<forward_range _Range, typename _Proj = identity,
indirect_unary_predicate<projected<iterator_t<_Range>, _Proj>>
_Pred>
requires permutable<iterator_t<_Range>>
constexpr borrowed_subrange_t<_Range>
operator()(_Range&& __r, _Pred __pred, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__pred), std::move(__proj));
}
};
inline constexpr __partition_fn partition{};
#if _GLIBCXX_HOSTED
namespace __detail
{
// Like find_if_not(), but uses and updates a count of the
// remaining range length instead of comparing against an end
// iterator.
template<typename _Iter, typename _Pred, typename _Distance>
constexpr _Iter
__find_if_not_n(_Iter __first, _Distance& __len, _Pred __pred)
{
for (; __len; --__len, (void) ++__first)
if (!__pred(*__first))
break;
return __first;
}
template<typename _Iter, typename _Sent, typename _Pointer,
typename _Pred, typename _Distance>
constexpr subrange<_Iter>
__stable_partition_adaptive(_Iter __first, _Sent __last,
_Pred __pred, _Distance __len,
_Pointer __buffer,
_Distance __buffer_size)
{
if (__len == 1)
return {__first, ranges::next(__first, 1)};
if (__len <= __buffer_size)
{
_Iter __result1 = __first;
_Pointer __result2 = __buffer;
// The precondition guarantees that !__pred(__first), so
// move that element to the buffer before starting the loop.
// This ensures that we only call __pred once per element.
*__result2 = ranges::iter_move(__first);
++__result2;
++__first;
for (; __first != __last; ++__first)
if (__pred(*__first))
{
*__result1 = ranges::iter_move(__first);
++__result1;
}
else
{
*__result2 = ranges::iter_move(__first);
++__result2;
}
ranges::move(__buffer, __result2, __result1);
return {__result1, __first};
}
_Iter __middle = __first;
ranges::advance(__middle, __len / 2);
_Iter __left_split
= __detail::__stable_partition_adaptive(__first, __middle, __pred,
__len / 2, __buffer,
__buffer_size).begin();
// Advance past true-predicate values to satisfy this
// function's preconditions.
_Distance __right_len = __len - __len / 2;
_Iter __right_split = __detail::__find_if_not_n(__middle, __right_len, __pred);
if (__right_len)
__right_split
= __detail::__stable_partition_adaptive(__right_split, __last, __pred,
__right_len, __buffer, __buffer_size).begin();
return ranges::rotate(__left_split, __middle, __right_split);
}
} // namespace __detail
struct __stable_partition_fn
{
template<bidirectional_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
indirect_unary_predicate<projected<_Iter, _Proj>> _Pred>
requires permutable<_Iter>
_GLIBCXX26_CONSTEXPR
subrange<_Iter>
operator()(_Iter __first, _Sent __last,
_Pred __pred, _Proj __proj = {}) const
{
__first = ranges::find_if_not(__first, __last, __pred, __proj);
if (__first == __last)
return {__first, __first};
using _DistanceType = iter_difference_t<_Iter>;
const _DistanceType __len = ranges::distance(__first, __last);
auto __pred_proj = __detail::__make_pred_proj(__pred, __proj);
#if __glibcxx_constexpr_algorithms >= 202306L // >= C++26
if consteval {
// Simulate a _Temporary_buffer of length 1:
iter_value_t<_Iter> __buf = ranges::iter_move(__first);
*__first = std::move(__buf);
return __detail::__stable_partition_adaptive(__first, __last,
__pred_proj,
__len, &__buf,
_DistanceType(1));
}
#endif
_Temporary_buffer<_Iter, iter_value_t<_Iter>> __buf(__first, ptrdiff_t(__len));
return __detail::__stable_partition_adaptive(__first, __last,
__pred_proj,
__len, __buf.begin(),
_DistanceType(__buf.size()));
}
template<bidirectional_range _Range, typename _Proj = identity,
indirect_unary_predicate<projected<iterator_t<_Range>, _Proj>>
_Pred>
requires permutable<iterator_t<_Range>>
_GLIBCXX26_CONSTEXPR
borrowed_subrange_t<_Range>
operator()(_Range&& __r, _Pred __pred, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__pred), std::move(__proj));
}
};
inline constexpr __stable_partition_fn stable_partition{};
#endif
template<typename _Iter, typename _Out1, typename _Out2>
struct in_out_out_result
{
[[no_unique_address]] _Iter in;
[[no_unique_address]] _Out1 out1;
[[no_unique_address]] _Out2 out2;
template<typename _IIter, typename _OOut1, typename _OOut2>
requires convertible_to<const _Iter&, _IIter>
&& convertible_to<const _Out1&, _OOut1>
&& convertible_to<const _Out2&, _OOut2>
constexpr
operator in_out_out_result<_IIter, _OOut1, _OOut2>() const &
{ return {in, out1, out2}; }
template<typename _IIter, typename _OOut1, typename _OOut2>
requires convertible_to<_Iter, _IIter>
&& convertible_to<_Out1, _OOut1>
&& convertible_to<_Out2, _OOut2>
constexpr
operator in_out_out_result<_IIter, _OOut1, _OOut2>() &&
{ return {std::move(in), std::move(out1), std::move(out2)}; }
};
template<typename _Iter, typename _Out1, typename _Out2>
using partition_copy_result = in_out_out_result<_Iter, _Out1, _Out2>;
struct __partition_copy_fn
{
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
weakly_incrementable _Out1, weakly_incrementable _Out2,
typename _Proj = identity,
indirect_unary_predicate<projected<_Iter, _Proj>> _Pred>
requires indirectly_copyable<_Iter, _Out1>
&& indirectly_copyable<_Iter, _Out2>
constexpr partition_copy_result<_Iter, _Out1, _Out2>
operator()(_Iter __first, _Sent __last,
_Out1 __out_true, _Out2 __out_false,
_Pred __pred, _Proj __proj = {}) const
{
for (; __first != __last; ++__first)
if (std::__invoke(__pred, std::__invoke(__proj, *__first)))
{
*__out_true = *__first;
++__out_true;
}
else
{
*__out_false = *__first;
++__out_false;
}
return {std::move(__first),
std::move(__out_true), std::move(__out_false)};
}
template<input_range _Range, weakly_incrementable _Out1,
weakly_incrementable _Out2,
typename _Proj = identity,
indirect_unary_predicate<projected<iterator_t<_Range>, _Proj>>
_Pred>
requires indirectly_copyable<iterator_t<_Range>, _Out1>
&& indirectly_copyable<iterator_t<_Range>, _Out2>
constexpr partition_copy_result<borrowed_iterator_t<_Range>, _Out1, _Out2>
operator()(_Range&& __r, _Out1 __out_true, _Out2 __out_false,
_Pred __pred, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__out_true), std::move(__out_false),
std::move(__pred), std::move(__proj));
}
};
inline constexpr __partition_copy_fn partition_copy{};
struct __partition_point_fn
{
template<forward_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
indirect_unary_predicate<projected<_Iter, _Proj>> _Pred>
[[nodiscard]] constexpr _Iter
operator()(_Iter __first, _Sent __last,
_Pred __pred, _Proj __proj = {}) const
{
auto __len = ranges::distance(__first, __last);
while (__len > 0)
{
auto __half = __len / 2;
auto __middle = __first;
ranges::advance(__middle, __half);
if (std::__invoke(__pred, std::__invoke(__proj, *__middle)))
{
__first = __middle;
++__first;
__len = __len - __half - 1;
}
else
__len = __half;
}
return __first;
}
template<forward_range _Range, typename _Proj = identity,
indirect_unary_predicate<projected<iterator_t<_Range>, _Proj>>
_Pred>
[[nodiscard]] constexpr borrowed_iterator_t<_Range>
operator()(_Range&& __r, _Pred __pred, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__pred), std::move(__proj));
}
};
inline constexpr __partition_point_fn partition_point{};
template<typename _Iter1, typename _Iter2, typename _Out>
using merge_result = in_in_out_result<_Iter1, _Iter2, _Out>;
struct __merge_fn
{
template<input_iterator _Iter1, sentinel_for<_Iter1> _Sent1,
input_iterator _Iter2, sentinel_for<_Iter2> _Sent2,
weakly_incrementable _Out, typename _Comp = ranges::less,
typename _Proj1 = identity, typename _Proj2 = identity>
requires mergeable<_Iter1, _Iter2, _Out, _Comp, _Proj1, _Proj2>
constexpr merge_result<_Iter1, _Iter2, _Out>
operator()(_Iter1 __first1, _Sent1 __last1,
_Iter2 __first2, _Sent2 __last2, _Out __result,
_Comp __comp = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
while (__first1 != __last1 && __first2 != __last2)
{
if (std::__invoke(__comp,
std::__invoke(__proj2, *__first2),
std::__invoke(__proj1, *__first1)))
{
*__result = *__first2;
++__first2;
}
else
{
*__result = *__first1;
++__first1;
}
++__result;
}
auto __copy1 = ranges::copy(std::move(__first1), std::move(__last1),
std::move(__result));
auto __copy2 = ranges::copy(std::move(__first2), std::move(__last2),
std::move(__copy1.out));
return { std::move(__copy1.in), std::move(__copy2.in),
std::move(__copy2.out) };
}
template<input_range _Range1, input_range _Range2, weakly_incrementable _Out,
typename _Comp = ranges::less,
typename _Proj1 = identity, typename _Proj2 = identity>
requires mergeable<iterator_t<_Range1>, iterator_t<_Range2>, _Out,
_Comp, _Proj1, _Proj2>
constexpr merge_result<borrowed_iterator_t<_Range1>,
borrowed_iterator_t<_Range2>,
_Out>
operator()(_Range1&& __r1, _Range2&& __r2, _Out __result,
_Comp __comp = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
return (*this)(ranges::begin(__r1), ranges::end(__r1),
ranges::begin(__r2), ranges::end(__r2),
std::move(__result), std::move(__comp),
std::move(__proj1), std::move(__proj2));
}
};
inline constexpr __merge_fn merge{};
namespace __detail
{
template<typename _Iter1, typename _Iter2, typename _Out, typename _Comp>
void
__move_merge_adaptive(_Iter1 __first1, _Iter1 __last1,
_Iter2 __first2, _Iter2 __last2,
_Out __result, _Comp __comp)
{
while (__first1 != __last1 && __first2 != __last2)
{
if (__comp(*__first2, *__first1))
{
*__result = ranges::iter_move(__first2);
++__first2;
}
else
{
*__result = ranges::iter_move(__first1);
++__first1;
}
++__result;
}
if (__first1 != __last1)
ranges::move(__first1, __last1, __result);
}
template<typename _Iter1, typename _Iter2, typename _Iter3, typename _Comp>
void
__move_merge_adaptive_backward(_Iter1 __first1, _Iter1 __last1,
_Iter2 __first2, _Iter2 __last2,
_Iter3 __result, _Comp __comp)
{
if (__first1 == __last1)
{
ranges::move_backward(__first2, __last2, __result);
return;
}
else if (__first2 == __last2)
return;
--__last1;
--__last2;
while (true)
{
if (__comp(*__last2, *__last1))
{
*--__result = ranges::iter_move(__last1);
if (__first1 == __last1)
{
ranges::move_backward(__first2, ++__last2, __result);
return;
}
--__last1;
}
else
{
*--__result = ranges::iter_move(__last2);
if (__first2 == __last2)
return;
--__last2;
}
}
}
template<typename _Iter1, typename _Iter2>
_Iter1
__rotate_adaptive(_Iter1 __first, _Iter1 __middle, _Iter1 __last,
iter_difference_t<_Iter1> __len1,
iter_difference_t<_Iter1> __len2,
_Iter2 __buffer,
iter_difference_t<_Iter1> __buffer_size)
{
_Iter2 __buffer_end;
if (__len1 > __len2 && __len2 <= __buffer_size)
{
if (__len2)
{
__buffer_end = ranges::move(__middle, __last, __buffer).out;
ranges::move_backward(__first, __middle, __last);
return ranges::move(__buffer, __buffer_end, __first).out;
}
else
return __first;
}
else if (__len1 <= __buffer_size)
{
if (__len1)
{
__buffer_end = ranges::move(__first, __middle, __buffer).out;
ranges::move(__middle, __last, __first);
return ranges::move_backward(__buffer, __buffer_end, __last).out;
}
else
return __last;
}
else
return ranges::rotate(__first, __middle, __last).begin();
}
template<typename _Iter, typename _Pointer, typename _Comp>
void
__merge_adaptive(_Iter __first, _Iter __middle, _Iter __last,
iter_difference_t<_Iter> __len1,
iter_difference_t<_Iter> __len2,
_Pointer __buffer, _Comp __comp)
{
if (__len1 <= __len2)
{
_Pointer __buffer_end = ranges::move(__first, __middle, __buffer).out;
__detail::__move_merge_adaptive(__buffer, __buffer_end, __middle, __last,
__first, __comp);
}
else
{
_Pointer __buffer_end = ranges::move(__middle, __last, __buffer).out;
__detail::__move_merge_adaptive_backward(__first, __middle, __buffer,
__buffer_end, __last, __comp);
}
}
template<typename _Iter, typename _Distance, typename _Pointer, typename _Comp>
void
__merge_adaptive_resize(_Iter __first, _Iter __middle, _Iter __last,
_Distance __len1, _Distance __len2,
_Pointer __buffer, _Distance __buffer_size,
_Comp __comp)
{
if (__len1 <= __buffer_size || __len2 <= __buffer_size)
__detail::__merge_adaptive(__first, __middle, __last,
__len1, __len2, __buffer, __comp);
else
{
_Iter __first_cut = __first;
_Iter __second_cut = __middle;
_Distance __len11 = 0;
_Distance __len22 = 0;
if (__len1 > __len2)
{
__len11 = __len1 / 2;
ranges::advance(__first_cut, __len11);
__second_cut = ranges::lower_bound(__middle, __last, *__first_cut,
__comp);
__len22 = ranges::distance(__middle, __second_cut);
}
else
{
__len22 = __len2 / 2;
ranges::advance(__second_cut, __len22);
__first_cut = ranges::upper_bound(__first, __middle, *__second_cut,
__comp);
__len11 = ranges::distance(__first, __first_cut);
}
_Iter __new_middle
= __detail::__rotate_adaptive(__first_cut, __middle, __second_cut,
_Distance(__len1 - __len11), __len22,
__buffer, __buffer_size);
__detail::__merge_adaptive_resize(__first, __first_cut, __new_middle,
__len11, __len22,
__buffer, __buffer_size, __comp);
__detail::__merge_adaptive_resize(__new_middle, __second_cut, __last,
_Distance(__len1 - __len11),
_Distance(__len2 - __len22),
__buffer, __buffer_size, __comp);
}
}
template<typename _Iter, typename _Distance, typename _Comp>
constexpr void
__merge_without_buffer(_Iter __first, _Iter __middle, _Iter __last,
_Distance __len1, _Distance __len2, _Comp __comp)
{
if (__len1 == 0 || __len2 == 0)
return;
if (__len1 + __len2 == 2)
{
if (__comp(*__middle, *__first))
ranges::iter_swap(__first, __middle);
return;
}
_Iter __first_cut = __first;
_Iter __second_cut = __middle;
_Distance __len11 = 0;
_Distance __len22 = 0;
if (__len1 > __len2)
{
__len11 = __len1 / 2;
ranges::advance(__first_cut, __len11);
__second_cut = ranges::lower_bound(__middle, __last, *__first_cut, __comp);
__len22 = ranges::distance(__middle, __second_cut);
}
else
{
__len22 = __len2 / 2;
ranges::advance(__second_cut, __len22);
__first_cut = ranges::upper_bound(__first, __middle, *__second_cut, __comp);
__len11 = ranges::distance(__first, __first_cut);
}
_Iter __new_middle = ranges::rotate(__first_cut, __middle, __second_cut).begin();
__detail::__merge_without_buffer(__first, __first_cut, __new_middle,
__len11, __len22, __comp);
__detail::__merge_without_buffer(__new_middle, __second_cut, __last,
__len1 - __len11, __len2 - __len22, __comp);
}
} // namespace __detail
struct __inplace_merge_fn
{
template<bidirectional_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Comp = ranges::less,
typename _Proj = identity>
requires sortable<_Iter, _Comp, _Proj>
_GLIBCXX26_CONSTEXPR
_Iter
operator()(_Iter __first, _Iter __middle, _Sent __last,
_Comp __comp = {}, _Proj __proj = {}) const
{
if constexpr (!same_as<_Iter, _Sent>)
return (*this)(__first, __middle, ranges::next(__middle, __last),
std::move(__comp), std::move(__proj));
else
{
using _DistanceType = iter_difference_t<_Iter>;
if (__first == __middle || __middle == __last)
return __last;
const _DistanceType __len1 = ranges::distance(__first, __middle);
const _DistanceType __len2 = ranges::distance(__middle, __last);
auto __comp_proj = __detail::__make_comp_proj(__comp, __proj);
#if _GLIBCXX_HOSTED
# if __glibcxx_constexpr_algorithms >= 202306L // >= C++26
if consteval {
__detail::__merge_without_buffer(__first, __middle, __last,
__len1, __len2, __comp_proj);
return __last;
}
# endif
using _TmpBuf = _Temporary_buffer<_Iter, iter_value_t<_Iter>>;
// __merge_adaptive will use a buffer for the smaller of
// [first,middle) and [middle,last).
_TmpBuf __buf(__first, ptrdiff_t(ranges::min(__len1, __len2)));
if (__buf.size() == __buf._M_requested_size()) [[likely]]
__detail::__merge_adaptive
(__first, __middle, __last, __len1, __len2, __buf.begin(), __comp_proj);
else if (__buf.begin() == 0) [[unlikely]]
__detail::__merge_without_buffer
(__first, __middle, __last, __len1, __len2, __comp_proj);
else
__detail::__merge_adaptive_resize
(__first, __middle, __last, __len1, __len2, __buf.begin(),
_DistanceType(__buf.size()), __comp_proj);
#else
__detail::__merge_without_buffer
(__first, __middle, __last, __len1, __len2, __comp_proj);
#endif
return __last;
}
}
template<bidirectional_range _Range,
typename _Comp = ranges::less, typename _Proj = identity>
requires sortable<iterator_t<_Range>, _Comp, _Proj>
_GLIBCXX26_CONSTEXPR
borrowed_iterator_t<_Range>
operator()(_Range&& __r, iterator_t<_Range> __middle,
_Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), std::move(__middle),
ranges::end(__r),
std::move(__comp), std::move(__proj));
}
};
inline constexpr __inplace_merge_fn inplace_merge{};
struct __includes_fn
{
template<input_iterator _Iter1, sentinel_for<_Iter1> _Sent1,
input_iterator _Iter2, sentinel_for<_Iter2> _Sent2,
typename _Proj1 = identity, typename _Proj2 = identity,
indirect_strict_weak_order<projected<_Iter1, _Proj1>,
projected<_Iter2, _Proj2>>
_Comp = ranges::less>
[[nodiscard]] constexpr bool
operator()(_Iter1 __first1, _Sent1 __last1,
_Iter2 __first2, _Sent2 __last2,
_Comp __comp = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
while (__first1 != __last1 && __first2 != __last2)
if (std::__invoke(__comp,
std::__invoke(__proj2, *__first2),
std::__invoke(__proj1, *__first1)))
return false;
else if (std::__invoke(__comp,
std::__invoke(__proj1, *__first1),
std::__invoke(__proj2, *__first2)))
++__first1;
else
{
++__first1;
++__first2;
}
return __first2 == __last2;
}
template<input_range _Range1, input_range _Range2,
typename _Proj1 = identity, typename _Proj2 = identity,
indirect_strict_weak_order<projected<iterator_t<_Range1>, _Proj1>,
projected<iterator_t<_Range2>, _Proj2>>
_Comp = ranges::less>
[[nodiscard]] constexpr bool
operator()(_Range1&& __r1, _Range2&& __r2, _Comp __comp = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
return (*this)(ranges::begin(__r1), ranges::end(__r1),
ranges::begin(__r2), ranges::end(__r2),
std::move(__comp),
std::move(__proj1), std::move(__proj2));
}
};
inline constexpr __includes_fn includes{};
template<typename _Iter1, typename _Iter2, typename _Out>
using set_union_result = in_in_out_result<_Iter1, _Iter2, _Out>;
struct __set_union_fn
{
template<input_iterator _Iter1, sentinel_for<_Iter1> _Sent1,
input_iterator _Iter2, sentinel_for<_Iter2> _Sent2,
weakly_incrementable _Out, typename _Comp = ranges::less,
typename _Proj1 = identity, typename _Proj2 = identity>
requires mergeable<_Iter1, _Iter2, _Out, _Comp, _Proj1, _Proj2>
constexpr set_union_result<_Iter1, _Iter2, _Out>
operator()(_Iter1 __first1, _Sent1 __last1,
_Iter2 __first2, _Sent2 __last2,
_Out __result, _Comp __comp = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
while (__first1 != __last1 && __first2 != __last2)
{
if (std::__invoke(__comp,
std::__invoke(__proj1, *__first1),
std::__invoke(__proj2, *__first2)))
{
*__result = *__first1;
++__first1;
}
else if (std::__invoke(__comp,
std::__invoke(__proj2, *__first2),
std::__invoke(__proj1, *__first1)))
{
*__result = *__first2;
++__first2;
}
else
{
*__result = *__first1;
++__first1;
++__first2;
}
++__result;
}
auto __copy1 = ranges::copy(std::move(__first1), std::move(__last1),
std::move(__result));
auto __copy2 = ranges::copy(std::move(__first2), std::move(__last2),
std::move(__copy1.out));
return {std::move(__copy1.in), std::move(__copy2.in),
std::move(__copy2.out)};
}
template<input_range _Range1, input_range _Range2, weakly_incrementable _Out,
typename _Comp = ranges::less,
typename _Proj1 = identity, typename _Proj2 = identity>
requires mergeable<iterator_t<_Range1>, iterator_t<_Range2>, _Out,
_Comp, _Proj1, _Proj2>
constexpr set_union_result<borrowed_iterator_t<_Range1>,
borrowed_iterator_t<_Range2>, _Out>
operator()(_Range1&& __r1, _Range2&& __r2,
_Out __result, _Comp __comp = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
return (*this)(ranges::begin(__r1), ranges::end(__r1),
ranges::begin(__r2), ranges::end(__r2),
std::move(__result), std::move(__comp),
std::move(__proj1), std::move(__proj2));
}
};
inline constexpr __set_union_fn set_union{};
template<typename _Iter1, typename _Iter2, typename _Out>
using set_intersection_result = in_in_out_result<_Iter1, _Iter2, _Out>;
struct __set_intersection_fn
{
template<input_iterator _Iter1, sentinel_for<_Iter1> _Sent1,
input_iterator _Iter2, sentinel_for<_Iter2> _Sent2,
weakly_incrementable _Out, typename _Comp = ranges::less,
typename _Proj1 = identity, typename _Proj2 = identity>
requires mergeable<_Iter1, _Iter2, _Out, _Comp, _Proj1, _Proj2>
constexpr set_intersection_result<_Iter1, _Iter2, _Out>
operator()(_Iter1 __first1, _Sent1 __last1,
_Iter2 __first2, _Sent2 __last2, _Out __result,
_Comp __comp = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
while (__first1 != __last1 && __first2 != __last2)
if (std::__invoke(__comp,
std::__invoke(__proj1, *__first1),
std::__invoke(__proj2, *__first2)))
++__first1;
else if (std::__invoke(__comp,
std::__invoke(__proj2, *__first2),
std::__invoke(__proj1, *__first1)))
++__first2;
else
{
*__result = *__first1;
++__first1;
++__first2;
++__result;
}
// TODO: Eliminating these variables triggers an ICE.
auto __last1i = ranges::next(std::move(__first1), std::move(__last1));
auto __last2i = ranges::next(std::move(__first2), std::move(__last2));
return {std::move(__last1i), std::move(__last2i), std::move(__result)};
}
template<input_range _Range1, input_range _Range2, weakly_incrementable _Out,
typename _Comp = ranges::less,
typename _Proj1 = identity, typename _Proj2 = identity>
requires mergeable<iterator_t<_Range1>, iterator_t<_Range2>, _Out,
_Comp, _Proj1, _Proj2>
constexpr set_intersection_result<borrowed_iterator_t<_Range1>,
borrowed_iterator_t<_Range2>, _Out>
operator()(_Range1&& __r1, _Range2&& __r2, _Out __result,
_Comp __comp = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
return (*this)(ranges::begin(__r1), ranges::end(__r1),
ranges::begin(__r2), ranges::end(__r2),
std::move(__result), std::move(__comp),
std::move(__proj1), std::move(__proj2));
}
};
inline constexpr __set_intersection_fn set_intersection{};
template<typename _Iter, typename _Out>
using set_difference_result = in_out_result<_Iter, _Out>;
struct __set_difference_fn
{
template<input_iterator _Iter1, sentinel_for<_Iter1> _Sent1,
input_iterator _Iter2, sentinel_for<_Iter2> _Sent2,
weakly_incrementable _Out, typename _Comp = ranges::less,
typename _Proj1 = identity, typename _Proj2 = identity>
requires mergeable<_Iter1, _Iter2, _Out, _Comp, _Proj1, _Proj2>
constexpr set_difference_result<_Iter1, _Out>
operator()(_Iter1 __first1, _Sent1 __last1,
_Iter2 __first2, _Sent2 __last2, _Out __result,
_Comp __comp = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
while (__first1 != __last1 && __first2 != __last2)
if (std::__invoke(__comp,
std::__invoke(__proj1, *__first1),
std::__invoke(__proj2, *__first2)))
{
*__result = *__first1;
++__first1;
++__result;
}
else if (std::__invoke(__comp,
std::__invoke(__proj2, *__first2),
std::__invoke(__proj1, *__first1)))
++__first2;
else
{
++__first1;
++__first2;
}
return ranges::copy(std::move(__first1), std::move(__last1),
std::move(__result));
}
template<input_range _Range1, input_range _Range2, weakly_incrementable _Out,
typename _Comp = ranges::less,
typename _Proj1 = identity, typename _Proj2 = identity>
requires mergeable<iterator_t<_Range1>, iterator_t<_Range2>, _Out,
_Comp, _Proj1, _Proj2>
constexpr set_difference_result<borrowed_iterator_t<_Range1>, _Out>
operator()(_Range1&& __r1, _Range2&& __r2, _Out __result,
_Comp __comp = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
return (*this)(ranges::begin(__r1), ranges::end(__r1),
ranges::begin(__r2), ranges::end(__r2),
std::move(__result), std::move(__comp),
std::move(__proj1), std::move(__proj2));
}
};
inline constexpr __set_difference_fn set_difference{};
template<typename _Iter1, typename _Iter2, typename _Out>
using set_symmetric_difference_result
= in_in_out_result<_Iter1, _Iter2, _Out>;
struct __set_symmetric_difference_fn
{
template<input_iterator _Iter1, sentinel_for<_Iter1> _Sent1,
input_iterator _Iter2, sentinel_for<_Iter2> _Sent2,
weakly_incrementable _Out, typename _Comp = ranges::less,
typename _Proj1 = identity, typename _Proj2 = identity>
requires mergeable<_Iter1, _Iter2, _Out, _Comp, _Proj1, _Proj2>
constexpr set_symmetric_difference_result<_Iter1, _Iter2, _Out>
operator()(_Iter1 __first1, _Sent1 __last1,
_Iter2 __first2, _Sent2 __last2,
_Out __result, _Comp __comp = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
while (__first1 != __last1 && __first2 != __last2)
if (std::__invoke(__comp,
std::__invoke(__proj1, *__first1),
std::__invoke(__proj2, *__first2)))
{
*__result = *__first1;
++__first1;
++__result;
}
else if (std::__invoke(__comp,
std::__invoke(__proj2, *__first2),
std::__invoke(__proj1, *__first1)))
{
*__result = *__first2;
++__first2;
++__result;
}
else
{
++__first1;
++__first2;
}
auto __copy1 = ranges::copy(std::move(__first1), std::move(__last1),
std::move(__result));
auto __copy2 = ranges::copy(std::move(__first2), std::move(__last2),
std::move(__copy1.out));
return {std::move(__copy1.in), std::move(__copy2.in),
std::move(__copy2.out)};
}
template<input_range _Range1, input_range _Range2, weakly_incrementable _Out,
typename _Comp = ranges::less,
typename _Proj1 = identity, typename _Proj2 = identity>
requires mergeable<iterator_t<_Range1>, iterator_t<_Range2>, _Out,
_Comp, _Proj1, _Proj2>
constexpr set_symmetric_difference_result<borrowed_iterator_t<_Range1>,
borrowed_iterator_t<_Range2>,
_Out>
operator()(_Range1&& __r1, _Range2&& __r2, _Out __result,
_Comp __comp = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
return (*this)(ranges::begin(__r1), ranges::end(__r1),
ranges::begin(__r2), ranges::end(__r2),
std::move(__result), std::move(__comp),
std::move(__proj1), std::move(__proj2));
}
};
inline constexpr __set_symmetric_difference_fn set_symmetric_difference{};
// min is defined in <bits/ranges_util.h>.
struct __max_fn
{
template<typename _Tp, typename _Proj = identity,
indirect_strict_weak_order<projected<const _Tp*, _Proj>>
_Comp = ranges::less>
[[nodiscard]] constexpr const _Tp&
operator()(const _Tp& __a, const _Tp& __b,
_Comp __comp = {}, _Proj __proj = {}) const
{
if (std::__invoke(__comp,
std::__invoke(__proj, __a),
std::__invoke(__proj, __b)))
return __b;
else
return __a;
}
template<input_range _Range, typename _Proj = identity,
indirect_strict_weak_order<projected<iterator_t<_Range>, _Proj>>
_Comp = ranges::less>
requires indirectly_copyable_storable<iterator_t<_Range>,
range_value_t<_Range>*>
[[nodiscard]] constexpr range_value_t<_Range>
operator()(_Range&& __r, _Comp __comp = {}, _Proj __proj = {}) const
{
auto __first = ranges::begin(__r);
auto __last = ranges::end(__r);
__glibcxx_assert(__first != __last);
range_value_t<_Range> __result(*__first);
while (++__first != __last)
{
auto&& __tmp = *__first;
if (std::__invoke(__comp,
std::__invoke(__proj, __result),
std::__invoke(__proj, __tmp)))
__result = std::forward<decltype(__tmp)>(__tmp);
}
return __result;
}
template<copyable _Tp, typename _Proj = identity,
indirect_strict_weak_order<projected<const _Tp*, _Proj>>
_Comp = ranges::less>
[[nodiscard]] constexpr _Tp
operator()(initializer_list<_Tp> __r,
_Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::subrange(__r),
std::move(__comp), std::move(__proj));
}
};
inline constexpr __max_fn max{};
struct __clamp_fn
{
template<typename _Tp, typename _Proj = identity,
indirect_strict_weak_order<projected<const _Tp*, _Proj>> _Comp
= ranges::less>
[[nodiscard]] constexpr const _Tp&
operator()(const _Tp& __val, const _Tp& __lo, const _Tp& __hi,
_Comp __comp = {}, _Proj __proj = {}) const
{
__glibcxx_assert(!(std::__invoke(__comp,
std::__invoke(__proj, __hi),
std::__invoke(__proj, __lo))));
auto&& __proj_val = std::__invoke(__proj, __val);
if (std::__invoke(__comp,
std::forward<decltype(__proj_val)>(__proj_val),
std::__invoke(__proj, __lo)))
return __lo;
else if (std::__invoke(__comp,
std::__invoke(__proj, __hi),
std::forward<decltype(__proj_val)>(__proj_val)))
return __hi;
else
return __val;
}
};
inline constexpr __clamp_fn clamp{};
template<typename _Tp>
struct min_max_result
{
[[no_unique_address]] _Tp min;
[[no_unique_address]] _Tp max;
template<typename _Tp2>
requires convertible_to<const _Tp&, _Tp2>
constexpr
operator min_max_result<_Tp2>() const &
{ return {min, max}; }
template<typename _Tp2>
requires convertible_to<_Tp, _Tp2>
constexpr
operator min_max_result<_Tp2>() &&
{ return {std::move(min), std::move(max)}; }
};
template<typename _Tp>
using minmax_result = min_max_result<_Tp>;
struct __minmax_fn
{
template<typename _Tp, typename _Proj = identity,
indirect_strict_weak_order<projected<const _Tp*, _Proj>>
_Comp = ranges::less>
[[nodiscard]] constexpr minmax_result<const _Tp&>
operator()(const _Tp& __a, const _Tp& __b,
_Comp __comp = {}, _Proj __proj = {}) const
{
if (std::__invoke(__comp,
std::__invoke(__proj, __b),
std::__invoke(__proj, __a)))
return {__b, __a};
else
return {__a, __b};
}
template<input_range _Range, typename _Proj = identity,
indirect_strict_weak_order<projected<iterator_t<_Range>, _Proj>>
_Comp = ranges::less>
requires indirectly_copyable_storable<iterator_t<_Range>, range_value_t<_Range>*>
[[nodiscard]] constexpr minmax_result<range_value_t<_Range>>
operator()(_Range&& __r, _Comp __comp = {}, _Proj __proj = {}) const
{
auto __first = ranges::begin(__r);
auto __last = ranges::end(__r);
__glibcxx_assert(__first != __last);
auto __comp_proj = __detail::__make_comp_proj(__comp, __proj);
minmax_result<range_value_t<_Range>> __result = {*__first, __result.min};
if (++__first == __last)
return __result;
else
{
// At this point __result.min == __result.max, so a single
// comparison with the next element suffices.
auto&& __val = *__first;
if (__comp_proj(__val, __result.min))
__result.min = std::forward<decltype(__val)>(__val);
else
__result.max = std::forward<decltype(__val)>(__val);
}
while (++__first != __last)
{
// Now process two elements at a time so that we perform at most
// 1 + 3*(N-2)/2 comparisons in total (each of the (N-2)/2
// iterations of this loop performs three comparisons).
range_value_t<_Range> __val1 = *__first;
if (++__first == __last)
{
// N is odd; in this final iteration, we perform at most two
// comparisons, for a total of 1 + 3*(N-3)/2 + 2 comparisons,
// which is not more than 3*N/2, as required.
if (__comp_proj(__val1, __result.min))
__result.min = std::move(__val1);
else if (!__comp_proj(__val1, __result.max))
__result.max = std::move(__val1);
break;
}
auto&& __val2 = *__first;
if (!__comp_proj(__val2, __val1))
{
if (__comp_proj(__val1, __result.min))
__result.min = std::move(__val1);
if (!__comp_proj(__val2, __result.max))
__result.max = std::forward<decltype(__val2)>(__val2);
}
else
{
if (__comp_proj(__val2, __result.min))
__result.min = std::forward<decltype(__val2)>(__val2);
if (!__comp_proj(__val1, __result.max))
__result.max = std::move(__val1);
}
}
return __result;
}
template<copyable _Tp, typename _Proj = identity,
indirect_strict_weak_order<projected<const _Tp*, _Proj>>
_Comp = ranges::less>
[[nodiscard]] constexpr minmax_result<_Tp>
operator()(initializer_list<_Tp> __r,
_Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::subrange(__r),
std::move(__comp), std::move(__proj));
}
};
inline constexpr __minmax_fn minmax{};
struct __min_element_fn
{
template<forward_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
indirect_strict_weak_order<projected<_Iter, _Proj>>
_Comp = ranges::less>
[[nodiscard]] constexpr _Iter
operator()(_Iter __first, _Sent __last,
_Comp __comp = {}, _Proj __proj = {}) const
{
if (__first == __last)
return __first;
auto __i = __first;
while (++__i != __last)
{
if (std::__invoke(__comp,
std::__invoke(__proj, *__i),
std::__invoke(__proj, *__first)))
__first = __i;
}
return __first;
}
template<forward_range _Range, typename _Proj = identity,
indirect_strict_weak_order<projected<iterator_t<_Range>, _Proj>>
_Comp = ranges::less>
[[nodiscard]] constexpr borrowed_iterator_t<_Range>
operator()(_Range&& __r, _Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__comp), std::move(__proj));
}
};
inline constexpr __min_element_fn min_element{};
struct __max_element_fn
{
template<forward_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
indirect_strict_weak_order<projected<_Iter, _Proj>>
_Comp = ranges::less>
[[nodiscard]] constexpr _Iter
operator()(_Iter __first, _Sent __last,
_Comp __comp = {}, _Proj __proj = {}) const
{
if (__first == __last)
return __first;
auto __i = __first;
while (++__i != __last)
{
if (std::__invoke(__comp,
std::__invoke(__proj, *__first),
std::__invoke(__proj, *__i)))
__first = __i;
}
return __first;
}
template<forward_range _Range, typename _Proj = identity,
indirect_strict_weak_order<projected<iterator_t<_Range>, _Proj>>
_Comp = ranges::less>
[[nodiscard]] constexpr borrowed_iterator_t<_Range>
operator()(_Range&& __r, _Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__comp), std::move(__proj));
}
};
inline constexpr __max_element_fn max_element{};
template<typename _Iter>
using minmax_element_result = min_max_result<_Iter>;
struct __minmax_element_fn
{
template<forward_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
indirect_strict_weak_order<projected<_Iter, _Proj>>
_Comp = ranges::less>
[[nodiscard]] constexpr minmax_element_result<_Iter>
operator()(_Iter __first, _Sent __last,
_Comp __comp = {}, _Proj __proj = {}) const
{
auto __comp_proj = __detail::__make_comp_proj(__comp, __proj);
minmax_element_result<_Iter> __result = {__first, __first};
if (__first == __last || ++__first == __last)
return __result;
else
{
// At this point __result.min == __result.max, so a single
// comparison with the next element suffices.
if (__comp_proj(*__first, *__result.min))
__result.min = __first;
else
__result.max = __first;
}
while (++__first != __last)
{
// Now process two elements at a time so that we perform at most
// 1 + 3*(N-2)/2 comparisons in total (each of the (N-2)/2
// iterations of this loop performs three comparisons).
auto __prev = __first;
if (++__first == __last)
{
// N is odd; in this final iteration, we perform at most two
// comparisons, for a total of 1 + 3*(N-3)/2 + 2 comparisons,
// which is not more than 3*N/2, as required.
if (__comp_proj(*__prev, *__result.min))
__result.min = __prev;
else if (!__comp_proj(*__prev, *__result.max))
__result.max = __prev;
break;
}
if (!__comp_proj(*__first, *__prev))
{
if (__comp_proj(*__prev, *__result.min))
__result.min = __prev;
if (!__comp_proj(*__first, *__result.max))
__result.max = __first;
}
else
{
if (__comp_proj(*__first, *__result.min))
__result.min = __first;
if (!__comp_proj(*__prev, *__result.max))
__result.max = __prev;
}
}
return __result;
}
template<forward_range _Range, typename _Proj = identity,
indirect_strict_weak_order<projected<iterator_t<_Range>, _Proj>>
_Comp = ranges::less>
[[nodiscard]] constexpr minmax_element_result<borrowed_iterator_t<_Range>>
operator()(_Range&& __r, _Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__comp), std::move(__proj));
}
};
inline constexpr __minmax_element_fn minmax_element{};
struct __lexicographical_compare_fn
{
template<input_iterator _Iter1, sentinel_for<_Iter1> _Sent1,
input_iterator _Iter2, sentinel_for<_Iter2> _Sent2,
typename _Proj1 = identity, typename _Proj2 = identity,
indirect_strict_weak_order<projected<_Iter1, _Proj1>,
projected<_Iter2, _Proj2>>
_Comp = ranges::less>
[[nodiscard]] constexpr bool
operator()(_Iter1 __first1, _Sent1 __last1,
_Iter2 __first2, _Sent2 __last2,
_Comp __comp = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
if constexpr (__detail::__is_normal_iterator<_Iter1>
&& same_as<_Iter1, _Sent1>)
return (*this)(__first1.base(), __last1.base(),
std::move(__first2), std::move(__last2),
std::move(__comp),
std::move(__proj1), std::move(__proj2));
else if constexpr (__detail::__is_normal_iterator<_Iter2>
&& same_as<_Iter2, _Sent2>)
return (*this)(std::move(__first1), std::move(__last1),
__first2.base(), __last2.base(),
std::move(__comp),
std::move(__proj1), std::move(__proj2));
else
{
constexpr bool __sized_iters
= (sized_sentinel_for<_Sent1, _Iter1>
&& sized_sentinel_for<_Sent2, _Iter2>);
if constexpr (__sized_iters)
{
using _ValueType1 = iter_value_t<_Iter1>;
using _ValueType2 = iter_value_t<_Iter2>;
// This condition is consistent with the one in
// __lexicographical_compare_aux in <bits/stl_algobase.h>.
constexpr bool __use_memcmp
= (__is_memcmp_ordered_with<_ValueType1, _ValueType2>::__value
&& __ptr_to_nonvolatile<_Iter1>
&& __ptr_to_nonvolatile<_Iter2>
&& (is_same_v<_Comp, ranges::less>
|| is_same_v<_Comp, ranges::greater>)
&& is_same_v<_Proj1, identity>
&& is_same_v<_Proj2, identity>);
if constexpr (__use_memcmp)
{
const auto __d1 = __last1 - __first1;
const auto __d2 = __last2 - __first2;
if (const auto __len = std::min(__d1, __d2))
{
const auto __c
= std::__memcmp(__first1, __first2, __len);
if constexpr (is_same_v<_Comp, ranges::less>)
{
if (__c < 0)
return true;
if (__c > 0)
return false;
}
else if constexpr (is_same_v<_Comp, ranges::greater>)
{
if (__c > 0)
return true;
if (__c < 0)
return false;
}
}
return __d1 < __d2;
}
}
for (; __first1 != __last1 && __first2 != __last2;
++__first1, (void) ++__first2)
{
if (std::__invoke(__comp,
std::__invoke(__proj1, *__first1),
std::__invoke(__proj2, *__first2)))
return true;
if (std::__invoke(__comp,
std::__invoke(__proj2, *__first2),
std::__invoke(__proj1, *__first1)))
return false;
}
return __first1 == __last1 && __first2 != __last2;
}
}
template<input_range _Range1, input_range _Range2,
typename _Proj1 = identity, typename _Proj2 = identity,
indirect_strict_weak_order<projected<iterator_t<_Range1>, _Proj1>,
projected<iterator_t<_Range2>, _Proj2>>
_Comp = ranges::less>
[[nodiscard]] constexpr bool
operator()(_Range1&& __r1, _Range2&& __r2, _Comp __comp = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
return (*this)(ranges::begin(__r1), ranges::end(__r1),
ranges::begin(__r2), ranges::end(__r2),
std::move(__comp),
std::move(__proj1), std::move(__proj2));
}
private:
template<typename _Iter, typename _Ref = iter_reference_t<_Iter>>
static constexpr bool __ptr_to_nonvolatile
= is_pointer_v<_Iter> && !is_volatile_v<remove_reference_t<_Ref>>;
};
inline constexpr __lexicographical_compare_fn lexicographical_compare;
template<typename _Iter>
struct in_found_result
{
[[no_unique_address]] _Iter in;
bool found;
template<typename _Iter2>
requires convertible_to<const _Iter&, _Iter2>
constexpr
operator in_found_result<_Iter2>() const &
{ return {in, found}; }
template<typename _Iter2>
requires convertible_to<_Iter, _Iter2>
constexpr
operator in_found_result<_Iter2>() &&
{ return {std::move(in), found}; }
};
template<typename _Iter>
using next_permutation_result = in_found_result<_Iter>;
struct __next_permutation_fn
{
template<bidirectional_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Comp = ranges::less, typename _Proj = identity>
requires sortable<_Iter, _Comp, _Proj>
constexpr next_permutation_result<_Iter>
operator()(_Iter __first, _Sent __last,
_Comp __comp = {}, _Proj __proj = {}) const
{
if (__first == __last)
return {std::move(__first), false};
auto __i = __first;
++__i;
if (__i == __last)
return {std::move(__i), false};
auto __lasti = ranges::next(__first, __last);
__i = __lasti;
--__i;
for (;;)
{
auto __ii = __i;
--__i;
if (std::__invoke(__comp,
std::__invoke(__proj, *__i),
std::__invoke(__proj, *__ii)))
{
auto __j = __lasti;
while (!(bool)std::__invoke(__comp,
std::__invoke(__proj, *__i),
std::__invoke(__proj, *--__j)))
;
ranges::iter_swap(__i, __j);
ranges::reverse(__ii, __last);
return {std::move(__lasti), true};
}
if (__i == __first)
{
ranges::reverse(__first, __last);
return {std::move(__lasti), false};
}
}
}
template<bidirectional_range _Range, typename _Comp = ranges::less,
typename _Proj = identity>
requires sortable<iterator_t<_Range>, _Comp, _Proj>
constexpr next_permutation_result<borrowed_iterator_t<_Range>>
operator()(_Range&& __r, _Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__comp), std::move(__proj));
}
};
inline constexpr __next_permutation_fn next_permutation{};
template<typename _Iter>
using prev_permutation_result = in_found_result<_Iter>;
struct __prev_permutation_fn
{
template<bidirectional_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Comp = ranges::less, typename _Proj = identity>
requires sortable<_Iter, _Comp, _Proj>
constexpr prev_permutation_result<_Iter>
operator()(_Iter __first, _Sent __last,
_Comp __comp = {}, _Proj __proj = {}) const
{
if (__first == __last)
return {std::move(__first), false};
auto __i = __first;
++__i;
if (__i == __last)
return {std::move(__i), false};
auto __lasti = ranges::next(__first, __last);
__i = __lasti;
--__i;
for (;;)
{
auto __ii = __i;
--__i;
if (std::__invoke(__comp,
std::__invoke(__proj, *__ii),
std::__invoke(__proj, *__i)))
{
auto __j = __lasti;
while (!(bool)std::__invoke(__comp,
std::__invoke(__proj, *--__j),
std::__invoke(__proj, *__i)))
;
ranges::iter_swap(__i, __j);
ranges::reverse(__ii, __last);
return {std::move(__lasti), true};
}
if (__i == __first)
{
ranges::reverse(__first, __last);
return {std::move(__lasti), false};
}
}
}
template<bidirectional_range _Range, typename _Comp = ranges::less,
typename _Proj = identity>
requires sortable<iterator_t<_Range>, _Comp, _Proj>
constexpr prev_permutation_result<borrowed_iterator_t<_Range>>
operator()(_Range&& __r, _Comp __comp = {}, _Proj __proj = {}) const
{
return (*this)(ranges::begin(__r), ranges::end(__r),
std::move(__comp), std::move(__proj));
}
};
inline constexpr __prev_permutation_fn prev_permutation{};
#if __glibcxx_ranges_contains >= 202207L // C++ >= 23
struct __contains_fn
{
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
typename _Tp _GLIBCXX26_RANGE_ALGO_DEF_VAL_T(_Iter, _Proj)>
requires indirect_binary_predicate<ranges::equal_to,
projected<_Iter, _Proj>, const _Tp*>
constexpr bool
operator()(_Iter __first, _Sent __last, const _Tp& __value, _Proj __proj = {}) const
{ return ranges::find(std::move(__first), __last, __value, std::move(__proj)) != __last; }
template<input_range _Range,
typename _Proj = identity,
typename _Tp
_GLIBCXX26_RANGE_ALGO_DEF_VAL_T(iterator_t<_Range>, _Proj)>
requires indirect_binary_predicate<ranges::equal_to,
projected<iterator_t<_Range>, _Proj>, const _Tp*>
constexpr bool
operator()(_Range&& __r, const _Tp& __value, _Proj __proj = {}) const
{ return (*this)(ranges::begin(__r), ranges::end(__r), __value, std::move(__proj)); }
};
inline constexpr __contains_fn contains{};
struct __contains_subrange_fn
{
template<forward_iterator _Iter1, sentinel_for<_Iter1> _Sent1,
forward_iterator _Iter2, sentinel_for<_Iter2> _Sent2,
typename _Pred = ranges::equal_to,
typename _Proj1 = identity, typename _Proj2 = identity>
requires indirectly_comparable<_Iter1, _Iter2, _Pred, _Proj1, _Proj2>
constexpr bool
operator()(_Iter1 __first1, _Sent1 __last1, _Iter2 __first2, _Sent2 __last2,
_Pred __pred = {}, _Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
return __first2 == __last2
|| !ranges::search(__first1, __last1, __first2, __last2,
std::move(__pred), std::move(__proj1), std::move(__proj2)).empty();
}
template<forward_range _Range1, forward_range _Range2,
typename _Pred = ranges::equal_to,
typename _Proj1 = identity, typename _Proj2 = identity>
requires indirectly_comparable<iterator_t<_Range1>, iterator_t<_Range2>,
_Pred, _Proj1, _Proj2>
constexpr bool
operator()(_Range1&& __r1, _Range2&& __r2, _Pred __pred = {},
_Proj1 __proj1 = {}, _Proj2 __proj2 = {}) const
{
return (*this)(ranges::begin(__r1), ranges::end(__r1),
ranges::begin(__r2), ranges::end(__r2),
std::move(__pred), std::move(__proj1), std::move(__proj2));
}
};
inline constexpr __contains_subrange_fn contains_subrange{};
#endif // __glibcxx_ranges_contains
#if __glibcxx_ranges_find_last >= 202207L // C++ >= 23
struct __find_last_fn
{
template<forward_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Proj = identity,
typename _Tp _GLIBCXX26_RANGE_ALGO_DEF_VAL_T(_Iter, _Proj)>
requires indirect_binary_predicate<ranges::equal_to, projected<_Iter, _Proj>, const _Tp*>
[[nodiscard]] constexpr subrange<_Iter>
operator()(_Iter __first, _Sent __last, const _Tp& __value, _Proj __proj = {}) const
{
if constexpr (same_as<_Iter, _Sent> && bidirectional_iterator<_Iter>)
{
_Iter __found = ranges::find(reverse_iterator<_Iter>{__last},
reverse_iterator<_Iter>{__first},
__value, std::move(__proj)).base();
if (__found == __first)
return {__last, __last};
else
return {ranges::prev(__found), __last};
}
else
{
_Iter __found = ranges::find(__first, __last, __value, __proj);
if (__found == __last)
return {__found, __found};
__first = __found;
for (;;)
{
__first = ranges::find(ranges::next(__first), __last, __value, __proj);
if (__first == __last)
return {__found, __first};
__found = __first;
}
}
}
template<forward_range _Range, typename _Proj = identity,
typename _Tp
_GLIBCXX26_RANGE_ALGO_DEF_VAL_T(iterator_t<_Range>, _Proj)>
requires indirect_binary_predicate<ranges::equal_to, projected<iterator_t<_Range>, _Proj>, const _Tp*>
[[nodiscard]] constexpr borrowed_subrange_t<_Range>
operator()(_Range&& __r, const _Tp& __value, _Proj __proj = {}) const
{ return (*this)(ranges::begin(__r), ranges::end(__r), __value, std::move(__proj)); }
};
inline constexpr __find_last_fn find_last{};
struct __find_last_if_fn
{
template<forward_iterator _Iter, sentinel_for<_Iter> _Sent, typename _Proj = identity,
indirect_unary_predicate<projected<_Iter, _Proj>> _Pred>
[[nodiscard]] constexpr subrange<_Iter>
operator()(_Iter __first, _Sent __last, _Pred __pred, _Proj __proj = {}) const
{
if constexpr (same_as<_Iter, _Sent> && bidirectional_iterator<_Iter>)
{
_Iter __found = ranges::find_if(reverse_iterator<_Iter>{__last},
reverse_iterator<_Iter>{__first},
std::move(__pred), std::move(__proj)).base();
if (__found == __first)
return {__last, __last};
else
return {ranges::prev(__found), __last};
}
else
{
_Iter __found = ranges::find_if(__first, __last, __pred, __proj);
if (__found == __last)
return {__found, __found};
__first = __found;
for (;;)
{
__first = ranges::find_if(ranges::next(__first), __last, __pred, __proj);
if (__first == __last)
return {__found, __first};
__found = __first;
}
}
}
template<forward_range _Range, typename _Proj = identity,
indirect_unary_predicate<projected<iterator_t<_Range>, _Proj>> _Pred>
[[nodiscard]] constexpr borrowed_subrange_t<_Range>
operator()(_Range&& __r, _Pred __pred, _Proj __proj = {}) const
{ return (*this)(ranges::begin(__r), ranges::end(__r), std::move(__pred), std::move(__proj)); }
};
inline constexpr __find_last_if_fn find_last_if{};
struct __find_last_if_not_fn
{
template<forward_iterator _Iter, sentinel_for<_Iter> _Sent, typename _Proj = identity,
indirect_unary_predicate<projected<_Iter, _Proj>> _Pred>
[[nodiscard]] constexpr subrange<_Iter>
operator()(_Iter __first, _Sent __last, _Pred __pred, _Proj __proj = {}) const
{
if constexpr (same_as<_Iter, _Sent> && bidirectional_iterator<_Iter>)
{
_Iter __found = ranges::find_if_not(reverse_iterator<_Iter>{__last},
reverse_iterator<_Iter>{__first},
std::move(__pred), std::move(__proj)).base();
if (__found == __first)
return {__last, __last};
else
return {ranges::prev(__found), __last};
}
else
{
_Iter __found = ranges::find_if_not(__first, __last, __pred, __proj);
if (__found == __last)
return {__found, __found};
__first = __found;
for (;;)
{
__first = ranges::find_if_not(ranges::next(__first), __last, __pred, __proj);
if (__first == __last)
return {__found, __first};
__found = __first;
}
}
}
template<forward_range _Range, typename _Proj = identity,
indirect_unary_predicate<projected<iterator_t<_Range>, _Proj>> _Pred>
[[nodiscard]] constexpr borrowed_subrange_t<_Range>
operator()(_Range&& __r, _Pred __pred, _Proj __proj = {}) const
{ return (*this)(ranges::begin(__r), ranges::end(__r), std::move(__pred), std::move(__proj)); }
};
inline constexpr __find_last_if_not_fn find_last_if_not{};
#endif // __glibcxx_ranges_find_last
#if __glibcxx_ranges_fold >= 202207L // C++ >= 23
template<typename _Iter, typename _Tp>
struct in_value_result
{
[[no_unique_address]] _Iter in;
[[no_unique_address]] _Tp value;
template<typename _Iter2, typename _Tp2>
requires convertible_to<const _Iter&, _Iter2>
&& convertible_to<const _Tp&, _Tp2>
constexpr
operator in_value_result<_Iter2, _Tp2>() const &
{ return {in, value}; }
template<typename _Iter2, typename _Tp2>
requires convertible_to<_Iter, _Iter2>
&& convertible_to<_Tp, _Tp2>
constexpr
operator in_value_result<_Iter2, _Tp2>() &&
{ return {std::move(in), std::move(value)}; }
};
namespace __detail
{
template<typename _Fp>
class __flipped
{
_Fp _M_f;
public:
template<typename _Tp, typename _Up>
requires invocable<_Fp&, _Up, _Tp>
invoke_result_t<_Fp&, _Up, _Tp>
operator()(_Tp&&, _Up&&); // not defined
};
template<typename _Fp, typename _Tp, typename _Iter, typename _Up>
concept __indirectly_binary_left_foldable_impl = movable<_Tp> && movable<_Up>
&& convertible_to<_Tp, _Up>
&& invocable<_Fp&, _Up, iter_reference_t<_Iter>>
&& assignable_from<_Up&, invoke_result_t<_Fp&, _Up, iter_reference_t<_Iter>>>;
template<typename _Fp, typename _Tp, typename _Iter>
concept __indirectly_binary_left_foldable = copy_constructible<_Fp>
&& indirectly_readable<_Iter>
&& invocable<_Fp&, _Tp, iter_reference_t<_Iter>>
&& convertible_to<invoke_result_t<_Fp&, _Tp, iter_reference_t<_Iter>>,
decay_t<invoke_result_t<_Fp&, _Tp, iter_reference_t<_Iter>>>>
&& __indirectly_binary_left_foldable_impl
<_Fp, _Tp, _Iter, decay_t<invoke_result_t<_Fp&, _Tp, iter_reference_t<_Iter>>>>;
template <typename _Fp, typename _Tp, typename _Iter>
concept __indirectly_binary_right_foldable
= __indirectly_binary_left_foldable<__flipped<_Fp>, _Tp, _Iter>;
} // namespace __detail
template<typename _Iter, typename _Tp>
using fold_left_with_iter_result = in_value_result<_Iter, _Tp>;
struct __fold_left_with_iter_fn
{
template<typename _Ret_iter,
typename _Iter, typename _Sent, typename _Tp, typename _Fp>
static constexpr auto
_S_impl(_Iter __first, _Sent __last, _Tp __init, _Fp __f)
{
using _Up = decay_t<invoke_result_t<_Fp&, _Tp, iter_reference_t<_Iter>>>;
using _Ret = fold_left_with_iter_result<_Ret_iter, _Up>;
if (__first == __last)
return _Ret{std::move(__first), _Up(std::move(__init))};
_Up __accum = std::__invoke(__f, std::move(__init), *__first);
for (++__first; __first != __last; ++__first)
__accum = std::__invoke(__f, std::move(__accum), *__first);
return _Ret{std::move(__first), std::move(__accum)};
}
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Tp _GLIBCXX26_DEF_VAL_T(iter_value_t<_Iter>),
__detail::__indirectly_binary_left_foldable<_Tp, _Iter> _Fp>
constexpr auto
operator()(_Iter __first, _Sent __last, _Tp __init, _Fp __f) const
{
using _Ret_iter = _Iter;
return _S_impl<_Ret_iter>(std::move(__first), __last,
std::move(__init), std::move(__f));
}
template<input_range _Range,
typename _Tp _GLIBCXX26_DEF_VAL_T(range_value_t<_Range>),
__detail::__indirectly_binary_left_foldable<_Tp, iterator_t<_Range>> _Fp>
constexpr auto
operator()(_Range&& __r, _Tp __init, _Fp __f) const
{
using _Ret_iter = borrowed_iterator_t<_Range>;
return _S_impl<_Ret_iter>(ranges::begin(__r), ranges::end(__r),
std::move(__init), std::move(__f));
}
};
inline constexpr __fold_left_with_iter_fn fold_left_with_iter{};
struct __fold_left_fn
{
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Tp _GLIBCXX26_DEF_VAL_T(iter_value_t<_Iter>),
__detail::__indirectly_binary_left_foldable<_Tp, _Iter> _Fp>
constexpr auto
operator()(_Iter __first, _Sent __last, _Tp __init, _Fp __f) const
{
return ranges::fold_left_with_iter(std::move(__first), __last,
std::move(__init), std::move(__f)).value;
}
template<input_range _Range,
typename _Tp _GLIBCXX26_DEF_VAL_T(range_value_t<_Range>),
__detail::__indirectly_binary_left_foldable<_Tp, iterator_t<_Range>> _Fp>
constexpr auto
operator()(_Range&& __r, _Tp __init, _Fp __f) const
{ return (*this)(ranges::begin(__r), ranges::end(__r), std::move(__init), std::move(__f)); }
};
inline constexpr __fold_left_fn fold_left{};
template<typename _Iter, typename _Tp>
using fold_left_first_with_iter_result = in_value_result<_Iter, _Tp>;
struct __fold_left_first_with_iter_fn
{
template<typename _Ret_iter, typename _Iter, typename _Sent, typename _Fp>
static constexpr auto
_S_impl(_Iter __first, _Sent __last, _Fp __f)
{
using _Up = decltype(ranges::fold_left(std::move(__first), __last,
iter_value_t<_Iter>(*__first), __f));
using _Ret = fold_left_first_with_iter_result<_Ret_iter, optional<_Up>>;
if (__first == __last)
return _Ret{std::move(__first), optional<_Up>()};
optional<_Up> __init(in_place, *__first);
for (++__first; __first != __last; ++__first)
*__init = std::__invoke(__f, std::move(*__init), *__first);
return _Ret{std::move(__first), std::move(__init)};
}
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
__detail::__indirectly_binary_left_foldable<iter_value_t<_Iter>, _Iter> _Fp>
requires constructible_from<iter_value_t<_Iter>, iter_reference_t<_Iter>>
constexpr auto
operator()(_Iter __first, _Sent __last, _Fp __f) const
{
using _Ret_iter = _Iter;
return _S_impl<_Ret_iter>(std::move(__first), __last, std::move(__f));
}
template<input_range _Range,
__detail::__indirectly_binary_left_foldable<range_value_t<_Range>, iterator_t<_Range>> _Fp>
requires constructible_from<range_value_t<_Range>, range_reference_t<_Range>>
constexpr auto
operator()(_Range&& __r, _Fp __f) const
{
using _Ret_iter = borrowed_iterator_t<_Range>;
return _S_impl<_Ret_iter>(ranges::begin(__r), ranges::end(__r), std::move(__f));
}
};
inline constexpr __fold_left_first_with_iter_fn fold_left_first_with_iter{};
struct __fold_left_first_fn
{
template<input_iterator _Iter, sentinel_for<_Iter> _Sent,
__detail::__indirectly_binary_left_foldable<iter_value_t<_Iter>, _Iter> _Fp>
requires constructible_from<iter_value_t<_Iter>, iter_reference_t<_Iter>>
constexpr auto
operator()(_Iter __first, _Sent __last, _Fp __f) const
{
return ranges::fold_left_first_with_iter(std::move(__first), __last,
std::move(__f)).value;
}
template<input_range _Range,
__detail::__indirectly_binary_left_foldable<range_value_t<_Range>, iterator_t<_Range>> _Fp>
requires constructible_from<range_value_t<_Range>, range_reference_t<_Range>>
constexpr auto
operator()(_Range&& __r, _Fp __f) const
{ return (*this)(ranges::begin(__r), ranges::end(__r), std::move(__f)); }
};
inline constexpr __fold_left_first_fn fold_left_first{};
struct __fold_right_fn
{
template<bidirectional_iterator _Iter, sentinel_for<_Iter> _Sent,
typename _Tp _GLIBCXX26_DEF_VAL_T(iter_value_t<_Iter>),
__detail::__indirectly_binary_right_foldable<_Tp, _Iter> _Fp>
constexpr auto
operator()(_Iter __first, _Sent __last, _Tp __init, _Fp __f) const
{
using _Up = decay_t<invoke_result_t<_Fp&, iter_reference_t<_Iter>, _Tp>>;
if (__first == __last)
return _Up(std::move(__init));
_Iter __tail = ranges::next(__first, __last);
_Up __accum = std::__invoke(__f, *--__tail, std::move(__init));
while (__first != __tail)
__accum = std::__invoke(__f, *--__tail, std::move(__accum));
return __accum;
}
template<bidirectional_range _Range,
typename _Tp _GLIBCXX26_DEF_VAL_T(range_value_t<_Range>),
__detail::__indirectly_binary_right_foldable<_Tp, iterator_t<_Range>> _Fp>
constexpr auto
operator()(_Range&& __r, _Tp __init, _Fp __f) const
{ return (*this)(ranges::begin(__r), ranges::end(__r), std::move(__init), std::move(__f)); }
};
inline constexpr __fold_right_fn fold_right{};
struct __fold_right_last_fn
{
template<bidirectional_iterator _Iter, sentinel_for<_Iter> _Sent,
__detail::__indirectly_binary_right_foldable<iter_value_t<_Iter>, _Iter> _Fp>
requires constructible_from<iter_value_t<_Iter>, iter_reference_t<_Iter>>
constexpr auto
operator()(_Iter __first, _Sent __last, _Fp __f) const
{
using _Up = decltype(ranges::fold_right(__first, __last,
iter_value_t<_Iter>(*__first), __f));
if (__first == __last)
return optional<_Up>();
_Iter __tail = ranges::prev(ranges::next(__first, std::move(__last)));
return optional<_Up>(in_place,
ranges::fold_right(std::move(__first), __tail,
iter_value_t<_Iter>(*__tail),
std::move(__f)));
}
template<bidirectional_range _Range,
__detail::__indirectly_binary_right_foldable<range_value_t<_Range>, iterator_t<_Range>> _Fp>
requires constructible_from<range_value_t<_Range>, range_reference_t<_Range>>
constexpr auto
operator()(_Range&& __r, _Fp __f) const
{ return (*this)(ranges::begin(__r), ranges::end(__r), std::move(__f)); }
};
inline constexpr __fold_right_last_fn fold_right_last{};
#endif // __glibcxx_ranges_fold
} // namespace ranges
#if __glibcxx_shift >= 201806L // C++ >= 20
template<typename _ForwardIterator>
constexpr _ForwardIterator
shift_left(_ForwardIterator __first, _ForwardIterator __last,
typename iterator_traits<_ForwardIterator>::difference_type __n)
{
__glibcxx_assert(__n >= 0);
if (__n == 0)
return __last;
auto __mid = ranges::next(__first, __n, __last);
if (__mid == __last)
return __first;
return std::move(std::move(__mid), std::move(__last), std::move(__first));
}
template<typename _ForwardIterator>
constexpr _ForwardIterator
shift_right(_ForwardIterator __first, _ForwardIterator __last,
typename iterator_traits<_ForwardIterator>::difference_type __n)
{
__glibcxx_assert(__n >= 0);
if (__n == 0)
return __first;
using _Cat
= typename iterator_traits<_ForwardIterator>::iterator_category;
if constexpr (derived_from<_Cat, bidirectional_iterator_tag>)
{
auto __mid = ranges::next(__last, -__n, __first);
if (__mid == __first)
return __last;
return std::move_backward(std::move(__first), std::move(__mid),
std::move(__last));
}
else
{
auto __result = ranges::next(__first, __n, __last);
if (__result == __last)
return __last;
auto __dest_head = __first, __dest_tail = __result;
while (__dest_head != __result)
{
if (__dest_tail == __last)
{
// If we get here, then we must have
// 2*n >= distance(__first, __last)
// i.e. we are shifting out at least half of the range. In
// this case we can safely perform the shift with a single
// move.
std::move(std::move(__first), std::move(__dest_head), __result);
return __result;
}
++__dest_head;
++__dest_tail;
}
for (;;)
{
// At the start of each iteration of this outer loop, the range
// [__first, __result) contains those elements that after shifting
// the whole range right by __n, should end up in
// [__dest_head, __dest_tail) in order.
// The below inner loop swaps the elements of [__first, __result)
// and [__dest_head, __dest_tail), while simultaneously shifting
// the latter range by __n.
auto __cursor = __first;
while (__cursor != __result)
{
if (__dest_tail == __last)
{
// At this point the ranges [__first, result) and
// [__dest_head, dest_tail) are disjoint, so we can safely
// move the remaining elements.
__dest_head = std::move(__cursor, __result,
std::move(__dest_head));
std::move(std::move(__first), std::move(__cursor),
std::move(__dest_head));
return __result;
}
std::iter_swap(__cursor, __dest_head);
++__dest_head;
++__dest_tail;
++__cursor;
}
}
}
}
#endif
namespace ranges
{
#if __glibcxx_shift >= 202202L // C++ >= 23
struct __shift_left_fn
{
template<permutable _Iter, sentinel_for<_Iter> _Sent>
constexpr subrange<_Iter>
operator()(_Iter __first, _Sent __last, iter_difference_t<_Iter> __n) const
{
__glibcxx_assert(__n >= 0);
if (__n == 0)
return {__first, ranges::next(__first, __last)};
auto __mid = ranges::next(__first, __n, __last);
if (__mid == __last)
return {__first, __first};
return {__first, ranges::move(__mid, __last, __first).out};
}
template<forward_range _Range>
requires permutable<iterator_t<_Range>>
constexpr borrowed_subrange_t<_Range>
operator()(_Range&& __r, range_difference_t<_Range> __n) const
{ return (*this)(ranges::begin(__r), ranges::end(__r), __n); }
};
inline constexpr __shift_left_fn shift_left{};
struct __shift_right_fn
{
template<permutable _Iter, sentinel_for<_Iter> _Sent>
constexpr subrange<_Iter>
operator()(_Iter __first, _Sent __last, iter_difference_t<_Iter> __n) const
{
__glibcxx_assert(__n >= 0);
if (__n == 0)
return {__first, ranges::next(__first, __last)};
if constexpr (bidirectional_iterator<_Iter> && same_as<_Iter, _Sent>)
{
auto __mid = ranges::next(__last, -__n, __first);
if (__mid == __first)
return {__last, __last};
return {ranges::move_backward(__first, __mid, __last).out, __last};
}
else
{
auto __result = ranges::next(__first, __n, __last);
if (__result == __last)
return {__result, __result};
auto __dest_head = __first, __dest_tail = __result;
while (__dest_head != __result)
{
if (__dest_tail == __last)
{
// If we get here, then we must have
// 2*n >= distance(__first, __last)
// i.e. we are shifting out at least half of the range. In
// this case we can safely perform the shift with a single
// move.
auto __lasti = ranges::move(__first, __dest_head, __result).out;
// __glibcxx_assert(__lasti == __last);
return {__result, __lasti};
}
++__dest_head;
++__dest_tail;
}
for (;;)
{
// At the start of each iteration of this outer loop, the range
// [__first, __result) contains those elements that after shifting
// the whole range right by __n, should end up in
// [__dest_head, __dest_tail) in order.
// The below inner loop swaps the elements of [__first, __result)
// and [__dest_head, __dest_tail), while simultaneously shifting
// the latter range by __n.
auto __cursor = __first;
while (__cursor != __result)
{
if (__dest_tail == __last)
{
// At this point the ranges [__first, result) and
// [__dest_head, dest_tail) are disjoint, so we can safely
// move the remaining elements.
__dest_head = ranges::move(__cursor, __result, __dest_head).out;
auto __lasti = ranges::move(__first, __cursor, __dest_head).out;
// __glibcxx_assert(__lasti == __last);
return {__result, __lasti};
}
ranges::iter_swap(__cursor, __dest_head);
++__dest_head;
++__dest_tail;
++__cursor;
}
}
}
}
template<forward_range _Range>
requires permutable<iterator_t<_Range>>
constexpr borrowed_subrange_t<_Range>
operator()(_Range&& __r, range_difference_t<_Range> __n) const
{ return (*this)(ranges::begin(__r), ranges::end(__r), __n); }
};
inline constexpr __shift_right_fn shift_right{};
#endif // C++23
} // namespace ranges
_GLIBCXX_END_NAMESPACE_VERSION
} // namespace std
#endif // concepts
#endif // C++20
#endif // _RANGES_ALGO_H
Reference in New Issue View Git Blame Copy Permalink