|
|
A081671
|
|
Expansion of e.g.f. exp(4x) * I_0(2x).
|
|
30
|
|
|
1, 4, 18, 88, 454, 2424, 13236, 73392, 411462, 2325976, 13233628, 75682512, 434662684, 2505229744, 14482673832, 83940771168, 487610895942, 2838118247064, 16547996212044, 96635257790352, 565107853947444, 3308820294176016, 19395905063796312, 113814537122646432
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
COMMENTS
|
Largest coefficient of (1+4x+x^2)^n. - Paul Barry, Dec 15 2003
Also number of paths from (0,0) to (n,0) using steps U=(1,1), H=(1,0) and D=(1,-1), the H steps come in 4 colors. - N-E. Fahssi, Feb 05 2008
Diagonal of rational function 1/(1 - (x^2 + 4*x*y + y^2)). - Gheorghe Coserea, Aug 01 2018
|
|
LINKS
|
|
|
FORMULA
|
a(n) = Sum_{m=0..n} Sum_{k=0..m} C(n, m)*C(m, k)*C(2k, k).
a(n) = Sum_{k=0..n} 2^(n-k) * C(n, k) * C(2*k, k). - Paul Barry, Jan 27 2005
a(n) = Sum_{k=0..n} 6^(n-k) * (-1)^k * C(n,k) * C(2*k,k). - Paul D. Hanna, Dec 09 2018
D-finite with recurrence: n*a(n) = 4*(2*n-1)*a(n-1) - 12*(n-1)*a(n-2). - Vaclav Kotesovec, Oct 13 2012
a(n) = 2^n*hypergeom([-n,1/2], [1], -2). - Peter Luschny, Apr 26 2016
a(n) = Sum_{k=0..floor(n/2)} 4^(n-2*k) * binomial(n,2*k) * binomial(2*k,k). - Seiichi Manyama, May 04 2019
a(n) = (1/Pi) * Integral_{x = -1..1} (2 + 4*x^2)^n/sqrt(1 - x^2) dx = (1/Pi) * Integral_{x = -1..1} (6 - 4*x^2)^n/sqrt(1 - x^2) dx . - Peter Bala, Jan 27 2020
exp(Sum_{n >= 1} a(n)*x^n/n) = 1 + 4*x + 17*x^2 + 76*x^3 + 354*x^4 + ... is the o.g.f. of A005572.
The Gauss congruences a(n*p^k) == a(n*p^(k-1)) (mod p^k) hold for prime p and positive integers n and k. (End)
|
|
MAPLE
|
seq(simplify(2^n*hypergeom([-n, 1/2], [1], -2)), n=0..23); # Peter Luschny, Apr 26 2016
seq(simplify(GegenbauerC(n, -n, -2)), n=0..23); # Peter Luschny, May 09 2016
|
|
MATHEMATICA
|
Table[SeriesCoefficient[1/Sqrt[(1-2*x)*(1-6*x)], {x, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Oct 13 2012 *)
|
|
PROG
|
(Maxima) a(n):=coeff(expand((1+4*x+x^2)^n), x^n);
(PARI) x='x+O('x^66); Vec(1/sqrt((1-2*x)*(1-6*x))) \\ Joerg Arndt, May 07 2013
(PARI) {a(n) = sum(k=0, n\2, 4^(n-2*k)*binomial(n, 2*k)*binomial(2*k, k))} \\ Seiichi Manyama, May 04 2019
|
|
CROSSREFS
|
m-th binomial transforms of A000984: A126869 (m = -2), A002426 (m = -1 and m = -3 for signed version), A000984 (m = 0 and m = -4 for signed version), A026375 (m = 1 and m = -5 for signed version), A081671 (m = 2 and m = -6 for signed version), A098409 (m = 3 and m = -7 for signed version), A098410 (m = 4 and m = -8 for signed version), A104454 (m = 5 and m = -9 for signed version).
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|