%I#26 2022年9月8日08:45:34
%S 0,1,4,8,9,16,17,25,36,40,49,52,64,72,73,80,81,89,97100116121136,
%电话:14414514816918019319622523223324124425626028289292,
%电话:30531332433736036136938840044409424449457
%N整数k,使k和k+1都是2个平方的和。
%C等价地,非负k,使得k*(k+1)是两个平方的和。
%另外,非负k,使得k*(k+1)/2是两个平方的和。这很容易从“两个平方和定理”得到:x是两个(非负)平方的和,如果它的素因式分解不包含p^e,其中p==3(mod 4),e是奇数_罗伯特·伊斯雷尔,2018年3月26日
%一般来说,序列包括所有正方形。
%H David A.Corneth,<A href=“/A140612/b140612.txt”>n,A(n)表,n=1.-10000</a>
%e 40=6^2+2^2,41=5^2+4^2,所以40在序列中。
%t(*M6*)A1={};做[If[SquaresR[2,n(n+1)/2]>0,AppendTo[A1,n]],{n,0,1500}];A1类
%t连接[{0},展平[Position[Accumulate[Range[500]],_?(平方R[2,#]>0&)]](*哈维·P·戴尔,2015年6月7日*)
%t序列位置[表[If[SquaresR[2,n]>0,1,0],{n,0500}],{1,1}][[All,1]-1(*H arvey P.Dale_,2021年7月28日*)
%o(岩浆)[0..460]中的k:k |对于[0,1]中的所有{k+a:a |范数方程(1,k+a)eq true}];//_Marius A.Burtea,2019年10月8日
%o(Python)
%o从itertools导入count、islice、starmap
%o来自sympy进口保理商
%o定义A140612_gen(startvalue=0):#术语生成器>=startvalue
%o表示计数中的k(最大值(startvalue,0)):
%o if all(星图(lambda d,e:e%2==0或d%4!=3,因子(k*(k+1)).items())):
%o产量k
%o A140612_list=list(岛屿(A140612_gen(),20))#_Chai Wah Wu_,2022年3月7日
%Y参见A000217、A000290、A000404、A001481、A002378、A050795、A073613。
%K nonn,简单
%氧1,3
%A _弗兰克林·T·亚当斯-沃特斯,2008年5月19日
%E a(1)=0由_Max Alekseyev_于2019年10月8日预先编制和编辑