OFFSET
1,2
EXAMPLE
6111+sigma(6111)+sigma(sigma(6111))=T(286).
MATHEMATICA
triangularQ[n_] := IntegerQ[Sqrt[8n + 1]]; Reap[For[n = 1, n < 12000, n++, s = DivisorSigma[1, n]; If[triangularQ[n + s + DivisorSigma[1, s]], Sow[n]]]][[2, 1]] (* Jean-François Alcover, Oct 15 2013 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Giovanni Resta, Feb 13 2006
STATUS
approved