|
kmax=25；清除[P]；Do[P[k]=乘积[P=素数[n]；如果[Mod[p，4]==3，1-2/（p*（p-1）^2）//N[#，40]&，1]，{N，1，2^k}]；打印[“P（”，k，“）=”，P[k]]，{k，10，kmax}]；β=加泰罗尼亚语^2*（2/3）*P[kmax]；RealDigits[beta，10，16]//第一个
(* -------------------------------------------------------------------------- *)
$MaxExtraPrecision=1000；数字=121；
f[p]：=（1-2/（p*（p-1）^2））；
coefs=Rest[CoefficientList[Series[Log[f[1/x]]，{x，0，1000}]，x]]；
S[m_，n_，S_]：=（t=1；总和=0；difs=1；当[Abs[difs]>10^（-数字-5）||difs==0，difs=（MoebiusMu[t]/t）*Log[If[S*t==1，DirichletL[m，n，S*t]，总和[Zeta[S*t，j/m]*Dirichlet字符[m，n，j]^t，{j，1，m}]/m^（S*t）]]；总和=总和+difs；t++]；总额）；
P[m_，n_，s_]：=1/EulerPhi[m]*和[Conjugate[DirichletCharacter[m，r，n]]*s[m，r，s]，{r，1，EulerPhi[m]}]+和[If[GCD[P，m]>1&&Mod[P，m]==n，1/P^s，0]，{P，1，m}]；
m=2；集水坑=0；difp=1；当[Abs[difp]>10^（-数字-5）||difp==0时，difp=coefs[[m]]*P[4，3，m]；集水坑=集水坑+difp；m++]；
RealDigits[Chop[N[2*Catalan^2/3*Exp[simp]，数字]]，10，数字-1][[1](*瓦茨拉夫·科特索维奇2021年1月16日*）
| 