A082953 a(n) = A000252(n) / A070732(n).
1, 2, 4, 8, 16, 8, 36, 32, 36, 32, 100, 32, 144, 72, 64, 128, 256, 72, 324, 128, 144, 200, 484, 128, 400, 288, 324, 288, 784, 128, 900, 512, 400, 512, 576, 288, 1296, 648, 576, 512, 1600, 288, 1764, 800, 576, 968, 2116
Offset: 1
Links
- Michael De Vlieger, Table of n, a(n) for n = 1..10000
Crossrefs
Programs
-
Maple
A082953 := proc(n) numtheory[phi](n)*numtheory[phi](2*n) ; end proc: seq(A082953(n),n=1..100) ; # R. J. Mathar, Jan 07 2011
-
Mathematica
Array[Times @@ Map[EulerPhi, {#, 2 #}] &, 47] (* Michael De Vlieger, Apr 21 2019 *)
-
PARI
a(n) = eulerphi(n)*eulerphi(2*n); \\ Michel Marcus, Jun 04 2025
Formula
Multiplicative with a(2^e) = 2^(2e-1) and a(p^e) = (p-1)^2*p^(2e-2) for p > 2. - R. J. Mathar, Apr 14 2011
a(n) = phi(n)^2 if n odd; 2*phi(n)^2 if n even, where phi(n) = A000010(n). - Jianing Song, Apr 20 2019
Sum_{k=1..n} a(k) ~ c * n^3, where c = (2/5) * Product_{p prime} (1 - (2*p-1)/p^3) = (2/5) * A065464 = 0.171299... . - Amiram Eldar, Oct 30 2022
Comments