OFFSET
1,1
LINKS
Harry J. Smith, Table of n, a(n) for n = 1..1000
EXAMPLE
gcd(sigma(12), sigma(13)) = gcd(28, 14) = 14 > 12, so 12 is in the sequence.
MATHEMATICA
Select[Range[1, 10^5], GCD[DivisorSigma[1, # ], DivisorSigma[ 1, # + 1]] > # &]
PROG
(PARI) { n=0; s=1; for (m=1, 10^9, t=s; s=sigma(m + 1); if (gcd(t, s) > m, write("b066025.txt", n++, " ", m); if (n==1000, return)) ) } \\ Harry J. Smith, Nov 07 2009
CROSSREFS
KEYWORD
nonn
AUTHOR
Joseph L. Pe, Dec 11 2001
EXTENSIONS
Definition and Example edited by Harry J. Smith, Nov 07 2009
STATUS
approved