%I #18 Aug 09 2024 11:16:01

%S 0,0,0,0,0,1,0,0,0,1,0,1,0,1,1,0,0,1,0,1,1,1,0,1,0,1,0,1,0,6,0,0,1,2,

%T 1,1,0,2,1,1,0,6,0,1,1,1,0,1,0,1,1,1,0,1,1,1,1,1,0,7,0,1,1,0,1,7,0,2,

%U 1,7,0,1,0,2,0,2,1,7,0,1,0,1,0,8,1,1,2

%N a(n) = A010846(n) - A373738(n).

%H Michael De Vlieger, <a href="/A372192/b372192.txt">Table of n, a(n) for n = 1..10000</a>

%H Michael De Vlieger, <a href="/A372192/a372192_1.png">Log log scatterplot of a(n)</a>, n = 6..2^20, ignoring a(n) = 0, with a color code where red indicates omega(n) = 2, gold omega(n) = 3, lime green = omega(n) = 4, spring green = omega(n) = 5, and blue = omega(n) = 6, where omega = A001221.

%t {0}~Join~Table[Count[Range[n], _?(Divisible[n, Times @@ FactorInteger[#][[All, 1]]] &)] - Floor[(1/PrimeNu[n]!)*Times @@ Map[Log[#, n] + 1 &, FactorInteger[n][[All, 1]]]], {n, 2, 120}]

%Y Cf. A010846, A373738.

%K nonn

%O 1,30

%A _Michael De Vlieger_, Aug 06 2024