|
|
A247789
|
|
Primes p such that sigma(2p-1) is a prime q.
|
|
2
|
|
|
5, 13, 1201, 1741, 2521, 5101, 8581, 14281, 41761, 139921, 353641, 595141, 697381, 1460341, 1639861, 3723721, 3889261, 6093541, 7956061, 11143921, 14199121, 33874681, 46938361, 51521401, 56063461, 62395621, 67523821, 90706981
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Corresponding values of primes q for a(n) are in A247836.
|
|
LINKS
|
|
|
EXAMPLE
|
Prime 13 is in sequence because sigma(2*13-1) = sigma(25) = 31 (prime).
|
|
PROG
|
(Magma) [p: p in PrimesUpTo(50000000) | IsPrime(SumOfDivisors(2*p-1))]
(PARI) forprime(p=1, 10^9, if(isprime(sigma(2*p-1)), print1(p, ", "))) \\ Derek Orr, Sep 25 2014
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|