OFFSET
1,2
COMMENTS
Each of the corresponding numbers of the form 16^n * (16^n - 1) - 1 is almost a repdigit in base 16 (i.e., a repdigit with only one change).
n is in the sequence iff 4*n is in A098845. - Robert Israel, Dec 14 2015
LINKS
Brady Haran and Simon Pampena, Glitch Primes and Cyclops Numbers, Numberphile video (2015)
MATHEMATICA
Select[Range[1000], PrimeQ[16^# (16^# - 1) - 1] &] (* _Shivam Patel_, Dec 13 2015 *)
PROG
(PARI) is(n)=ispseudoprime(16^n * (16^n - 1) - 1) \\ Anders Hellström, Dec 13 2015
(Magma) [n: n in [1..200] | IsPrime(256^n-16^n-1)]; // Vincenzo Librandi, Dec 14 2015
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Shivam N. Patel, Dec 13 2015
EXTENSIONS
a(12)-a(16) computed from A098845 by Ray Chandler, Sep 25 2019
STATUS
approved