OFFSET
1,2
COMMENTS
The primes arising from this construction (e.g., 999998999999) are among the primes counted in A266148. In particular, it follows that A266148(a(n)) > 0. - David A. Corneth, May 19 2016
a(16) > 188484. - Ben Meekins, Sep 08 2018
LINKS
Brady Haran and Simon Pampena, Glitch Primes and Cyclops Numbers, Numberphile video (2015).
M. Kamada, Near-repdigit-related prime numbers
EXAMPLE
6 is in the sequence because 10^12 - 10^6 - 1 = 999998999999 is prime.
MATHEMATICA
Select[Range[15000], PrimeQ[10^# (10^# - 1) - 1] &] (* Vincenzo Librandi, Dec 08 2015 *)
PROG
(PARI) for(n=1, 9999, if(ispseudoprime(10^n*(10^n-1)-1), print1(n", ")))
(Magma) [n: n in [0..200] | IsPrime(10^n*(10^n-1)-1)]; // Vincenzo Librandi, Dec 08 2015
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Serge Batalov, Dec 07 2015
EXTENSIONS
a(11) from Kazuyoshi Asao, Feb 11 2002
a(12) from Serge Batalov, Dec 25 2015
a(13) from Ben Meekins, Feb 16 2016
a(14) from Ben Meekins, Dec 17 2016
a(15) from Ben Meekins, Sep 08 2018
STATUS
approved