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-4 of 4 results.

A318146 Coefficients of the Omega polynomials of order 2, triangle T(n,k) read by rows with 0<=k<=n.

Original entry on oeis.org

1, 0, 1, 0, -2, 3, 0, 16, -30, 15, 0, -272, 588, -420, 105, 0, 7936, -18960, 16380, -6300, 945, 0, -353792, 911328, -893640, 429660, -103950, 10395, 0, 22368256, -61152000, 65825760, -36636600, 11351340, -1891890, 135135
Offset: 0

Views

Author

Peter Luschny, Aug 22 2018

Keywords

Comments

The name 'Omega polynomial' is not a standard name. The Omega numbers are the coefficients of the Omega polynomials, the associated Omega numbers are the weights of P(m, k) in the recurrence formula given below.
The signed Euler secant numbers appear as values at x=-1 and the signed Euler tangent numbers as the coefficients of x.

Examples

			Row n in the triangle below is the coefficient list of OmegaPolynomial(2, n). For other cases than m = 2 see the cross-references.
[0] [1]
[1] [0,        1]
[2] [0,       -2,         3]
[3] [0,       16,       -30,       15]
[4] [0,     -272,       588,     -420,       105]
[5] [0,     7936,    -18960,    16380,     -6300,      945]
[6] [0,  -353792,    911328,  -893640,    429660,  -103950,    10395]
[7] [0, 22368256, -61152000, 65825760, -36636600, 11351340, -1891890, 135135]

Links

Crossrefs

Variant is A088874 (unsigned).
T(n,1) = A000182(n), T(n,n) = A001147(n).
All row sums are 1, alternating row sums are A028296 (A000364).
A023531 (m=1), this seq (m=2), A318147 (m=3), A318148 (m=4).
Associated Omega numbers: A318254 (m=2), A318255 (m=3).
Coefficients of x for Omega polynomials of all orders are in A318253.

Programs

  • Maple
    OmegaPolynomial := proc(m, n) local Omega;
Omega := m -> hypergeom([], [seq(i/m, i=1..m-1)], (z/m)^m):
series(Omega(m)^x, z, m*(n+1)):
sort(expand((m*n)!*coeff(%, z, n*m)), [x], ascending) end:
CL := p -> PolynomialTools:-CoefficientList(p, x):
FL := p -> ListTools:-Flatten(p):
FL([seq(CL(OmegaPolynomial(2, n)), n=0..8)]);
# Alternative:
ser := series(sech(z)^(-x), z, 24): row := n -> n!*coeff(ser, z, n):
seq(seq(coeff(row(2*n), x, k), k=0..n), n=0..6); # Peter Luschny, Jul 01 2019
  • Mathematica
    OmegaPolynomial[m_,n_] :=  Module [{ },
S = Series[MittagLefflerE[m,z]^x, {z,0,10}];
Expand[(m n)! Coefficient[S,z,n]] ]
Table[CoefficientList[OmegaPolynomial[2,n],x], {n,0,7}] // Flatten
(* Second program: *)
T[n_, k_] := (2n)! SeriesCoefficient[Sech[z]^-x, {z, 0, 2n}, {x, 0, k}];
Table[T[n, k], {n, 0, 7}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jul 23 2019, after Peter Luschny *)
  • Sage
    def OmegaPolynomial(m, n):
    R = ZZ[x]; z = var('z')
    f = [i/m for i in (1..m-1)]
    h = lambda z: hypergeometric([], f, (z/m)^m)
    return R(factorial(m*n)*taylor(h(z)^x, z, 0, m*n + 1).coefficient(z, m*n))
[list(OmegaPolynomial(2, n)) for n in (0..6)]
# Recursion over the polynomials, returns a list of the first len polynomials:
def OmegaPolynomials(m, len, coeffs=true):
    R = ZZ[x]; B = [0]*len; L = [R(1)]*len
    for k in (1..len-1):
        s = x*sum(binomial(m*k-1, m*(k-j))*B[j]*L[k-j] for j in (1..k-1))
        B[k] = c = 1 - s.subs(x=1)
        L[k] = R(expand(s + c*x))
    return [list(l) for l in L] if coeffs else L
print(OmegaPolynomials(2, 6))

Formula

OmegaPolynomial(m, n) = (m*n)! [z^n] E(m, z)^x where E(m, z) is the Mittag-Leffler function.
OmegaPolynomial(m, n) = (m*n)!*[z^(n*m)] H(m, z)^x where H(m, z) = hypergeom([], [seq(i/m, i=1..m-1)], (z/m)^m).
The Omega polynomials can be computed by the recurrence P(m, 0) = 1 and for n >= 1 P(m, n) = x * Sum_{k=0..n-1} binomial(m*n-1, m*k)*T(m, n-k)*P(m, k) where T(m, n) are the generalized tangent numbers A318253. A separate computation of the T(m, n) can be avoided, see the Sage implementation below for the details.
T(n, k) = [x^k] (2*n)! [z^(2*n)] sech(z)^(-x). - Peter Luschny, Jul 01 2019

A318254 Associated Omega numbers of order 2, triangle T(n,k) read by rows for n >= 0 and 0 <= k <= n.

Original entry on oeis.org

