proposed
approved
proposed
approved
editing
proposed
Essentially the same as A001679. - Eric W. Weisstein, Mar 25 2022
approved
editing
editing
approved
a(n) ~ c * d^n / n^(3/2), where d = A246403 = 2.18946198566085056388702757711... and c = 0.4213018528699249210965028.421301852869924921096502830935802411658488216342994235732491571594804013 . - Vaclav Kotesovec, Jun 26 2014
approved
editing
proposed
approved
editing
proposed
max = 36; G001678[x_] := Sum[a[k]*x^k, {k, 0, max}]; a[0] = a[1] = a[3] = 0; a[2] = a[4] = a[5] = 1; coes = CoefficientList[ Series[G001678[x] - x^2/((1 + x)*Product[(1 - x^k)^a[k + 1], {k, 1, max}]), {x, 0, max}], x]; sol = Solve[Thread[coes == 0]][[1]]; f[x_] = G001678[x] /. sol; se = Series[ (1 + ((1 + x)/x)*f[x] - (f[x]^2 + f[x^2])/(2*x)), {x, 0, max}]; A059123 = Join[{0}, CoefficientList[se, x][[2 ;; -2]]] (* Jean-François Alcover, May 25 2012, from g.f. *)
terms = 36; (* F = G001678 *) F[_] = 0; Do[F[x_] = (x^2/(1 + x))*Exp[Sum[ F[x^k]/(k*x^k), {k, 1, j}]] + O[x]^j // Normal, {j, 1, terms + 1}];
G[x_] = 1 + ((1 + x)/x)*F[x] - (F[x]^2 + F[x^2])/(2*x) + O[x]^terms;
CoefficientList[G[x] - 1, x] (* Jean-François Alcover, May 25 2012, updated Jan 12 2018 *)
approved
editing
editing
approved
Wolfdieter Lang , Jan 09 2001