OFFSET
1,1
EXAMPLE
2 is in the sequence because (2^2 - 2 - 1)*2^2 + 1 = 5 is prime.
MATHEMATICA
lst={}; Do[If[PrimeQ[(2^n - n-1)*2^n+1], AppendTo[lst, n]], {n, 10000}]; lst
PROG
(PARI) is(n)=ispseudoprime((2^n-n-1)<<n+1) \\ Charles R Greathouse IV, Feb 17 2017
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
Michel Lagneau, Nov 30 2011
EXTENSIONS
a(20)-a(23) from Michael S. Branicky, Jul 12 2023
STATUS
approved