|
MATHEMATICA
|
terms = 17;
s = (1 - 1/(1 + Sum[(n+1)^(n-1)*t^n, {n, 1, terms}]))/t + O[t]^(terms-1);
A141313 = 2^Range[terms-1]*CoefficientList[s, t];
did[m_, n_] := If[Mod[m, n] == 0, 1, 0];
EulerTransform[seq_] := Module[{coeff, final = {}}, coeff = Table[Sum[d* did[i, d]*seq[[d]], {d, 1, i}], {i, 1, Length[seq]}]; For[i = 1, i <= Length[seq], i++, AppendTo[final, (coeff[[i]] + Sum[coeff[[d]]*final[[i - d]], {d, 1, i-1}])/i]]; final];
|