OFFSET
1,2
PROG
(PARI) fourkp1(n) = { forstep(x=1, n, 4, for(y=1, floor(sqrt(x)), z=x-y^2; if(issquare(z), print1(floor(sqrt(z))" "); break; ) ) ) }
CROSSREFS
KEYWORD
nonn
AUTHOR
Cino Hilliard, Oct 06 2003
STATUS
approved