OFFSET
1,1
MATHEMATICA
nn = 41; f[n_] := n; h[n_] := DivisorSum[n, MoebiusMu[#] # &]; A = Accumulate[Table[Table[h[GCD[n, k]], {k, 1, nn}], {n, 1, nn}]]; B = -Abs[A]; B[[All, 1]] = Table[f[n], {n, 1, nn}]; b = 1 + Total[Sign[1 + Sign[Accumulate[Transpose[B]]]]]; d = Table[Sum[Abs[A[[n, k]]]/k, {k, 2, b[[n]]}], {n, 1, nn}]; Denominator[d]
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Mats Granvik, Jun 09 2020
STATUS
approved