OFFSET
1,2
COMMENTS
See A357843 for more details.
LINKS
László Tóth, Alternating Sums Concerning Multiplicative Arithmetic Functions, Journal of Integer Sequences, Vol. 20 (2017), Article 17.2.1.
FORMULA
a(n) = denominator(Sum_{k=1..n} (-1)^(k+1)/d(k)), where d(k) = A000005(k).
MATHEMATICA
Denominator[Accumulate[Array[(-1)^(# + 1)/DivisorSigma[0, #] &, 60]]]
PROG
(PARI) lista(nmax) = {my(s = 0); for(k = 1, nmax, s += (-1)^(k+1) / numdiv(k); print1(denominator(s), ", "))};
(Python)
from fractions import Fraction
from sympy import divisor_count
def A357844(n): return sum(Fraction(1 if k&1 else -1, divisor_count(k)) for k in range(1, n+1)).denominator # Chai Wah Wu, Oct 16 2022
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Amiram Eldar, Oct 16 2022
STATUS
approved