A247791 Primes p such that there is a prime q for which sigma(sigma(2*q-1)) = p.
7, 131071, 524287
Offset: 1
Examples
Prime 7 is in sequence because there is prime 2 such that sigma(sigma(2*2-1)) = sigma(sigma(3)) = sigma(4) = 7.
Programs
-
Magma
[SumOfDivisors(SumOfDivisors(2*n-1)): n in [A247790(n)]];
-
Magma
[SumOfDivisors(SumOfDivisors(2*n-1)): n in[1..1000000] | IsPrime(SumOfDivisors(SumOfDivisors(2*n-1))) and IsPrime(n)];
-
PARI
forprime(p=1,10^7,if(ispseudoprime(sigma(sigma(2*p-1))),print1(sigma(sigma(2*p-1)),", "))) \\ Derek Orr, Sep 29 2014
Comments