OFFSET
1,1
MATHEMATICA
r = 1456; lst = Table[0, {r}]; lim = Floor[Sqrt[8*r - 7]]; Do[num = (i^2 + i)/2 + (j^2 + j)/2; If[num <= r, lst[[num]]++], {i, lim}, {j, i - 1}]; Flatten@Position[lst, n_ /; n > 2]
Module[{nn=60, trnos}, trnos=Accumulate[Range[nn]]; Select[Sort[Flatten[ Table[ PositionIndex[Counts[Total/@Subsets[trnos, {2}]]][i], {i, 3, nn}]]], #<= Last[trnos]&]] (* The program uses the PositionIndex and Counts functions from Mathematica version 10 *) (* Harvey P. Dale, Dec 26 2015 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Arkadiusz Wesolowski, Dec 02 2015
STATUS
approved