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