OFFSET
1,3
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..10000
MATHEMATICA
Table[GCD[DivisorSigma[0, n], DivisorSigma[1, n]], {n, 120}] (* Harvey P. Dale, Dec 05 2017 *)
PROG
(PARI) A009205(n) = gcd(numdiv(n), sigma(n)); \\ Antti Karttunen, May 22 2017
(Python)
from math import prod, gcd
from sympy import factorint
def A009205(n):
f = factorint(n).items()
return gcd(prod(e+1 for p, e in f), prod((p**(e+1)-1)//(p-1) for p, e in f)) # Chai Wah Wu, Jul 27 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
Data section extended to 120 terms by Antti Karttunen, May 22 2017
STATUS
approved