OFFSET
2,2
PROG
(PARI) a353387(nmax)={my(v=vectorsmall(2*nmax^2)); for(n=1, nmax, my(dfirst=0);
for(k=0, n, my(s=n^2+k^2); if(!v[s], if(!dfirst, print1(s, ", "); dfirst=1); v[s]++)))};
a353387(52)
CROSSREFS
KEYWORD
nonn
AUTHOR
Hugo Pfoertner, Apr 16 2022
STATUS
approved