1, 1, 1, 1, 3, -2, 1, 5, -20, 16, 1, 7, -70, 336, -272, 1, 9, -168, 2016, -9792, 7936, 1, 11, -330, 7392, -89760, 436480, -353792, 1, 13, -572, 20592, -466752, 5674240, -27595776, 22368256, 1, 15, -910, 48048, -1750320, 39719680, -482926080, 2348666880, -1903757312
Offset: 0

Views

Author

Peter Luschny, Aug 26 2018

Keywords

Comments

The Omega polynomials A318146 are defined by the recurrence P(m, 0) = 1 and for n>=1 P(m, n) = x * Sum_{k=0..n-1} binomial(m*n-1, m*k)*t(m, n-k)*P(m, k) where t(m, n) are the generalized tangent numbers A318253. The Omega numbers are the coefficients of the Omega polynomials. The associated Omega numbers are the weights of P(m, k) in the recurrence formula.

Examples

			Triangle starts:
[0] [1]
[1] [1,  1]
[2] [1,  3,   -2]
[3] [1,  5,  -20,    16]
[4] [1,  7,  -70,   336,    -272]
[5] [1,  9, -168,  2016,   -9792,    7936]
[6] [1, 11, -330,  7392,  -89760,  436480,   -353792]
[7] [1, 13, -572, 20592, -466752, 5674240, -27595776, 22368256]

Crossrefs

Even-indexed rows of A220901 (up to signs).
T(n, 0) = A005408, T(n, n) = A220901 (up to signs), row sums are A040000.
Cf. A318146, A318253, A318255 (m=3).

Programs

  • Maple
    # The function TNum is defined in A318253.
T := (m, n, k) -> `if`(k=0, 1, binomial(m*n-1, m*(n-k))*TNum(m, k)):
for n from 0 to 6 do seq(T(2, n, k), k=0..n) od;
  • Sage
    def AssociatedOmegaNumberTriangle(m, len):
    R = ZZ[x]; B = [1]*len; L = [R(1)]*len; T = [[1]]
    for k in (1..len-1):
        s = x*sum(binomial(m*k-1, m*(k-j))*B[j]*L[k-j] for j in (1..k-1))
        B[k] = c = 1 - s.subs(x=1); L[k] = R(expand(s + c*x))
        T.append([1] + [binomial(m*k-1, m*(k-j))*B[j] for j in (1..k)])
    return T
A318254Triangle = lambda dim: AssociatedOmegaNumberTriangle(2, dim)
print(A318254Triangle(8))

Formula

T(m, n, k) = binomial(m*n-1, m*(n-k))*A318253(m, k) for k>0 and 1 for k=0. We consider here the case m=2.

A318258 a(n) = [x] OmegaPolynomial(5, n). OmegaPolynomials are defined in A318146.

Original entry on oeis.org

0, 1, -125, 249250, -2886735625, 122209131374375, -14455143383196875000, 4006210678487307667578125, -2297417123000769120910212890625, 2485076260705905645263720799941406250, -4719878705811419698488114573981055908203125
Offset: 0

Views

Author

Peter Luschny, Aug 22 2018

Keywords

Crossrefs

Cf. A318253 (case n=5), A318146.
Cf. A000182 (m=2), A293951 (m=3), A273352 (m=4), this seq (m=5).

Programs

  • Maple
    # The function OmegaPolynomial is defined in A318146.
seq(coeff(OmegaPolynomial(5, n), x, 1), n=0..11);
  • Mathematica
    LMlist[m_, len_] := Table[(m n)!, {n, 0, len}]*
CoefficientList[Series[Log[MittagLefflerE[m, z]], {z, 0, len}], z];
LMlist[5, 13]

A318255 Associated Omega numbers of order 3, triangle T(n,k) read by rows for n >= 0 and 0 <= k <= n.

Original entry on oeis.org

1, 1, 1, 1, 10, -9, 1, 28, -504, 477, 1, 55, -4158, 78705, -74601, 1, 91, -18018, 1432431, -27154764, 25740261, 1, 136, -55692, 11595870, -923261976, 17503377480, -16591655817, 1, 190, -139536, 60087690, -12529983960, 997692516360, -18914487631380, 17929265150637
Offset: 0

Views

Author

Peter Luschny, Aug 26 2018

Keywords

Comments

See the comments in A318254.

Examples

			Triangle starts:
[0] 1
[1] 1,   1
[2] 1,  10,     -9
[3] 1,  28,   -504,      477
[4] 1,  55,  -4158,    78705,     -74601
[5] 1,  91, -18018,  1432431,  -27154764,    25740261
[6] 1, 136, -55692, 11595870, -923261976, 17503377480, -16591655817

Crossrefs

T(n, 0) = A060544, T(n, n) = A293951(n+1) (up to signs), row sums are A040000.
Cf. A318146, A318253, A318254 (m=2).

Programs

  • Maple
    # The function TNum is defined in A318253.
T := (m, n, k) -> `if`(k=0, 1, binomial(m*n-1, m*(n-k))*TNum(m, k)):
for n from 0 to 6 do seq(T(3, n, k), k=0..n) od;
  • Sage
    # uses[AssociatedOmegaNumberTriangle from A318254]
A318255Triangle = lambda dim: AssociatedOmegaNumberTriangle(3, dim)
print(A318255Triangle(8))

Formula

T(m, n, k) = binomial(m*n-1, m*(n-k))*A318253(m, k) for k>0 and 1 for k=0. We consider here the case m=3.
Showing 1-4 of 4 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.