A009191 a(n) = gcd(n, d(n)), where d(n) is the number of divisors of n (A000005).
1, 2, 1, 1, 1, 2, 1, 4, 3, 2, 1, 6, 1, 2, 1, 1, 1, 6, 1, 2, 1, 2, 1, 8, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 9, 1, 2, 1, 8, 1, 2, 1, 2, 3, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 8, 1, 2, 1, 12, 1, 2, 3, 1, 1, 2, 1, 2, 1, 2, 1, 12, 1, 2, 3, 2, 1, 2, 1, 10, 1, 2, 1, 12, 1, 2, 1, 8, 1, 6, 1, 2, 1, 2, 1, 12, 1, 2, 3, 1, 1, 2, 1, 8, 1
Offset: 1
Keywords
Links
- Antti Karttunen, Table of n, a(n) for n = 1..65537 (terms 1..1000 from T. D. Noe)
Crossrefs
Programs
-
Haskell
a009191 n = gcd n $ a000005 n -- Reinhard Zumkeller, May 09 2013, Aug 14 2011
-
Mathematica
f[n_] := GCD[n, DivisorSigma[0, n]]; Array[f, 105] (* Robert G. Wilson v, Mar 27 2013 *)
-
PARI
a(n)=gcd(numdiv(n),n) \\ Charles R Greathouse IV, Mar 26 2013
Comments