std::riemann_zeta, std::riemann_zetaf, std::riemann_zetal
来自cppreference.com
| 在标头 <cmath> 定义
|
||
| (1) | (C++17 起) (C++23 前) |
|
| (C++23 起) | ||
| |
(2) | (C++17 起) |
| |
(3) | (C++17 起) |
| 在标头 <cmath> 定义
|
||
| |
(A) | (C++17 起) |
A) 为所有整数类型提供额外重载,将它们当做
double。参数
| num | - | 浮点数或整数 |
返回值
如果没有发生错误,那么返回 num 的黎曼 Zeta 函数值,ζ(num),定义于整个实轴:
- 对于 num>1,是 Σ∞
n=1n-num - 对于 0≤num≤1,是
Σ∞1 1-21-num
n=1(-1)n-1
n-num - 对于 num<0,是 2num
πnum-1
sin(
)Γ(1−num)ζ(1−num)πnum 2
错误处理
可能报告 math_errhandling 中指定的错误
- 如果实参是 NaN,那么返回 NaN 且不报告定义域错误
注解
不支持 C++17,但支持 ISO 29124:2010 的实现在定义了 __STDCPP_MATH_SPEC_FUNCS__ 为至少 201003L 的值,且用户在包含任何标准库头文件前定义了 __STDCPP_WANT_MATH_SPEC_FUNCS__ 时也会提供此函数。
不支持 ISO 29124:2010 但支持 TR 19768:2007 (TR1) 的实现也会在标头 tr1/cmath 及命名空间 std::tr1 中提供此函数。
此函数的一种实现参考 boost.math。
额外重载不需要严格以 (A) 的形式提供。它们只需要能够对它们的整数类型实参 num 确保 std::riemann_zeta(num) 和 std::riemann_zeta(static_cast<double>(num)) 的效果相同即可。
示例
运行此代码
#include <cmath>
#include <format>
#include <iostream>
#include <numbers>
int main()
{
constexpr auto π = std::numbers::pi;
// 广为人知的值的点检查
for (const double x : {-1.0, 0.0, 1.0, 0.5, 2.0})
std::cout << std::format("ζ({})\t= {:+.5f}\n", x, std::riemann_zeta(x));
std::cout << std::format("π²/6\t= {:+.5f}\n", π * π / 6);
}
输出:
ζ(-1) = -0.08333
ζ(0) = -0.50000
ζ(1) = +inf
ζ(0.5) = -1.46035
ζ(2) = +1.64493
π²/6 = +1.64493
外部链接
| Weisstein, Eric W. “黎曼 Zeta 函数”来自 MathWorld--A Wolfram Web Resource。 |