%I#19 2021年4月15日05:25:01
%S 0,1,2,2,3,4,5,5,6,7,8,9,10,10,11,11,12,13,14,15,16,16,17,17,
%电话:18,19,19,20,21,22,23,23,24,25,26,26,27,28,29,30,31,32,32,32,
%U 33,34,35,35,36,36,37,38,39,39,40,41,41,42,43,44,45
%N a(N)是不超过N的非光滑数(“锯齿”数)的数量。这是A064052的计数函数。
%C此序列不同于移位的A072490,在22个术语后。
%D Steven R.Finch,《数学常数》,剑桥大学出版社,2003年,第2.21章,第166页。
%D Daniel H.Greene和Donald E.Knuth,《算法分析数学》,第三版,Birkhäuser,1990年,第95-98页。
%F自2019年11月7日起(开始)
%F a(n)=总和{i=1..floor(sqrt(n))}(pi(floor(n/i))-pi(i))。
%F a(n)=和{p<=平方(n)}(p-1)+和{sqrt(n)<p<=n}楼层(n/p),其中p是素数。
%F a(n)=n-A064775(n)。(结束)
%F a(n)~log(2)*n-A153810*n/log(n)-A242610*n/log(n)^2+O(n/log(n)^3)(Greene和Knuth,1990)。-_Amiram Eldar_,2021年4月15日
%t jaggedQ[n_]:=jaggedQ[n]=(f=系数整数[n][[All,1]];s=Sqrt[n];计数[f,p_/;p>s]>0);a[n]:=(对于[cnt=0;j=2,j<=n,j++,If[jaggedQ[j],cnt++]];cnt);表[a[n],{n,1100}]
%Y参见A064052、A064775、A072490、A153810、A242610。
%K nonn公司
%氧1,3
%A _Jean-François Alcover,2014年5月16日
|