%I #11 Jan 12 2022 15:13:14
%S 1,13,1468,2701,2959,3735,8686,11920
%N Numbers k such that (2^k + k + 1)*2^k - 1 is prime.
%e 13 is in the sequence because (2^13 + 13 + 1)*2^13 - 1 = 67223551 is prime.
%t lst={};Do[If[PrimeQ[(2^n + n+1)*2^n-1],AppendTo[lst,n]],{n,10000}];lst
%o (PARI) is(n)=isprime((2^n+n+1)<<n-1) \\ _Charles R Greathouse IV_, Feb 17 2017
%o (Python)
%o from sympy import isprime
%o def afind(limit, startk=1):
%o pow2 = 2**startk
%o for k in range(startk, limit+1):
%o if isprime((pow2 + k + 1)*pow2 - 1):
%o print(k, end=", ")
%o pow2 *= 2
%o afind(1500) # _Michael S. Branicky_, Jan 12 2022
%Y Cf. A201356, A201358, A201359, A201360, A201361, A201362, A201363.
%K nonn,hard,more
%O 1,2
%A _Michel Lagneau_, Nov 30 2011
%E a(8) from _Michael S. Branicky_, Jan 12 2022
|