A378086 Number of nonsquarefree numbers < prime(n).
0, 0, 1, 1, 3, 4, 5, 6, 7, 11, 11, 13, 14, 14, 16, 20, 22, 23, 25, 26, 27, 29, 31, 33, 36, 39, 39, 40, 41, 42, 49, 50, 53, 53, 57, 58, 61, 63, 64, 68, 70, 71, 74, 75, 76, 77, 81, 84, 86, 87, 88, 90, 91, 97, 99, 101, 103, 104, 107, 109, 109, 113, 119, 120, 121
Offset: 1
Keywords
Examples
The nonsquarefree numbers counted under each term begin: n=1: n=2: n=3: n=4: n=5: n=6: n=7: n=8: n=9: n=10: n=11: n=12: -------------------------------------------------------------- . . 4 4 9 12 16 18 20 28 28 36 8 9 12 16 18 27 27 32 4 8 9 12 16 25 25 28 4 8 9 12 24 24 27 4 8 9 20 20 25 4 8 18 18 24 4 16 16 20 12 12 18 9 9 16 8 8 12 4 4 9 8 4
Crossrefs
Programs
-
Mathematica
Table[Length[Select[Range[Prime[n]],!SquareFreeQ[#]&]],{n,100}]
-
Python
from math import isqrt from sympy import prime, mobius def A378086(n): return (p:=prime(n))-sum(mobius(k)*(p//k**2) for k in range(1,isqrt(p)+1)) # Chai Wah Wu, Dec 05 2024
Formula
a(n) = A057627(prime(n)).