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.10, pp. 30-31.
FORMULA
a(n) = denominator(Sum_{k=1..n} (-1)^(k+1)/A047994(k)).
MATHEMATICA
uphi[n_] := Times @@ (-1 + Power @@@ FactorInteger[n]); uphi[1] = 1; Denominator[Accumulate[Table[(-1)^(n+1)/uphi[n], {n, 1, 50}]]]
PROG
(PARI) uphi(n) = {my(f = factor(n)); prod(i = 1, #f~, -1 + f[i, 1]^f[i, 2]); }
list(nmax) = {my(s = 0); for(k = 1, nmax, s += (-1)^(k+1) / uphi(k); print1(denominator(s), ", "))};
CROSSREFS
KEYWORD
nonn,easy,frac,new
AUTHOR
Amiram Eldar, Dec 24 2024
STATUS
approved