OFFSET
1,1
COMMENTS
Also primes p that are the sum of the divisors of some n where 2^sigma(n) - 1 is a Mersenne prime (A000668).
If a(6) exists, it must be greater than A000043(48) = 57885161, and also not equal to any of the Mersenne prime exponents 74207281, 77232917, 82589933, 136279841. - Gord Palameta, Oct 22 2024
EXAMPLE
Mersenne exponent 7 is in the sequence because sigma(4) = 7.
Mersenne exponent 31 is in the sequence because there are two numbers n (16 and 25) with sigma(n) = 31.
PROG
(Magma) Set(Sort([SumOfDivisors(n): n in[1..10000] | IsPrime((2^SumOfDivisors(n))- 1)]));
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Jaroslav Krizek, Jan 16 2015
STATUS
approved