c – 多次执行实施

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

惯用的相应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算法(或者自己编写一个).