A201637 → (θ+1)²

nonnegative coefficients

Proved for all n all roots real, in (−∞, 0]coefficients of one signHadamard kernelinterlacing (weak)

Recurrence

\(\left(-5 + 2\,n\right)\,P_{n}(x) = \left(-72\,x - 870\,x^{2} + 72\,n x + 1178\,n x^{2} - 16\,n^{2} x - 512\,n^{2} x^{2} + 72\,n^{3} x^{2}\right)\,P_{n-2}(x) + \left(6\,x + 50\,x^{2} - 600\,x^{3} - 4\,n x - 56\,n x^{2} + 740\,n x^{3} + 16\,n^{2} x^{2} - 280\,n^{2} x^{3} + 32\,n^{3} x^{3}\right)\,P_{n-2}'(x) + \left(3\,x^{2} + 62\,x^{3} - 65\,x^{4} - 2\,n x^{2} - 64\,n x^{3} + 66\,n x^{4} + 16\,n^{2} x^{3} - 16\,n^{2} x^{4}\right)\,P_{n-2}''(x) + \left(-5 + 80\,x + 2\,n - 72\,n x + 16\,n^{2} x\right)\,P_{n-1}(x) + \left(-8\,x - 60\,x^{2} + 4\,n x + 64\,n x^{2} - 16\,n^{2} x^{2}\right)\,P_{n-1}'(x)\)
\(P_{1}(x) = 1\)
\(P_{2}(x) = 1\)

Rows ?

nx0x1x2x3x4x5x6x7x8x9x10x11x12x13Pn(1)
111
211
3189
41325487
5188522384995
6120829527104300013265
7145613068704009270025920202545
8196050490513920145350012234242469603489255
9119761795503132800161005002827094416715664258048066981915
101401660534017000000144137500447842304540288504238970880293932801418241825
1118104196963285218560111912500055928658241170634029610379962368359146224036288000032839832025
121162886255342403902464785924300059341834944196132059312297098492928203105978496568356480004829932800825613363575
131326641953563418380220165128706940055976258822427460285404166475547713536746818300814440788301824009471540038406897623040022397626926675
1416542460297264811614841631674544440048356431280643369828227281611636407310745620653910334734418830796618240084449461852800166083963494401052366515200652180214711025
15113095218457189235037276928187540130570039048354446976374106775320208181119424423116846570442506319526472837513987200480243894475276818073301785804803060577420992001708694507520020289055572410625
Technical details
Node id
n-741874bf2346879a · revision 0 · JSON
Definition source
derived:theta_plus_one_squared
Fingerprint
f6-3c9a530453f07c86

Observed on the stored rows ?

  • Degrees: 0, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 (irregular)
  • Palindromic rows: no
  • γ-positive (gamma_positive): no
  • Nonnegative coefficients (nonn): yes
  • Only rational roots (rr_by_factors: every nonzero row is a product of linear factors over ℚ): no

Matches ?

No OEIS matches.

Submission text ?

Atlas rows n = 1, …, ascending powers of x
%e Example
Triangle begins:
   1: 1;
   2: 1;
   3: 1,    8;
   4: 1,   32,     54;
   5: 1,   88,    522,      384;
   6: 1,  208,   2952,     7104,      3000;
   7: 1,  456,  13068,    70400,     92700,     25920;
   8: 1,  960,  50490,   513920,   1453500,   1223424,    246960;
   9: 1, 1976, 179550,  3132800,  16100500,  28270944,  16715664,   2580480;
  10: 1, 4016, 605340, 17000000, 144137500, 447842304, 540288504, 238970880, 29393280;
%F Formula
The row polynomials satisfy (2*n-5)*P(n,x) = -((16*n^2-72*n+72)*x-(72*n^3-512*n^2+1178*n-870)*x^2)*P(n-2,x) - ((4*n-6)*x-(16*n^2-56*n+50)*x^2-(32*n^3-280*n^2+740*n-600)*x^3)*d/dx P(n-2,x) - ((2*n-3)*x^2-(16*n^2-64*n+62)*x^3+(16*n^2-66*n+65)*x^4)*d^2/dx^2 P(n-2,x) + (2*n-5+(16*n^2-72*n+80)*x)*P(n-1,x) + ((4*n-8)*x-(16*n^2-64*n+60)*x^2)*d/dx P(n-1,x) for n >= 3, with P(1,x) = 1, P(2,x) = 1.
%t Mathematica
P[1, x_] := 1;
P[2, x_] := 1;
P[n_, x_] := P[n, x] = Expand[Together[(-((16*n^2 - 72*n + 72)*x - (72*n^3 - 512*n^2 + 1178*n - 870)*x^2)*P[n - 2, x] - ((4*n - 6)*x - (16*n^2 - 56*n + 50)*x^2 - (32*n^3 - 280*n^2 + 740*n - 600)*x^3)*D[P[n - 2, x], x] - ((2*n - 3)*x^2 - (16*n^2 - 64*n + 62)*x^3 + (16*n^2 - 66*n + 65)*x^4)*D[P[n - 2, x], {x, 2}] + (2*n - 5 + (16*n^2 - 72*n + 80)*x)*P[n - 1, x] + ((4*n - 8)*x - (16*n^2 - 64*n + 60)*x^2)*D[P[n - 1, x], x])/(2*n - 5)]];
Table[CoefficientList[P[n, x], x], {n, 1, 10}] // Flatten
%o PARI/GP
P(N) = {my(v = vector(N)); v[1] = 1; v[2] = 1; for(n = 3, N, my(i = n); v[i] = (-((16*n^2 - 72*n + 72)*x - (72*n^3 - 512*n^2 + 1178*n - 870)*x^2)*v[i - 2] - ((4*n - 6)*x - (16*n^2 - 56*n + 50)*x^2 - (32*n^3 - 280*n^2 + 740*n - 600)*x^3)*deriv(v[i - 2], x) - ((2*n - 3)*x^2 - (16*n^2 - 64*n + 62)*x^3 + (16*n^2 - 66*n + 65)*x^4)*deriv(deriv(v[i - 2], x), x) + (2*n - 5 + (16*n^2 - 72*n + 80)*x)*v[i - 1] + ((4*n - 8)*x - (16*n^2 - 64*n + 60)*x^2)*deriv(v[i - 1], x))/(2*n - 5)); v};
row(n, p) = {my(c = Vecrev(p)); c};
T(N) = {my(v = P(N)); concat(vector(#v, i, row(i, v[i])))};
T(10)

polytool c48629374597

Properties

Proof chains ? JSON

Same fingerprint ?

References

  1. [Bra15] Petter Brändén. Unimodality, log-concavity, real-rootedness and beyond. Handbook of enumerative combinatorics, :437–483, March 2015. doi:10.1201/b18255-10
  2. [LW07] Lily L. Liu and Yi Wang. A unified approach to polynomial sequences with only real zeros. Advances in Applied Mathematics, 38(4):542–560, May 2007. doi:10.1016/j.aam.2006.02.003
  3. [SP14] J Schur and G Polya. Über zwei Arten von Faktorenfolgen in der Theorie der algebraischen Gleichungen. J. Reine Angew. Math., 1914.