%I#50 2021年9月28日09:16:20
%S 9,28,35,65,72,911261331521892172243280344351370407,
%电话:4685135205395595767728730737567938548559451001,
%电话:1008102710641072112512161241133213391343135813951456151215471674
%N个数,是两个不同的正立方体的和。
%这个序列不包含素数,因为x^3+y^3=(x^2-x*y+y^2)*(x+y)_M.F.Hasler,2008年4月12日
%C没有术语==3、4、5或6 mod 9。-_罗伯特·伊斯雷尔,2014年10月7日
%C a(n)第2版:{1,0,1,1,0,1,0,0,1,1,1,0,1,0,1,0
%H Chai Wah Wu,n表,n=1..100000的a(n)(M.F.Hasler的术语1..902)
%H<a href=“/index/Su#ssq”>与立方体和相关的序列索引</a>
%e9的顺序是2^3+1^3=9。
%e 35的顺序是3^3+2^3=35。
%p N:=10000:#获取所有术语<=N
%p S:=选择(`<=`,{seq(seq(i^3+j^3,j=1..i-1),i=2..楼层(N^(1/3)))},N);
%p#如果使用Maple 11或更早版本,请取消注释下一行
%p#排序(转换(S,列表));
%p#_Robert Israel_,2014年10月7日
%t lst={};Do[Do[x=a^3;Do[y=b^3;If[x+y==n,AppendTo[lst,n]],{b,Floor[(n-x)^(1/3)],a+1,-1}],{a,Floor[n^(1/3)],1,-1}],{n,6!}];第1期(*_Vladimir Joseph Stephan Orlovsky,2009年1月22日*)
%t选择[Range@1700,Total@Boole@Map[And[!MemberQ[#,0],UnsameQ@@#]&,PowersRepresentations[#,2,3]]>0&](*_Michael De Vlieger_,2017年5月13日*)
%o(PARI)是A024670(n)=表示(i=ceil(sqrtn(n\2+1,3)),sqrtn(n-.5,3),isA000578(n-i^3)&return(1))/*也可以使用“表示(i=2,sqartn(n\2-1,3)…)”,但这要慢得多,因为[n/2,n]中的立方体比[1,n/2]中的要少。将这里的-1替换为+.5将得到A003325,并考虑到a(n)=x^3+x^3。将-1替换为0可能会由于舍入错误而错过此形式的一些a(n)。-_M.F.Hasler,2008年4月12日*/
%o(Python)
%o从itertools导入计数,takewhile
%o定义缺陷(极限):
%o cbs=列表(takewhile(λx:x<=极限,(i**3表示计数(1)中的i))
%o sms=集合(c+d代表i,c代表枚举(cbs)代表cbs中的d[i+1:])
%o如果s<=限制,则返回已排序的(s代表短信中的s)
%o印刷(aupto(1674))#_Michael S.Branicky_,2021年9月28日
%Y另请参见:2个正立方体的和(不一定是不同的):A003325。3个不同正数立方体的总和:A024975。不同正立方体的总和:A003997。两个不同非负立方体的和:A114090。2个非负立方体的和:A004999。两个不同正方形的和:A004431。立方体:A000578。
%K非n
%O 1,1号机组
%百灵鸟金伯利_
%E Name编辑人:Zak Seidov,2011年5月31日
|