数字=106;精度=数字+10;
tmax=500;(*被积函数在tmax之外可忽略不计*)
kmax=300;(*f(k)在kmax之外可忽略不计*)
InLogZeta[k_]:=NIntegrate[(t-2k)Log[Zeta[t]],{t,2k,tmax},工作精度->精度,最大递归->20,精度目标->精度];
f[k_]:=与[{mu=MoebiusMu[k]},如果[mu==0,0,(mu/k^3)*InLogZeta[k]]];
s=0;
Do[s=s+f[k];打印[k,“”,s],{k,1,kmax}];
实数字[s][[1]][[1;;数字]](*让-弗朗索瓦·奥尔科弗2022年6月21日之后瓦茨拉夫·科特索维奇*)