OFFSET
1,1
COMMENTS
See A237607 for the numbers themselves and for more comments.
The ratio r(n) = a(n+1)/a(n) seems to converge to a constant value of about 4.245...
LINKS
Stanislav Sykora, PARI/GP scripts for genetic threads
EXAMPLE
One-digit numbers k with moebius(k)=0 are 4, 8 and 9. Hence a(1) = 3.
PROG
(PARI) \\ See the link (use GT_DivMod0 with "IsMoebiusZero" for the "property" argument).
(Python)
from sympy.ntheory.factor_ import core
def P(n): return core(n) < n
reach = {4, 8, 9}
for d in range(1, 10001):
print(len(reach), end=", ")
newreach = set()
for q in reach:
for digit in range(10):
t = 10*q + digit
if P(t):
newreach.add(t)
reach = newreach
# Michael S. Branicky, Jun 03 2026
CROSSREFS
KEYWORD
nonn,hard,more,base
AUTHOR
Stanislav Sykora, Feb 15 2014
EXTENSIONS
a(12)-a(13) from Lars Blomberg, Jan 02 2016
STATUS
approved
