%I#19 2024年8月13日02:13:13
%S 1,2,3,3,4,5,5,6,7,7,7,7,8,9,8,9,9,10,11,12,13,13,14,15,14,
%电话:15,13,13,13,13,14,15,16,15,16,17,18,19,19,20,19,20,21,20,21,21,
%U 22,23,23,232,24,25,25,23,24,25,25,26,27,27,28,29,29,30
%N 2n和第N个无平方数之间的差值:a(N)=2n-A005117(N)。
%H Antti Karttunen,n的表,n=1..10001的a(n)</a>
%F a(n)=2n-A005117(n)。
%F a(n)=n-A243348(n)。
%F a(n)~c*n,其中c=2-Pi^2/6(A152416)。-阿米拉姆·埃尔达尔,2024年3月4日
%t模块[{sf=Select[Range[150],SquareFreeQ]},表[2n-sf[[n]],{n,长度[sf]}](*_哈维P.戴尔,2021年6月26日*)
%o(Python)
%o从数学导入isqrt
%o来自sympy import mobius
%o定义A243351(n):
%o定义f(x):对于范围(1,isqrt(x)+1)中的k,返回n+x-sum(mobius(k)*(x//k**2))
%o m,k=n,f(n)
%哦,而我!=k:
%o m,k=k,f(k)
%o返回(n<<1)-m#Chai Wah Wu,2024年8月12日
%Y参考A005117、A152416、A243348。
%K nonn公司
%O 1,2号机组
%2014年6月4日,安蒂·卡图宁