OFFSET
0,2
FORMULA
G.f.: B(x) + B(x)^2/2 + B(x^2)/2 - B(x)*(B(x)-x), where B(x) is g.f. for A029856.
MATHEMATICA
max = 30; B[_] = 1; Do[B[x_] = x*Exp[Sum[(B[x^k] + x^k)/k + O[x]^n, {k, 1, n}]] // Normal, {n, 1, max}]; A[x_] = 1 + x + x^2 + B[x] - B[x]^2/2 + B[x^2]/2; CoefficientList[A[x] + O[x]^max, x] (* Jean-François Alcover, Jan 28 2019 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Christian G. Bower, Nov 15 1998
STATUS
approved