A101842: Triangle read by rows: T(n,k), for k=-n..n-1, is the scaled (by 2^n n!) probability that the sum of n uniform [-1, 1] variables is between k and k+1.

A101842 ↗ palindromicγ-positivenonnegative coefficients

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

Recurrence from A101842

\(P_{n}(x) = \left(1 + x + 2\,n x^{2}\right)\,P_{n-1}(x) + \left(x - x^{3}\right)\,P_{n-1}'(x)\)
\(P_{0}(x) = 1 + x\)

Rows ?

nx0x1x2x3x4x5x6x7x8x9x10x11x12x13x14x15Pn(1)
0112
113318
21716167148
31156111511561151384
4131206626105610566262063113840
51636592989755411774117747554298965963146080
6112720521330847349105099154624154624105099473491330820521271645120
712556297569352746778243311660957233750723375071660957824331274677569356297255110321920
Technical details
Node id
n-b3bdd177c3ec2fb9 · revision 0 · JSON
Definition source
oeis_catalog:A101842
Fingerprint
f6-0214d827b805d73a

Observed on the stored rows ?

  • Degrees: 1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29 (affine: slope 2, intercept 1)
  • Palindromic rows: yes
  • γ-positive (gamma_positive): yes
  • Nonnegative coefficients (nonn): yes
  • Only rational roots (rr_by_factors: every nonzero row is a product of linear factors over ℚ): no

Matches ?

  • A101842 reversed triangle 61 terms, 7 rows · tabf · nonn · catalog b-file?
  • A101842 triangle 61 terms, 7 rows · tabf · nonn · catalog b-file?

Row sums ? A000165, A052576

Submission text ?

Atlas rows n = 0, …, ascending powers of x
%e Example
Triangle begins:
  0: 1,    1;
  1: 1,    3,     3,      1;
  2: 1,    7,    16,     16,       7,        1;
  3: 1,   15,    61,    115,     115,       61,        15,         1;
  4: 1,   31,   206,    626,    1056,     1056,       626,       206,        31,         1;
  5: 1,   63,   659,   2989,    7554,    11774,     11774,      7554,      2989,       659,        63,         1;
  6: 1,  127,  2052,  13308,   47349,   105099,    154624,    154624,    105099,     47349,     13308,      2052,       127,         1;
  7: 1,  255,  6297,  56935,  274677,   824331,   1660957,   2337507,   2337507,   1660957,    824331,    274677,     56935,      6297,      255,       1;
  8: 1,  511, 19162, 237862, 1518478,  5960818,  15747154,  29428654,  39984640,  39984640,  29428654,  15747154,   5960818,   1518478,   237862,   19162,    511,     1;
  9: 1, 1023, 58015, 979297, 8136844, 40851316, 137449588, 328667020, 578256262, 763546234, 763546234, 578256262, 328667020, 137449588, 40851316, 8136844, 979297, 58015, 1023, 1;
%F Formula
The row polynomials satisfy P(n,x) = (1+x+2*n*x^2)*P(n-1,x) + (x-x^3)*d/dx P(n-1,x) for n >= 1, with P(0,x) = 1 + x.
%t Mathematica
P[0, x_] := 1 + x;
P[n_, x_] := P[n, x] = Expand[(1 + x + 2*n*x^2)*P[n - 1, x] + (x - x^3)*D[P[n - 1, x], x]];
Table[CoefficientList[P[n, x], x], {n, 0, 9}] // Flatten
%o PARI/GP
P(N) = {my(v = vector(N + 1)); v[1] = 1 + x; for(n = 1, N, my(i = n + 1); v[i] = (1 + x + 2*n*x^2)*v[i - 1] + (x - x^3)*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 - 1, v[i])))};
T(9)

polytool c48629374597

