std::inner_product

template< class InputIt1, class InputIt2, class T >
T inner_product( InputIt1 first1, InputIt1 last1,
                 InputIt2 first2, T init );

template< class InputIt1, class InputIt2, class T,
          class BinaryOp1, class BinaryOp2 >
T inner_product( InputIt1 first1, InputIt1 last1,
                 InputIt2 first2, T init,
                 BinaryOp1 op1, BinaryOp2 op2 );

Вычисляет внутреннее произведение (т.е. сумму произведений) или выполняет упорядоченную операцию map/reduce над диапазоном [first1, last1) и диапазоном из std::distance(first1, last1) элементов, начиная с first2.

Если last2 — это std::distance(first1, last1)-й следующий итератор после first2, и выполняется любое из следующих условий, поведение не определено:

• T не является CopyConstructible.
• T не является CopyAssignable.
• op1 или op2 изменяет любой элемент [first1, last1) или [first2, last2).
• op1 или op2 делает недействительным любой итератор или поддиапазон в [first1, last1] или [first2, last2].

Параметры

Возвращаемое значение

acc после всех модификаций.

Возможная реализация

template<class InputIt1, class InputIt2, class T>
constexpr // since C++20
T inner_product(InputIt1 first1, InputIt1 last1, InputIt2 first2, T init)
{
    while (first1 != last1)
    {
        init = std::move(init) + (*first1) * (*first2); // std::move since C++20
        ++first1;
        ++first2;
    }
    
    return init;
}

template<class InputIt1, class InputIt2, class T,
         class BinaryOp1, class BinaryOp2>
constexpr // since C++20
T inner_product(InputIt1 first1, InputIt1 last1, InputIt2 first2, T init,
                BinaryOp1 op1, BinaryOp2 op2)
{
    while (first1 != last1)
    {
        init = op1(std::move(init), op2(*first1, *first2)); // std::move since C++20
        ++first1;
        ++first2;
    }
    
    return init;
}

Примечания

Параллелизуемая версия этого алгоритма, std::transform_reduce, требует, чтобы op1 и op2 были коммутативны и ассоциативны, но std::inner_product не накладывает такого требования и всегда выполняет операции в заданном порядке.

Пример

#include <functional>
#include <iostream>
#include <numeric>
#include <vector>

int main()
{
    std::vector<int> a{0, 1, 2, 3, 4};
    std::vector<int> b{5, 4, 2, 3, 1};
    
    int r1 = std::inner_product(a.begin(), a.end(), b.begin(), 0);
    std::cout << "Inner product of a and b: " << r1 << '\n';
    
    int r2 = std::inner_product(a.begin(), a.end(), b.begin(), 0,
                                std::plus<>(), std::equal_to<>());
    std::cout << "Number of pairwise matches between a and b: " <<  r2 << '\n';
}

Inner product of a and b: 21
Number of pairwise matches between a and b: 2

Отчёты о дефектах

Следующие отчёты о дефектах, изменяющие поведение, были применены ретроактивно к ранее опубликованным стандартам C++.

Смотрите также