login
A385975
a(n) = Sum_{k=0..n} binomial(n,k) * binomial(2*k,n-k).
4
1, 1, 5, 16, 61, 231, 896, 3515, 13917, 55501, 222595, 896930, 3628120, 14724022, 59922175, 244456581, 999393021, 4093381925, 16793794625, 69001682216, 283889914171, 1169402930621, 4822302149450, 19905773689763, 82243524966936, 340087268899656, 1407397396315006
OFFSET
0,3
LINKS
FORMULA
a(n) = [x^n] (1 + x * (1 + x)^2)^n.
The g.f. exp( Sum_{k>=1} a(k) * x^k/k ) has integer coefficients and equals (1/x) * Series_Reversion( x / (1 + x * (1 + x)^2) ). See A161634.
MATHEMATICA
Table[Sum[Binomial[n, k]Binomial[2*k, n-k], {k, 0, n}], {n, 0, 30}] (* Vincenzo Librandi, Oct 08 2025 *)
PROG
(PARI) a(n) = sum(k=0, n, binomial(n, k)*binomial(2*k, n-k));
(Magma) [&+[Binomial(n, k) * Binomial(2*k, n-k) : k in [0..n] ]: n in [0..40]]; // Vincenzo Librandi, Oct 08 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Sep 25 2025
STATUS
approved