登录
OEIS由OEIS基金会的许多慷慨捐赠者.

 

标志
提示
(来自的问候整数序列在线百科全书!)
A024670号 两个不同的正立方体的和。 39

%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日

查找|欢迎光临|维基|寄存器|音乐|地块2|演示|索引|浏览|更多|网络摄像头
贡献新序列。或评论|格式|样式表|变换|超级搜索|最近
OEIS社区|维护人OEIS基金会。

许可协议、使用条款、隐私政策。.

上次修改时间:美国东部夏令时2024年5月8日02:29。包含372317个序列。(在oeis4上运行。)