|
|
A020894
|
|
Nonnegative numbers that are sums of two nonzero cubes.
|
|
2
|
|
|
0, 2, 7, 9, 16, 19, 26, 28, 35, 37, 54, 56, 61, 63, 65, 72, 91, 98, 117, 124, 126, 127, 128, 133, 152, 169, 189, 208, 215, 217, 218, 224, 243, 250, 271, 279, 280, 296, 316, 331, 335, 341, 342, 344, 351, 370, 386, 387, 397, 407, 432, 448, 468, 469
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
When calculating terms, there is no need to search beyond a value x defined by x^3 - (x-1)^3 = n. The positive solution is given by x = 1/2 + (sqrt(12n-3))/6.
There are no cubes in this sequence, but the numbers before and after a cube are all included. (End)
|
|
LINKS
|
|
|
EXAMPLE
|
7 is in the sequence because 2^3 + (-1)^3 = 7
8 is not in the sequence because the only solutions to x^3 + y^3 = 8 have either x=0 or y=0. (End)
|
|
MATHEMATICA
|
Reap[For[n = 0, n < 500, n++, fi = FindInstance[x > 0 && y != 0 && n == x^3 + y^3, {x, y}, Integers, 1]; If[fi =!= {}, Print[n, " = ", Hold[x^3 + y^3] /. fi[[1]]]; Sow[n]]]][[2, 1]] (* Jean-François Alcover, Nov 05 2016 *)
|
|
PROG
|
(PARI) isA020894(n) = {r=0; x=1.0/2+sqrt(12*n-3.0)/6; for(i=1, floor(x), if(ispower(n-i^3, 3) & (n != i^3), r++)); r>0}; \\ Michael B. Porter, Oct 16 2009
(PARI) T=thueinit('z^3+1);
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|