OFFSET
1,1
COMMENTS
No term can be written as x^2 + y^2 + z^9.
LINKS
Colin Barker, Table of n, a(n) for n = 1..1000
W. Jagy and I. Kaplansky, Sums of Squares, Cubes and Higher Powers, Experimental Mathematics, vol. 4 (1995), p. 171 (see Theorem under Higher powers).
MAPLE
P := select(p -> isprime(p), [seq(n, n=5..1000, 4)]):
seq((6*p)^3, p in P); # Peter Luschny, Jun 22 2018
MATHEMATICA
P = Select[Range[5, 300, 4], PrimeQ];
A305728 = (6P)^3 (* Jean-François Alcover, Jun 22 2018 *)
PROG
(Magma) [216*p^3: p in PrimesUpTo(300) | IsOne(p mod 4)];
(PARI) first(n) = my(res=List()); forprime(p=5, oo, if(p%4 == 1, listput(res, (6*p)^3); n--; if(n==0, return(res)))) \\ David A. Corneth, Jun 27 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Bruno Berselli, Jun 22 2018
STATUS
approved