OFFSET
0,2
COMMENTS
Partial sums of A213282.
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..500
FORMULA
G.f.: Sum_{k>=0} ( binomial(3*k,k) / (2*k + 1) ) * x^k / (1 - x)^(3*k+1).
a(0) = 1; a(n) = a(n-1) + Sum_{i=0..n-1} Sum_{j=0..n-i-1} a(i) * a(j) * a(n-i-j-1).
a(n) ~ 2^(n - 1/2) / (sqrt(3*Pi*(2 - (2 - sqrt(2))^(1/3)/2^(2/3) - 1/(2*(2 - sqrt(2)))^(1/3))) * n^(3/2) * (2 - 3/(sqrt(2) - 1)^(1/3) + 3*(sqrt(2) - 1)^(1/3))^n). - Vaclav Kotesovec, Nov 04 2021
a(n) = (1/n) * Sum_{k=0..floor((n-1)/2)} 2^(n-k) * binomial(n,k) * binomial(2*n-k,n-1-2*k) for n > 0. - Seiichi Manyama, Apr 01 2024
MATHEMATICA
nmax = 23; A[_] = 0; Do[A[x_] = (1 + x A[x]^3)/(1 - x) + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
nmax = 23; CoefficientList[Series[Sum[(Binomial[3 k, k]/(2 k + 1)) x^k/(1 - x)^(3 k + 1), {k, 0, nmax}], {x, 0, nmax}], x]
a[0] = 1; a[n_] := a[n] = a[n - 1] + Sum[Sum[a[i] a[j] a[n - i - j - 1], {j, 0, n - i - 1}], {i, 0, n - 1}]; Table[a[n], {n, 0, 23}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Jul 25 2021
STATUS
approved