OFFSET
1,1
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..20000
MAPLE
a:= n-> ilcm(seq(d+n/d, d=numtheory[divisors](n))):
seq(a(n), n=1..60); # Alois P. Heinz, Mar 09 2016
MATHEMATICA
eq[n_] := And @@ Table[y = n/x; Mod[a[x*y], x + y] == 0, {x, Divisors[n]}]; r[n_] := Reduce[ eq[n], a[n], Integers] /. C[1] -> 1; Table[a[n] /. ToRules[r[n]], {n, 1, 51}] (* Jean-François Alcover, Aug 03 2012 *)
CROSSREFS
KEYWORD
nonn,nice
AUTHOR
Andrew Weimholt, Jul 05 2010
STATUS
approved