%I#21 2022年6月22日09:35:50

%S 6,3,7,0,5,6,1,8,4,0,7,4,6,7,6,6,4,3,5,9,6,8,5,0,4,7,5,2,7，

%T 6,9,4,5,7,9,8,9,6,0,7,7,1,9,9,5,3,6,7,0,9,0,6,1,3,7,10,75,5,8,8，

%U 3,1,6,0,4,3,3,2,7,1,5,1,6,1,8,3,6,7,5,3,8,5,5,6,6,12,3,1,3,1,8,1,7,5

%N和{p=素数}1/（p*log p）^2的十进制展开式。

%C通过将arXiv:0811.4739的形式扩展为Riemann zeta函数上的二重积分得到。

%H R.J.Mathar，<a href=“https://arxiv.org/abs/0811.4739（网址：https://arxiv.org/abs/0811.4739）“>素数zeta函数某些积分的二十位数</a>，arXiv:0811.4739（2008-2009）。

%e 1/A016627^2+1/A016650^2+1/8.047189^2+…=0.637056184074676....

%t位数=106；精度=数字+10；

%t tmax=500；（*被积函数在tmax之外可忽略不计*）

%t kmax=300；（*f（k）在kmax之外可忽略不计*）

%t InLogZeta[k_]：=NIntegrate[（t-2k）Log[Zeta[t]]，{t，2k，tmax}，工作精度->精度，MaxRecursion->20，AccuracyGoal->precision]；

%tf[k_]：=与[{mu=MoebiusMu[k]}，如果[mu==0，0，（mu/k^3）*InLogZeta[k]]]；

%t s=0；

%t做[s=s+f[k]；打印[k，“”，s]，{k，1，kmax}]；

%t真实数字[s][[1]][[1；；数字]]（*_Jean-François Alcover_，2022年6月21日，在_Vaclav Kotesovec_*之后）

%o（PARI）默认值（realprecision，200）；s=0；对于（k=1300，s=s+moebius（k）/k^3*intnum（x=2*k，[1]，1]，（x-2*k）*log（zeta（x））；打印）；\\_Vaclav Kotesovec_，2022年6月12日

%Y参见A137245、A115563、A221711、A319231。

%K nonn，cons公司

%0、1

%A R.J.Mathar，2018年9月14日

%E更多条款，来自_Vaclav Kotesovec_，2022年6月12日