OFFSET
1,1
LINKS
Donovan Johnson, Table of n, a(n) for n = 1..1000
EXAMPLE
4186 is in sequence because it is triangular (4186 = 91*92/2) and the sum of its triangular divisors, 4186+91+1 = 4278 is also triangular (4278 = 92*93/2).
MATHEMATICA
triQ[n_] := n > 1 && IntegerQ[Sqrt[8*n+1]]; q[n_] := triQ[1 + DivisorSum[n, #&, triQ[#] &]]; Select[Accumulate[Range[22000]], q] (* Amiram Eldar, Aug 12 2023 *)
PROG
(PARI) istriangular(n)=issquare(8*n+1)
{t=0; for(n=1, 10^8, if(istriangular(n), k=sumdiv(n, d, istriangular(d)*d) ; if(istriangular(k)&&k>>1, t+=1; write("b209310.txt", t, " ", n))))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Antonio Roldán, Jan 18 2013
STATUS
approved