OFFSET
0,3
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..1000
FORMULA
G.f.: 1+A'(x)*(x*A(x)-x^2)/A(x)^2, where A(x) is g.f. of A109081.
Recurrence: 2*(n-1)*(2*n - 1)*(38*n^2 - 162*n + 163)*a(n) = 2*(380*n^4 - 2380*n^3 + 5200*n^2 - 4676*n + 1431)*a(n-1) + 2*(n-2)*(76*n^3 - 362*n^2 + 502*n - 189)*a(n-2) + 3*(n-3)*(n-2)*(38*n^2 - 86*n + 39)*a(n-3). - Vaclav Kotesovec, Sep 23 2015
a(n) = n*hypergeom([1 - n, 1 - n, n + 1], [1, 3/2], 1/4) for n >= 1. - Peter Luschny, Mar 07 2022
MATHEMATICA
Join[{1}, Table[Sum[Binomial[n-1, n-k] Binomial[ n+k-1, n-k], {k, n}], {n, 25}]] (* Vincenzo Librandi, Sep 23 2015 *)
PROG
(Maxima)
a(n):=sum(binomial(n, k)*binomial(n+k-2, n-k-1), k, 0, n-1)/n;
A(x):=sum(a(n)*x^n, n, 1, 30);
taylor(x*diff(A(x), x)/A(x)-x^2*diff(1/x-1/A(x), x), x, 0, 10);
(Magma) [&+[Binomial(n-1, n-k)*Binomial(n+k-1, n-k): k in [0..n]]: n in [0..25]]; // Vincenzo Librandi, Sep 23 2015
(PARI) a(n) = sum(k=0, n, (binomial(n-1, n-k)*binomial(n+k-1, n-k))) \\ Anders Hellström, Sep 23 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Kruchinin, Sep 23 2015
STATUS
approved