OFFSET
1,1
COMMENTS
This sequence was inspired by e-mail from Ray Chandler, Nov 07 2007
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = A004215(n)^3. [From Ray Chandler, Jan 29 2009]
MAPLE
N:= 10^10: # to get all terms <= N
sort([seq(seq(4^(3*i) * (8*j + 7)^3, j = 0 .. floor((N^(1/3)/4^i - 7)/8)), i = 0 .. floor(log[4](N^(1/3))))]); # Robert Israel, Dec 26 2017
MATHEMATICA
b = Table[x^3, {x, 1, 300}]; a = {}; Do[Do[Do[AppendTo[a, (x^2 + y^2 + z^2)^3], {x, 0, 30}], {y, 0, 30}], {z, 0, 30}]; Union[a]; Complement[b, a] (*Artur Jasinski*)
Select[Range[200]^3, SquaresR[3, #]==0&] (* Harvey P. Dale, Feb 03 2015 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Nov 07 2007
EXTENSIONS
Extended by Ray Chandler, Jan 29 2009
STATUS
approved