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”).

Table of coefficients of a polynomial sequence related to the Springer numbers.
6

%I #13 Apr 15 2015 11:00:27

%S 1,1,2,3,4,4,11,26,12,8,57,120,136,32,16,361,970,760,560,80,32,2763,

%T 7052,8860,3680,2000,192,64,24611,72530,72884,58520,15120,6496,448,

%U 128,250737,716528,976464,538048,314720,55552,19712,1024,256

%N Table of coefficients of a polynomial sequence related to the Springer numbers.

%C Define a polynomial sequence S(n,x) recursively by

%C (1)... S(n+1,x) = x*S(n,x-1)+(x+1)*S(n,x+1) with S(0,x) = 1.

%C This table lists the coefficients of these polynomials (for n>=1) in ascending powers of x.

%C The first few polynomials are

%C S(0,x) = 1

%C S(1,x) = 2*x+1

%C S(2,x) = 4*x^2+4*x+3

%C S(3,x) = 8*x^3+12*x^2+26*x+11.

%C The sequence [1,1,3,11,57,...] of constant terms of the polynomials is the sequence of Springer numbers A001586. The zeros of the polynomials S(n,-x) lie on the vertical line Re x = 1/2 in the complex plane.

%C Compare the recurrence (1) with the recurrence relation satisfied by the coefficients T(n,k) of the polynomials of A104035, namely

%C (2)... T(n+1,k) = k*T(n,k-1)+(k+1)*T(n,k+1).

%H Peter Bala, <a href="/A185417/a185417.txt">The zeros of the row polynomials of A185417 </a>

%F E.g.f: F(x,t) = 1/(cos(t)-sin(t))*(tan(2*t)+sec(2*t))^x

%F = (cos(t)+sin(t))^x/(cos(t)-sin(t))^(x+1)

%F = 1 + (2*x+1)*t + (4*x^2+4*x+3)*t^2/2! + ....

%F Note that (tan(t)+sec(t))^x is the e.g.f for table A147309.

%F ROW POLYNOMIALS

%F The easily checked identity d/dt F(x,t) = x*F(x-1,t)+(x+1)*F(x+1,t) shows that the row generating polynomials of this table are the polynomials S(n,x) described in the Comments section above.

%F The polynomials S(n,-x) satisfy a Riemann hypothesis: that is, the zeros of S(n,-x) lie on the vertical line Re(x) = 1/2 in the complex plane - see the link.

%F RELATION WITH OTHER SEQUENCES

%F 1st column [1,1,3,11,57,...] is A001586.

%F Row sums sequence [1,3,11,57,...] is also A001586.

%F For n>=1, the values 1/2^n*P(2*n,-1/2) = [1,7,139,5473,...] appear to be A126156.

%e Table begin

%e n\k|.....0.....1.....2.....3.....4.....5......6

%e ===============================================

%e 0..|.....1

%e 1..|.....1.....2

%e 2..|.....3.....4.....4

%e 3..|....11....26....12.....8

%e 4..|....57...120...136....32...16

%e 5..|...361...970...760...560...80.....32

%e 6..|..2763..7052..8860..3680..2000...192....64

%e ...

%p #A185417

%p S := proc(n,x) option remember;

%p description 'polynomials S(n,x)'

%p if n = 0 return 1 else return x*S(n-1,x-1)+(x+1)*S(n-1,x+1)

%p end proc:

%p with(PolynomialTools):

%p for n from 1 to 10 CoefficientList(S(n,x),x); end do;

%t S[0, _] = 1; S[n_, x_] := S[n, x] = x*S[n-1, x-1] + (x+1)*S[n-1, x+1]; Table[ CoefficientList[S[n, x], x], {n, 0, 8}] // Flatten (* _Jean-François Alcover_, Apr 15 2015 *)

%Y Cf A001586 (1st column and row sums), A104035, A126156, A147309, A185415, A185418, A185419

%K nonn,easy,tabl

%O 1,3

%A _Peter Bala_, Jan 28 2011