OFFSET
1,2
LINKS
Eric Weisstein's World of Mathematics, Centered Square Number
MATHEMATICA
anmax = 1000; nmax = Floor[Sqrt[anmax/2]] + 1; Select[Union[Flatten[Table[Sum[k^2 + (k + 1)^2, {k, i, j}], {i, 0, nmax}, {j, i, nmax}]]], # <= anmax &] (* Vaclav Kotesovec, Dec 21 2018 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Dec 20 2018
STATUS
approved