OFFSET
1,2
COMMENTS
See A357820 for more details.
FORMULA
a(n) = denominator(Sum_{k=1..n} (-1)^(k+1)/psi(k)).
MATHEMATICA
psi[n_] := n * Times @@ (1 + 1/Transpose[FactorInteger[n]][[1]]); psi[1] = 1; Denominator[Accumulate[1/Array[(-1)^(# + 1)*psi[#] &, 50]]]
PROG
(PARI) f(n) = n * sumdivmult(n, d, issquarefree(d)/d); \\ A001615
a(n) = denominator(sum(k=1, n, (-1)^(k+1)/f(k))); \\ Michel Marcus, Oct 15 2022
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Amiram Eldar, Oct 14 2022
STATUS
approved