OFFSET
1,2
FORMULA
MATHEMATICA
f1[p_, e_] := e/p; f2[p_, e_] := (p + 1)^e; a[1] = 1; a[n_] := GCD[n*(1 + Plus @@ f1 @@@ (f = FactorInteger[n])), Times @@ f2 @@@ f]; Array[a, 100] (* Amiram Eldar, Nov 06 2021 *)
PROG
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Nov 06 2021
STATUS
approved