OFFSET
1,1
COMMENTS
The corresponding primes are given in A046842.
a(27) > 50000. - Tyler NeSmith, Apr 16 2022
a(30) > 10^5. - Michael S. Branicky, Sep 11 2024
MAPLE
A046843:=n->`if`(isprime(n*2^n-n-1), n, NULL): seq(A046843(n), n=1..1000); # Wesley Ivan Hurt, Oct 12 2014
MATHEMATICA
Select[Range[10000], PrimeQ[# 2^#-#-1]&] (* Harvey P. Dale, Oct 10 2011 *)
PROG
(PARI) is(n)=ispseudoprime(n*2^n-n-1) \\ Charles R Greathouse IV, May 22 2017
CROSSREFS
KEYWORD
nonn,more
AUTHOR
EXTENSIONS
Sign in definition corrected by R. J. Mathar, Oct 10 2010
a(10)-a(21) from Harvey P. Dale, Oct 10 2011
a(22)-a(26) from Tyler NeSmith, Apr 16 2022
a(27)-a(28) from Michael S. Branicky, May 17 2023
a(29) from Michael S. Branicky, Sep 11 2024
STATUS
approved