login
A374812
10^n-th nonsquarefree number.
1
4, 27, 260, 2550, 25506, 255051, 2550519, 25505464, 255054546, 2550546152, 25505460948, 255054609189, 2550546096325, 25505460966984, 255054609673752, 2550546096737236, 25505460967313208, 255054609673067871, 2550546096730376472, 25505460967304302850, 255054609673042999696, 2550546096730430494608
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