cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-5 of 5 results.

A291447 Triangle read by rows, numerators of coefficients (in rising powers) of rational polynomials P(n, x) such that Integral_{x=0..1} P'(n, x) = BernoulliMedian(n).

Original entry on oeis.org

0, 1, 0, 0, 0, 1, 0, 0, 0, 1, -1, 4, 0, 0, 0, 1, -3, 48, -12, 36, 0, 0, 0, 1, -7, 268, -176, 1968, -216, 64, 0, 0, 0, 1, -15, 240, -1580, 37140, -9900, 10400, -5760, 14400, 0, 0, 0, 1, -31, 4924, -11680, 488640, -238680, 496320, -639360, 5486400, -216000, 518400
Offset: 0

Views

Author

Peter Luschny, Aug 24 2017

Keywords

Comments

The Bernoulli median numbers are A212196/A181131. See A290694 for further comments.

Examples

			Triangle starts:
[0, 1]
[0, 0, 0, 1]
[0, 0, 0, 1, -1, 4]
[0, 0, 0, 1, -3, 48, -12, 36]
[0, 0, 0, 1, -7, 268, -176, 1968, -216, 64]
[0, 0, 0, 1, -15, 240, -1580, 37140, -9900, 10400, -5760, 14400]
The first few polynomials are:
P_0(x) = x.
P_1(x) = (1/3)*x^3.
P_2(x) = (4/5)*x^5 - x^4 + (1/3)*x^3.
P_3(x) = (36/7)*x^7 - 12*x^6 + (48/5)*x^5 - 3*x^4 + (1/3)*x^3.
P_4(x) = 64*x^9 - 216*x^8 + (1968/7)*x^7 - 176*x^6 + (268/5)*x^5 - 7*x^4 +(1/3)*x^3.
Evaluated at x = 1 this gives a decomposition of the Bernoulli median numbers:
BM(0) = 1     =    1.
BM(1) = 1/3   =  1/3.
BM(2) = 2/15  =  4/5 -   1 +    1/3.
BM(3) = 8/105 = 36/7 -  12 +   48/5 -   3 +   1/3.
BM(4) = 8/105 =   64 - 216 + 1968/7 - 176 + 268/5 - 7 + 1/3.

Links

Crossrefs

Cf. A164555/A027642, A212196/A181131, A291449/A291450, A290694/A290695, A291447/A291448.

Programs

  • Maple
    # The function BG_row is defined in A290694.
seq(BG_row(2, n, "num", "val"), n=0..12);        # A212196
seq(BG_row(2, n, "den", "val"), n=0..12);        # A181131
seq(print(BG_row(2, n, "num", "poly")), n=0..7); # A291447 (this seq.)
seq(print(BG_row(2, n, "den", "poly")), n=0..9); # A291448
  • Mathematica
    T[n_] := Integrate[Sum[(-1)^(n-j+1) StirlingS2[n, j] j! x^j, {j, 0, n}]^2, x];
Trow[n_] := CoefficientList[T[n], x] // Numerator;
Table[Trow[r], {r, 0, 6}] // Flatten

Formula

T(n,k) = Numerator([x^k] Integral(Sum_{j=0..n}(-1)^(n-j)*Stirling2(n,j)*j!*x^j)^m) for m = 2, n >= 0 and k = 0..m*n+1.

A291448 Triangle read by rows, denominators of coefficients (in rising powers) of rational polynomials P(n,x) such that Integral_{x=0..1} P'(n,x) = BernoulliMedian(n).

Original entry on oeis.org

1, 1, 1, 1, 1, 3, 1, 1, 1, 3, 1, 5, 1, 1, 1, 3, 1, 5, 1, 7, 1, 1, 1, 3, 1, 5, 1, 7, 1, 1, 1, 1, 1, 3, 1, 1, 1, 7, 1, 1, 1, 11, 1, 1, 1, 3, 1, 5, 1, 7, 1, 1, 1, 11, 1, 13, 1, 1, 1, 3, 1, 5, 1, 1, 1, 1, 1, 11, 1, 13, 1, 1, 1, 1, 1, 3, 1, 5, 1, 1, 1, 1, 1, 11, 1, 13, 1, 1
Offset: 0

Views

Author

Peter Luschny, Aug 24 2017

Keywords

Comments

See A291447 and A290694 for comments.

Examples

			Triangle starts:
[1, 1]
[1, 1, 1, 3]
[1, 1, 1, 3, 1, 5]
[1, 1, 1, 3, 1, 5, 1, 7]
[1, 1, 1, 3, 1, 5, 1, 7, 1, 1]
[1, 1, 1, 3, 1, 1, 1, 7, 1, 1, 1, 11]
[1, 1, 1, 3, 1, 5, 1, 7, 1, 1, 1, 11, 1, 13]

Crossrefs

Cf. A164555/A027642, A212196/A181131, A291449/A291450, A290694/A290695, A291447/A291448.

Programs

  • Maple
    # See A291447.
  • Mathematica
    T[n_] := Integrate[Sum[(-1)^(n-j+1) StirlingS2[n, j] j! x^j, {j,0,n}]^2, x];
Trow[n_] := CoefficientList[T[n], x] // Denominator;
Table[Trow[r], {r, 0, 7}] // Flatten

Formula

