OFFSET
1,2
FORMULA
MATHEMATICA
nn = 32; Denominator[Table[Sum[Sum[If[Mod[n, k] == 0, MoebiusMu[n/k]*HarmonicNumber[k], 0], {k, 1, n}], {n, 1, m}], {m, 1, nn}]]
PROG
(PARI) a(n) = denominator(sum(m=1, n, sumdiv(m, d, moebius(m/d)*sum(i=1, d, 1/i)))); \\ Michel Marcus, Apr 23 2020
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Mats Granvik, Apr 23 2020
STATUS
approved