对数N（s，N）={my（B=getlocalbitprec（））；B+=细胞（abs（实数）-1）*指数（2*N））；localbitprec（B）；如果（位精度<B，s=位精度，B）；log（zeta（s）*prodeuler（p=2，N，1-p^（-s））；}/*计算和{p素数}log（log（p））/（p^s log（p））*/总和Eulerloglog（s，N=最大值（2,30/abs（s））={my（B=getlocalbitprec（），S=0，LN，T，lim，E）；localbitprec（32）；LN=对数（N）；lim=天花板（B*log（2）/LN）；本地比特币（B+32）；对于素数（p=2，N，S+=log（log（p））/（p^S*log（p））；LN=位精度（LN，B+32）；T=整数（0，[oo，1]）；E=欧拉（）；无石英（K=1，lim，我的（[k]=k，m=moebius（k）/k，a=1/（k*LN））；/*Tk=整数（0，[oo，1/a]）*/my（Tk=向量（#T，i，如果（i==1，T[1]，T[i]*a））；S-=m*intnum（t=0，oo，（log（t）+E）*LogzetaN（k*（t+s），N），Tk）；返回（S）；}SumEulerloglog（1）