OFFSET
1,2
FORMULA
G.f. A(x) satisfies: A(x) = x * exp(2 * Sum_{i>=1} Sum_{j>=1} A(x^(i*(2*j-1)))/(2*j - 1)).
MATHEMATICA
a[n_] := a[n] = SeriesCoefficient[x Product[Product[((1 + x^(j k))/(1 - x^(j k)))^a[j], {k, 1, n - 1}], {j, 1, n - 1}], {x, 0, n}]; Table[a[n], {n, 1, 25}]
terms = 25; A[_] = 0; Do[A[x_] = x Exp[2 Sum[Sum[A[x^(i (2 j - 1))]/(2 j - 1), {j, 1, terms}], {i, 1, terms}]] + O[x]^(terms + 1) // Normal, terms + 1]; Rest[CoefficientList[A[x], x]]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, May 14 2019
STATUS
approved