OFFSET
1,1
COMMENTS
a(n) is a divisor of A020696(n). - Ivan Neretin, May 27 2015
LINKS
Carl R. White, Table of n, a(n) for n = 1..10000
EXAMPLE
Since the positive divisors of 6 are 1, 2, 3 and 6, a(6) = LCM(1+1,2+1,3+1,6+1) = LCM(2,3,4,7) = 84.
MAPLE
f:= n -> ilcm(op(map(`+`, numtheory:-divisors(n), 1)));
seq(f(n), n=1..100); # Robert Israel, Jul 24 2014
MATHEMATICA
a057643[n_Integer] := Apply[LCM, Map[# + 1 &, Divisors[n]]]; Table[a057643[n], {n, 10000}] (* Michael De Vlieger, Jul 19 2014 *)
PROG
(PARI) a(n)=lcm(apply(d->d+1, divisors(n))) \\ Charles R Greathouse IV, Feb 14 2013
(Python)
from math import lcm
from sympy import divisors
def A057643(n): return lcm(*(d+1 for d in divisors(n, generator=True))) # Chai Wah Wu, Jun 30 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, Oct 11 2000
STATUS
approved