OFFSET
1,3
MATHEMATICA
Do[ If[ PrimeQ[ 4*9^n + 1], Print[ n ]], {n, 0, 10^4}] (* Ryan Propper, Mar 25 2007 *)
Select[Range[0, 5000], PrimeQ[4 9^#+1]&] (* Harvey P. Dale, Dec 07 2011 *)
PROG
(PARI) is(n)=ispseudoprime(4*9^n+1) \\ Charles R Greathouse IV, Jun 13 2017
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Robert G. Wilson v, Aug 22 2000
EXTENSIONS
a(6)-a(8) from Ryan Propper, Mar 25 2007
a(9) from Michael S. Branicky, Aug 11 2024
STATUS
approved