%I#12 2023年8月25日17:20:51

%S 1,101431212121249274815563963190945117182662083338618333559，

%电话：11290828893752821828633312591566094562158526435665650141135，

%电话：888640364800590255288276580897412861251176745390297425986625434820160690743301503751949108731388102309925375

%N Yu。V.Matiyasevich的黎曼假设检验。

%C Python程序定义的序列是无限的，当且仅当黎曼假设为真，在这种情况下，程序永远不会停止；否则返回负数并停止。

%H Yuri Matiyasevich，<a href=“https://doi.org/10.1016/j.tcs.2019.07.028“>计算机科学中的黎曼假设，《理论计算机科学》，第807卷，2020年，第257-265页。

%t表[Factorial2[2n+3]/15+PrimePi[n]^2（Factorial2[2n]-Factorial 2[2n-2]总和[（-1）^（k+1）/k，{k，n-1}]楼层[Log[2，LCM@@Range[n]]）/2，{n，20}]（*_Eric W.Weisstein_，2023年8月25日*）

%o（Python）

%o从数学导入gcd

%o定义Riemann测试（停止）：

%o d=m=p=0

%o f0=f1=f3=n=q=r=1

%o b=正确

%o#仅添加条件“n<=stop”

%o#保存测试人员的资源。

%当r>=0且n<=停止时：

%o打印（r，end=“，”）

%o d*=2*n

%o d+=-f1，如果（b:=非b）否则f1

%o n+=1

%o g=gcd（n，q）

%o q=（n*q）//g

%o如果g==1:p+=1

%o m=0

%o q2=q

%当q2>1时：

%o q2//=2

%o m+=d

%o f1=2*f0

%o f0*=2*n

%o f3*=（2*n+3）

%o r=f3-p*p*（m-f0）

%o如果n<停止：

%o打印（“运气不好，伯恩哈德！”）

%o返回r

%o黎曼测试（22）

%Y参考A000720。

%K nonn公司

%O 1,2号机组

%A _彼得·卢什尼，2022年8月8日