A154227 → θ+1

nonnegative coefficients

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

Recurrence

\(\left(8 + 4\,n + 4\,n^{2}\right)\,P_{n}(x) = \left(4\,n^{2} x^{2} + 3\,n^{3} x^{2} - 3\,n^{4} x^{2} - 3\,n^{5} x^{2} - n^{6} x^{2}\right)\,P_{n-4}(x) + \left(-16\,n x - 12\,n x^{2} - 28\,n^{2} x - 20\,n^{2} x^{2} - 16\,n^{3} x - 12\,n^{3} x^{2} - 4\,n^{4} x - 4\,n^{4} x^{2}\right)\,P_{n-3}(x) + \left(-16 - 16\,x - 8\,x^{2} - 12\,n + 8\,n x - 4\,n x^{2} - 4\,n^{2} + 20\,n^{2} x - 4\,n^{2} x^{2} + 16\,n^{3} x + 4\,n^{4} x\right)\,P_{n-2}(x) + \left(24 + 16\,x + 16\,n + 8\,n x + 8\,n^{2} + 8\,n^{2} x\right)\,P_{n-1}(x)\)
\(P_{0}(x) = 1\)
\(P_{1}(x) = 1 + 2\,x\)
\(P_{2}(x) = 1 + 16\,x + 3\,x^{2}\)
\(P_{3}(x) = 1 + 38\,x + 57\,x^{2} + 4\,x^{3}\)

Rows ?

nx0x1x2x3x4x5x6x7x8x9x10x11x12x13x14Pn(1)
011
1123
2116320
3138574100
41704741405690
511141776236828564550
611724887224328145516736160
7124611301114196142745226028618291960
8133823280422484154729063372654320135292682800
91450440821282264992637011911644224396211755220251025528360
10158478219338968046384130146368656649377826779360234657292011268176000
111742131745807850017564197011311378441319660818281027152181766254391504081122934298640
12192621257417744732571316195640705520418178793948854274027210283691514436183077943855561334791379840
1311138330828364724641653800235292235228581642180775761876778029444383528428733076004701002992761323312739714430054522800
1411380499215709254084361100615113477615220109529991613728684473473281408242749319189129358700959442135321277622421632659660155688838883840
Technical details
Node id
n-97f32c3deac53582 · revision 0 · JSON
Definition source
derived:theta_plus_one
Fingerprint
f6-1bf11212c793d813

Observed on the stored rows ?

  • Degrees: 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14 (affine: slope 1, intercept 0)
  • 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 = 0, …, ascending powers of x
%e Example
Triangle begins:
  0: 1;
  1: 1,   2;
  2: 1,  16,     3;
  3: 1,  38,    57,       4;
  4: 1,  70,   474,     140,       5;
  5: 1, 114,  1776,    2368,     285,        6;
  6: 1, 172,  4887,   22432,    8145,      516,       7;
  7: 1, 246, 11301,  114196,  142745,    22602,     861,      8;
  8: 1, 338, 23280,  422484, 1547290,   633726,   54320,   1352,    9;
  9: 1, 450, 44082, 1282264, 9926370, 11911644, 2243962, 117552, 2025, 10;
%F Formula
The row polynomials satisfy (4*n^2+4*n+8)*P(n,x) = -(n^6+3*n^5+3*n^4-3*n^3-4*n^2)*x^2*P(n-4,x) - ((4*n^4+16*n^3+28*n^2+16*n)*x+(4*n^4+12*n^3+20*n^2+12*n)*x^2)*P(n-3,x) - (4*n^2+12*n+16-(4*n^4+16*n^3+20*n^2+8*n-16)*x+(4*n^2+4*n+8)*x^2)*P(n-2,x) + (8*n^2+16*n+24+(8*n^2+8*n+16)*x)*P(n-1,x) for n >= 4, with P(0,x) = 1, P(1,x) = 1 + 2*x, P(2,x) = 1 + 16*x + 3*x^2, P(3,x) = 1 + 38*x + 57*x^2 + 4*x^3.
%t Mathematica
P[0, x_] := 1;
P[1, x_] := 1 + 2*x;
P[2, x_] := 1 + 16*x + 3*x^2;
P[3, x_] := 1 + 38*x + 57*x^2 + 4*x^3;
P[n_, x_] := P[n, x] = Expand[Together[(-(n^6 + 3*n^5 + 3*n^4 - 3*n^3 - 4*n^2)*x^2*P[n - 4, x] - ((4*n^4 + 16*n^3 + 28*n^2 + 16*n)*x + (4*n^4 + 12*n^3 + 20*n^2 + 12*n)*x^2)*P[n - 3, x] - (4*n^2 + 12*n + 16 - (4*n^4 + 16*n^3 + 20*n^2 + 8*n - 16)*x + (4*n^2 + 4*n + 8)*x^2)*P[n - 2, x] + (8*n^2 + 16*n + 24 + (8*n^2 + 8*n + 16)*x)*P[n - 1, x])/(4*n^2 + 4*n + 8)]];
Table[CoefficientList[P[n, x], x], {n, 0, 9}] // Flatten
%o PARI/GP
P(N) = {my(v = vector(N + 1)); v[1] = 1; v[2] = 1 + 2*x; v[3] = 1 + 16*x + 3*x^2; v[4] = 1 + 38*x + 57*x^2 + 4*x^3; for(n = 4, N, my(i = n + 1); v[i] = (-(n^6 + 3*n^5 + 3*n^4 - 3*n^3 - 4*n^2)*x^2*v[i - 4] - ((4*n^4 + 16*n^3 + 28*n^2 + 16*n)*x + (4*n^4 + 12*n^3 + 20*n^2 + 12*n)*x^2)*v[i - 3] - (4*n^2 + 12*n + 16 - (4*n^4 + 16*n^3 + 20*n^2 + 8*n - 16)*x + (4*n^2 + 4*n + 8)*x^2)*v[i - 2] + (8*n^2 + 16*n + 24 + (8*n^2 + 8*n + 16)*x)*v[i - 1])/(4*n^2 + 4*n + 8)); v};
row(n, p) = {my(c = Vecrev(p)); c};
T(N) = {my(v = P(N)); concat(vector(#v, i, row(i - 1, v[i])))};
T(9)

polytool c48629374597

Properties

Proof chains ? JSON

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.