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