A295310 a(n) = gcd(n, A062401(n)), where A062401(n) = phi(sigma(n)).
1, 2, 1, 2, 1, 2, 1, 8, 3, 2, 1, 12, 1, 2, 1, 2, 1, 6, 1, 4, 1, 2, 1, 8, 5, 2, 1, 4, 1, 6, 1, 4, 1, 2, 1, 36, 1, 2, 3, 8, 1, 2, 1, 4, 3, 2, 1, 12, 1, 10, 3, 2, 1, 2, 1, 8, 1, 2, 1, 12, 1, 2, 3, 2, 1, 6, 1, 4, 1, 2, 1, 24, 1, 2, 15, 4, 1, 6, 1, 20, 1, 2, 1, 12, 1, 2, 1, 8, 1, 18, 1, 4, 1, 2, 1, 24, 1, 2, 3, 20, 1, 6, 1, 8, 1
Offset: 1
Keywords
Links
- Antti Karttunen, Table of n, a(n) for n = 1..16384
Programs
-
Mathematica
Array[GCD[#, EulerPhi[DivisorSigma[1, #]]] &, 105] (* Michael De Vlieger, Nov 23 2017 *)
-
PARI
a(n) = gcd(n, eulerphi(sigma(n))); \\ Michel Marcus, Nov 23 2017