OFFSET
1,3
PROG
(PARI) ok(n)={my(i=sqrtint(n)); i=i-(i%2==0); while(i>0, my(a=i^2, j=i); while(j>0 && a<=n, if(a==n, return(1)); j-=2; a=a+j^2); i-=2); 0}
concat([0], select(ok, [1..1200])) \\ Antonio Roldán, Mar 12 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Dec 21 2018
STATUS
approved