%I#10 2022年9月8日08:45:35

%S 9,28,35,65,72,911261331521892172243280344351370407，

%电话：468513520539559576637728855100110081027106411251216，

%电话：133213391343135813951456151215471674172917361755179218431853

%A031980的N互补。

%C A024670的子序列；A141806给出了A024670中不在此序列中的术语。

%C不是A001235的超序列；7094269是A001235中的最小数字，但不在此序列中（参见下面的第三个示例），下一个数字是11261376。

%H Klaus Brockhaus，n的表格，n=1..24834的a（n）</a>

%H<a href=“/index/Su#ssq”>与平方和和立方体和相关的序列索引</a>

%e9是两个完全不同的非零立方体的和：9=1^3+2^3。9不在A031980中，因为1和2是A031980。因此，9是这个序列的一个术语。

%e 1729是两个不同的非零立方体的和，有两种方式：1729=9^3+10^3=1^3+12^3。1729不在A031980中，因为1和12是A031980。因此1729是该序列的一个术语。

%e 7094269是两个不同的非零立方体的和，有两种方式：7094269=70^3+189^3=133^3+168^3。7094269在A031980中，因为它不是A031980.两个早期项的立方体之和；在第一种情况下189和第二种情况下133不是A031980的项。因此，7094269不是该序列的术语。

%t最大值=2000；A031980={1}；Do[m=天花板[（n-1）^（1/3）]；s=选择[A031980，#<=m&]；ls=长度[s]；sumOfCubes=并集[扁平[表[s[i]]^3+s[[j]]^3，{i，1，ls}，{j，i+1，ls}]]；如果[FreeQ[sumOfCubes，n]，AppendTo[A031980，n]]，{n，2，max}]；补体[Range[max]，A031980]（*Jean-François Alcover_，2013年9月3日*）

%o（岩浆）m:1853；a： =[]；a2:={}；对于[1..m]中的n，做p:=1；u： =a2连接{x:x在a}中；而u中的p表示p:=p+1；结束while；如果p gt m，则断裂；结束条件：；a2:=在|x^3+p^3 lem}中的a2连接{x^3+p2:3:x；追加（~a，p）；结束；打印a2；

%Y参见A141806，A031980（未出现在前面的最小数，不是两个不同的早期项的立方体之和），A024670（两个不同正整数的立方体之和）和A001235（两个立方体以多种方式的总和）。

%K nonn公司

%O 1,1号机组

%A _Klaus Brockhaus，2008年7月16日