OFFSET
1,2
LINKS
Robert Israel, Table of n, a(n) for n = 1..8446
MAPLE
N:= 10000: # for members of A003072 up to N
S:= {}:
for x from 1 while 3*x^3 <= N do
for y from x while x^3+2*y^3<=N do
for z from y do
v:= x^3+y^3+z^3;
if v > N then break fi;
S:= S union {v}
od od od:
S:= sort(convert(S, list)):
select(i -> isprime(S[i]), [$1..nops(S)]); # Robert Israel, Aug 27 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
More terms from Sean A. Irvine, Jul 30 2019
STATUS
approved