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