对数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））；}/*计算sum_{pprime}1/（p^slog（p））*/SumEulerlog（s，N=最大值（2,30/abs（s））={my（B=getlocalbitprec（），S=0，LN，T，lim）；localbitprec（32）；LN=对数（N）；lim=天花板（B*log（2）/LN）；本地比特币（B+32）；对于素数（p=2，N，S+=1/（p^S*log（p）））；LN=位精度（LN，B+32）；T=整数（0，[oo，1]）；无石英（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，对数N（k*（t+S），N），Tk）；返回（S）；}