OFFSET
1,2
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.6, pp. 24-26.
FORMULA
a(n) = denominator(Sum_{k=1..n} (-1)^(k+1)/A007947(k)).
MATHEMATICA
rad[n_] := Times @@ FactorInteger[n][[;; , 1]]; Denominator[Accumulate[Table[(-1)^(n+1)/rad[n], {n, 1, 50}]]]
PROG
(PARI) rad(n) = vecprod(factor(n)[, 1]);
list(nmax) = {my(s = 0); for(k = 1, nmax, s += (-1)^(k+1) / rad(k); print1(denominator(s), ", "))};
CROSSREFS
KEYWORD
nonn,easy,frac
AUTHOR
Amiram Eldar, Dec 21 2024
STATUS
approved