c – 多次执行实施

有人知道实现的多指数算法吗?我正在寻找给定矢量A,B将使用一些快速算法计算A [i] ^ B [i]的乘积的东西.

谢谢!

最佳答案 以下假设您的数据是浮点数.如果您有多精度整数,请指定您的要求.

干净的数字方式当然是先记录日志.实际上,即使结果是有限的,部分产品也可以轻易地溢出/溢出.

惯用的相应C程序是:

#include <cmath>
#include <functional>
#include <numeric>

double f(double x, double y)
{
    return y * std::log(x);
}

template <typename I, typename J>
double multi_exponentiation(I a0, I an, J b0)
{
    return std::exp(std::inner_product(a0, an, b0, 0., std::plus<double>(), f));
}

// Example program
int main()
{
    std::vector<double> a, b;
    ...
    double e = multi_exponentiation(a.begin(), a.end(), b.begin());
}

使用inner_product而不是自己编写循环具有以下优点:一旦您知道性能是一个问题,您可以使用第三方库提供的parallel_inner_product算法替换inner_product算法(或者自己编写一个).

点赞