%I#50 2021年9月28日09:16:20

%第9、28、35、65、72、911261331521892172242432803413443513370407页，

%电话：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替换为+0.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 nonn公司

%O 1，1

%百灵鸟金伯利_

%E Name编辑人：Zak Seidov，2011年5月31日