OFFSET
0,3
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..386 (terms 0..100 from Vincenzo Librandi)
FORMULA
a(n) = [x^n] (1+n*x+x^2)^n.
a(n) = Sum_{k=0..floor(n/2)} binomial(n, k)*binomial(n-k, n-2*k)*n^(n-2*k).
a(n) ~ BesselI(0,2) * n^n. - Vaclav Kotesovec, Apr 17 2014
a(n) = GegenbauerPoly(n,-n,-n/2). - Emanuele Munarini, Oct 20 2016
From Ilya Gutkovskiy, Sep 20 2017: (Start)
a(n) = [x^n] 1/sqrt((1 + 2*x - n*x)*(1 - 2*x - n*x)).
a(n) = n! * [x^n] exp(n*x)*BesselI(0,2*x). (End)
From Seiichi Manyama, May 01 2019: (Start)
a(n) = Sum_{k=0..n} (n-2)^(n-k) * binomial(n,k) * binomial(2*k,k).
a(n) = Sum_{k=0..n} (n+2)^(n-k) * (-1)^k * binomial(n,k) * binomial(2*k,k). (End)
MATHEMATICA
Flatten[{1, Table[Sum[Binomial[n, k]*Binomial[n-k, n-2*k]*n^(n-2*k), {k, 0, Floor[n/2]}], {n, 1, 20}]}] (* Vaclav Kotesovec, Apr 17 2014 *)
Table[GegenbauerC[n, -n, -n/2] + KroneckerDelta[n, 0], {n, 0, 100}] (* Emanuele Munarini, Oct 20 2016 *)
Table[SeriesCoefficient[(1 + n*x + x^2)^n, {x, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Feb 13 2023 *)
PROG
(Maxima) a(n):=coeff(expand((1+n*x+x^2)^n), x, n);
(Maxima) makelist(ultraspherical(n, -n, -n/2), n, 0, 12); /* Emanuele Munarini, Oct 20 2016 */
makelist(a(n), n, 0, 20);
(Magma) P<x>:=PolynomialRing(Integers()); [ Coefficients((1+n*x+x^2)^n)[n+1]: n in [0..22] ]; // Klaus Brockhaus, Mar 02 2011
(PARI) {a(n) = sum(k=0, n, (n-2)^(n-k)*binomial(n, k)*binomial(2*k, k))} \\ Seiichi Manyama, May 01 2019
(PARI) a(n) = polcoef((1+n*x+x^2)^n, n); \\ Michel Marcus, May 01 2019
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Emanuele Munarini, Mar 02 2011
STATUS
approved