%I#18 2022年1月1日09:52:55

%S 2,16,54,72128250432520576686756102414581944200026623456，

%电话41604394460854860486686750756081929000982616641371814040，

%电话：1555215750160001852219656197102041221296243246962764831250328323280

%N是两个立方除数之和的数字（不一定是不同的）。

%C所有条款都是平等的。-_阿洛伊斯·海因茨，2020年8月30日

%C形式为a^3+b^3的数字，其中a除以b^3，b除以a^3。-_罗伯特·伊斯雷尔，2020年11月1日

%H Robert Israel，n的表，a（n）表示n=1..10000</a>

%e16在序列中，因为2是16的除数，2^3+2^3=16。

%pq:=n->（s->ormap（x->n-x in s，s））（映射（x->x^3，numtheory[除数]（n））：

%p选择（q，[2*i$i=1..17000]）[]；#_阿洛伊斯·海因茨，2020年8月30日

%p N:=50000:#获得项<=N

%p R：=空：

%p表示从1到地板（N^（1/3））do

%pBs:=选择（b->b^3+a^3<=N和b^3 moda=0，数字理论：-除数（a^3））；

%pR:=R联合映射（b->a^3+b^3，Bs）；

%日期：

%p排序（转换（R，列表））；#_罗伯特·伊斯雷尔，2020年11月1日

%t M=40000；（*获取术语<=M*）

%t R={}；

%t对于[a=1，a<=楼层[M^（1/3）]，a++，

%t Bs=选择[Divisors[a^3]，#^3+a^3<=M&&Mod[#^3，a]==0&]；

%t R=联合[R，a^3+#^3&/@Bs]]；

%t排序[R]（*Jean-François Alcover_，2022年1月1日，在_Robert Israel_*之后）

%Y参考A094147（两个平方除数的和）。

%K nonn公司

%O 1,1号机组

%A_Wesley Ivan Hurt_，2020年8月30日

%E a（12）-a（44）摘自_Alois P.Heinz，2020年8月30日