黎曼zeta函数的快捷c++实现

最新推荐文章于 2024-07-14 18:15:43 发布
最新推荐文章于 2024-07-14 18:15:43 发布
阅读量368 收藏
点赞数
文章标签： 算法
复制或抄袭本人的文章一律举报，我会每天全网查
本文链接：https://blog.csdn.net/PIPE111/article/details/109647896
版权 
const double bn[130]={
   0,0.16666666666666665741,-0.033333333333333332871,0.023809523809523808202,-0.033333333333333332871,0.075757575757575759678,-0.25311355311355310249,1.1666666666666667407,-7.0921568627450977118,54.971177944862155584,-529.12424242424242493,6192.1231884057970092,-86580.253113553117146,1425517.1666666667443,-27298231.067816093564,601580873.90064239502,-15116315767.092157364,429614643061.16668701,-13711655205088.332031,488332318973593.1875,-19296579341940068.,841693047573682560.,-40338071854059454464.,2.1150748638081992622e+021,-1.2086626522296526202e+023,7.500866746076964166e+024,-5.0387781014810688499e+026,3.6528776484818122276e+028,-2.8498769302450882361e+030,2.3865427499683627448e+032,-2.1399949257225334859e+034,2.050097572347809739e+036,-2.0938005911346379301e+038,2.2752696488463514863e+040,-2.6257710286239577207e+042,3.2125082102718031743e+044,-4.1598278166794711978e+046,5.6920695482035283174e+048,-8.2183629419784577665e+050,1.2502904327166994004e+053,-2.0015583233248370052e+055,3.3674982915364375556e+057,-5.9470970503135450205e+059,1.1011910323627976762e+062,-2.1355259545253502079e+064,4.3328896986641193847e+066,-9.1885528241669331811e+068,2.0346896776329073708e+071,-4.7003833958035730158e+073,1.1318043445484249411e+076,-2.8382249570693707354e+078,7.4064248979678852935e+080,-2.0096454802756605262e+083,5.6657170050805942089e+085,-1.6584511154136215904e+088,5.036885995049237839e+090,-1.586146823765818563e+093,5.1756743617545625189e+095,-1.7488921840217115846e+098,6.1160519994952182359e+100,-2.2122776912707833194e+103,8.2722776798770968766e+105,-3.1958925111415708478e+108,1.275008222338779314e+111,-5.2500923086774131347e+113,2.2301817894241626634e+116,-9.7684521930955207422e+118,4.4098361978452949821e+121,-2.0508570886464088957e+124,9.8214433279791276581e+126,-4.8412600798208880569e+129,2.455308880148098232e+132,-1.2806926804084747554e+135,6.8676167104668579435e+137,-3.7846468581969105852e+140,2.1426101250665290651e+143,-1.2456727137183695399e+146,7.4345787551000155753e+148,-4.5535795304641703729e+151,2.861211281685887046e+154,-1.8437723552033869952e+157,1.2181154536221046544e+160,-8.2482187185314121679e+162,5.7225877937832942168e+165,-4.066853052505910459e+168,2.9596092064642052191e+171,-2.2049522565189457143e+174,1.6812597072889599375e+177,-1.3116736213556958091e+180,1.0467894009478038522e+183,-8.5432893578833710248e+185,7.1287821322486546855e+188,-6.0802931455535904967e+191,5.2996776424849921167e+194,-4.719425916874586047e+197,4.2928413791402982607e+200,-3.9876744968232205368e+203,3.7819780419358879766e+206,-3.6614233683681191899e+209,3.6176090272372859917e+212,-3.6470772645191356311e+215,3.7508755436454406323e+218,-3.934586729643903006e+221,
最低0.47元/天 解锁文章
C++的boost
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
黎曼zeta函数的c++实现
c++的boost的博客
11-10 419
#include<cmath> using namespace std; template<class T> T Zeta(T s,double prec){ if(s<0&&(int)s==s&&(int)-s%2==0) return T(0); T sc=static_cast<T>(1 - static_cast<T>(s)); T sum=0,mult=0.5,change; T powers[10000
Zeta(2)
jiaboming4的博客
06-18 1335
An elegant proof of the Basel's Problem. i.e. the value of zeta(2).
riemann zeta 函数
01-04
riemann zeta 函数 riemann zeta 函数 riemann zeta 函数 riemann zeta 函数
黎曼 zeta 函数与黎曼猜想
06-16 1万+
ζ(s)\zeta (s), is a function of a complex variable ss that analytically continues the sum of the infinite series： ζ(s)=∑n=1∞1ns \zeta(s)=\sum_{n=1}^\infty\frac1{n^s} 1. 基本性质 （1）ss 的实部如果大于 1，则级数收敛； （2）
zeta.js:打包为 NPM 模块的 Riemann Zeta 函数的实现
07-03
泽塔 javascript 中的，精度误差为 1.3% 例子 > var zeta = require('riemann-zeta') > zeta(0) -0.5 > zeta(10) 1.00071 方法 var zeta = require('riemann-zeta') zeta(s) 返回s的 riemann zeta 函数。 安装 用做： npm install riemann-zeta 执照 麻省理工学院 学分 近似遵循欧拉-麦克劳林求和公式（Ralston and Rabinowitz, 2001，不考虑伯努利数对计算的影响，并以积分下限具体为 2 的情况为例。那么，近似遵循并运行良好适用于现代浏览器和数学应用程序。
常对象只能调用常函数成员 c++
Zetaa的博客
08-18 9805
关于 常对象只能调用常函数成员，不能调用非const 函数成员问题，我有点想法 就是我想到的一个角度来解释为什么。 你看，常函数成员的特性就是要在函数后面加一个const 修饰关键字 以  int getX() const;为例，，将之作解析， 即, int getX(const Coordinate * this); 这次我们重点看的不是 const 修饰谁的问题， 而是 t
黎曼的zeta函数(0)
qtxzh的博客
09-19 1万+
Gamma函数 考虑如下等式 ∑k=0∞xk=1+x+x2+⋯=11−x\sum^{\infin}_{k=0}x^k = 1 + x +x^2 + \cdots = \frac{1}{1-x}k=0∑∞​xk=1+x+x2+⋯=1−x1​ 而恰有 ∫0∞e−(1−x)tdt=[−11−xe−(1−x)t]0∞=11−x \int_0^{\infin}e^{-(1-x)t}\rm dt = [-\frac{1}{1-x}e^{-(1-x)t}]_0^{\infin} = \frac{1}{1-x}∫0∞​e−
一维黎曼问题的精确解
03-21
主要采用fortran编译了一维黎曼问题的精确解，可用于计算空气动力学计算
论文研究 - 关于零的黎曼Zeta函数，贝塞尔函数和高斯函数的共同性质
05-22
研究了复数域中Riemann Xi函数（到zeta函数），修改后的Bessel函数和高斯（bell）函数的有限泰勒级数逼近中零的行为。 可以看出，在逼近高阶逼近的过程中，有限逼近中的零点如何逼近真正的零点，在高斯（贝尔）函数...
zetaRS 逼近黎曼 zeta(0.5 + i*t) 以获得大 t：使用黎曼西格尔公式快速计算临界带上的黎曼 zeta 函数。-matlab开发
05-31
对于大 t，该例程比 Matlab 的符号工具箱附带的 zeta 函数要快得多。 大 t 的误差非常小，但 zetaRS 对于小 t 不是很精确。 示例：>> tic; zetaRS(10000000), toc 答案 = 11.4580 - 8.6434i 经过的时间是 0.009604 ...
通过局部因子的黎曼zeta矩阵模型
04-19
我们提出了构造黎曼zeta函数的random随机矩阵（UMM）的集合。 从这个意义上说，我们考虑zeta函数的Euler乘积表示中的每个因素，首先要为每个素数p构造一个UMM，从这个意义上说，我们的方法是“ p-iecemeal”。 我们...
黎曼函数的性质及其证明.pdf
06-02
黎曼函数（Riemann function）是一个特殊函数，由德国数学家黎曼发现提出，黎曼函数定义在[0,1]上，其基本定义是：R(x)=1/q，当x=p/q(p,q都属于正整数，p/q为既约真分数)；R(x)=0，当x=0，1和(0,1)内的无理数。 [1] ...
黎曼假设
06-03
这项工作代表了[1] [2]中工作的发展和细节，该工作致力于促进Abe​​ls通过修改zeta函数获得的结果。 研究了zeta函数的性质； 这些属性导致zeta函数的新规律。 选择特殊类型的修改后的zeta函数可以估计Riemann的zeta...
en笔记音标
weixin_33873846的博客
05-02 958
清辅音和浊辅音区别 开音节和闭音节区别 1 2 3 4 5 6 7...
C++ array使用方法详细介绍
热门推荐
终曲
11-05 4万+
文章目录
36. Adam 算法详解
qq_41238579的博客
07-10 945
Adam（Adaptive Moment Estimation）是一种结合动量法和自适应学习率的优化算法，自2014年提出以来，迅速成为深度学习中最流行和常用的优化算法之一。Adam算法的核心思想是利用梯度的一阶动量和二阶动量来动态调整学习率，既保持了动量法的优点，又结合了RMSProp的自适应学习率特性。
LC53最大子数组和、lc5最长回文子串、lc283移动零
最新发布
m0_47536537的博客
07-14 55
请你找出一个具有最大和的连续子数组（子数组最少包含一个元素），返回其最大和。连续子数组 [4,-1,2,1] 的和最大，为 6。移动到数组的末尾，同时保持非零元素的相对顺序。，必须在不复制数组的情况下原地对数组进行操作。"aba" 同样是符合题意的答案。是数组中的一个连续部分。，编写一个函数将所有。LC53最大子数组和。
算法014：找到字符串中所有字母异位词
m0_62319039的博客
07-14 180
算法014：找到字符串中所有字母异位词
“素数的频率紧密相关于一个精心构造的所谓黎曼zeta函数&zeta;(s)的性态。复平面上使黎曼&zeta; 函数取值为零的点被称为黎曼&zeta;函数的零点。s=-2n （n 为正整数）是黎曼&zeta; 函数的零点，这些零点分布有序、 性质简单，被称为黎曼&zeta; 函数的平凡零点 (trivial zero)”证明这个数据观点，给出严谨的推理流程
05-05
首先，我们需要了解黎曼 zeta 函数的定义： $$\zeta(s)=\sum_{n=1}^{\infty}\frac{1}{n^s}$$ 其中，$s$ 是一个复数。当 $s$ 的实部大于 $1$ 时，黎曼 zeta 函数是收敛的。当 $s=1$ 时，黎曼 zeta 函数的值为无穷大。 黎曼 zeta 函数的一个重要性质是，它可以被解析延拓到整个复平面，除了 $s=1$ 这个点存在一个极点。具体来说，黎曼 zeta 函数可以被写成以下形式： $$\zeta(s)=\frac{1}{\Gamma(s)}\int_{0}^{\infty}\frac{x^{s-1}}{e^x-1}\mathrm{d}x$$ 其中，$\Gamma(s)$ 是欧拉伽玛函数。 接下来，我们来证明素数的频率与黎曼 zeta 函数的零点相关。为此，我们需要引入另一个函数 $\psi(x)$，它被定义为： $$\psi(x)=\sum_{n\leq x}\Lambda(n)$$ 其中，$\Lambda(n)$ 表示 n 的 von Mangoldt 函数，即： $$\Lambda(n)=\begin{cases}\ln p, & \text{if }n=p^k\text{ for some prime }p\text{ and integer }k\geq 1\\0, & \text{otherwise}\end{cases}$$ 通过分部积分，可以得到： $$\psi(x)=x-\sum_{\rho}\frac{x^{\rho}}{\rho}-\ln 2\pi-\frac{1}{2}\ln(1-x^{-2})$$ 其中，$\rho$ 是黎曼 zeta 函数的零点。 接下来，我们需要证明的是，当 $x$ 趋近于正无穷时，$\psi(x)$ 与素数的个数 $\pi(x)$ 之间的关系是： $$\psi(x)\sim x\ \text{as}\ x\rightarrow\infty$$ 其中，“$\sim$”表示“渐进等于”。 这个结论可以通过黎曼-底格尔公式得到。黎曼-底格尔公式是一个重要的数学公式，它描述了黎曼 zeta 函数与素数分布之间的关系。具体来说，黎曼-底格尔公式可以写成以下形式： $$\pi(x)=\text{li}(x)+O\left(\frac{x}{\ln x}\right)$$ 其中，$\text{li}(x)$ 是对数积分函数，$O\left(\frac{x}{\ln x}\right)$ 是渐进符号，表示当 $x$ 趋近于正无穷时，剩余的误差可以被一个与 $\frac{x}{\ln x}$ 同阶的函数所控制。 我们可以对黎曼-底格尔公式进行微调，得到： $$\psi(x)=\text{li}(x)-\sum_{p}\text{li}(x^{1/p})-\ln 2-\frac{1}{2}\ln(1-x^{-2})+O\left(\frac{x}{\ln x}\right)$$ 其中，$p$ 是素数。这个公式的证明可以参考数论中的相关文献。 接下来，我们证明当 $x$ 趋近于正无穷时，$\psi(x)$ 与素数的个数 $\pi(x)$ 之间的关系是： $$\psi(x)\sim x\ \text{as}\ x\rightarrow\infty$$ 我们可以先证明当 $x$ 趋近于正无穷时，$\text{li}(x)$ 与 $x$ 的差距是比较小的。具体来说，根据定义，$\text{li}(x)$ 可以写成以下积分的形式： $$\text{li}(x)=\int_{2}^{x}\frac{\mathrm{d}t}{\ln t}$$ 通过分部积分，可以得到： $$\text{li}(x)=\frac{x}{\ln x}-\int_{2}^{x}\frac{\mathrm{d}t}{\ln^2 t}+\frac{2}{\ln 2}$$ 因此，当 $x$ 趋近于正无穷时，$\text{li}(x)$ 与 $x$ 的差距是 $O\left(\frac{x}{\ln x}\right)$ 级别的。 接下来，我们证明当 $x$ 趋近于正无穷时，$\sum_{p}\text{li}(x^{1/p})$ 与 $x$ 的差距也是比较小的。具体来说，我们可以写出： $$\sum_{p}\text{li}(x^{1/p})\leq\sum_{n}\text{li}(x^{1/n})=\sum_{n}\int_{2}^{x^{1/n}}\frac{\mathrm{d}t}{\ln t}=\sum_{n}\frac{x^{1/n}}{n\ln x}=O\left(\frac{x}{\ln x}\right)$$ 其中，第一个等式是因为 $\text{li}(x^{1/n})$ 可以看成是 $\text{li}(y)$，其中 $y$ 是满足 $y^n=x$ 的最小整数；第二个等式是通过换元积分得到的；第三个等式是通过级数展开得到的。 因此，当 $x$ 趋近于正无穷时，$\psi(x)$ 与 $x$ 的差距是 $O\left(\frac{x}{\ln x}\right)$ 级别的。这意味着，当 $x$ 趋近于正无穷时，$\psi(x)$ 与 $x$ 是同阶的，即： $$\psi(x)\sim x\ \text{as}\ x\rightarrow\infty$$ 因此，我们证明了素数的频率与黎曼 zeta 函数的零点分布相关的结论。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值