OFFSET
0,2
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..300
FORMULA
a(n) = Sum_{k=0..n} (-9)^k * binomial(-1/3,k) * binomial(n,k).
n*a(n) = (11*n-7)*a(n-1) - 10*(n-1)*a(n-2) for n > 1.
a(n) ~ 10^(n + 2/3) / (Gamma(1/3) * 3^(4/3) * n^(2/3)). - Vaclav Kotesovec, May 02 2025
a(n) = hypergeom([1/3, -n], [1], -9). - Stefano Spezia, May 04 2025
MATHEMATICA
Table[Sum[(-9)^k *Binomial[-1/3, k]* Binomial[n, k], {k, 0, n}], {n, 0, 25}] (* Vincenzo Librandi, May 04 2025 *)
PROG
(PARI) a(n) = sum(k=0, n, (-9)^k*binomial(-1/3, k)*binomial(n, k));
(Magma) I:=[4, 25]; [1] cat [n le 2 select I[n] else ((11*n-7)*Self(n-1) - 10*(n-1) *Self(n-2))/n : n in [1..30]]; // Vincenzo Librandi, May 04 2025
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Seiichi Manyama, May 01 2025
STATUS
approved
