A321084 Primes prime(n) such that 1 + Sum_{k=1..n} 2^(prime(k)-1) is prime.
2, 3, 5, 19, 2039, 2879
Offset: 1
Examples
a(3) = 5 since 1 + 2^(2-1) + 2^(3-1) + 2^(5-1) = 10111_2 = 23 is prime. Note that prime(3) = 5 and A080355(3+1) = 23 prime.
Programs
-
Mathematica
Prime@ Select[Range[10^3], PrimeQ[1 + Total@ Array[2^(Prime[#] - 1) &, #]] &] (* Michael De Vlieger, Oct 31 2018 *)
-
PARI
isok(p) = isprime(p) && isprime(1 + sum(k=1, primepi(p), 2^(prime(k)-1))); \\ Michel Marcus, Oct 27 2018
Extensions
a(5)-a(6) from Robert Israel, Oct 27 2018
Comments