login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A006442
Expansion of 1/sqrt(1 - 10*x + x^2).
12
1, 5, 37, 305, 2641, 23525, 213445, 1961825, 18205345, 170195525, 1600472677, 15122515985, 143457011569, 1365435096485, 13033485491077, 124715953657025, 1195966908404545, 11490534389896325, 110584004488276645, 1065853221648055025
OFFSET
0,2
COMMENTS
Number of Delannoy paths from (0,0) to (n,n) with steps U(0,1), H(1,0) and D(1,1) where H can choose from two colors. - Paul Barry, May 25 2005
Number of lattice paths from (0,0) to (n,n) using steps (0,1), (1,1), and two kinds of steps (1,0). - Joerg Arndt, Jul 01 2011
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. - Peter Bala, Jan 09 2022
LINKS
Hacène Belbachir and Abdelghani Mehdaoui, Recurrence relation associated with the sums of square binomial coefficients, Quaestiones Mathematicae (2021) Vol. 44, Issue 5, 615-624.
Hacène Belbachir, Abdelghani Mehdaoui and László Szalay, Diagonal Sums in the Pascal Pyramid, II: Applications, J. Int. Seq., Vol. 22 (2019), Article 19.3.5.
Ömür Deveci and Anthony G. Shannon, Some aspects of Neyman triangles and Delannoy arrays, Mathematica Montisnigri (2021) Vol. L, 36-43.
Tony D. Noe, On the Divisibility of Generalized Central Trinomial Coefficients, Journal of Integer Sequences, Vol. 9 (2006), Article 06.2.7.
FORMULA
Legendre polynomial evaluated at 5. - Michael Somos, Dec 04 2001
G.f.: 1/sqrt(1 - 10*x + x^2).
a(n) equals the central coefficient of (1 + 5*x + 6*x^2)^n. - Paul D. Hanna, Jun 03 2003
a(n) equals the (n+1)-th term of the binomial transform of 1/(1-2x)^(n+1). - Paul D. Hanna, Sep 29 2003
a(n) = Sum_{k=0..n} 2^k*binomial(n, k)*binomial(n+k, k). - Benoit Cloitre, Apr 13 2004
a(n) = Sum_{k=0..n} binomial(n,k)^2 * 2^k * 3^(n-k). - Paul D. Hanna, Feb 04 2012
E.g.f.: exp(5*x)*Bessel_I(0, 2*sqrt(6)*x); - Paul Barry, May 25 2005
D-finite with recurrence: n*a(n) - 5*(2n-1)*a(n-1) + (n-1)*a(n-2) = 0 [Eq (4) in the T. D. Noe article]. R. J. Mathar, Jun 26 2012
a(n) ~ (5 + 2*sqrt(6))^n/(2*sqrt(Pi*n)*sqrt(5*sqrt(6) - 12)). - Vaclav Kotesovec, Oct 05 2012
a(n) = hypergeom([-n, n+1], [1], -2). - Peter Luschny, May 23 2014
a(n) = Sum_{k=0..n} 2^k * C(2*k, k) * C(n+k, n-k). - Paul D. Hanna, Aug 17 2014
a(n) = Sum_{k=0..n} (k+1) * 3^k * (-1)^(n-k) * binomial(n,k) * binomial(n+k+1,n) / (n+k+1). - Vladimir Kruchinin, Nov 23 2014
From Peter Bala, Nov 28 2021: (Start)
a(n) = (1/3)*(1/2)^n*Sum_{k >= n} binomial(k,n)^2*(2/3)^k.
a(n) = (1/3)^(n+1)*hypergeom([n+1, n+1], [1], 2/3).
a(n) = (2^n)*hypergeom([-n, -n], [1], 3/2).
a(n) = [x^n] ((x - 1)*(3 - 2*x))^n
a(n) = (1/2)^n*A098270(n). (End)
MAPLE
seq(orthopoly[P](n, 5), n = 0 .. 20); # Robert Israel, Aug 18 2014
MATHEMATICA
Table[LegendreP[n, 5], {n, 0, 19}] (* Arkadiusz Wesolowski, Aug 13 2012 *)
CoefficientList[Series[1 / Sqrt[1 - 10 x + x^2], {x, 0, 20}], x] (* Vincenzo Librandi, Nov 23 2014 *)
PROG
(PARI) a(n)=subst(pollegendre(n), x, 5)
(PARI) /* as lattice paths: same as in A092566 but use */
steps=[[1, 0], [1, 0], [0, 1], [1, 1]]; /* note the double [1, 0] */
/* Joerg Arndt, Jul 01 2011 */
(PARI) {a(n)=sum(k=0, n, binomial(n, k)^2*2^k*3^(n-k))} /* Paul D. Hanna */
(PARI) {a(n) = sum(k=0, n, 2^k * binomial(2*k, k) * binomial(n+k, n-k) )}
for(n=0, 25, print1(a(n), ", ")) \\ Paul D. Hanna, Aug 17 2014
(Magma) [Evaluate(LegendrePolynomial(n), 5): n in [0..40]]; // G. C. Greubel, May 21 2023
(SageMath) [gen_legendre_P(n, 0, 5) for n in range(41)] # G. C. Greubel, May 21 2023
CROSSREFS
Sequences of the form LegendreP(n, 2*m+1): A000012 (m=0), A001850 (m=1), this sequence (m=2), A084768 (m=3), A084769 (m=4).
Cf. A098270, A243943 (a(n)^2).
Sequence in context: A091126 A066381 A078253 * A208675 A084212 A358883
KEYWORD
nonn,easy
STATUS
approved