A166972 → (θ+1)²

nonnegative coefficients

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

Recurrence

\(\left(n\right)\,P_{n}(x) = \left(-8\,x - 7\,n x^{2} + 8\,n^{2} x - 2\,n^{2} x^{2} + 9\,n^{3} x^{2}\right)\,P_{n-2}(x) + \left(16\,x - 10\,x^{2} + 16\,n x + 5\,n x^{2} - 13\,n x^{3} + 19\,n^{2} x^{2} - 14\,n^{2} x^{3} + 4\,n^{3} x^{3}\right)\,P_{n-2}'(x) + \left(26\,x^{2} - 2\,x^{3} + 26\,n x^{2} + 31\,n x^{3} - 3\,n x^{4} + 43\,n^{2} x^{3} - 4\,n^{2} x^{4}\right)\,P_{n-2}''(x) + \left(6\,x^{3} + 6\,n x^{3} + 9\,n x^{4} + 12\,n^{2} x^{4}\right)\,P_{n-2}'''(x) + \left(n + 4\,n^{2} x\right)\,P_{n-1}(x) + \left(-2\,x + 5\,n x - 4\,n x^{2} - 4\,n^{2} x^{2}\right)\,P_{n-1}'(x) + \left(3\,n x^{2}\right)\,P_{n-1}''(x)\)
\(P_{0}(x) = 1\)
\(P_{1}(x) = 1 + 4\,x\)

Rows ?

nx0x1x2x3x4x5x6x7x8x9x10x11x12x13x14Pn(1)
011
1145
2140950
3133236916718
412672999019522514640
5121396233838131360772536394356
61171192515130369172481117875257044913383372
71136956411048845532247928011589467572968047648964557605332
8110956544234182818814077714720100525232501368046512399683202110728127891248688
918765238849375609740592251339584785392586750203710632552127640397921937592325525011001643776172640
101701219144103866276539724392284477440575282671654502634201267451230839813194261007872204808599923091392200121112487558226480
1115609753196218276955040779923099833439524044240921722250310871814978088862529595232550637878836972160703733573205357497270034103851448830984461358464
1214487802561645852044697721840080447253830560276830019052720575345026597797205568112783537364149092115669537906306563980420691855754476219878000141400921868176032169787148653020871248
131359024204980963014109321321016121845122274435521862020822676793987536710270435454993180187571667565569565323066947637337501440177396020914355919371039137562750026486675644946538194559681927225319678957541021266609552
14128721936398962022434644154730756469282349908629056012380229560438252779253792316960269137679480294440399145639888521173577603681134232985667085343930267008595235295418454129940031479824217472961148043670030720198759288247448122642258844884940597899488416
Technical details
Node id
n-ee924b5c34b83af8 · revision 0 · JSON
Definition source
derived:theta_plus_one_squared
Fingerprint
f6-ffe5be3080cea858

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,        4;
  2: 1,       40,           9;
  3: 1,      332,         369,           16;
  4: 1,     2672,        9990,         1952,           25;
  5: 1,    21396,      233838,       131360,         7725,           36;
  6: 1,   171192,     5151303,      6917248,      1117875,        25704,          49;
  7: 1,  1369564,   110488455,    322479280,    115894675,      7296804,       76489,        64;
  8: 1, 10956544,  2341828188,  14077714720,  10052523250,   1368046512,    39968320,    211072,     81;
  9: 1, 87652388, 49375609740, 592251339584, 785392586750, 203710632552, 12764039792, 193759232, 552501, 100;
%F Formula
The row polynomials satisfy n*P(n,x) = ((8*n^2-8)*x+(9*n^3-2*n^2-7*n)*x^2)*P(n-2,x) + ((16*n+16)*x+(19*n^2+5*n-10)*x^2+(4*n^3-14*n^2-13*n)*x^3)*d/dx P(n-2,x) + ((26*n+26)*x^2+(43*n^2+31*n-2)*x^3-(4*n^2+3*n)*x^4)*d^2/dx^2 P(n-2,x) + ((6*n+6)*x^3+(12*n^2+9*n)*x^4)*d^3/dx^3 P(n-2,x) + (n+4*n^2*x)*P(n-1,x) + ((5*n-2)*x-(4*n^2+4*n)*x^2)*d/dx P(n-1,x) + 3*n*x^2*d^2/dx^2 P(n-1,x) for n >= 2, with P(0,x) = 1, P(1,x) = 1 + 4*x.
%t Mathematica
P[0, x_] := 1;
P[1, x_] := 1 + 4*x;
P[n_, x_] := P[n, x] = Expand[Together[(((8*n^2 - 8)*x + (9*n^3 - 2*n^2 - 7*n)*x^2)*P[n - 2, x] + ((16*n + 16)*x + (19*n^2 + 5*n - 10)*x^2 + (4*n^3 - 14*n^2 - 13*n)*x^3)*D[P[n - 2, x], x] + ((26*n + 26)*x^2 + (43*n^2 + 31*n - 2)*x^3 - (4*n^2 + 3*n)*x^4)*D[P[n - 2, x], {x, 2}] + ((6*n + 6)*x^3 + (12*n^2 + 9*n)*x^4)*D[P[n - 2, x], {x, 3}] + (n + 4*n^2*x)*P[n - 1, x] + ((5*n - 2)*x - (4*n^2 + 4*n)*x^2)*D[P[n - 1, x], x] + 3*n*x^2*D[P[n - 1, x], {x, 2}])/(n)]];
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 + 4*x; for(n = 2, N, my(i = n + 1); v[i] = (((8*n^2 - 8)*x + (9*n^3 - 2*n^2 - 7*n)*x^2)*v[i - 2] + ((16*n + 16)*x + (19*n^2 + 5*n - 10)*x^2 + (4*n^3 - 14*n^2 - 13*n)*x^3)*deriv(v[i - 2], x) + ((26*n + 26)*x^2 + (43*n^2 + 31*n - 2)*x^3 - (4*n^2 + 3*n)*x^4)*deriv(deriv(v[i - 2], x), x) + ((6*n + 6)*x^3 + (12*n^2 + 9*n)*x^4)*deriv(deriv(deriv(v[i - 2], x), x), x) + (n + 4*n^2*x)*v[i - 1] + ((5*n - 2)*x - (4*n^2 + 4*n)*x^2)*deriv(v[i - 1], x) + 3*n*x^2*deriv(deriv(v[i - 1], x), x))/(n)); 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. [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. [SP14] J Schur and G Polya. Über zwei Arten von Faktorenfolgen in der Theorie der algebraischen Gleichungen. J. Reine Angew. Math., 1914.