OFFSET
1,2
COMMENTS
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
MATHEMATICA
Table[GCD[DivisorSigma[0, n^3], DivisorSigma[0, n]], {n, 1, 500}]
PROG
(PARI) a(n)=my(f=factor(n)[, 2]); gcd(prod(i=1, #f, 3*f[i]+1), prod(i=1, #f, f[i]+1)) \\ Charles R Greathouse IV, Oct 16 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Apr 08 2002
STATUS
approved