editing
approved
editing
approved
(tail $ concat $ transpose [a073711_list, a073712_list])
(PARI) {a(n)=local(A=1); for(i=0, #binary(n), A=subst(A, x, x^2+x*O(x^n))+x*subst(A, x, x^2+x*O(x^n))^2); polcoeff(A, n)} \\ _Paul D. Hanna_, Dec 21 2012
for(n=0, 65, print1(a(n), ", ")) \\ _Paul D. Hanna_, Dec 21 2012
approved
editing
reviewed
approved
proposed
reviewed
editing
proposed
max = 64; f[x_] := Sum[ a[k]*x^k, {k, 0, max}]; a[0] = a[1] = a[2] = 1; coes = CoefficientList[ Series[ f[x] - f[x^2] - x*f[x^2]^2, {x, 0, max}], x]; Table[a[k] , {k, 0, max}] /. Solve[Thread[coes == 0]] // First (* Jean-François Alcover, Mar 06 2013 *)
For[A = 1; n = 1, n <= 65, n++, A = (Normal[A] /. x -> x^2) + x*(Normal[A] /. x -> x^2)^2 + O[x]^n]; CoefficientList[A, x] (* Jean-François Alcover, Mar 06 2013, updated Apr 23 2016 *)
approved
editing
_Reinhard Zumkeller_, , <a href="/A073711/b073711.txt">Table of n, a(n) for n = 0..10000</a>
proposed
approved
editing
proposed
max = 64; f[x_] := Sum[ a[k]*x^k, {k, 0, max}]; a[0] = a[1] = a[2] = 1; coes = CoefficientList[ Series[ f[x] - f[x^2] - x*f[x^2]^2, {x, 0, max}], x]; Table[a[k] , {k, 0, max}] /. Solve[Thread[coes == 0]] // First (* Jean-François Alcover, Mar 06 2013 *)
approved
editing