OFFSET
1,2
COMMENTS
For each integer j in A059609, 2^(j-1)*(2^j - 7) is in the sequence. E.g., for j = A059609(1) = 39 we get 151115727449904501489664. - M. F. Hasler and Farideh Firoozbakht, Dec 03 2013
LINKS
Farideh Firoozbakht and M. F. Hasler, Variations on Euclid's formula for perfect numbers, Journal of Integer Sequences 13 (2010), 18 pp. Article ID 10.3.1.
EXAMPLE
Sigma(25) = 31 = 1*25 + 6, so 31 mod 25 = 6.
MATHEMATICA
Select[Range[1000000], Mod[DivisorSigma[1, #] - 6, #] == 0 &] (* T. D. Noe, Dec 03 2013 *)
PROG
(PARI) isok(n) = Mod(sigma(n), n) == 6; \\ Michel Marcus, Jan 03 2023
CROSSREFS
KEYWORD
nonn,more,hard
AUTHOR
Labos Elemer, Oct 29 2003
EXTENSIONS
Terms corrected by Charles R Greathouse IV and Farideh Firoozbakht, Dec 03 2013
STATUS
approved
