OFFSET
1,1
LINKS
Michael S. Branicky, Table of n, a(n) for n = 1..344 (using factordb.com)
factordb.com, Status of (10^(2*n+2)-1)/99.
FORMULA
a(n) = 3 if n mod 3 == 2.
a(n) = A020639((10^(2*n+2)-1)//99).
EXAMPLE
a(1) = 101 since 101 is prime.
a(2) = 3 since 10101 = 3*7*13*37.
a(3) = 73 since 1010101 = 73*101*137.
a(4) = 41 since 101010101 = 41*271*9091.
PROG
(Python)
from sympy import factorint
def a(n): return min(factorint((10**(2*n+2)-1)//99))
print([a(n) for n in range(1, 46)]) # Michael S. Branicky, Apr 18 2026
(Julia)
function a(n)
num = (big(10)^(2n + 2) - 1) รท 99
for i in 3:2:num
if num % i == 0
return i
end
end
end # Hoang Nguyen, Apr 21 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Hoang Nguyen, Apr 16 2026
EXTENSIONS
a(9) onward from Michael S. Branicky, Apr 18 2026
STATUS
approved
