OFFSET
0,2
COMMENTS
Self-convolution of A351970.
FORMULA
G.f. A(x) satisfies: A(x) = A(x^3/(1 - x)^3) / (1 - x)^2.
E.g.f.: d/dx ( exp(x) * Sum_{n>=0} a(n) * x^(3*n+1) / (3*n+1)! ).
MATHEMATICA
a[0] = 1; a[n_] := a[n] = Sum[Binomial[n + 1, 3 k + 1] a[k], {k, 0, Floor[n/3]}]; Table[a[n], {n, 0, 32}]
nmax = 32; A[_] = 1; Do[A[x_] = A[x^3/(1 - x)^3]/(1 - x)^2 + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Apr 07 2022
STATUS
approved