|
MATHEMATICA
|
A222626 = {1, 2, 3, 7, 14, 46, 164, 996, 10195, 208168, 8016530, 542165050};
mob[m_, n_] := If[Mod[m, n] == 0, MoebiusMu[m/n], 0];
EULERi[b_] := Module[{a, c, i, d}, c = {}; For[i = 1, i <= Length[b], i++, c = Append[c, i b[[i]] - Sum[c[[d]] b[[i - d]], {d, 1, i - 1}]]]; a = {}; For[i = 1, i <= Length[b], i++, a = Append[a, (1/i)*Sum[mob[i, d] c[[d]], {d, 1, i}]]]; Return[a]];
|