%I #24 Feb 11 2019 02:53:51
%S 1,1,-2,1,-8,12,1,-2,60,-120,1,-128,-168,0,1680,1,2638,7320,-5040,
%T -25200,-30240,1,-98408,-300828,52800,1053360,1330560,665280,1,
%U 5307118,17914260,2522520,-56456400,-90810720,-60540480,-17297280
%N Triangle of the RBS1 polynomial coefficients.
%C In A160480 we defined the BS1 matrix by BS1[2*m-1,n=1] = 2*beta(2*m) and the recurrence relation BS1 [2*m-1,n] = (2*n-3)/(2*n-2)*(BS1[2*m-1,n-1]- BS1[2*m-3,n-1]/(2*n-3)^2), for positive and negative values of m and n= 1, 2, .. . As usual beta(m) = sum((-1)^k/(1+2*k)^m, k=0..infinity). It is well-known that BS1[1-2*m,n=1] = euler(2*m-2) for m = 1, 2, .., with euler(2*m-2) the Euler numbers A000364. These values together with the recurrence relation lead to BS1[ -1,n] = 1 for n = 1, 2, .. .
%C We discovered that the n-th term of the row coefficients BS1[1-2*m,n] for m = 1, 2, .., can be generated with the rather simple polynomials RBS1(1-2*m,n). Our discovery was enabled by the recurrence relation for the RBS1(1-2*m,n) polynomials which we derived from the recurrence relation for the BS[2*m-1,n] coefficients and the fact that RBS1(-1,n) = 1.
%C The RBS1 polynomials and the polynomials defined by sequence A083061 are related by a shift of +-1/2 and scaling by a power of 2 (see arXiv link). - _Richard P. Brent_, Jul 15 2014
%D B. C. Berndt, Ramanujan's Notebooks Part II, Springer-Verlag, Chapter 10, p. 21.
%H P. Bala, <a href="/A245244/a245244.pdf">A245244 and A160485 and some hypergeometric series evaluations of Ramanujan</a>
%H R. P. Brent, <a href="http://arxiv.org/abs/1407.3533">Generalising Tuenter's binomial sums</a>, arXiv:1407.3533 [math.CO], 2014.
%H R. B. Brent, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL18/Brent/brent5.html">Generalizing Tuenter's Binomial Sums</a>, J. Int. Seq. 18 (2015) # 15.3.2.
%F RBS1(1-2*m,n) = (2*n-1)^2*RBS1(3-2*m,n)-(2*n)*(2*n-1)*RBS1(3-2*m,n+1) for m = 2, 3, .., with RBS1(-1,n) =1 for n = 1, 2, .. .
%F From _Peter Bala_, Jan 22 2019: (Start)
%F The row polynomials RBS1 of the triangle are related to the polynomials Qbar(r,n), r = 0,1,2,..., introduced by Brent by Qbar(r,n) = RBS1(-2*r-1,-n).
%F Recurrence: Qbar(r+1,n) = (2*n + 1)^2*Qbar(r,n) - 2*n(2*n + 1)*Qbar(r,n-1) with Qbar(0,n) = 1 (Brent, equation 19, p.7).
%F Qbar(r,n) = binomial(2*n + 2,n + 1)/(2^(2*n + 1)) * Sum_{k = 0..n} binomial(n,k)/binomial(n + k + 1,k)*(2*k + 1)^(2*r); this follows easily from the above recurrence. Two examples are given below.
%F Qbar(r,n) = 1/4^n * Sum_{k = 0..n} binomial(2*n + 1,n - k)*(2*k + 1)^(2*r) For related polynomial sequences see A036970, A083061 and A245244. (End)
%e The first few rows of the triangle are:
%e [1]
%e [1, -2]
%e [1, -8, 12]
%e [1, -2, 60, -120]
%e [1, -128, -168, 0, 1680]
%e The first few RBS1(1-2*m,n) polynomials are:
%e RBS1(-1,n) = 1
%e RBS1(-3,n) = 1 - 2*n
%e RBS1(-5,n) = 1 - 8*n + 12*n^2
%e RBS1(-7,n) = 1 - 2*n + 60*n^2 - 120*n^3
%e From _Peter Bala_, Jan 22 2019: (Start)
%e Qbar(r,n) = binomial(2*n+2,n+1)/(2^(2*n+1)) * Sum_{k = 0..n} binomial(n,k)/binomial(n+k+1,k)*(2*k + 1)^(2*r):
%e Case r = 2: Qbar(2,n) = binomial(2*n+2,n+1)/2^(2*n+1) * ( 1 + 3^4*n/(n+2) + 5^4*n*(n-1)/((n+2)*(n+3)) + 7^4*n*(n-1)*(n-2)/((n+2)*(n+3)*(n+4)) + ... ) = 12*n^2 + 8*n + 1, valid for n a nonnegative integer (when the series terminates). The identity is also valid for complex n with real part greater than 1 (provided the factor binomial(2*n,n) is replaced with the appropriate expression involving the gamma function).
%e Case r = 3: Qbar(3,n) = binomial(2*n+2,n+1)/(2^(2*n+1)) * ( 1 + 3^6*n/(n+2) + 5^6*n*(n-1)/((n+2)*(n+3)) + 7^6*n*(n-1)*(n-2)/((n+2)*(n+3)*(n+4)) + ... ) = 120*n^3 + 60*n^2 + 2*n + 1, valid for n a nonnegative integer. The identity is also valid for complex n with real part greater than 2.
%e Note, the case r = 0 is equivalent to the identity 1 = binomial(2*n,n)/2^(2*n-1) * ( 1 + (n-1)/(n+1) + (n-1)*(n-2)/((n+1)*(n+2)) + (n-1)*(n-2)*(n-3)/((n+1)*(n+2)*(n+3)) + ... ), which is valid for complex n with real part greater than 0. This identity was found by Ramanujan. See Example 6, Chapter 10 in Berndt. (End)
%p nmax := 8; mmax := nmax: A(1, 1) := 1: RBS1(n, 2) := (2*n-1)^2*1-(2*n)*(2*n-1)*1: for m from 3 to mmax do for k from 0 to m-1 do A(m-1, k+1) := coeff(RBS1(n, m-1), n, k) od; RBS1(n+1, m-1) := 0: for k from 0 to m-1 do RBS1(n+1, m-1) := RBS1(n+1, m-1) + A(m-1, k+1)*(n+1)^k od: RBS1(n, m) := (2*n-1)^2*RBS1(n, m-1)-(2*n)*(2*n-1) * RBS1(n+1, m-1) od: for k from 0 to nmax-1 do A(nmax, k+1) := coeff(RBS1(n, nmax), n, k) od: seq(seq(A(n, m), m=1..n), n=1..nmax);
%Y A160480 is the Beta triangle.
%Y A009389(2*n) equals the second left hand column divided by 2.
%Y A001813 equals the first right hand column.
%Y The absolute values of the row sums equal the Euler numbers A000364.
%Y Cf. A083061, A036970, A245244.
%K easy,sign,tabl
%O 1,3
%A _Johannes W. Meijer_, May 24 2009, Jul 06 2009, Sep 19 2012