%I#38 2024年1月6日13:10:35
%S 0,1,4,5,9,13,14,16,25,29,30,36,41,49,50,54,55,61,64,77,81,85,86,90,
%电话:91100110113121126135139140144149169174181190194,
%电话:196199203204221225230245255256265271280284289294302
%N是连续平方和的数字。
%C此外,任何一对方形金字塔数的差异(A000330)。这些可以被称为“截断的正方形金字塔数”_富兰克林·亚当斯·沃特斯,2006年11月29日
%C如果n是截至m^2的d个连续平方和,n=A000330(m)-A000330(m-d)=d*(m^2-(d-1)m+(d-1。此外,A000330(x)>=x^3/3,因此m和d不能大于(3n)^(1/3)_M.F.Hasler,2024年1月2日
%H T.D.Noe,n的表格,n=1..10000的a(n)</a>
%H<a href=“/index/Su#ssq”>与平方和相关的序列的索引项</a>
%e所有方块(A000290:0,1,4,9,…)都在这个序列中,因为定义中的“连续”是指没有中断的子序列,所以一个术语合格。
%e5=1^2+2^2=A000330(2。
%t nMax=1000;t={0};Do[k=n;s=0;While[s=s+k^2;s<=nMax,AppendTo[t,s];k++],{n,平方[nMax]}];t=工会[t](*_t.D.Noe_,2012年10月23日*)
%o(哈斯克尔)
%o导入数据。Set(deleteFindMin,union,fromList);导入数据。列表(inits)
%o a034705 n=a034705_列表!!(n-1)
%o a034705_list=f 0(尾部$inits$a000290_list)(fromList[0]),其中
%o f x vss’@(vs:vss)s
%o | y<x=y:f x vss
%o|否则=f w vss(并集s$fromList$scanl1(+)ws)
%o其中ws@(w:_)=反向vs
%o(y,s')=删除查找最小值
%o--_Reinhard Zumkeller_,2015年5月12日
%o(PARI){is_A034705(n)=表示(d=1,sqrtnint(n*3,3),my(b=(d-1)/2,s=n/d-(d-1*(d*2-1)/6+b^2);分母(s)==分母(b)^2&&issquare(s,&s)&return(b+s)));!n}\\返回和的最大平方的索引(如果n在序列中,则返回1表示n=0),否则返回0。-_M.F.Hasler,2024年1月2日
%o(Python)
%o导入heapq
%o从itertools导入islice
%o def agen():#术语生成器
%o m=0;h=[(m,0,0)];nextcount=1;v1=无
%o为True时:
%o(v,s,l)=堆q堆p(h)
%如果v!=,则为ov1:产量v;v1=v
%o如果v>=m:
%o m+=下一计数*下一计数
%o堆堆(h,(m,1,nextcount))
%o下一个计数+=1
%o v-=s*s;s+=1;l+=1;v+=l*l
%o堆堆(h,(v,s,l))
%o打印(列表(islice(agen),60))#_Michael S.Branicky_,2024年1月6日
%Y参考A000290、A000330、A034706。
%Y参见A217843-A217850(3到10的连续幂之和)。
%Y参考A368570(此序列中每对连续整数的第一个)。
%K nonn公司
%氧1,3
%A _弗里德曼_
%E术语a(1..10^4)由M.F.Hasler_于2024年1月2日用独立代码进行双重检查
|