%I M3787 N1543#22 2022年2月1日00:53:22
%S 5,9,21,37,69,69,891371774214816577498851085130513531489,
%电话:1861261726933125424957617129810994659465107171240112401,
%电话:162371623724833307253523746701474414744161493677976780567805
%N设A(N)=#{(i,j):i^2+j^2<=N},V(N)=Pi*N,P(N)=1(N)-V(N;A000099给出n的值,其中|P(n)|设置新记录;序列给出A(A000099(n))。
%D N.J.A.Sloane,《整数序列手册》,学术出版社,1973年(包括该序列)。
%D N.J.A.Sloane和Simon Plouffe,《整数序列百科全书》,学术出版社,1995年(包括该序列)。
%H David W.Wilson,n的表,n=1..200的a(n)</a>
%H W.C.Mitchell,<a href=“http://dx.doi.org/10.1090/S0025-5718-1966-0195834-3“>k维超球面中晶格点的数量,《数学与计算机》,20(1966),300-310。
%t nmax=3*10^4;A[n_]:=1+4*楼层[Sqrt[n]]+4*楼层[平方[n/2]]^2+8*总和[楼层[Squart[n-j^2],{j,楼层[n/2]]+1,楼层[Sqrt[n]]}];V[n_]:=Pi*n;P[n]:=A[n]-V[n];记录=0;A000099=收获[对于[k=0;n=1,n<=nmax,n++,p=Abs[p[n]];如果[p>记录,记录=p;k++;母猪[an=A[n]];打印[“a(”,k,“)=”,an];]]][[2,1]](*_Jean-François Alcover_,2016年2月7日*)
%Y参考A000099、A000036、A000092、A000413、A000223。
%K nonn公司
%O 1,1
%A _N.J.A.斯隆_
%E条目于2005年6月28日修订
|