A335417 a(n) = denominator of Sum_{k=2..A335138(n)} abs(A309229(n, k))/k.
6, 15, 20, 15, 4, 35, 360, 630, 420, 1386, 280, 5005, 280, 1287, 3080, 3465, 2520, 770, 27720, 9009, 5720, 16380, 360360, 10010, 360360, 3465, 40040, 25740, 360360, 323323, 1750320, 765765, 123760, 3423420, 720720, 4849845, 360360, 45045
Offset: 1
Programs
-
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]