OFFSET
1,2
COMMENTS
For k > 1, numbers k such that the digits 66 followed by k-1 occurrences of the digit 3 followed by the digit 1 is prime (see Example section).
a(28) > 3*10^5.
LINKS
Makoto Kamada, Factorization of near-repdigit-related numbers.
Makoto Kamada, Search for 663w1.
EXAMPLE
4 is in this sequence because (199*10^4 - 7)/3 = 663331 is prime.
Initial terms and associated primes:
a(1) = 1, 661;
a(2) = 4, 663331;
a(3) = 16, 663333333333333331;
a(4) = 25, 663333333333333333333333331;
a(5) = 27, 66333333333333333333333333331; etc.
MATHEMATICA
Select[Range[0, 100000], PrimeQ[(199*10^# - 7)/3] &]
PROG
(PARI) isok(k) = isprime((199*10^k - 7)/3); \\ Michel Marcus, Feb 04 2017
CROSSREFS
KEYWORD
nonn,more,hard
AUTHOR
Robert Price, Feb 04 2017
EXTENSIONS
a(26)-a(27) from Robert Price, May 21 2020
STATUS
approved