OFFSET
0,2
COMMENTS
What a(n)/n is converging to?
Limit_{n->oo} a(n)/10^n = 1/zeta(2). [Max Alekseyev, Oct 18 2008]
LINKS
Lucas A. Brown, Python program.
FORMULA
a(n) = Sum_{j=0..-1+10^n} abs(mu(10^n + j)).
EXAMPLE
n=10: among numbers {10,...,19} seven are squarefree [10,11,13,14,15,17,19], so a(1)=7.
MAPLE
with(numtheory): for n from 0 to 5 do ct:=0: for k from 10^n to 2*10^n-1 do if abs(mobius(k))>0 then ct:=ct+1 else ct:=ct: fi: od: a[n]:=ct: od: seq(a[n], n=0..5); # Emeric Deutsch, Mar 28 2005
MATHEMATICA
Table[Apply[Plus, Table[Abs[MoebiusMu[10^w+j]], {j, 0, -1+10^(w-1)}]], {w, 0, 6}]
PROG
(PARI) { a(n) = sum(m=1, sqrtint(2*10^n-1), moebius(m) * ((2*10^n-1)\m^2 - (10^n-1)\m^2) ) } \\ Max Alekseyev, Oct 18 2008
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Nov 14 2002
EXTENSIONS
6079284 from Emeric Deutsch, Mar 28 2005
a(8)-a(15) from Max Alekseyev, Oct 18 2008
a(16)-a(21) from Lucas A. Brown, Feb 25 2024
STATUS
approved