A103406: Triangle read by rows: n-th row = unsigned coefficients of the characteristic polynomials of an n X n matrix with 2's on the diagonal and 1's elsewhere.
A007318 → θ+1
A103406 nonnegative coefficients only rational roots
Proved for all n: all roots real, in (−∞, 0] · coefficients of one sign · interlacing (weak) · Hadamard kernel what this means
Rows ?
| n | x0 | x1 | x2 | x3 | x4 | x5 | x6 | x7 | x8 | x9 | x10 | x11 | x12 | x13 | x14 | Pn(1) |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | 1 | 1 | ||||||||||||||
| 1 | 1 | 2 | 3 | |||||||||||||
| 2 | 1 | 4 | 3 | 8 | ||||||||||||
| 3 | 1 | 6 | 9 | 4 | 20 | |||||||||||
| 4 | 1 | 8 | 18 | 16 | 5 | 48 | ||||||||||
| 5 | 1 | 10 | 30 | 40 | 25 | 6 | 112 | |||||||||
| 6 | 1 | 12 | 45 | 80 | 75 | 36 | 7 | 256 | ||||||||
| 7 | 1 | 14 | 63 | 140 | 175 | 126 | 49 | 8 | 576 | |||||||
| 8 | 1 | 16 | 84 | 224 | 350 | 336 | 196 | 64 | 9 | 1280 | ||||||
| 9 | 1 | 18 | 108 | 336 | 630 | 756 | 588 | 288 | 81 | 10 | 2816 | |||||
| 10 | 1 | 20 | 135 | 480 | 1050 | 1512 | 1470 | 960 | 405 | 100 | 11 | 6144 | ||||
| 11 | 1 | 22 | 165 | 660 | 1650 | 2772 | 3234 | 2640 | 1485 | 550 | 121 | 12 | 13312 | |||
| 12 | 1 | 24 | 198 | 880 | 2475 | 4752 | 6468 | 6336 | 4455 | 2200 | 726 | 144 | 13 | 28672 | ||
| 13 | 1 | 26 | 234 | 1144 | 3575 | 7722 | 12012 | 13728 | 11583 | 7150 | 3146 | 936 | 169 | 14 | 61440 | |
| 14 | 1 | 28 | 273 | 1456 | 5005 | 12012 | 21021 | 27456 | 27027 | 20020 | 11011 | 4368 | 1183 | 196 | 15 | 131072 |
Generated b-files: A001792 (row sums), A093375 (reversed triangle), A103406 ?
Computed in your browser ?
Other recurrences
Image under θ+1 · proved equal for all n (proof step)
Reversed rows · proved equal for all n (proof step)
Technical details
- Node id
n-801895951aa20f29· revision 3 · JSON- Definition source
oeis_catalog:A103406derived:theta_plus_one(equal bye-2526de52e9387fae)derived:reversal(equal bye-3dbd2ded62b9bb8b)- Fingerprint
f6-e99c3824c4864868
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 ℚ): yes
OEIS matches ?
- A001792: row sums, 15 terms agree OEIS nonn (row-sum match: finite agreement only) · view at the OEIS · generated b-file 200 terms from n = 0 ?
- A103406: triangle, 69 terms agree, 11 full rows row-aligned (tabl) OEIS nonn in polytool catalog · view at the OEIS · generated b-file 990 terms from n = 0 ?
- A093375: reversed triangle, 67 terms agree, 11 full rows row-aligned (tabl) OEIS nonn in polytool catalog · view at the OEIS · generated b-file 990 terms from n = 1 ?
- A049610: row sums, 15 terms agree OEIS nonn (row-sum match: finite agreement only) · view at the OEIS
A generated b-file is computed from the node's definition along the row alignment of its OEIS link and starts at the entry's offset. It is an identification on a finite prefix, not a proof that the node and the OEIS entry agree for every n.
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.
The atlas's own rows: n = 0, …, ascending powers of x
%e ExampleTriangle begins: 0: 1; 1: 1, 2; 2: 1, 4, 3; 3: 1, 6, 9, 4; 4: 1, 8, 18, 16, 5; 5: 1, 10, 30, 40, 25, 6; 6: 1, 12, 45, 80, 75, 36, 7; 7: 1, 14, 63, 140, 175, 126, 49, 8; 8: 1, 16, 84, 224, 350, 336, 196, 64, 9; 9: 1, 18, 108, 336, 630, 756, 588, 288, 81, 10;
%F FormulaThe row polynomials satisfy P(n,x) = (2+2*x)*P(n-1,x) - (1+2*x+x^2)*P(n-2,x) for n >= 2, with P(0,x) = 1, P(1,x) = 1 + 2*x.
%t MathematicaP[0, x_] := 1;
P[1, x_] := 1 + 2*x;
P[n_, x_] := P[n, x] = Expand[(2 + 2*x)*P[n - 1, x] - (1 + 2*x + x^2)*P[n - 2, x]];
Table[CoefficientList[P[n, x], x], {n, 0, 9}] // Flatten%o PARI/GPP(N) = {my(v = vector(N + 1)); v[1] = 1; v[2] = 1 + 2*x; for(n = 2, N, my(i = n + 1); v[i] = (2 + 2*x)*v[i - 1] - (1 + 2*x + x^2)*v[i - 2]); 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)Generated with polytool c48629374597.
A093375: reversed triangle, offset 1, from the row P0
%e ExampleTriangle begins: 1: 1; 2: 2, 1; 3: 3, 4, 1; 4: 4, 9, 6, 1; 5: 5, 16, 18, 8, 1; 6: 6, 25, 40, 30, 10, 1; 7: 7, 36, 75, 80, 45, 12, 1; 8: 8, 49, 126, 175, 140, 63, 14, 1; 9: 9, 64, 196, 336, 350, 224, 84, 16, 1; 10: 10, 81, 288, 588, 756, 630, 336, 108, 18, 1;
%F FormulaRow n lists the coefficients of P(n,x) in descending powers of x, where P(n,x) = (2+2*x)*P(n-1,x) - (1+2*x+x^2)*P(n-2,x) for n >= 3, with P(1,x) = 1, P(2,x) = 1 + 2*x.
%t MathematicaP[1, x_] := 1;
P[2, x_] := 1 + 2*x;
P[n_, x_] := P[n, x] = Expand[(2 + 2*x)*P[n - 1, x] - (1 + 2*x + x^2)*P[n - 2, x]];
Table[Reverse[PadRight[CoefficientList[P[n, x], x], n]], {n, 1, 10}] // Flatten%o PARI/GPP(N) = {my(v = vector(N)); v[1] = 1; v[2] = 1 + 2*x; for(n = 3, N, my(i = n); v[i] = (2 + 2*x)*v[i - 1] - (1 + 2*x + x^2)*v[i - 2]); v};
row(n, p) = {my(c = Vecrev(p)); c = concat(c, vector(n - #c)); c = Vecrev(c); c};
T(N) = {my(v = P(N)); concat(vector(#v, i, row(i, v[i])))};
T(10)Generated with polytool c48629374597.
A103406: triangle, offset 0, from the row P0
%e ExampleTriangle begins: 0: 1; 1: 1, 2; 2: 1, 4, 3; 3: 1, 6, 9, 4; 4: 1, 8, 18, 16, 5; 5: 1, 10, 30, 40, 25, 6; 6: 1, 12, 45, 80, 75, 36, 7; 7: 1, 14, 63, 140, 175, 126, 49, 8; 8: 1, 16, 84, 224, 350, 336, 196, 64, 9; 9: 1, 18, 108, 336, 630, 756, 588, 288, 81, 10;
%F FormulaThe row polynomials satisfy P(n,x) = (2+2*x)*P(n-1,x) - (1+2*x+x^2)*P(n-2,x) for n >= 2, with P(0,x) = 1, P(1,x) = 1 + 2*x.
%t MathematicaP[0, x_] := 1;
P[1, x_] := 1 + 2*x;
P[n_, x_] := P[n, x] = Expand[(2 + 2*x)*P[n - 1, x] - (1 + 2*x + x^2)*P[n - 2, x]];
Table[PadRight[CoefficientList[P[n, x], x], n + 1], {n, 0, 9}] // Flatten%o PARI/GPP(N) = {my(v = vector(N + 1)); v[1] = 1; v[2] = 1 + 2*x; for(n = 2, N, my(i = n + 1); v[i] = (2 + 2*x)*v[i - 1] - (1 + 2*x + x^2)*v[i - 2]); v};
row(n, p) = {my(c = Vecrev(p)); c = concat(c, vector(n + 1 - #c)); c};
T(N) = {my(v = P(N)); concat(vector(#v, i, row(i - 1, v[i])))};
T(9)Generated with polytool c48629374597.
Proved properties
-
Real-rooted. Every row P_n, n >= 0, is zero or real-rooted with all roots in (-inf, 0]
Every polynomial of the sequence has only real roots, inside the stated interval.
-
Coefficients of one sign. Every nonzero row, n >= 0: the coefficients of P_n have one sign
All coefficients of each polynomial have the same sign (after x -> -x when alternating).derived from the cited steps by a stated theorem:
e-a112169ce63438eb,e-ac2f606fc4a9c0c0,e-becc9d9cf02b45c7 -
Hadamard kernel (Malo). Rows P_n, n >= 0, are Hadamard kernels (preserve real-rootedness)
Multiplying the coefficients of any real-rooted polynomial termwise by those of a row (Hadamard product) gives a real-rooted polynomial again.derived from the cited steps by a stated theorem:
e-a112169ce63438eb,e-ac2f606fc4a9c0c0,e-becc9d9cf02b45c7 -
Interlacing of consecutive rows. Rows P_(m-1), P_m interlace (weak, unordered) for every m >= 1
The roots of consecutive polynomials alternate on the real line.by the cited steps:
e-a112169ce63438eb,e-becc9d9cf02b45c7
Proof chains
5 incoming proof steps. Each step is expanded back to the roots of the proof DAG, the theorem checks of a node's own rule. JSON
-
definition_equality
e-2526de52e9387faeroot: theorem checkOre-algebra identity of two definitions polytool ?
-
definition_equality
e-3dbd2ded62b9bb8broot: theorem checkOre-algebra identity of two definitions polytool ?
-
preserver_image euler_polynomial, shift 0
e-a112169ce63438ebLaguerre’s theorem on multiplier sequences [SP14; Bra15, Thm 40] polytool ?
-
premise (sequence):
binomial
uses: Every row P_n, n >= 0, is zero or real-rooted with all roots in (-inf, 0]
-
theorem_check liu_wang
e-317eb4ccca60c6acroot: theorem checkLiu–Wang theorem [LW07, Thm 2.1, Cor 2.4] ?
-
preserver_image reversal, shift 0
e-033cf842aad3bdf7Affine, reversal and Möbius substitutions (elementary) polytool ?
-
premise (sequence):
binomial
uses: Every row P_n, n >= 0, is zero or real-rooted with all roots in (-inf, 0]
-
↵ theorem_check liu_wang
e-317eb4ccca60c6ac(expanded above) -
basis_map A007318
e-095baa2a9dfe188fGeometric basis map (elementary) polytool ?
-
premise (sequence):
a(n) = 1 if n is of the form m(m+3)/2, otherwise 0.
uses: Every row P_n, n >= 0, is zero or real-rooted with all roots in (-inf, 0]
-
veronese_section r=2, j=0, shift 0
e-292206ad4ae86d2dVeronese sections of Pólya frequency sequences [ASW52; Edr53; AW24] polytool ?
-
premise (sequence):
basis-map coefficients of Irregular triangle read by rows formed by taking every other row of the Delannoy array (A0082…
uses: Every row P_n, n >= 0, is zero or real-rooted with all roots in (-inf, 0]
-
parent_provers periodic_product
e-1991a148b9538b93root: theorem checkParent-row provers of polytool polytool ?
-
preserver_image mobius, shift 0
e-12303ad0e0b9f8ebAffine, reversal and Möbius substitutions (elementary) polytool ?
-
premise (sequence):
basis-map coefficients of Irregular triangle read by rows formed by taking every other row of the Delannoy array (A0082…
uses: Every row P_n, n >= 0, is zero or real-rooted with all roots in (-inf, 0]
-
premise (sequence):
basis-map coefficients of Irregular triangle read by rows formed by taking every other row of the Delannoy array (A0082…
-
-
premise (sequence):
basis-map coefficients of Irregular triangle read by rows formed by taking every other row of the Delannoy array (A0082…
-
parent_provers periodic_product
e-b7eca039219a7f53root: theorem checkParent-row provers of polytool polytool ?
-
preserver_image mobius, shift 0
e-7e44b66e04c497d6Affine, reversal and Möbius substitutions (elementary) polytool ?
-
premise (sequence):
a(n) = 1 if n is of the form m(m+3)/2, otherwise 0.
uses: Every row P_n, n >= 0, is zero or real-rooted with all roots in (-inf, 0]
-
premise (sequence):
a(n) = 1 if n is of the form m(m+3)/2, otherwise 0.
-
-
premise (sequence):
a(n) = 1 if n is of the form m(m+3)/2, otherwise 0.
-
parent_provers periodic_product
e-13983862a411779croot: theorem checkParent-row provers of polytool polytool ?
-
preserver_image reversal, shift 0
e-989c95e5c893fce9Affine, reversal and Möbius substitutions (elementary) polytool ?
-
premise (sequence):
Riordan array (1-x,x/(1-x)).
uses: Every row P_n, n >= 0, is zero or real-rooted with all roots in (-inf, 0]
-
parent_provers periodic_product
e-b62074ba45057e82root: theorem checkParent-row provers of polytool polytool ?
-
-
premise (sequence):
Riordan array (1-x,x/(1-x)).
-
gamma_polynomial scale 1, shift 0
e-c76adc2a99cd3deaGamma criterion RealRooted IsGammaExpansion.isRealRooted_and_hasRootsNonpos_iff ?
-
premise (sequence):
Triangle read by rows, T(n,k) = GegenbauerC(m,-n,-3/2) where m = k if k<n else 2*n-k, for n>=0 and 0<=k<=2n.
uses: Every row P_n, n >= 0, is zero or real-rooted with all roots in (-inf, 0]
-
preserver_image reversal, shift 0
e-0e058c428dfdfeb6Affine, reversal and Möbius substitutions (elementary) polytool ?
-
premise (sequence):
Triangle read by rows, T(n,k) = GegenbauerC(m,-n,-3/2) where m = k if k<n else 2*n-k, for n>=0 and 0<=k<=2n.
uses: Every row P_n, n >= 0, is zero or real-rooted with all roots in (-inf, 0]
-
parent_provers periodic_product
e-d82c9b48e4b31a6broot: theorem checkParent-row provers of polytool polytool ?
-
-
premise (sequence):
Triangle read by rows, T(n,k) = GegenbauerC(m,-n,-3/2) where m = k if k<n else 2*n-k, for n>=0 and 0<=k<=2n.
-
↵ parent_provers periodic_product
e-d82c9b48e4b31a6b(expanded above)
-
-
premise (sequence):
Triangle read by rows, T(n,k) = GegenbauerC(m,-n,-3/2) where m = k if k<n else 2*n-k, for n>=0 and 0<=k<=2n.
-
preserver_image affine, shift 0
e-3a52de7742ed9314Affine, reversal and Möbius substitutions (elementary) polytool ?
-
premise (sequence):
a(n) = 1 if n is of the form m(m+3)/2, otherwise 0.
uses: Every row P_n, n >= 0, is zero or real-rooted with all roots in (-inf, 0]
-
premise (sequence):
a(n) = 1 if n is of the form m(m+3)/2, otherwise 0.
-
↵ theorem_check liu_wang
-
premise (sequence):
binomial
-
↵ basis_map A007318
e-095baa2a9dfe188f(expanded above) -
↵ parent_provers periodic_product
e-13983862a411779c(expanded above) -
↵ preserver_image reversal, shift 0
e-989c95e5c893fce9(expanded above) -
↵ gamma_polynomial scale 1, shift 0
e-c76adc2a99cd3dea(expanded above) -
↵ preserver_image affine, shift 0
e-3a52de7742ed9314(expanded above)
-
-
premise (sequence):
binomial
-
preserver_image reversal, shift 0
e-ac2f606fc4a9c0c0Affine, reversal and Möbius substitutions (elementary) polytool ?
-
premise (sequence):
Array T(m,n) read by ascending antidiagonals: T(m,n) = m*binomial(n+m-2, n-1) for m, n >= 1.
uses: Every row P_n, n >= 0, is zero or real-rooted with all roots in (-inf, 0]
-
preserver_image derivative_polynomial, shift 0
e-5c65a1a70531c253Laguerre’s theorem on differential operators [SP14] polytool ?
-
premise (sequence):
binomial
uses: Every row P_n, n >= 0, is zero or real-rooted with all roots in (-inf, 0]
-
↵ theorem_check liu_wang
e-317eb4ccca60c6ac(expanded above) -
↵ preserver_image reversal, shift 0
e-033cf842aad3bdf7(expanded above) -
↵ basis_map A007318
e-095baa2a9dfe188f(expanded above) -
↵ parent_provers periodic_product
e-13983862a411779c(expanded above) -
↵ preserver_image reversal, shift 0
e-989c95e5c893fce9(expanded above) -
↵ gamma_polynomial scale 1, shift 0
e-c76adc2a99cd3dea(expanded above) -
↵ preserver_image affine, shift 0
e-3a52de7742ed9314(expanded above)
-
↵ theorem_check liu_wang
-
premise (sequence):
binomial
-
preserver_image reversal, shift 0
e-c2af0ae48d7f11c8Affine, reversal and Möbius substitutions (elementary) polytool ?
-
premise (sequence):
Triangle read by rows: n-th row = unsigned coefficients of the characteristic polynomials of an n X n matrix with 2's o…
uses: Every row P_n, n >= 0, is zero or real-rooted with all roots in (-inf, 0]
-
↵ preserver_image euler_polynomial, shift 0
e-a112169ce63438eb(expanded above) -
preserver_image euler_polynomial, shift 0
e-becc9d9cf02b45c7Laguerre’s theorem on multiplier sequences [SP14; Bra15, Thm 40] polytool ?
-
premise (sequence):
binomial
uses: Every row P_n, n >= 0, is zero or real-rooted with all roots in (-inf, 0]
-
↵ theorem_check liu_wang
e-317eb4ccca60c6ac(expanded above) -
↵ preserver_image reversal, shift 0
e-033cf842aad3bdf7(expanded above) -
↵ basis_map A007318
e-095baa2a9dfe188f(expanded above) -
↵ parent_provers periodic_product
e-13983862a411779c(expanded above) -
↵ preserver_image reversal, shift 0
e-989c95e5c893fce9(expanded above) -
↵ gamma_polynomial scale 1, shift 0
e-c76adc2a99cd3dea(expanded above) -
↵ preserver_image affine, shift 0
e-3a52de7742ed9314(expanded above)
-
↵ theorem_check liu_wang
-
premise (sequence):
binomial
-
↵ preserver_image euler_polynomial, shift 0
-
premise (sequence):
Triangle read by rows: n-th row = unsigned coefficients of the characteristic polynomials of an n X n matrix with 2's o…
-
preserver_image derivative_polynomial, shift 0
e-d2751b1a8ff0ac45Laguerre’s theorem on differential operators [SP14] polytool ?
-
premise (sequence):
binomial
uses: Every row P_n, n >= 0, is zero or real-rooted with all roots in (-inf, 0]
-
↵ theorem_check liu_wang
e-317eb4ccca60c6ac(expanded above) -
↵ preserver_image reversal, shift 0
e-033cf842aad3bdf7(expanded above) -
↵ basis_map A007318
e-095baa2a9dfe188f(expanded above) -
↵ parent_provers periodic_product
e-13983862a411779c(expanded above) -
↵ preserver_image reversal, shift 0
e-989c95e5c893fce9(expanded above) -
↵ gamma_polynomial scale 1, shift 0
e-c76adc2a99cd3dea(expanded above) -
↵ preserver_image affine, shift 0
e-3a52de7742ed9314(expanded above)
-
↵ theorem_check liu_wang
-
premise (sequence):
binomial
-
-
premise (sequence):
Array T(m,n) read by ascending antidiagonals: T(m,n) = m*binomial(n+m-2, n-1) for m, n >= 1.
-
↵ preserver_image euler_polynomial, shift 0
e-becc9d9cf02b45c7(expanded above)
Used in 9 proof steps
e-115dd329a45c57a0preserver_image (affine, shift 0) → shift_x_plus_one of Triangle read by rows: n-th row = unsigned coefficients of the charac…e-6a0b4c5cd2e9e2b4preserver_image (derivative_polynomial, shift 0) → one_plus_derivative of Triangle read by rows: n-th row = unsigned coefficients of the cha…e-85c5d5448d617672preserver_image (euler_polynomial, shift 0) → Triangle read by rows: T(n,k) = (k+1)^2*binomial(n,k) (0 <= k <= n).e-96faac355005fea1preserver_image (euler_polynomial, shift 0) → theta_plus_one_squared of Triangle read by rows: n-th row = unsigned coefficients of the …e-c2af0ae48d7f11c8preserver_image (reversal, shift 0) → Array T(m,n) read by ascending antidiagonals: T(m,n) = m*binomial(n+m-2, n-1) for m, n >=…e-c370da7cbe42718epreserver_image (derivative_polynomial, shift 0) → derivative of Triangle read by rows: n-th row = unsigned coefficients of the characterist…e-e68d80e1d57e8e24preserver_image (mobius, shift 0) → mobius_x_over_x_plus_one of Triangle read by rows: n-th row = unsigned coefficients of th…e-274dee66a9035b75veronese_section (r=2, j=1, shift 0) → veronese_section_2_1 of Triangle read by rows: n-th row = unsigned coefficients of the ch…e-f6f2e37f93f1c1e6veronese_section (r=2, j=0, shift 0) → veronese_section_2_0 of Triangle read by rows: n-th row = unsigned coefficients of the ch…
Related sequences ?
Same fingerprint (27)
Equal up to content, sign, reversal and a small index shift.
- A135838 → θ+1
- A007318 → ×(−(1/2)n+(1/2)n²) → D+1
- basis-map coefficients of A113214, shifted → θ+1
- A028326 → θ+1
- A344565 → D+1
- A178820 → D+1
- A028326 → D+1
- A007318 → ×(−(1/2)n+(1/2)n²) → θ+1
- A178821 → D+1
- A135838 → D+1
- A178822 → D+1
- basis-map coefficients of A113214, shifted → D+1
- A038208 → θ+1
- A007318 → ×(1+(3/2)n+(1/2)n²) → θ+1
- A178820 → θ+1
- A085880 → θ+1
- A178821 → θ+1
- A007318 → D → D+1
- A038208 → D+1
- A344565 → θ+1
- A007318 → ×(1+(3/2)n+(1/2)n²) → D+1
- A085880 → D+1
- A007318 → D → θ+1
- A007318 → ×(1+2n+n²) → D+1
- A093375: Array T(m,n) read by ascending antidiagonals: T(m,n) = m*binomial(n+m-2, n-1) for m, n >= 1.
- A007318 → ×(1+2n+n²) → θ+1
- A178822 → θ+1
References
- [ASW52] Michael Aissen, I. J. Schoenberg and A. M. Whitney. On the generating functions of totally positive sequences I. Journal d’Analyse Mathématique, 2(1):93–103, 1952. doi:10.1007/BF02786970
- [AW24] Christos A. Athanasiadis and David G. Wagner. Veronese sections and interlacing matrices of polynomials and formal power series. arXiv:2404.12989. arXiv
- [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
- [Edr53] Albert Edrei. Proof of a Conjecture of Schoenberg on the Generating Function of a Totally Positive Sequence. Canadian Journal of Mathematics, 5:86–94, 1953. doi:10.4153/CJM-1953-010-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
- [SP14] J Schur and G Polya. Über zwei Arten von Faktorenfolgen in der Theorie der algebraischen Gleichungen. J. Reine Angew. Math., 1914.