/*============================================================================= Copyright (c) 2011-2017 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_DETAIL_ALGORITHM_OVERLAP_COUNT_HPP #define SPROUT_DETAIL_ALGORITHM_OVERLAP_COUNT_HPP #include #include #include #include #include #include #include HDR_FUNCTIONAL_SSCRISK_CEL_OR_SPROUT namespace sprout { namespace detail { template inline SPROUT_CONSTEXPR typename std::iterator_traits::difference_type overlap_count_impl_ra( RandomAccessIterator const& first, RandomAccessIterator const& last, BinaryPredicate pred, typename std::iterator_traits::difference_type pivot ) { return pivot == 0 ? (pred(*first, *last) ? 1 : 0) : sprout::detail::overlap_count_impl_ra( first, sprout::next(first, pivot), pred, pivot / 2 ) + sprout::detail::overlap_count_impl_ra( sprout::next(first, pivot), last, pred, (sprout::distance(first, last) - pivot) / 2 ) ; } template inline SPROUT_CONSTEXPR typename std::enable_if< sprout::is_constant_distance_iterator::value, typename std::iterator_traits::difference_type >::type overlap_count( RandomAccessIterator const& first, RandomAccessIterator const& last, BinaryPredicate pred, std::random_access_iterator_tag* ) { return first == last || sprout::distance(first, last) == 1 ? 0 : sprout::detail::overlap_count_impl_ra( first, sprout::next(first, sprout::distance(first, last) - 1), pred, (sprout::distance(first, last) - 1) / 2 ) ; } template inline SPROUT_CONSTEXPR sprout::tuples::tuple::value_type, typename std::iterator_traits::difference_type> overlap_count_impl_1( sprout::tuples::tuple::value_type, typename std::iterator_traits::difference_type> const& current, InputIterator const& last, BinaryPredicate pred, typename std::iterator_traits::difference_type n ) { typedef sprout::tuples::tuple::value_type, typename std::iterator_traits::difference_type> type; return sprout::tuples::get<0>(current) == last ? current : n == 1 ? type( sprout::next(sprout::tuples::get<0>(current)), *sprout::tuples::get<0>(current), sprout::tuples::get<2>(current) + (pred(*sprout::tuples::get<0>(current), sprout::tuples::get<1>(current)) ? 1 : 0) ) : sprout::detail::overlap_count_impl_1( sprout::detail::overlap_count_impl_1( current, last, pred, n / 2 ), last, pred, n - n / 2 ) ; } template inline SPROUT_CONSTEXPR sprout::tuples::tuple::value_type, typename std::iterator_traits::difference_type> overlap_count_impl( sprout::tuples::tuple::value_type, typename std::iterator_traits::difference_type> const& current, InputIterator const& last, BinaryPredicate pred, typename std::iterator_traits::difference_type n ) { return sprout::tuples::get<0>(current) == last ? current : sprout::detail::overlap_count_impl( sprout::detail::overlap_count_impl_1( current, last, pred, n ), last, pred, n * 2 ) ; } template inline SPROUT_CONSTEXPR typename std::iterator_traits::difference_type overlap_count( InputIterator const& first, InputIterator const& last, BinaryPredicate pred, std::input_iterator_tag* ) { typedef sprout::tuples::tuple::value_type, typename std::iterator_traits::difference_type> type; return first == last ? 0 : sprout::tuples::get<2>( sprout::detail::overlap_count_impl(type(sprout::next(first), *first, 0), last, pred, 1) ) ; } // // overlap_count // // recursion depth: // O(log N) // template inline SPROUT_CONSTEXPR typename std::iterator_traits::difference_type overlap_count(InputIterator const& first, InputIterator const& last, BinaryPredicate pred) { typedef typename std::iterator_traits::iterator_category* category; return sprout::detail::overlap_count(first, last, pred, category()); } template inline SPROUT_CONSTEXPR typename std::iterator_traits::difference_type overlap_count(InputIterator const& first, InputIterator const& last) { return sprout::detail::overlap_count( first, last, NS_SSCRISK_CEL_OR_SPROUT::equal_to::value_type>() ); } } // namespace detail } // namespace sprout #endif // #ifndef SPROUT_DETAIL_ALGORITHM_OVERLAP_COUNT_HPP