cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A245244 Triangle of coefficients of the Pbar polynomials, read by rows.

Original entry on oeis.org

1, -3, 4, 25, -56, 32, -427, 1228, -1184, 384, 12465, -41840, 52416, -29184, 6144, -555731, 2079892, -3076288, 2258688, -829440, 122880, 35135945, -142843304, 237829600, -208562688, 102279168, -26787840, 2949120, -2990414715, 12987478876, -23672564832, 23581133952, -13947525120, 4929576960, -970260480
Offset: 0

Views

Author

Richard P. Brent, Jul 14 2014

Keywords

Comments

Pbar(r,n) is a polynomial of degree r defined by the recurrence
Pbar(r+1,n) = (2*n-1)^2 * Pbar(r,n) - 4*(n-1)^2 * Pbar(r,n-1)
with initial condition Pbar(0,n) = 1.

Examples

			Pbar(1,n) = 4*n-3, Pbar(2,n) = 32*n^2 - 56*n + 25.
Triangle begins:
1,
-3, 4,
25, -56, 32,
-427, 1228, -1184, 384,
12465, -41840, 52416, -29184, 6144,
...
From _Peter Bala_, Jan 22 2018: (Start)
The polynomials Pbar(r,n) as hypergeometric series:
r = 0: n*Pbar(0,n) = n = 1 + 3*(n-1)/(n+1) + 5*(n-1)*(n-2)/((n+1)*(n+2)) + 7*(n-1)*(n-2)*(n-3)/((n+1)*(n+2)*(n+3)) + ..., for n a positive integer (when the series terminates). The identity is also valid for complex n with real part greater than 1/2.
r = 1: n*Pbar(1,n) = n*(4*n - 3) = 1 + 3^3*(n-1)/(n+1) + 5^3*(n-1)*(n-2)/((n+1)*(n+2)) + 7^3*(n-1)*(n-2)*(n-3)/((n+1)*(n+2)*(n+3)) + ..., for n a positive integer (when the series terminates). The identity is also valid for complex n with real part greater than 3/2.
r = 2: n*Pbar(2,n) = n*(32*n^2 - 56*n + 25) = 1 + 3^5*(n-1)/(n+1) + 5^5*(n-1)*(n-2)/((n+1)*(n+2)) + 7^5*(n-1)*(n-2)*(n-3)/((n+1)*(n+2)*(n+3)) + ..., for n a positive integer (when the series terminates). The identity is also valid for complex n with real part greater than 5/2.
The above identities when r = 0 and r = 1 were found by Ramanujan. See Example 5 and Example 13 in Chapter 10 of Berndt. (End)
		

References

  • B. C. Berndt, Ramanujan's Notebooks Part II, Springer-Verlag, Chapter 10, p. 20 and p. 23.

Crossrefs

(-1)^r Pbar(r,0) is sequence A009843. The leading coefficient of Pbar(r,n) is sequence A047053. Cf. also A036970, A083061, A160485 for analogous moments of Bernoulli random walks.

Programs

  • PARI
    N=10; P=vector(N+2); P[1]=1;
    Pbar(r)=P[r+1];
    for (r=0, N, P[r+2] = (2*n-1)^2 * Pbar(r) - 4*(n-1)^2 * subst(Pbar(r),n,n-1) );
    seq=[];  for(r=1,N, seq=concat(seq, Vecrev(P[r])); );  seq
    \\ Joerg Arndt, Jan 27 2015

Formula

In terms of Dumont-Foata polynomials F(r,x,y,z),
Pbar(r,n) = (-4)^r F(r+1,1/2-n,1/2,1/2).
In terms of odd absolute moments of a symmetric Bernoulli random walk with an odd number of steps,
n*C(2*n,n)*Pbar(r,n) = Sum_{k} C(2*n-1,k) * |2*n-1-2*k|^(2*r+1).
In terms of the Pochhammer symbol or ascending factorial (x)_k,
Pbar(r,n) = Sum_{1 <= j <= k <= r+1} (-1)^(j+1)*(1-n)_{k-1}*(2j-1)^(2r+1)/((k-j)!(k)_j).
n*Pbar(r,n) = 1 + 3^(2*r+1)*(n-1)/(n+1) + 5^(2*r+1)*(n-1)*(n-2)/((n+1)*(n+2)) + 7^(2*r+1)*(n-1)*(n-2)*(n-3)/((n+1)*(n+2)*(n+3)) + ... = Sum_{k = 0..n-1} binomial(n-1,k)/binomial(n+k,k)*(2*k + 1)^(2*r+1); this follows easily from the above recurrence. Examples are given below. - Peter Bala, Jan 22 2018

Extensions

More terms from Joerg Arndt, Jan 27 2015