OFFSET
0,2
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..1000
FORMULA
G.f.: B(x)^(1/3), where B(x) is the g.f. of A016137.
a(n) = (-3)^n * Sum_{k=0..n} 2^k * binomial(-1/3,k) * binomial(-1/3,n-k).
a(n) ~ 2^(n + 1/3) * 3^n / (Gamma(1/3) * n^(2/3)). - Vaclav Kotesovec, Aug 18 2025
D-finite with recurrence: (-9*n-3)*a(n)+(18*n-6)*a(n-1)+(n+1)*a(n+1) = 0. - Georg Fischer, Aug 29 2025
MATHEMATICA
CoefficientList[Series[1/((1-3*x)*(1-6*x))^(1/3), {x, 0, 33}], x] (* Vincenzo Librandi, Aug 28 2025 *)
PROG
(PARI) a(n) = (-3)^n*sum(k=0, n, 2^k*binomial(-1/3, k)*binomial(-1/3, n-k));
(Magma) R<x> := PowerSeriesRing(Rationals(), 34); f := 1 / ( (1-3*x) * (1-6*x) )^(1/3); coeffs := [ Coefficient(f, n) : n in [0..33] ]; coeffs; // Vincenzo Librandi, Aug 28 2025
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Seiichi Manyama, Aug 18 2025
STATUS
approved
