OFFSET
1,1
COMMENTS
Let W = 2^p - 9 and s = (W+7)/(2*p), then 3^s == 4 (mod W) for terms 1..6.
MATHEMATICA
Select[Prime[Range[3, 565]], PrimeQ[2^#-9]&] (* The program generates the first five terms of the sequence. *) (* Harvey P. Dale, Aug 24 2024 *)
PROG
(PARI)
forprime(p=5, 10^5, W= 2^p-9; if(ispseudoprime(W), print1(p, ", ")))
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Dmitry Ezhov, Mar 07 2017
STATUS
approved