A101842 triangle, offset 1, from P0
%e Example
Triangle begins:
   1: 1,    1;
   2: 1,    3,     3,      1;
   3: 1,    7,    16,     16,       7,        1;
   4: 1,   15,    61,    115,     115,       61,        15,         1;
   5: 1,   31,   206,    626,    1056,     1056,       626,       206,        31,         1;
   6: 1,   63,   659,   2989,    7554,    11774,     11774,      7554,      2989,       659,        63,         1;
   7: 1,  127,  2052,  13308,   47349,   105099,    154624,    154624,    105099,     47349,     13308,      2052,       127,         1;
   8: 1,  255,  6297,  56935,  274677,   824331,   1660957,   2337507,   2337507,   1660957,    824331,    274677,     56935,      6297,      255,       1;
   9: 1,  511, 19162, 237862, 1518478,  5960818,  15747154,  29428654,  39984640,  39984640,  29428654,  15747154,   5960818,   1518478,   237862,   19162,    511,     1;
  10: 1, 1023, 58015, 979297, 8136844, 40851316, 137449588, 328667020, 578256262, 763546234, 763546234, 578256262, 328667020, 137449588, 40851316, 8136844, 979297, 58015, 1023, 1;
%F Formula
The row polynomials satisfy P(n,x) = (1+x+(2*n-2)*x^2)*P(n-1,x) + (x-x^3)*d/dx P(n-1,x) for n >= 2, with P(1,x) = 1 + x.
%t Mathematica
P[1, x_] := 1 + x;
P[n_, x_] := P[n, x] = Expand[(1 + x + (2*n - 2)*x^2)*P[n - 1, x] + (x - x^3)*D[P[n - 1, x], x]];
Table[CoefficientList[P[n, x], x], {n, 1, 10}] // Flatten
%o PARI/GP
P(N) = {my(v = vector(N)); v[1] = 1 + x; for(n = 2, N, my(i = n); v[i] = (1 + x + (2*n - 2)*x^2)*v[i - 1] + (x - x^3)*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)

polytool c48629374597

A101842 reversed triangle, offset 1, from P0
%e Example
Triangle begins:
   1: 1,    1;
   2: 1,    3,     3,      1;
   3: 1,    7,    16,     16,       7,        1;
   4: 1,   15,    61,    115,     115,       61,        15,         1;
   5: 1,   31,   206,    626,    1056,     1056,       626,       206,        31,         1;
   6: 1,   63,   659,   2989,    7554,    11774,     11774,      7554,      2989,       659,        63,         1;
   7: 1,  127,  2052,  13308,   47349,   105099,    154624,    154624,    105099,     47349,     13308,      2052,       127,         1;
   8: 1,  255,  6297,  56935,  274677,   824331,   1660957,   2337507,   2337507,   1660957,    824331,    274677,     56935,      6297,      255,       1;
   9: 1,  511, 19162, 237862, 1518478,  5960818,  15747154,  29428654,  39984640,  39984640,  29428654,  15747154,   5960818,   1518478,   237862,   19162,    511,     1;
  10: 1, 1023, 58015, 979297, 8136844, 40851316, 137449588, 328667020, 578256262, 763546234, 763546234, 578256262, 328667020, 137449588, 40851316, 8136844, 979297, 58015, 1023, 1;
%F Formula
Row n lists the coefficients of P(n,x) in descending powers of x, where P(n,x) = (1+x+(2*n-2)*x^2)*P(n-1,x) + (x-x^3)*d/dx P(n-1,x) for n >= 2, with P(1,x) = 1 + x.
%t Mathematica
P[1, x_] := 1 + x;
P[n_, x_] := P[n, x] = Expand[(1 + x + (2*n - 2)*x^2)*P[n - 1, x] + (x - x^3)*D[P[n - 1, x], x]];
Table[Reverse[CoefficientList[P[n, x], x]], {n, 1, 10}] // Flatten
%o PARI/GP
P(N) = {my(v = vector(N)); v[1] = 1 + x; for(n = 2, N, my(i = n); v[i] = (1 + x + (2*n - 2)*x^2)*v[i - 1] + (x - x^3)*deriv(v[i - 1], x)); v};
row(n, p) = {my(c = Vecrev(p)); c = Vecrev(c); c};
T(N) = {my(v = P(N)); concat(vector(#v, i, row(i, v[i])))};
T(10)

polytool c48629374597

Properties

Proof chains ? JSON

References

  1. [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