OFFSET
1,2
EXAMPLE
7 is in the sequence because x^2-7*y^2=84 has integer solutions, including (x,y) = (14,4) and (28,10).
MATHEMATICA
Select[Range[300], Length[FullSimplify[Solve[x^2-#*y^2==#*(#+1)*(#+2)/6, {x, y}, Integers]/.C[1]->1]]>0&] (* Vaclav Kotesovec, Oct 05 2014, Mathematica version >= 8 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Colin Barker, Oct 04 2014
EXTENSIONS
More terms from Vaclav Kotesovec, Oct 05 2014
STATUS
approved