A357819 Denominators of the partial sums of the reciprocals of the Dedekind psi function (A001615).
1, 3, 12, 4, 12, 1, 8, 24, 24, 72, 72, 36, 252, 504, 126, 504, 504, 504, 2520, 840, 3360, 10080, 10080, 10080, 2016, 2016, 2016, 2016, 10080, 10080, 5040, 1260, 5040, 15120, 7560, 3780, 71820, 17955, 143640, 17955, 35910, 574560, 6320160, 6320160, 6320160, 6320160
Offset: 1
Programs
-
Mathematica
psi[n_] := n * Times @@ (1 + 1/Transpose[FactorInteger[n]][[1]]); psi[1] = 1; Denominator[Accumulate[1/Array[psi[#] &, 50]]]
-
PARI
f(n) = n * sumdivmult(n, d, issquarefree(d)/d); \\ A001615 a(n) = denominator(sum(k=1, n, 1/f(k))); \\ Michel Marcus, Oct 15 2022
Formula
a(n) = denominator(Sum_{k=1..n} 1/psi(k)).
Comments