login
Primes p which equal prime(2^k) - prime(2^(k-1)) + 1 for some k.
1

%I #26 Jul 15 2021 02:49:30

%S 2,5,13,79,181,409,2053,21011,96487,1235199066139957,

%T 180572921869744273,747909478827063599

%N Primes p which equal prime(2^k) - prime(2^(k-1)) + 1 for some k.

%e Prime(2^12)-prime(2^11) + 1 = 21011 is prime.

%t Select[#[[2]]-#[[1]]+1&/@Partition[Prime[2^Range[0,30]],2,1],PrimeQ] (* _Harvey P. Dale_, Jul 08 2021 *)

%o (PARI) f(n)= for(x=1,n,y=prime(2^x)-prime(2^(x-1))+1;if(isprime(y),print1(y",")))

%Y Cf. A345982 (values of k), A033844.

%K nonn,more,hard

%O 1,1

%A _Cino Hilliard_, Feb 20 2004

%E Corrected and edited by _N. J. A. Sloane_, following a suggestion from _Harvey P. Dale_, Jul 07 2021

%E a(10)-a(12) from _Mohammed Yaseen_, Jul 08 2021