OFFSET
1,1
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
MATHEMATICA
Module[{nn=100, pent}, pent=PolygonalNumber[5, Range[nn]]; Select[Flatten[ Outer[ Plus, pent, Range[Sqrt[Max[pent]]]^2]]//Union, IntegerQ[ (1+Sqrt[ 1+24#])/6]&]] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Apr 27 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Jonathan Vos Post, Jan 21 2008
STATUS
approved