OFFSET
1,1
COMMENTS
For k > 1, numbers k such that the digit 1 followed by k-2 occurrences of the digit 6 followed by the digits 33 is prime (see Example section).
a(25) > 3*10^5.
LINKS
Makoto Kamada, Factorization of near-repdigit-related numbers.
Makoto Kamada, Search for 16w33.
EXAMPLE
4 is in this sequence because (5*10^4 - 101)/3 = 16633 is prime.
Initial terms and associated primes:
a(1) = 4, 16633;
a(2) = 19, 16666666666666666633;
a(3) = 29, 166666666666666666666666666633;
a(4) = 42, 1666666666666666666666666666666666666666633;
a(5) = 48, 1666666666666666666666666666666666666666666666633; etc.
MATHEMATICA
Select[Range[2, 100000], PrimeQ[(5*10^# - 101)/3] &]
PROG
(PARI) isok(k) = ispseudoprime((5*10^k-101)/3); \\ Altug Alkan, Apr 19 2018
CROSSREFS
KEYWORD
nonn,more,hard
AUTHOR
Robert Price, Feb 16 2017
EXTENSIONS
a(22) from Robert Price, Apr 19 2018
a(23)-a(24) from Robert Price, Oct 25 2023
STATUS
approved