OFFSET
1,2
LINKS
MATHEMATICA
Select[Range[2000], PrimeQ[# 16^# + 1] &]
PROG
(Magma) [n: n in [0..2000] | IsPrime(n*16^n+1)];
(Sage) [n for n in (1..2000) if is_prime(n*16^n + 1)] # Bruno Berselli, May 09 2014
(PARI) is(n)=ispseudoprime(n*16^n+1) \\ Charles R Greathouse IV, Jun 06 2017
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Vincenzo Librandi, May 08 2014
EXTENSIONS
a(7)-a(8) from Loeh's list (see Links)
STATUS
approved