|
%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
|