T(n,k) = Denominator([x^k] Integral(Sum_{j=0..n}(-1)^(n-j)*Stirling2(n,j)*j!*x^j)^m) for m = 2, n >= 0 and k = 0..m*n+1.

A290695 Triangle read by rows, denominators of coefficients (in rising powers) of rational polynomials P(n, x) such that Integral_{x=0..1} P'(n, x) = Bernoulli(n, 1).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 2, 3, 1, 1, 2, 1, 2, 1, 1, 2, 3, 1, 5, 1, 1, 2, 1, 2, 1, 1, 1, 1, 2, 3, 1, 1, 1, 7, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 2, 3, 1, 5, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 2, 3, 1, 1, 1, 1, 1, 1, 1, 11, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1
Offset: 0

Views

Author

Peter Luschny, Aug 24 2017

Keywords

Comments

See A290694 for comments.

Examples

			Triangle starts:
[1, 1]
[1, 1, 2]
[1, 1, 2, 3]
[1, 1, 2, 1, 2]
[1, 1, 2, 3, 1, 5]
[1, 1, 2, 1, 2, 1, 1]
[1, 1, 2, 3, 1, 1, 1, 7]
[1, 1, 2, 1, 2, 1, 1, 1, 1]

Crossrefs

Cf. A164555/A027642, A212196/A181131, A291449/A291450, A290694/A290695, A291447/A291448.

Programs

  • Maple
    T_row := n -> denom(PolynomialTools:-CoefficientList(add((-1)^(n-j+1)*Stirling2(n, j-1)*(j-1)!*x^j/j, j=1..n+1), x)): for n from 0 to 7 do T_row(n) od;
  • Mathematica
    T[n_] := Denominator[CoefficientList[Sum[(-1)^(n-j+1) StirlingS2[n, j-1] (j-1)! x^j/j, {j, 1, n+1}], x]];
Table[T[n], {n, 0, 7}] (* Jean-François Alcover, Jun 15 2019, from Maple *)

Formula

T(n, k) = Denominator([x^k] Integral (Sum_{j=0..n} (-1)^(n-j)*Stirling2(n,j)*j!* x^j)^m) for m = 1 and k = 0..n+1.

A291449 Numerators of Integral_{x=0..1} P(n, x)^3 with P(n, x) = Sum_{k=0..n} (-1)^(n-k)* Stirling2(n, k)*k!*x^k.

Original entry on oeis.org

1, 1, 13, 1, 43, -61, 728877, 81739, -1779449713, -2112052153, 730622680308569, 113221320488699, -3660430816956396309, -3021604582205161, 21842539561810574341396283, 66747470298418575790593659, -124586733960451680357554181608419, -28471605423890788373026535240299
Offset: 0

Views

Author

Peter Luschny, Aug 24 2017

Keywords

Comments

Consider a family of integrals I(m, n) = Integral_{x=0..1} P(n, x)^m with P(n, x) = Sum_{k=0..n} (-1)^(n-k)*Stirling2(n, k)*k!*x^k. I(1, n) are the Bernoulli numbers A164555/A027642, I(2, n) are the Bernoulli median numbers A212196/A181131, I(3, n) are the numbers A291449/A291450. The coefficients of the polynomials P(n, x)^m are for m = 1 A290694/A290695, for m = 2 A291447/A291448. (See A290694 for further comments.)

Crossrefs

Cf. A164555/A027642, A212196/A181131, A291449/A291450, A290694/A290695, A291447/A291448.

Programs

  • Maple
    # Function BG_row is defined in A290694.
seq(BG_row(3, n, "num", "val"), n=0..17);
  • Mathematica
    P[n_, x_] := Sum[(-1)^(n-k)*StirlingS2[n, k]*k!*x^k, {k, 0, n}];
a[n_] := Integrate[P[n, x]^3, {x, 0, 1}] // Numerator;
Table[a[n], {n, 0, 17}] (* Jean-François Alcover, Jun 15 2019 *)

A291450 Denominators of Integral_{x=0..1} P(n, x)^3 with P(n, x) = Sum_{k=0..n}(-1)^(n-k)* Stirling2(n, k)*k!*x^k.

Original entry on oeis.org

1, 4, 140, 28, 20020, 4004, 6466460, 184756, 148728580, 29745716, 133706993420, 2431036244, 449741705140, 31885268, 670910837521540, 134182167504308, 409926521725660940, 4822664961478364, 1278006214791766460, 1921813856829724, 242081282475556183660, 4401477863191930612
Offset: 0

Views

Author

Peter Luschny, Aug 24 2017

Keywords

Comments

See A291449 and A290694 for comments.

Crossrefs

Cf. A164555/A027642, A212196/A181131, A291449/A291450, A290694/A290695, A291447/A291448.

Programs

  • Maple
    # Function BG_row is defined in A290694.
seq(BG_row(3, n, "den", "val"), n=0..20);
  • Mathematica
    P[n_, x_] := Sum[(-1)^(n-k)*StirlingS2[n, k]*k!*x^k, {k, 0, n}];
a[n_] := Integrate[P[n, x]^3, {x, 0, 1}] // Denominator;
Table[a[n], {n, 0, 21}] (* Jean-François Alcover, Jun 15 2019 *)
Showing 1-5 of 5 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.