OFFSET
1,3
MATHEMATICA
nn = 23; f[x_] := Product[(1 + x^n/n!)^a[n], {n, 1, nn}]; sol = SolveAlways[0 == Series[f[x] - Exp[-x]/(1 - x), {x, 0, nn}], x]; Table[a[n], {n, 1, nn}] /. sol // Flatten
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, May 08 2022
STATUS
approved