|
|
A066025
|
|
Numbers k such that gcd(sigma(k), sigma(k+1)) > k.
|
|
1
|
|
|
5, 12, 14, 70, 88, 125, 180, 204, 206, 220, 957, 1253, 1334, 1364, 1634, 1673, 1750, 1888, 1919, 2685, 2759, 2958, 2974, 3127, 4364, 5191, 7615, 8142, 8632, 9114, 11219, 12000, 12035, 14664, 14841, 18414, 18762, 18873, 19358, 20118, 20145, 20712, 24957, 25194
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
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
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|