OFFSET
0,1
COMMENTS
a(n)/10^n converges to Pi^2/(Pi^2-6) = zeta(2)/(zeta(2)-1) as n -> oo. - Chai Wah Wu, Jul 21 2024
LINKS
Chai Wah Wu, Algorithms for complementary sequences, arXiv:2409.05844 [math.NT], 2024.
FORMULA
a(n) = A013929(10^n).
PROG
(Python)
from functools import lru_cache
from math import isqrt
from sympy import mobius
def A374812(n):
r = 10**n
@lru_cache(maxsize=None)
def g(x): return mobius(x)
def f(x): return r+sum(g(k)*(x//k**2) for k in range(1, isqrt(x)+1))
m, k = r, f(r)
while m != k:
m, k = k, f(k)
return m
CROSSREFS
KEYWORD
nonn
AUTHOR
Chai Wah Wu, Jul 20 2024
EXTENSIONS
a(15) from Chai Wah Wu, Jul 21 2024
a(16)-a(21) from Chai Wah Wu, Aug 08 2024
STATUS
approved