std::transform_reduce返回值错误，如何修正以获取正确计算结果？

问题背景 #

现有向量定义：

vector<double> X {0.95, 3.75};
vector<double> p {3, 0};

需要计算的目标运算为：

(0.95^3)*(3.75^0) // 预期结果：0.8573749999999999

当前使用的代码返回值为1，不符合预期：

std::transform_reduce(
  std::execution::par,
  X.begin(),
  X.end(),
  p.end(),
  1.0,
  &pow,
  std::multiplies()
);

错误原因与修正方案 #

核心错误点 #

1. 第二个范围迭代器错误：双范围版本的transform_reduce需要传入第二个范围的起始迭代器，而非结束迭代器。传入p.end()会导致无法遍历到p的有效元素，运算逻辑完全偏离预期。
2. 操作参数顺序颠倒：transform_reduce的参数规则是：初始值之后先放归约操作（用来合并转换后的结果），再放转换操作（用来处理两个范围的对应元素）。你把计算幂次的pow和乘积归约的std::multiplies()顺序写反了。

修正后的代码 #

std::transform_reduce(
  std::execution::par,
  X.begin(),
  X.end(),
  p.begin(), // 改为p的起始迭代器
  1.0,
  std::multiplies<double>(), // 归约操作：将所有幂次结果相乘
  [](double x, double exp) { return std::pow(x, exp); } // 转换操作：计算x的exp次幂
);

额外注意事项 #

• 避免直接传递&pow：std::pow存在多个重载，直接传函数指针可能引发编译器匹配歧义，使用lambda包装更安全。
• 确保向量长度一致：X和p的元素数量必须相同，否则会触发越界访问的未定义行为。

内容的提问来源于stack exchange，提问作者dearn44