|
|
A335416
|
|
a(n) = numerator of Sum_{k=2..A335138(n)} abs(A309229(n, k))/k.
|
|
1
|
|
|
5, 16, 27, 17, 7, 37, 679, 1343, 911, 3163, 741, 10652, 793, 3310, 7027, 8263, 8381, 1777, 114419, 27952, 18301, 47987, 1480987, 27593, 1346549, 12821, 159847, 92503, 1612207, 861440, 7920089, 3251923, 544897, 15681619, 3047987, 18092504, 1643309, 193472
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
MATHEMATICA
|
nn = 38; 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}]; Numerator[d]
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,frac
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|