检验过的
经核准的
提出
编辑
返回m**2#柴华武2024年8月18日
(Python)
从数学导入isqrt
从sympy导入primepi,primerange
定义A085986号(n) :
def f(x):返回int(n+x+(t:=primepi(s:=isqrt(x)))+(t*(t-1)>>1)-sum(素数范围(1,s+1)中k的primepi(x//k))
m、 k=n,f(n)
而m!=克:
m、 k=k,f(k)
使用(数字理论):P:=进程(n)
如果nops(因子集(n))=2且bigomega(n)=2,则n^2;fi;
结束:seq(P(i),i=1..141)#保罗·拉瓦,2019年7月18日
Sum_{n>=1}1/a(n)=(P(2)^2-P(4))/2=(A085548号^2 -A085964号)/2=0.063767…,其中P是素数zeta函数。 -阿米拉姆·埃尔达尔2020年7月6日
(MAGMA公司岩浆)[1..150]|IsSquarefree(k)和#PrimeDivisors(k)eq 2]中的[k^2:k//马吕斯·A·伯蒂2019年10月24日
OEIS服务器: https://oeis.org/edit/global/2944