%I#30 2022年6月23日12:32:43

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

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

%U 1,0,6,4,4,3,1,0,4,1,0,4,7,2,0,66,3,0,7,6,1,9,5,2,5,2,7,5,1,3,4,6,0

%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“>素数zeta函数某些积分的二十位数</a>，arXiv:0811.4739[math.NT]，2008-2018。

%e 1/（2*A253191）+1/（3*A175478）+1/1（5*2.59029…）+1/（7*3.7865）+…=1.52097043...

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

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

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

%t InLogZeta[k_]：=NIntegrate[（t-k）Log[Zeta[t]]，{t，k，tmax}，工作精度->精度，最大递归->20，精度目标->精度]；

%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=1500，s=s+moebius（k）/k^3*intnum（x=k，[1]，1]，（x-k）*log（zeta（x））；打印）；\\_Vaclav Kotesovec_，2022年6月12日

%Y参见A137245、A115563、A221711、A319232、A354917。

%K nonn，cons公司

%O 1,2号机组

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

%E更多数字来自_Vaclav Kotesovec_，2022年6月12日