A166343 → Pₙ + Pₙ₊₁

OEIS matches ?

No OEIS matches.

OEIS submission text ?

Fields to paste into an OEIS submission: the example %e, the formula %F and programs %t (Mathematica) and %o (PARI/GP) that list the same rows. For an OEIS link the row numbers, the offset and the orientation are those of the entry, as matched; the matched terms were compared with the published ones, which identifies the node with the entry on a finite prefix, not for every n. The PARI/GP program is not machine-checked here.

Triangle begins:
   1: 2,    39,      28,         1;
   2: 2,    85,     189,        59,          1;
   3: 2,   179,     880,       776,        122,          1;
   4: 2,   369,    3477,      6462,       2880,        249,          1;
   5: 2,   751,   12524,     42515,      39530,      10013,        504,          1;
   6: 2,  1517,   42581,    242445,     404481,     215439,      33319,       1015,          1;
   7: 2,  3051,  139376,   1258082,    3440304,    3278268,    1085224,     107574,       2038,         1;
   8: 2,  6121,  444573,   6114520,   25802420,   39943722,   23782234,    5168672,     340050,      4085,       1;
   9: 2, 12263, 1392892,  28352677,  176875852,  417206678,  403065376,  159208138,   23628578,   1058963,    8180,     1;
  10: 2, 24549, 4309485, 127001615, 1134491241, 3895517058, 5705033250, 3669327534, 1003635780, 104765929, 3262785, 16371, 1;

The row polynomials satisfy P(n,x) = (1+n*x)*P(n-2,x) + (x-x^2)*d/dx P(n-2,x) + (n+2)*x*P(n-1,x) + (x-x^2)*d/dx P(n-1,x) for n >= 3, with P(1,x) = 2 + 39*x + 28*x^2 + x^3, P(2,x) = 2 + 85*x + 189*x^2 + 59*x^3 + x^4.

P[1, x_] := 2 + 39*x + 28*x^2 + x^3;
P[2, x_] := 2 + 85*x + 189*x^2 + 59*x^3 + x^4;
P[n_, x_] := P[n, x] = Expand[(1 + n*x)*P[n - 2, x] + (x - x^2)*D[P[n - 2, x], x] + (n + 2)*x*P[n - 1, x] + (x - x^2)*D[P[n - 1, x], x]];
Table[CoefficientList[P[n, x], x], {n, 1, 10}] // Flatten

P(N) = {my(v = vector(N)); v[1] = 2 + 39*x + 28*x^2 + x^3; v[2] = 2 + 85*x + 189*x^2 + 59*x^3 + x^4; for(n = 3, N, my(i = n); v[i] = (1 + n*x)*v[i - 2] + (x - x^2)*deriv(v[i - 2], x) + (n + 2)*x*v[i - 1] + (x - x^2)*deriv(v[i - 1], x)); 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)

Proved properties

• Real-rooted. Every row P_n, n >= 1, is zero or real-rooted with all roots in (-inf, +inf)
  Every polynomial of the sequence has only real roots, inside the stated interval.
  by the cited steps: e-e8d26d2660ac53e4

Proof chains

1 incoming proof step. Each step is expanded back to the roots of the proof DAG, the theorem checks of a node's own rule. JSON

• obreschkoff_combination e-e8d26d2660ac53e4

  Hermite–Kakeya–Obreschkoff theorem [Obr63; Fis06, Thm 1.43] ?

  • premise (sequence): Triangle T(n, k) = coefficients of ( t(n, x) ) where t(n, x) = (1-x)^(n+1)*p(n, x)/x, p(n, x) = x*D( p(n-1, x) ), with …
    uses: Every row P_n, n >= 1, is zero or real-rooted with all roots in (-inf, +inf)
    • theorem_check liu_wang e-f953fa70efc27616 root: theorem check

      Liu–Wang theorem [LW07, Thm 2.1, Cor 2.4] ?

    • parent_provers root_domain e-0968a5ea4bacadd5 root: theorem check

      Parent-row provers of polytool polytool ?

  • premise (sequence): Triangle T(n, k) = coefficients of ( t(n, x) ) where t(n, x) = (1-x)^(n+1)*p(n, x)/x, p(n, x) = x*D( p(n-1, x) ), with …
    uses: Rows P_(m-1), P_m interlace (weak, unordered) for every m >= 2
    no recorded property covers this requirement

Related sequences ?

No candidate equalities or related sequences.

References

1. [Fis06] Steve Fisk. Polynomials, roots, and interlacing. arXiv:math/0612833. arXiv
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. [Obr63] Nikola Obreschkoff. Verteilung und berechnung der nullstellen reeller polynome. Hochschulbücher für mathematik, VEB Deutscher Verlag der Wissenschaften, 1963.