OFFSET
1,1
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
FORMULA
MAPLE
a := n -> nops(select(issqrfree, [$n..(2*n)])):
seq(a(n), n=1..75); # Peter Luschny, Mar 02 2017
MATHEMATICA
a[n_] := Select[Range[n, 2n], SquareFreeQ] // Length;
Array[a, 75] (* Jean-François Alcover, Jun 22 2018 *)
PROG
(PARI) f(n)=my(s); forfactored(k=1, sqrtint(n), s += n\k[1]^2*moebius(k)); s
a(n)=f(2*n)-f(n-1) \\ Charles R Greathouse IV, Nov 05 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Jul 05 2010
STATUS
approved