A062401 a(n) = phi(sigma(n)).
1, 2, 2, 6, 2, 4, 4, 8, 12, 6, 4, 12, 6, 8, 8, 30, 6, 24, 8, 12, 16, 12, 8, 16, 30, 12, 16, 24, 8, 24, 16, 36, 16, 18, 16, 72, 18, 16, 24, 24, 12, 32, 20, 24, 24, 24, 16, 60, 36, 60, 24, 42, 18, 32, 24, 32, 32, 24, 16, 48, 30, 32, 48, 126, 24, 48, 32, 36, 32, 48, 24, 96, 36, 36, 60
Offset: 1
Keywords
Examples
a(9) = 12 because sigma(9) = 13 and phi(13) = 12.
References
- D. S. Mitrinovic et al., Handbook of Number Theory, Kluwer, p. 14.
Links
- T. D. Noe, Table of n, a(n) for n = 1..10000
Crossrefs
Programs
-
Haskell
a062401 = a000010 . a000203 -- Reinhard Zumkeller, Jan 04 2013
-
Maple
with(numtheory); A062401:=n->phi(sigma(n)); seq(A062401(n), n=1..50); # Wesley Ivan Hurt, Apr 07 2014
-
Mathematica
Table[EulerPhi[DivisorSigma[1, n]], {n, 1, 80}] (* Carl Najafi, Aug 16 2011 *)
-
PARI
vector(150, n, eulerphi(sigma(n)))
-
PARI
for (n=1, 10000, write("b062401.txt", n, " ", eulerphi(sigma(n))) ) \\ Harry J. Smith, Aug 07 2009