OFFSET
1,1
LINKS
Zak Seidov, Table of n, a(n) for n = 1..1000
EXAMPLE
29 = 2^2 + 3^2 + 4^2 = 4 + 9 + 16.
89 = 2^2 + 6^2 + 7^2 = 4 + 36 + 49; also 89 = 3^2 + 4^2 + 8^2 = 9 + 16 + 64.
353 = 2^2 + 5^2 + 18^2 = 4 + 25 + 324; also 353 = 4^2 + 9^2 + 16^2 = 16 + 81 + 256.
PROG
(PARI) {m=22; p=m^2; v=vector(m, x, x^2); w=[]; for(i=1, m-2, for(j=i+1, m-1, for(k=j+1, m, if((n=v[i]+v[j]+v[k])<p&&isprime(n), w=concat(w, n))))); w=listsort(List(w), 1); for(j=1, #w-1, print1(w[j], ", "))} /* Klaus Brockhaus, Feb 11 2007 */
CROSSREFS
KEYWORD
nonn
AUTHOR
Tomas Xordan, Jan 21 2007
EXTENSIONS
Edited, corrected and extended by Klaus Brockhaus, Feb 11 2007
STATUS
approved