%I#10 2021年1月15日06:04:51

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

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

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

%乘积{素数p==5（mod 6）}1/（1+1/p^2）的十进制展开式。

%电子邮箱：0.94450093450470009867342910941914443425461104807808690667695735771。。。

%t位数=105；

%t精度=数字+5；

%t prodeuler[p_，a_，b_，expr_]：=乘积[If[a<=p<=b，expr，1]，{p，素数[Range[PrimePi[a]，PrimePi[b]]}]；

%t Lv3[s_]：=prodeuler[p，1，2^（精度/s），1/（1-KroneckerSymbol[-3，p]*p^-s）]//N[#，精度]&；

%t Lv4[s_]：=2*Im[PolyLog[s，Exp[2*I*Pi/3]]/Sqrt[3]；

%t Lv[s_]：=如果[s>=10000，Lv3[s]，Lv4[s]]；

%t gv[s_]：=（1-3^（-s））*Zeta[s]/Lv[s]；

%t pB=（3/4）*乘积[gv[2^n*2]^（2^-（n+1）），{n，0，11}]//n[#，精度]&；

%t pD=（45*pB*Lv[2]）/（4*Pi^2）；

%t RealDigits[pD，10，digits][[1]（*此代码大部分是由于_Artur Jasinski_*）

%t（*-----------------------------------------------------------------------*）

%tS[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++]；总额）；

%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}]；

%t Z[m_，n_，s_]：=（w=1；sumz=0；difz=1；当[Abs[difz]>10^（-数字-5）时，difz=P[m，n，s*w]/w；sumz=sumz+difz；w++]；实验[sumz]）；

%t$MaxExtraPrecision=1000；数字=121；真实数字[Chop[N[Z[6,5,4]/Z[6,15,2]，数字]]，10，数字-1][[1]（*_Vaclav Kotesovec_，2021年1月15日*）

%Y相似常数：A175646、A175647、A248930、A24893、A301429、A333240、A334826、A335963、A340576、A34057。

%K nonn，cons公司

%0、1

%A _Jean-François Alcover，2021年1月12日