A241001 G.f. satisfies: A(x)^3 = x + A(x*A(x)^3)^2.

1, 1, 4, 67, 1769, 62499, 2718712, 139232874, 8166885453, 538622419492, 39414368809017, 3167946276763685, 277436199521729856, 26298695327956996229, 2683220559729561789860, 293242531673596767855914, 34182433657089947169431226, 4234005174244450641927919376

Offset

0,3

Links

Table of n, a(n) for n=0..17.

Example

G.f.: A(x) = 1 + x + 4*x^2 + 67*x^3 + 1769*x^4 + 62499*x^5 + 2718712*x^6 +...
RELATED SERIES.
A(x)^2 = 1 + 2*x + 9*x^2 + 142*x^3 + 3688*x^4 + 129072*x^5 + 5581063*x^6 +...
A(x)^3 = 1 + 3*x + 15*x^2 + 226*x^3 + 5769*x^4 + 199968*x^5 + 8594032*x^6 +...
A(x*A(x)^3) = 1 + x + 7*x^2 + 106*x^3 + 2754*x^4 + 96488*x^5 + 4175632*x^6 +...
A(x*A(x)^3)^2 = 1 + 2*x + 15*x^2 + 226*x^3 + 5769*x^4 + 199968*x^5 + 8594032*x^6 +...

Prog

(PARI) {a(n)=local(A=[1, 1], Ax); for(i=1, n, A=concat(A, 0); Ax=Ser(A);
A[#A]=Vec(1+subst(Ax^2, x, x*Ax^3) - Ax^3)[#A]); A[n+1]}
for(n=0, 25, print1(a(n), ", "))

Crossrefs

Cf. A240997, A240996.

Sequence in context: A367121 A225772 A329478 * A308382 A083931 A133881

Adjacent sequences: A240998 A240999 A241000 * A241002 A241003 A241004

Keyword

nonn

Author

Paul D. Hanna, Aug 06 2014

Status

approved