std::riemann_zeta, std::riemann_zetaf, std::riemann_zetal
来自cppreference.com
< cpp | experimental | special functions
double riemann_zeta( double arg ); double riemann_zeta( float arg ); |
(1) | |
double riemann_zeta( IntegralType arg ); |
(2) | |
与所有特殊函数一样,仅当实现将 __STDCPP_MATH_SPEC_FUNCS__
定义为至少 201003L 的值,且用户在包含任何标准库头文件前定义了 __STDCPP_WANT_MATH_SPEC_FUNCS__
时,riemann_zeta
才保证在 <cmath>
中可用。
参数
n | - | 多项式的度数,无符号整数类型的值 |
arg | - | 浮点或整数类型的值 |
返回值
如果未发生错误,则返回 arg 的黎曼 zeta 函数的值,即 ζ(arg),在整个实轴上定义:
- 对于 arg > 1,Σ∞
n=1n-arg
。 - 对于 0 ≤ arg ≤ 1,
Σ∞1 1 - 21-arg
n=1(-1)n-1
n-arg
。 - 对于 arg < 0,2arg
πarg-1
sin(
)Γ(1 − arg)ζ(1 − arg)。πarg 2
错误处理
根据 math_errhandling 的规定进行错误报告。
- 如果实参为 NaN,则返回 NaN 但不报告定义域错误。
注解
不支持 TR 29124 但支持 TR 19768 的实现,在头文件 tr1/cmath
和命名空间 std::tr1
中提供此函数。
此函数的实现也在 boost.math 中可用。
示例
(以 gcc 6.0 运行)
运行此代码
#define __STDCPP_WANT_MATH_SPEC_FUNCS__ 1 #include <cmath> #include <iostream> int main() { // spot checks for well-known values std::cout << "ζ(-1) = " << std::riemann_zeta(-1) << '\n' << "ζ(0) = " << std::riemann_zeta(0) << '\n' << "ζ(1) = " << std::riemann_zeta(1) << '\n' << "ζ(0.5) = " << std::riemann_zeta(0.5) << '\n' << "ζ(2) = " << std::riemann_zeta(2) << ' ' << "(π²/6 = " << std::pow(std::acos(-1), 2) / 6 << ")\n"; }
输出:
ζ(-1) = -0.0833333 ζ(0) = -0.5 ζ(1) = inf ζ(0.5) = -1.46035 ζ(2) = 1.64493 (π²/6 = 1.64493)
外部链接
Weisstein, Eric W. "Riemann Zeta Function." From MathWorld--A Wolfram Web Resource.