%I #13 May 10 2024 09:25:53
%S 1,1,3,2,5,3,6,4,9,5,10,6,12,6,15,8,16,9,18,10,18,10,22,12,25,12,27,
%T 12,28,15,30,16,30,16,30,18,36,18,36,20,40,18,42,20,45,22,46,24,42,25,
%U 48,24,52,27,50,24,54,28,58,30,60,30,54,32,60,30,66,32,66,30,70,36,72,36,75,36,60,36
%N a(n) = phi(15 * n)/8.
%F G.f.: Sum_{k>=1} mu(15 * k) * x^k / (1 - x^k)^2, where mu() is the Moebius function (A008683).
%F Multiplicative with a(p^e) = p^e if p = 3 or 5, and (p-1)*p^(e-1) otherwise.
%F Sum_{k=1..n} a(k) ~ (225/(64*Pi^2)) * n^2. - _Amiram Eldar_, May 10 2024
%t a[n_] := EulerPhi[15 * n]/8; Array[a, 100] (* _Amiram Eldar_, May 10 2024 *)
%o (PARI) a(n) = eulerphi(15*n)/8;
%o (PARI) my(N=80, x='x+O('x^N)); Vec(sum(k=1, N, moebius(15*k)*x^k/(1-x^k)^2))
%Y Column k=15 of A372673.
%Y Cf. A008683.
%K nonn,mult,easy
%O 1,3
%A _Seiichi Manyama_, May 10 2024