/*============================================================================= Copyright (c) 2011-2016 Bolero MURAKAMI https://github.com/bolero-MURAKAMI/Sprout Distributed under the Boost Software License, Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt) =============================================================================*/ #ifndef SPROUT_ALGORITHM_MAX_ELEMENT_HPP #define SPROUT_ALGORITHM_MAX_ELEMENT_HPP #include #include #include #include #include #include #include #include HDR_FUNCTIONAL_SSCRISK_CEL_OR_SPROUT namespace sprout { namespace detail { template inline SPROUT_CONSTEXPR ForwardIterator iter_max(ForwardIterator const& a, ForwardIterator const& b, Compare comp) { return comp(*a, *b) ? b : a; } template inline SPROUT_CONSTEXPR RandomAccessIterator max_element_impl_ra( RandomAccessIterator const& first, RandomAccessIterator const& last, Compare comp, typename std::iterator_traits::difference_type pivot ) { return pivot == 0 ? first : sprout::detail::iter_max( sprout::detail::max_element_impl_ra( first, sprout::next(first, pivot), comp, pivot / 2 ), sprout::detail::max_element_impl_ra( sprout::next(first, pivot), last, comp, (sprout::distance(first, last) - pivot) / 2 ), comp ) ; } template inline SPROUT_CONSTEXPR typename std::enable_if< sprout::is_constant_distance_iterator::value, RandomAccessIterator >::type max_element( RandomAccessIterator const& first, RandomAccessIterator const& last, Compare comp, std::random_access_iterator_tag* ) { return first == last ? last : sprout::detail::max_element_impl_ra( first, last, comp, sprout::distance(first, last) / 2 ) ; } template inline SPROUT_CONSTEXPR sprout::pair max_element_impl_1( sprout::pair const& current, ForwardIterator const& last, Compare comp, typename std::iterator_traits::difference_type n ) { typedef sprout::pair type; return current.first == last ? current : n == 1 ? type(sprout::next(current.first), sprout::detail::iter_max(current.second, current.first, comp)) : sprout::detail::max_element_impl_1( sprout::detail::max_element_impl_1( current, last, comp, n / 2 ), last, comp, n - n / 2 ) ; } template inline SPROUT_CONSTEXPR sprout::pair max_element_impl( sprout::pair const& current, ForwardIterator const& last, Compare comp, typename std::iterator_traits::difference_type n ) { return current.first == last ? current : sprout::detail::max_element_impl( sprout::detail::max_element_impl_1( current, last, comp, n ), last, comp, n * 2 ) ; } template inline SPROUT_CONSTEXPR ForwardIterator max_element( ForwardIterator const& first, ForwardIterator const& last, Compare comp, std::forward_iterator_tag* ) { typedef sprout::pair type; return first == last ? last : sprout::detail::max_element_impl(type(sprout::next(first), first), last, comp, 1).second ; } } // namespace detail // 25.4.7 Minimum and maximum // // recursion depth: // O(log N) // template inline SPROUT_CONSTEXPR ForwardIterator max_element(ForwardIterator first, ForwardIterator last, Compare comp) { typedef typename std::iterator_traits::iterator_category* category; return sprout::detail::max_element(first, last, comp, category()); } template inline SPROUT_CONSTEXPR ForwardIterator max_element(ForwardIterator first, ForwardIterator last) { return sprout::max_element( first, last, NS_SSCRISK_CEL_OR_SPROUT::less::value_type>() ); } } // namespace sprout #endif // #ifndef SPROUT_ALGORITHM_MAX_ELEMENT_HPP