OFFSET
1,1
COMMENTS
Equivalently, numbers k such that 125*10^(k-3) - 1 is prime.
Searched up to k = 20000. - Michael B. Porter, Dec 07 2017
Searched up to k = 40000. - Michael S. Branicky, Aug 07 2024
EXAMPLE
1249, 1249999, and 12499999 are prime, so 4, 7, and 8 are in the sequence.
MATHEMATICA
Select[Range[2000], PrimeQ[10^#/8 - 1] &] (* Michael De Vlieger, Dec 08 2017 *)
PROG
(PARI) lista(nn) = for(n=4, nn, if(ispseudoprime(10^n/8 - 1), print1(n, ", "))) \\ Iain Fox, Dec 04 2017
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Patrick A. Thomas, Dec 03 2017
EXTENSIONS
a(12), a(13), a(15), a(16) from Jon E. Schoenfield, Dec 03 2017
a(14) and a(17)-a(21) from Michael B. Porter, Dec 07 2017
a(22) from Michael S. Branicky, Aug 07 2024
STATUS
approved