A346504 G.f. A(x) satisfies: A(x) = 1 + x + x^3 * A(x)^2 / (1 - x).

1, 1, 0, 1, 3, 4, 6, 14, 28, 49, 95, 196, 386, 754, 1524, 3102, 6258, 12700, 26032, 53440, 109772, 226457, 468863, 972300, 2020274, 4208530, 8784556, 18365322, 38461110, 80682740, 169501696, 356579216, 751138916, 1584281062, 3345404514, 7072055268, 14965933024, 31702754496

Offset

0,5

Links

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

Formula

a(0) = a(1) = 1, a(2) = 0; a(n) = a(n-1) + Sum_{k=0..n-3} a(k) * a(n-k-3).

G.f.: (1-x)*(1 - sqrt((1 - x - 4*x^3 - 4*x^4)/(1-x))) / (2*x^3). - Vaclav Kotesovec, Sep 27 2023

Mathematica

nmax = 37; A[_] = 0; Do[A[x_] = 1 + x + x^3 A[x]^2/(1 - x) + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
a[0] = a[1] = 1; a[2] = 0; a[n_] := a[n] = a[n - 1] + Sum[a[k] a[n - k - 3], {k, 0, n - 3}]; Table[a[n], {n, 0, 37}]
CoefficientList[Series[(1 - x)*(1 - Sqrt[(1 - x - 4*x^3 - 4*x^4)/(1 - x)]) / (2*x^3), {x, 0, 40}], x] (* Vaclav Kotesovec, Sep 27 2023 *)

Crossrefs

Cf. A023426, A023431, A025243, A090345, A248100, A257290, A346503.

Sequence in context: A271618 A137820 A049892 * A063477 A168219 A129827

Adjacent sequences: A346501 A346502 A346503 * A346505 A346506 A346507

Keyword

nonn

Author

Ilya Gutkovskiy, Jul 21 2021

Status

approved