%I#30 2024年1月3日10:46:42
%S 1,2,3,4,4,5,6,7,7,7,8,9,10,11,12,12,13,13,14,14,15,16,17,17,17,18,
%电话18,18,19,20,21,21,212,23,24,25,26,27,27,28,29,30,30,30-31,32,32,
%U 32,33,33,34,35,36,37,38,38,39,40,40,41,42,43,44,45,46,47
%长度为N的无平方间隙中平方素数的最小值。
%H Antti Karttunen,<a href=“/A107079/b107079.txt”>n,a(n)表,n=1..100000</a>
%H Louis Marmet,<a href=“http://www.marmet.org/louis/sqfgap网站/“>无平方间隙的首次出现及其计算算法</a>
%H L.Marmet,<a href=“http://arxiv.org/abs/1210.3829“>首次出现无平方间隙及其计算算法,arXiv预打印arXiv:1210.3829[math.NT],2012.-发件人:N.J.A.Sloane,1月1日
%F a(n)=和{k=0..n-1,moebius_mu(n-k-1)mod 2}。
%F a(n)=A013928(n+1)+A107078(n)。
%F From _Antti Karttune_,2016年10月7日:(开始)
%F a(n)=1+A013928(n)。[参考Charles R Greathouse IV的PARI计划。]
%F对于所有n>=1,a(A005117(n))=n。
%F(结束)
%t a[n_]:=和[Boole[SquareFreeQ[k]],{k,1,n-1}]+1;
%t数组[a,100](*_Jean-François Alcover_,2018年9月11日,摘自A013928*)
%o(PARI)A107079(n)=1+总和(k=1,n-1,比特和(moebius(k),1))
%o(Python)
%o从数学导入isqrt
%o来自sympy import mobius
%o def A107079(n):对于范围(1,isqrt(n-1)+1)中的k,返回1+总和(mobius(k)*((n-1,//k**2)#_Chai Wah Wu_,2024年1月3日
%Y一个大于A013928。A005117的左反转。
%Y参考A045882,A107078。
%K nonn公司
%O 1,2号机组
%A Paul Barry,2005年5月10日
%E来自_Charles R Greathouse IV_的新定义,2008年9月22日
|