有人知道实现的多指数算法吗?我正在寻找给定矢量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算法(或者自己编写一个).