%I #8 Oct 16 2023 09:30:14
%S 1,1,1,1,1,1,2,1,1,1,2,1,2,2,1,1,2,1,2,1,2,2,2,1,1,2,1,2,2,1,2,1,2,2,
%T 2,1,2,2,2,1,2,2,2,2,1,2,2,1,3,1,2,2,2,1,2,2,2,2,2,1,2,2,2,1,2,2,2,2,
%U 2,2,2,1,2,2,1,2,4,2,2,1,1,2,2,2,2,2,2
%N a(n) is the number of divisors of n that are coprime to the terms of A366642.
%H Amiram Eldar, <a href="/A366643/b366643.txt">Table of n, a(n) for n = 1..10000</a>
%F Multiplicative with a(p^e) = 1 if p is a term of A366642, and e+1 otherwise.
%F Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k)/A000005(k) = 1/2.
%F Dirichlet g.f.: zeta(s)^2 * Product_{p in A366642} (1 - 1/p^s).
%F Sum_{k=1..n} a(k) ~ c * n * log(n), where c = Product_{p in A366642} (1 - 1/p) = 0.26485234983834588444... .
%t seq[max_] := With[{ps = {2, 3, 5, 149, 10771}}, If[max <= Max[ps], f[p_, e_] := If[MemberQ[ps, p], 1, e + 1]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, max], Print["Add to ps more terms from A366642"]]]; seq[10^6]
%Y Cf. A000005, A001227, A035218, A366642.
%K nonn,easy,mult
%O 1,7
%A _Amiram Eldar_, Oct 15 2023