login
A153331
Numbers k such that (10^k - 1)*140/99 + 1 is prime.
0
10, 276, 478, 582, 1630, 6342, 14688
OFFSET
1,1
COMMENTS
These numbers are always even. If k is odd, then 10^k - 1 produces a number with an odd number of 9's which 99 does not divide. Also these numbers are palindromic.
a(8) > 50000. - Michael S. Branicky, Jan 15 2026
EXAMPLE
10 is a term, because (10^10 - 1)*140/99 + 1 = 14141414141 is prime.
PROG
(PARI) /* n=number of values to test, r=repeat digits, e.g., 14, 121, 177, 1234, etc. d = last digit appended to the end */
repr(n, r, d) = ln=length(Str(r)); for(x=0, n, y=(10^(ln*x)-1)*10*r/ (10^ln-1)+1; if(ispseudoprime(y), print1(ln*x", ")))
repr(10^3, 14, 1)
(Python)
from sympy import isprime
def afind(limit, startk=2):
k = startk + (startk%2)
t = int("1" + "41"*(k//2))
for k in range(startk, limit+1, 2):
if isprime(t): print(k, end=", ")
t *= 100
t += 41
afind(600) # Michael S. Branicky, Dec 13 2021
CROSSREFS
Sequence in context: A222998 A368770 A221526 * A290202 A203149 A117654
KEYWORD
nonn,base,hard,more
AUTHOR
Cino Hilliard, Dec 23 2008
EXTENSIONS
a(6) from Michael S. Branicky, Dec 13 2021
a(7) from Michael S. Branicky, May 14 2023
STATUS
approved