OFFSET
1,3
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..1000
László Tóth, Alternating Sums Concerning Multiplicative Arithmetic Functions, Journal of Integer Sequences, Vol. 20 (2017), Article 17.2.1. See section 4.14, p. 35.
FORMULA
a(n) = denominator(Sum_{k=1..n} (-1)^(k+1)/A206369(k)).
MATHEMATICA
f[p_, e_] := Sum[(-1)^(e-k)*p^k, {k, 0, e}]; beta[1] = 1; beta[n_] := Times @@ f @@@ FactorInteger[n]; Denominator[Accumulate[Table[(-1)^(n+1)/beta[n], {n, 1, 50}]]]
PROG
(PARI) beta(n) = {my(f = factor(n)); prod(i=1, #f~, p = f[i, 1]; e = f[i, 2]; sum(k = 0, e, (-1)^(e-k)*p^k)); }
list(nmax) = {my(s = 0); for(k = 1, nmax, s += (-1)^(k+1) / beta(k); print1(denominator(s), ", "))};
CROSSREFS
KEYWORD
nonn,easy,frac
AUTHOR
Amiram Eldar, Dec 27 2024
STATUS
approved
