OFFSET
1,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
EXAMPLE
Sums not divisible by 17 are indicated in asterisks:
....|...1....8...27...64...125...216...343...512...729..1000..1331
------------------------------------------------------------------
1...|...*....*....*....*.....*.....*.....*.....*.....*.....*.....*
8...|...*....*....*....*.....*.....*.....*.....*.....*.....*.....*
27..|...*....*....*....*.....*.....*.....*.....*.....*.....*.....*
64..|...*....*....*....*.....*.....*.....*.....*.....*.....*.....*
125.|...*....*....*....*.....*.....*.....*.....*.....*.....*.....*
216.|...*....*....*....*.....*.....*.....*.....*.....*.....*..1547
343.|...*....*....*....*.....*.....*.....*.....*.....*..1343.....*
512.|...*....*....*....*.....*.....*.....*.....*...1241....*.....*
729.|...*....*....*....*.....*.....*.....*..1241.....*.....*.....*
1000|...*....*....*....*.....*.....*..1343.....*.....*.....*.....*
1331|...*....*....*....*.....*..1547.....*.....*.....*.....*.....*
MATHEMATICA
upto[n_] := Block[{t}, Union@ Reap[ Do[If[Mod[t = x^3 + y^3, 17] == 0, Sow@t], {x, n^(1/3)}, {y, Min[x, (n - x^3)^(1/3)]}]][[2, 1]]]; upto[47000] (* Giovanni Resta, Jun 12 2020 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Jun Mizuki (suzuki32(AT)sanken.osaka-u.ac.jp), Nov 15 2004
STATUS
approved