A173018 → x ↦ x+1 → (θ+1)²

nonnegative coefficients

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

Recurrence

\(\left(4 - 4\,n + n^{2}\right)\,P_{n}(x) = \left(-4 - 60\,x + 132\,x^{2} + 12\,n + 128\,n x - 284\,n x^{2} - 13\,n^{2} - 96\,n^{2} x + 221\,n^{2} x^{2} + 6\,n^{3} + 32\,n^{3} x - 74\,n^{3} x^{2} - n^{4} - 4\,n^{4} x + 9\,n^{4} x^{2}\right)\,P_{n-2}(x) + \left(-11\,x + 21\,x^{2} + 124\,x^{3} + 20\,n x - 64\,n x^{2} - 244\,n x^{3} - 11\,n^{2} x + 68\,n^{2} x^{2} + 167\,n^{2} x^{3} + 2\,n^{3} x - 28\,n^{3} x^{2} - 46\,n^{3} x^{3} + 4\,n^{4} x^{2} + 4\,n^{4} x^{3}\right)\,P_{n-2}'(x) + \left(-3\,x^{2} + 17\,x^{3} + 20\,x^{4} + 4\,n x^{2} - 32\,n x^{3} - 36\,n x^{4} - n^{2} x^{2} + 20\,n^{2} x^{3} + 21\,n^{2} x^{4} - 4\,n^{3} x^{3} - 4\,n^{3} x^{4}\right)\,P_{n-2}''(x) + \left(-6 - 32\,x + 13\,n + 48\,n x - 9\,n^{2} - 24\,n^{2} x + 2\,n^{3} + 4\,n^{3} x\right)\,P_{n-1}(x) + \left(-7\,x + 16\,x^{2} + 8\,n x - 32\,n x^{2} - 2\,n^{2} x + 20\,n^{2} x^{2} - 4\,n^{3} x^{2}\right)\,P_{n-1}'(x)\)
\(P_{1}(x) = 1\)
\(P_{2}(x) = 1\)

Rows ?

nx0x1x2x3x4x5x6Pn(1)
111
211
3246
4624939
52414412616310
61209601350480252935
7720720014040864015503632186
850406048015120013440045150453649400855
Technical details
Node id
n-2ee019a0ed5b4470 · revision 0 · JSON
Definition source
derived:theta_plus_one_squared
Fingerprint
f6-856ad640314d7fe1

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:      2,       4;
   4:      6,      24,        9;
   5:     24,     144,      126,       16;
   6:    120,     960,     1350,      480,       25;
   7:    720,    7200,    14040,     8640,     1550,      36;
   8:   5040,   60480,   151200,   134400,    45150,    4536,     49;
   9:  40320,  564480,  1723680,  2016000,  1020600,  208656,  12446,    64;
  10: 362880, 5806080, 20956320, 30481920, 20853000, 6713280, 889350, 32640, 81;
%F Formula
The row polynomials satisfy (n^2-4*n+4)*P(n,x) = -(n^4-6*n^3+13*n^2-12*n+4+(4*n^4-32*n^3+96*n^2-128*n+60)*x-(9*n^4-74*n^3+221*n^2-284*n+132)*x^2)*P(n-2,x) + ((2*n^3-11*n^2+20*n-11)*x+(4*n^4-28*n^3+68*n^2-64*n+21)*x^2+(4*n^4-46*n^3+167*n^2-244*n+124)*x^3)*d/dx P(n-2,x) - ((n^2-4*n+3)*x^2+(4*n^3-20*n^2+32*n-17)*x^3+(4*n^3-21*n^2+36*n-20)*x^4)*d^2/dx^2 P(n-2,x) + (2*n^3-9*n^2+13*n-6+(4*n^3-24*n^2+48*n-32)*x)*P(n-1,x) - ((2*n^2-8*n+7)*x+(4*n^3-20*n^2+32*n-16)*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[(-(n^4 - 6*n^3 + 13*n^2 - 12*n + 4 + (4*n^4 - 32*n^3 + 96*n^2 - 128*n + 60)*x - (9*n^4 - 74*n^3 + 221*n^2 - 284*n + 132)*x^2)*P[n - 2, x] + ((2*n^3 - 11*n^2 + 20*n - 11)*x + (4*n^4 - 28*n^3 + 68*n^2 - 64*n + 21)*x^2 + (4*n^4 - 46*n^3 + 167*n^2 - 244*n + 124)*x^3)*D[P[n - 2, x], x] - ((n^2 - 4*n + 3)*x^2 + (4*n^3 - 20*n^2 + 32*n - 17)*x^3 + (4*n^3 - 21*n^2 + 36*n - 20)*x^4)*D[P[n - 2, x], {x, 2}] + (2*n^3 - 9*n^2 + 13*n - 6 + (4*n^3 - 24*n^2 + 48*n - 32)*x)*P[n - 1, x] - ((2*n^2 - 8*n + 7)*x + (4*n^3 - 20*n^2 + 32*n - 16)*x^2)*D[P[n - 1, x], x])/(n^2 - 4*n + 4)]];
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] = (-(n^4 - 6*n^3 + 13*n^2 - 12*n + 4 + (4*n^4 - 32*n^3 + 96*n^2 - 128*n + 60)*x - (9*n^4 - 74*n^3 + 221*n^2 - 284*n + 132)*x^2)*v[i - 2] + ((2*n^3 - 11*n^2 + 20*n - 11)*x + (4*n^4 - 28*n^3 + 68*n^2 - 64*n + 21)*x^2 + (4*n^4 - 46*n^3 + 167*n^2 - 244*n + 124)*x^3)*deriv(v[i - 2], x) - ((n^2 - 4*n + 3)*x^2 + (4*n^3 - 20*n^2 + 32*n - 17)*x^3 + (4*n^3 - 21*n^2 + 36*n - 20)*x^4)*deriv(deriv(v[i - 2], x), x) + (2*n^3 - 9*n^2 + 13*n - 6 + (4*n^3 - 24*n^2 + 48*n - 32)*x)*v[i - 1] - ((2*n^2 - 8*n + 7)*x + (4*n^3 - 20*n^2 + 32*n - 16)*x^2)*deriv(v[i - 1], x))/(n^2 - 4*n + 4)); 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. [Bra06] Petter Brändén. On linear transformations preserving the Pólya frequency property. Transactions of the American Mathematical Society, 358(8):3697–3716, February 2006. doi:10.1090/s0002-9947-06-03856-6
  2. [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
  3. [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
  4. [MW08] Shi-Mei Ma and Yi Wang. q-Eulerian polynomials and polynomials with only real zeros. The Electronic Journal of Combinatorics, 15(1):Research Paper 17, January 2008. doi:10.37236/741
  5. [SP14] J Schur and G Polya. Über zwei Arten von Faktorenfolgen in der Theorie der algebraischen Gleichungen. J. Reine Angew. Math., 1914.