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