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

A123125 Triangle of Eulerian numbers T(n,k), 0 <= k <= n, read by rows.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 4, 1, 0, 1, 11, 11, 1, 0, 1, 26, 66, 26, 1, 0, 1, 57, 302, 302, 57, 1, 0, 1, 120, 1191, 2416, 1191, 120, 1, 0, 1, 247, 4293, 15619, 15619, 4293, 247, 1, 0, 1, 502, 14608, 88234, 156190, 88234, 14608, 502, 1, 0, 1, 1013, 47840, 455192, 1310354, 1310354, 455192, 47840, 1013, 1
Offset: 0

Views

Author

Philippe Deléham, Sep 30 2006

Keywords

Comments

The beginning of this sequence does not quite agree with the usual version, which is A173018. - N. J. A. Sloane, Nov 21 2010
Each row of A123125 is the reverse of the corresponding row in A173018. - Michael Somos, Mar 17 2011
A008292 (subtriangle for k>=1 and n>=1) is the main entry for these numbers.
Triangle T(n,k), 0 <= k <= n, read by rows given by [0,1,0,2,0,3,0,4,0,5,0,...] DELTA [1,0,2,0,3,0,4,0,5,0,6,...] where DELTA is the operator defined in A084938.
Row sums are the factorials. - Roger L. Bagula and Gary W. Adamson, Aug 14 2008
If the initial zero column is deleted, the result is A008292. - Roger L. Bagula and Gary W. Adamson, Aug 14 2008
This result gives an alternative method of calculating the Eulerian numbers by an Umbral Calculus expansion from Comtet. - Roger L. Bagula, Nov 21 2009
This function seems to be equivalent to the PolyLog expansion. - Roger L. Bagula, Nov 21 2009
A raising operator formed from the e.g.f. of this entry is the generator of a sequence of polynomials p(n,x;t) defined in A046802 that specialize to those for A119879 as p(n,x;-1), A007318 as p(n,x;0), A073107 as p(n,x;1), and A046802 as p(n,0;t). See Copeland link for more associations. - Tom Copeland, Oct 20 2015
The Eulerian numbers in this setup count the permutation trees of power n and width k (see the Luschny link). For the associated combinatorial statistic over permutations see the Sage program below and the example section. - Peter Luschny, Dec 09 2015 [See Elder et al. link. Peter Luschny, Jul 13 2022]
From Wolfdieter Lang, Apr 03 2017: (Start)
The row polynomials R(n, x) = Sum_{k=0..n} T(n, k)*x^k are the numerator polynomials of the o.g.f. G(n, x) of n-powers {m^n}_{m>=0} (with 0^0 = 1): G(n, x) = R(n, x)/(1-x)^(n+1). See the Aug 14 2008 formula, where f(x,n) = R(n, x). The e.g.f. of R(n, t) is given in Copeland's Oct 14 2015 formula below.
The first nine column sequences are A000007, A000012, A000295, A000460, A000498, A000505, A000514, A001243, A001244. (End)
With all offsets 0, let A_n(x;y) = (y + E.(x))^n, an Appell sequence in y where E.(x)^k = E_k(x) are the Eulerian polynomials of this entry, A123125. Then the row polynomials of A046802 (the h-polynomials of the stellahedra) are given by h_n(x) = A_n(x;1); the row polynomials of A248727 (the face polynomials of the stellahedra), by f_n(x) = A_n(1 + x;1); the Swiss-knife polynomials of A119879, by Sw_n(x) = A_n(-1;1 + x); and the row polynomials of the Worpitsky triangle (A130850), by w_n(x) = A(1 + x;0). Other specializations of A_n(x;y) give A090582 (the f-polynomials of the permutohedra, cf. also A019538) and A028246 (another version of the Worpitsky triangle). - Tom Copeland, Jan 24 2020
Let b(n) = (1/(n+1))*Sum_{k=0..n-1} (-1)^(n-k+1)*T(n, k+1) / binomial(n, k+1). Then b(n) = Bernoulli(n, 1) = -n*Zeta(1 - n) = Integral_{x=0..1} F_n(x) for n >= 1. Here F_n(x) are the signed Fubini polynomials (A278075). (See also Rzadkowski and Urlinska, example 1.) - Peter Luschny, Feb 15 2021
Patrick J. Burchell (see link) describes the following method: To get the k-th row of the triangle write the nonnegative integers with a fixed exponent k as a sequence, 0^k, 1^k, 2^k, ..., and then apply the first differences to them k + 1 times. - Peter Luschny, Apr 02 2023

Examples

			The triangle T(n, k) begins:
  n\k 0 1    2     3      4       5       6      7     8    9 10...
  0:  1
  1:  0 1
  2:  0 1    1
  3:  0 1    4     1
  4:  0 1   11    11      1
  5:  0 1   26    66     26       1
  6:  0 1   57   302    302      57       1
  7:  0 1  120  1191   2416    1191     120      1
  8:  0 1  247  4293  15619   15619    4293    247     1
  9:  0 1  502 14608  88234  156190   88234  14608   502    1
 10:  0 1 1013 47840 455192 1310354 1310354 455192 47840 1013  1
...  Reformatted. - _Wolfdieter Lang_, Feb 14 2015
------------------------------------------------------------------
The width statistic over permutations, n=4.
  [1, 2, 3, 4] => 3; [1, 2, 4, 3] => 2; [1, 3, 2, 4] => 2; [1, 3, 4, 2] => 2;
  [1, 4, 2, 3] => 2; [1, 4, 3, 2] => 1; [2, 1, 3, 4] => 3; [2, 1, 4, 3] => 2;
  [2, 3, 1, 4] => 2; [2, 3, 4, 1] => 3; [2, 4, 1, 3] => 2; [2, 4, 3, 1] => 2;
  [3, 1, 2, 4] => 3; [3, 1, 4, 2] => 3; [3, 2, 1, 4] => 2; [3, 2, 4, 1] => 3;
  [3, 4, 1, 2] => 3; [3, 4, 2, 1] => 2; [4, 1, 2, 3] => 4; [4, 1, 3, 2] => 3;
  [4, 2, 1, 3] => 3; [4, 2, 3, 1] => 3; [4, 3, 1, 2] => 3; [4, 3, 2, 1] => 2;
Gives row(4) = [0, 1, 11, 11, 1]. - _Peter Luschny_, Dec 09 2015
------------------------------------------------------------------
From _Wolfdieter Lang_, Apr 03 2017: (Start)
Recurrence: T(5, 3) = (6-3)*T(4, 2) + 3*T(4, 3) = 3*11 + 3*11 = 66.
O.g.f. column k=2: (x/(1 - 2*x))*E_x*(x/(1-x)) = (x/(1-x))^2/(1-2*x).
E.g.f. column k=2: A(2, x) = x*A(1, x) + x*E(1, x) = x*1 + x*(exp(x)-1) = x*exp(x), hence E(2, x) = (1 + int(x*exp(-x),x ))*exp(2*x) = exp(x)*(exp(x) - (1+x)). See A000295. (End)
		

References

  • L. Comtet, Advanced Combinatorics, Reidel, Holland, 1978, page 245. [Roger L. Bagula, Nov 21 2009]
  • Ronald L. Graham, Donald E. Knuth and Oren Patashnik, Concrete Mathematics, 2nd ed.; Addison-Wesley, 1994, p. 268, Row reversed table 268. - Wolfdieter Lang, Apr 03 2017
  • Douglas C. Montgomery and Lynwood A. Johnson, Forecasting and Time Series Analysis, MaGraw-Hill, New York, 1976, page 91. - Roger L. Bagula and Gary W. Adamson, Aug 14 2008

Crossrefs

See A008292 (subtriangle for k>=1 and n>=1), which is the main entry for these numbers. Another version has the zeros at the ends of the rows, as in Concrete Mathematics: see A173018.
T(2n,n) gives A180056.

Programs

  • Haskell
    a123125 n k = a123125_tabl !! n !! k
    a123125_row n = a123125_tabl !! n
    a123125_tabl = [1] : zipWith (:) [0, 0 ..] a008292_tabl
    -- Reinhard Zumkeller, Nov 06 2013
    
  • Maple
    gf := 1/(1 - t*exp(x)): ser := series(gf, x, 12):
    cx := n -> (-1)^(n + 1)*factor(n!*coeff(ser, x, n)*(t - 1)^(n + 1)):
    seq(print(seq(coeff(cx(n), t, k), k = 0..n)), n = 0..9); # Peter Luschny, Feb 11 2021
    A123125 := proc(n, k) option remember; if k = n then 1 elif k <= 0 or k > n then 0 else k*procname(n-1, k) + (n-k+1)*procname(n-1, k-1) fi end:
    seq(print(seq(A123125(n, k), k=0..n)), n=0..10); # Peter Luschny, Mar 28 2021
    # Alternative (Patrick J. Burchell):
    t := a -> Statistics:-Difference([0, a]): Trow := k -> (t@@(k+1))([seq(n^k, n = 0..k)]):
    seq(print(Trow(n)), n = 0..6); # Peter Luschny, Apr 02 2023
  • Mathematica
    f[x_, n_] := f[x, n] = (1 - x)^(n + 1)*Sum[k^n*x^k, {k, 0, Infinity}];
    Table[CoefficientList[f[x, n], x], {n,0,9}] // Flatten (* Roger L. Bagula, Aug 14 2008 *)
    t[n_ /; n >= 0, 0] = 1; t[n_, k_] /; k<0 || k>n = 0; t[n_, k_] := t[n, k] = (n-k) t[n-1, k-1] + (k+1) t[n-1, k]; T[n_, k_] := t[n, n-k];
    Table[T[n, k], {n,0,10}, {k, 0, n}] // Flatten (* Jean-François Alcover, May 26 2019 *)
    A123125[n_, k_] := Sum[(-1)^j*(n - j - k + 1)^n * Binomial[n + 1, j], {j, 0, n - k}];
    Table[A123125[n, k], {n, 0, 9}, {k, 0, n}] // TableForm  (* Peter Luschny, Aug 12 2022 *)
  • Python
    from math import isqrt, comb
    def A123125(n):
        a = (m:=isqrt(k:=n+1<<1))+(k>m*(m+1))
        b = comb(a+1,2)-n
        return sum(-(b-j)**(a-1)*comb(a,j) if j&1 else (b-j)**(a-1)*comb(a,j) for j in range(b)) # Chai Wah Wu, Nov 13 2024
  • Sage
    def statistic_eulerian(pi):
        if not pi: return 0
        h, i, branch, next = 0, len(pi), [0], pi[0]
        while True:
            while next < branch[len(branch)-1]:
                del(branch[len(branch)-1])
            current = 0
            h += 1
            while next > current:
                i -= 1
                if i == 0: return h
                branch.append(next)
                current, next = next, pi[i]
    def A123125_row(n):
        L = [0]*(n+1)
        for p in Permutations(n):
            L[statistic_eulerian(p)] += 1
        return L
    [A123125_row(n) for n in range(7)] # Peter Luschny, Dec 09 2015
    

Formula

Sum_{k=0..n} T(n,k) = n! = A000142(n).
Sum_{k=0..n} 2^k*T(n,k) = A000629(n).
Sum_{k=0..n} 3^k*T(n,k) = abs(A009362(n+1)).
Sum_{k=0..n} 2^(n-k)*T(n,k) = A000670(n).
Sum_{k=0..n} T(n,k)*3^(n-k) = A122704(n). - Philippe Deléham, Nov 07 2007
G.f.: f(x,n) = (1 - x)^(n + 1)*Sum_{k>=0} k^n*x^k. - Roger L. Bagula and Gary W. Adamson, Aug 14 2008. f is not the g.f. of the triangle, it is the polynomial of row n. See an Apr 03 2017 comment above - Wolfdieter Lang, Apr 03 2017
Sum_{k=0..n} T(n,k)*x^k = A000007(n), A000142(n), A000629(n), A123227(n), A201355(n), A201368(n) for x = 0, 1, 2, 3, 4, 5 respectively. - Philippe Deléham, Dec 01 2011
E.g.f. (1-t)/(1-t*exp((1-t)x)). A123125 * A007318 = A130850 = unsigned A075263, related to reversed A028246. A007318 * A123125 = A046802. Evaluating the row polynomials at -1, giving the alternating-sign row sum, generates A009006. - Tom Copeland, Oct 14 2015
From Wolfdieter Lang, Apr 03 2017: (Start)
T(n, k) = A173018(n, n-k), 0 <= k <= n. Row reversed Euler's triangle. See Graham et al., p. 268.
Recurrence (from A173018): T(n, 0) = 1 if n=0 else 0; T(n, k) = 0 if n < k and T(n, k) = (n+1-k)*T(n-1, k-1) + k*T(n-1, k) else.
T(n, k) = Sum_{j=0..k} (-1)^(k-j)*binomial(n-j, k-j)*S2(n, j)*j!, 0 <= k <= n, else 0. For S2(n, k)*k! see A131689.
The recurrence for the o.g.f. of the sequence of column k is
G(k, x) = (x/(1 - k*x))*(E_x - (k-2))*G(k-1, x), with the Euler operator E_x = x*d_x, for k >= 1, with G(0, x) = 1. (Proof from the recurrence of T(n, k)).
The e.g.f of the sequence of column k is found from E(k, x) = (1 + int(A(k, x),x)*exp(-k*x))*exp(k*x), k >= 1, with the recurrence
A(k, x) = x*A(k-1, x) +(1 + (1-k)*(1-x))*E(k-1, x) for k >= 1, with A(0,x)= 0. (Proof from the recurrence of T(n, k)). (End)
T(n, k) = Sum_{j=0..n-k} (-1)^j*(n-j-k+1)^n*binomial(n + 1, j). - Peter Luschny, Aug 12 2022
G.f.: Sum_{m >= 0} x^m/(1/(1-x)-m*t). - Mamuka Jibladze, Mar 12 2025

A145905 Square array read by antidiagonals: Hilbert transform of triangle A060187.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 9, 5, 1, 1, 27, 25, 7, 1, 1, 81, 125, 49, 9, 1, 1, 243, 625, 343, 81, 11, 1, 1, 729, 3125, 2401, 729, 121, 13, 1, 1, 2187, 15625, 16807, 6561, 1331, 169, 15, 1, 1, 6561, 78125, 117649, 59049, 14641, 2197, 225, 17, 1, 1, 19683, 390625, 823543
Offset: 0

Views

Author

Peter Bala, Oct 27 2008

Keywords

Comments

Definition of the Hilbert transform of a triangular array:
For many square arrays in the database the entries in a row are polynomial in the column index, of degree d say and hence the row generating function has the form P(x)/(1-x)^(d+1), where P is some polynomial function. Often the array whose rows are formed from the coefficients of these P polynomials is of independent interest. This suggests the following definition.
Let [L(n,k)]n,k>=0 be a lower triangular array and let R(n,x) := sum {k = 0 .. n} L(n,k)*x^k, denote the n-th row generating polynomial of L. Then we define the Hilbert transform of L, denoted Hilb(L), to be the square array whose n-th row, n >= 0, has the generating function R(n,x)/(1-x)^(n+1).
In this particular case, L is the array A060187, the array of Eulerian numbers of type B, whose row polynomials are the h-polynomials for permutohedra of type B. The Hilbert transform is an infinite Vandermonde matrix V(1,3,5,...).
We illustrate the Hilbert transform with a few examples:
(1) The Delannoy number array A008288 is the Hilbert transform of Pascal's triangle A007318 (view as the array of coefficients of h-polynomials of n-dimensional cross polytopes).
(2) The transpose of the array of nexus numbers A047969 is the Hilbert transform of the triangle of Eulerian numbers A008292 (best viewed in this context as the coefficients of h-polynomials of n-dimensional permutohedra of type A).
(3) The sequence of Eulerian polynomials begins [1, x, x + x^2, x + 4*x^2 + x^3, ...]. The coefficients of these polynomials are recorded in triangle A123125, whose Hilbert transform is A004248 read as square array.
(4) A108625, the array of crystal ball sequences for the A_n lattices, is the Hilbert transform of A008459 (viewed as the triangle of coefficients of h-polynomials of n-dimensional associahedra of type B).
(5) A142992, the array of crystal ball sequences for the C_n lattices, is the Hilbert transform of A086645, the array of h-vectors for type C root polytopes.
(6) A108553, the array of crystal ball sequences for the D_n lattices, is the Hilbert transform of A108558, the array of h-vectors for type D root polytopes.
(7) A086764, read as a square array, is the Hilbert transform of the rencontres numbers A008290.
(8) A143409 is the Hilbert transform of triangle A073107.

Examples

			Triangle A060187 (with an offset of 0) begins
1;
1, 1;
1, 6, 1;
so the entries in the first three rows of the Hilbert transform of
A060187 come from the expansions:
Row 0: 1/(1-x) = 1 + x + x^2 + x^3 + ...;
Row 1: (1+x)/(1-x)^2 = 1 + 3*x + 5*x^2 + 7*x^3 + ...;
Row 2: (1+6*x+x^2)/(1-x)^3 = 1 + 9*x + 25*x^2 + 49*x^3 + ...;
The array begins
n\k|..0....1.....2.....3......4
================================
0..|..1....1.....1.....1......1
1..|..1....3.....5.....7......9
2..|..1....9....25....49.....81
3..|..1...27...125...343....729
4..|..1...81...625..2401...6561
5..|..1..243..3125.16807..59049
...
		

Crossrefs

Cf. A008292, A039755, A052750 (first superdiagonal), A060187, A114172, A145901.

Programs

  • Maple
    T:=(n,k) -> (2*k + 1)^n: seq(seq(T(n-k,k),k = 0..n),n = 0..10);

Formula

T(n,k) = (2*k + 1)^n, (see equation 4.10 in [Franssens]). This array is the infinite Vandermonde matrix V(1,3,5,7, ....) having a LDU factorization equal to A039755 * diag(2^n*n!) * transpose(A007318).

A093375 Array T(m,n) read by ascending antidiagonals: T(m,n) = m*binomial(n+m-2, n-1) for m, n >= 1.

Original entry on oeis.org

1, 2, 1, 3, 4, 1, 4, 9, 6, 1, 5, 16, 18, 8, 1, 6, 25, 40, 30, 10, 1, 7, 36, 75, 80, 45, 12, 1, 8, 49, 126, 175, 140, 63, 14, 1, 9, 64, 196, 336, 350, 224, 84, 16, 1, 10, 81, 288, 588, 756, 630, 336, 108, 18, 1, 11, 100, 405, 960, 1470, 1512, 1050, 480, 135, 20, 1, 12
Offset: 1

Views

Author

Ralf Stephan, Apr 28 2004

Keywords

Comments

Number of n-long m-ary words avoiding the pattern 1-1'2'.
T(n,n+1) = Sum_{i=1..n} T(n,i).
Exponential Riordan array [(1+x)e^x, x] as a number triangle. - Paul Barry, Feb 17 2009
From Peter Bala, Jul 22 2014: (Start)
Call this array M and for k = 0,1,2,... define M(k) to be the lower unit triangular block array
/I_k 0\
\ 0 M/
having the k X k identity matrix I_k as the upper left block; in particular, M(0) = M. The infinite matrix product M(0)*M(1)*M(2)*..., which is clearly well-defined, is equal to A059298. (End)

Examples

			Array T(m,n) (with rows m >= 1 and columns n >= 1) begins as follows:
   1   1   1   1   1   1 ...
   2   4   6   8  10  12 ...
   3   9  18  30  45  63 ...
   4  16  40  80 140 224 ...
   5  25  75 175 350 630 ...
   ...
Triangle S(n,k) = T(n-k+1, k+1) begins
.n\k.|....0....1....2....3....4....5....6
= = = = = = = = = = = = = = = = = = = = =
..0..|....1
..1..|....2....1
..2..|....3....4....1
..3..|....4....9....6....1
..4..|....5...16...18....8....1
..5..|....6...25...40...30...10....1
..6..|....7...36...75...80...45...12....1
...
		

Crossrefs

Rows include A045943. Columns include A002411, A027810.
Main diagonal is A037965. Subdiagonals include A002457.
Antidiagonal sums are A001792.
See A103283 for a signed version.
Cf. A103406, A059298, A073107 (unsigned inverse).

Programs

  • GAP
    nmax:=14;; T:=List([1..nmax],n->List([1..nmax],k->k*Binomial(n+k-2,n-1)));;
    b:=List([2..nmax],n->OrderedPartitions(n,2));;
    a:=Flat(List([1..Length(b)],i->List([1..Length(b[i])],j->T[b[i][j][1]][b[i][j][2]]))); # Muniru A Asiru, Aug 07 2018
    
  • Mathematica
    nmax = 10;
    T = Transpose[CoefficientList[# + O[z]^(nmax+1), z]& /@ CoefficientList[(1 - x z)/(1 - z - x z)^2 + O[x]^(nmax+1), x]];
    row[n_] := T[[n+1, 1 ;; n+1]];
    Table[row[n], {n, 0, nmax}] // Flatten (* Jean-François Alcover, Aug 07 2018 *)
  • Sage
    # uses[riordan_array from A256893]
    riordan_array((1+x)*exp(x), x, 8, exp=true) # Peter Luschny, Nov 02 2019

Formula

Triangle = P*M, the binomial transform of the infinite bidiagonal matrix M with (1,1,1,...) in the main diagonal and (1,2,3,...) in the subdiagonal, and zeros elsewhere. P = Pascal's triangle as an infinite lower triangular matrix. - Gary W. Adamson, Nov 05 2006
From Peter Bala, Sep 20 2012: (Start)
E.g.f. for triangle: (1 + z)*exp((1 + x)*z) = 1 + (2 + x)*z + (3 + 4*x + x^2)*z^2/2! + ....
O.g.f. for triangle: (1 - x*z)/(1 - z - x*z)^2 = 1 + (2 + x)*z + (3 + 4*x + x^2)*z^2 + ....
The n-th row polynomial R(n,x) of the triangle equals (1+x)^n + n*(1+x)^(n-1) for n >= 0 and satisfies d/dx(R(n,x)) = n*R(n-1,x), as well as R(n,x+y) = Sum_{k = 0..n} binomial(n,k)*R(k,x)*y^(n-k). The row polynomials are a Sheffer sequence of Appell type.
Matrix inverse of the triangle is a signed version of A073107. (End)
From Tom Copeland, Oct 20 2015: (Start)
With offset 0 and D = d/dx, the raising operator for the signed row polynomials P(n,x) is RP = x - d{log[e^D/(1-D)]}/dD = x - 1 - 1/(1-D) = x - 2 - D - D^2 + ..., i.e., RP P(n,x) = P(n+1,x).
The e.g.f. for the signed array is (1-t) * e^(-t) * e^(x*t).
From the Appell formalism, the row polynomials PI(n,x) of A073107 are the umbral inverse of this entry's row polynomials; that is, P(n,PI(.,x)) = x^n = PI(n,P(.,x)) under umbral composition. (End)
From Petros Hadjicostas, Nov 01 2019: (Start)
As a triangle, we let S(n,k) = T(n-k+1, k+1) = (n-k+1)*binomial(n, k) for n >= 0 and 0 <= k <= n. See the example below.
As stated above by Peter Bala, Sum_{n,k >= 0} S(n,k)*z^n*x^k = (1 - x*z)/(1 - z -x*z)^2.
Also, Sum_{n, k >= 0} S(n,k)*z^n*x^k/n! = (1+z)*exp((1+x)*z).
As he also states, the n-th row polynomial is R(n,x) = Sum_{k = 0..n} S(n, k)*x^k = (1 + x)^n + n*(1 + x)^(n-1).
If we define the signed triangle S*(n,k) = (-1)^(n+k) * S(n,k) = (-1)^(n+k) * T(n-k+1, k+1), as Tom Copeland states, Sum_{n,k >= 0} S^*(n,k)*t^n*x^k/n! = (1-t)*exp((1-x)*(-t)) = (1-t) * e^(-t) * e^(x*t).
Apparently, S*(n,k) = A103283(n,k).
As he says above, the signed n-th row polynomial is P(n,x) = (-1)^n*R(n,-x) = (x - 1)^n - n*(x - 1)^(n-1).
According to Gary W. Adamson, P(n,x) is "the monic characteristic polynomial of the n X n matrix with 2's on the diagonal and 1's elsewhere." (End)

A357479 a(n) = (n!/6) * Sum_{k=0..n-3} 1/k!.

Original entry on oeis.org

0, 0, 0, 1, 8, 50, 320, 2275, 18256, 164388, 1644000, 18084165, 217010200, 2821132886, 39495860768, 592437911975, 9479006592160, 161143112067400, 2900576017214016, 55110944327067273, 1102218886541346600, 23146596617368279930, 509225125582102160000
Offset: 0

Views

Author

Seiichi Manyama, Sep 30 2022

Keywords

Crossrefs

Column k=3 of A073107.

Programs

  • Mathematica
    Table[n!/6 Sum[1/k!,{k,0,n-3}],{n,0,30}] (* Harvey P. Dale, Apr 02 2023 *)
  • PARI
    a(n) = n!/6*sum(k=0, n-3, 1/k!);
    
  • PARI
    a(n) = n!*sum(k=0, n, binomial(k, 3)/k!);
    
  • PARI
    my(N=30, x='x+O('x^N)); concat([0, 0, 0], Vec(serlaplace(x^3/6*exp(x)/(1-x))))
    
  • PARI
    my(N=30, x='x+O('x^N)); concat([0, 0, 0], Vec(sum(k=3, N, k!*x^k/(1-x)^(k+1))/6))

Formula

a(n) = n! * Sum_{k=0..n} binomial(k,3)/k!.
a(0) = 0; a(n) = n * a(n-1) + binomial(n,3).
E.g.f.: x^3/6 * exp(x)/(1-x).
G.f.: (1/6) * Sum_{k>=3} k! * x^k/(1-x)^(k+1).

A357480 a(n) = (n!/24) * Sum_{k=0..n-4} 1/k!.

Original entry on oeis.org

0, 0, 0, 0, 1, 10, 75, 560, 4550, 41076, 410970, 4521000, 54252495, 705283150, 9873965101, 148109477880, 2369751647900, 40285778016680, 725144004303300, 13777736081766576, 275554721635336365, 5786649154342069650, 127306281395525539615, 2928044472097087420000
Offset: 0

Views

Author

Seiichi Manyama, Sep 30 2022

Keywords

Crossrefs

Column k=4 of A073107.

Programs

  • PARI
    a(n) = n!/24*sum(k=0, n-4, 1/k!);
    
  • PARI
    a(n) = n!*sum(k=0, n, binomial(k, 4)/k!);
    
  • PARI
    my(N=30, x='x+O('x^N)); concat([0, 0, 0, 0], Vec(serlaplace(x^4/24*exp(x)/(1-x))))
    
  • PARI
    my(N=30, x='x+O('x^N)); concat([0, 0, 0, 0], Vec(sum(k=4, N, k!*x^k/(1-x)^(k+1))/24))

Formula

a(n) = n! * Sum_{k=0..n} binomial(k,4)/k!.
a(0) = 0; a(n) = n * a(n-1) + binomial(n,4).
E.g.f.: x^4/24 * exp(x)/(1-x).
G.f.: (1/24) * Sum_{k>=4} k! * x^k/(1-x)^(k+1).

A073474 Triangle T(n,k) read by rows, where o.g.f. for T(n,k) is n!*Sum_{k=0..n} (1+x)^(n-k)/k!.

Original entry on oeis.org

1, 2, 1, 5, 6, 2, 16, 33, 24, 6, 65, 196, 228, 120, 24, 326, 1305, 2120, 1740, 720, 120, 1957, 9786, 20550, 23160, 14760, 5040, 720, 13700, 82201, 212352, 305970, 265440, 138600, 40320, 5040, 109601, 767208, 2356424, 4146576, 4571280, 3232320, 1431360, 362880, 40320
Offset: 0

Views

Author

Vladeta Jovovic, Aug 26 2002

Keywords

Comments

Row sums give A010844.

Examples

			Triangle begins:
    1;
    2,    1;
    5,    6,    2;
   16,   33,   24,    6;
   65,  196,  228,  120,  24;
  326, 1305, 2120, 1740, 720, 120;
  ...
		

Crossrefs

Cf. A000142, A000522, A073107, A010844 (row sums).

Programs

  • Maple
    G:=simplify(series(exp(x)/(1-x-x*y),x=0,13)): P[0]:=1: for n from 1 to 11 do P[n]:=sort(n!*coeff(G,x^n)) od: seq(seq(coeff(y*P[n],y^k),k=1..n+1),n=0..9);
    # second Maple program:
    b:= proc(n, k) option remember; `if`(k>n, 0, `if`(k=0, 1,
          n*(b(n-1, k-1)+b(n-1, k))))
        end:
    T:= (n, k)-> b(n+1, k+1)/(n+1):
    seq(seq(T(n, k), k=0..n), n=0..9);  # Alois P. Heinz, Sep 12 2019
  • Mathematica
    b[n_, k_] := b[n, k] = If[k>n, 0, If[k==0, 1, n (b[n-1, k-1]+b[n-1, k])]];
    T[n_, k_] := b[n+1, k+1]/(n+1);
    Table[T[n, k], {n, 0, 9}, {k, 0, n}] // Flatten (* Jean-François Alcover, Dec 09 2019, after Alois P. Heinz *)
    T[n_, k_] := Sum[Binomial[j, k] FactorialPower[n, j], {j, 0, n}]; (* Peter Luschny, Oct 16 2024 *)
  • SageMath
    def T(n, k): return sum(binomial(j, k) * falling_factorial(n, j) for j in range(n+1))
    for n in range(8): print([T(n, k) for k in range(n+1)])
    # Peter Luschny, Oct 16 2024

Formula

E.g.f.: exp(x)/(1-x-x*y). - Vladeta Jovovic, Oct 17 2003
T(n, k) = Sum_{j=0..n} binomial(j, k)*FallingFactorial(n, j). - Peter Luschny, Oct 16 2024

Extensions

Edited by Emeric Deutsch, Jun 10 2004

A377661 Triangle read by rows: T(n, k) = e*Gamma(n - k + 1, 1)*binomial(n, k)^2.

Original entry on oeis.org

1, 2, 1, 5, 8, 1, 16, 45, 18, 1, 65, 256, 180, 32, 1, 326, 1625, 1600, 500, 50, 1, 1957, 11736, 14625, 6400, 1125, 72, 1, 13700, 95893, 143766, 79625, 19600, 2205, 98, 1, 109601, 876800, 1534288, 1022336, 318500, 50176, 3920, 128, 1
Offset: 0

Views

Author

Peter Luschny, Nov 03 2024

Keywords

Examples

			Triangle starts:
[0]      1;
[1]      2,      1;
[2]      5,      8,       1;
[3]     16,     45,      18,       1;
[4]     65,    256,     180,      32,      1;
[5]    326,   1625,    1600,     500,     50,     1;
[6]   1957,  11736,   14625,    6400,   1125,    72,    1;
[7]  13700,  95893,  143766,   79625,  19600,  2205,   98,   1;
[8] 109601, 876800, 1534288, 1022336, 318500, 50176, 3920, 128, 1;
		

Crossrefs

Cf. A000522 (column 0), A001105 (subdiagonal), A377662 (row sums), A073107.

Programs

  • Maple
    T := (n, k) -> exp(1)*GAMMA(n - k + 1, 1)*binomial(n, k)^2:
    seq(seq(simplify(T(n, k)), k = 0..n), n=0..8);
    # Alternative:
    A377661 := (n, k) -> n!*binomial(n,k)*add(1/(k!*(j-k)!), j = k..n):
    for n from 0 to 8 do seq(A377661(n, k), k = 0..n) od;
    # Or:
    T := (n, k) -> binomial(n, k)^2 * KummerU(k - n, k - n, 1):
    for n from 0 to 8 do seq(simplify(T(n, k)), k= 0..n) od;
  • Mathematica
    T[n_, k_] := E Binomial[n, k]^2 Gamma[1 - k + n, 1];
    Table[T[n, k], {n, 0, 8}, {k, 0, n}] // Flatten
  • Python
    from math import comb, isqrt, factorial
    def A377661(n):
        a = (m:=isqrt(k:=n+1<<1))-(k<=m*(m+1))
        b = n-comb(a+1,2)
        fa, fb = factorial(a), factorial(b)
        return comb(a,b)*sum(fa//(fb*factorial(j-b)) for j in range(b,a+1)) # Chai Wah Wu, Nov 12 2024

Formula

T(n, k) = binomial(n, k)*Sum_{j=k..n} n!/(k!*(j-k)!).
T(n, k) = binomial(n, k)^2 * KummerU(k - n, k - n, 1).
T(n, k) = binomial(n, k) * A073107(n, k).

A377662 a(n) = Sum_{k=0..n} binomial(n, k) * Sum_{j=k..n} n!/(k!*(j - k)!). Row sums of A377661.

Original entry on oeis.org

1, 3, 14, 80, 534, 4102, 35916, 354888, 3915750, 47754938, 637840356, 9256590928, 144977618044, 2436460447020, 43719637179224, 834042701945520, 16852447379512710, 359468276129261730, 8070500634880125300, 190211302604157871680, 4695001374741310892820
Offset: 0

Views

Author

Peter Luschny, Nov 07 2024

Keywords

Crossrefs

Programs

  • Maple
    A377662 := n -> add(binomial(n, k)*add(n!/(k!*(j-k)!), j = k..n), k = 0..n):
    seq(A377662(n), n = 0..20);
    # Or:
    a := n -> add(binomial(n, k)^2 * KummerU(k-n, k-n, 1), k = 0..n):
    seq(simplify(a(n)), n = 0..20);
  • Mathematica
    a[n_] := E Sum[Gamma[n - k + 1, 1] Binomial[n, k]^2, {k, 0, n}];
    Table[a[n], {n, 0, 20}]

Formula

a(n) = e*Sum_{k=0..n} binomial(n, k)^2 * Gamma(n - k + 1, 1).
a(n) = Sum_{k=0..n} binomial(n, k)^2 * KummerU(k - n, k - n, 1).
a(n) = Sum_{k=0..n} binomial(n, k) * A073107(n, k).
From Vaclav Kotesovec, Nov 07 2024: (Start)
Recurrence: n*(n^2 - 9*n + 12)*a(n) = 2*(n^4 - 7*n^3 - 5*n^2 + 29*n - 12)*a(n-1) - (n-1)*(n^4 - 2*n^3 - 47*n^2 + 104*n - 48)*a(n-2) + 2*(n-2)^2*(2*n - 3)*(n^2 - 7*n + 4)*a(n-3).
a(n) ~ 2^(-1/2) * exp(2*sqrt(n) - n + 1/2) * n^(n + 1/4) * (1 + 31/(48*sqrt(n))). (End)

A073480 Triangle T(n,k) read by rows, where e.g.f. for T(n,k) is exp(x*y)*log(1+x)/(1-x).

Original entry on oeis.org

1, 1, 2, 5, 3, 3, 14, 20, 6, 4, 94, 70, 50, 10, 5, 444, 564, 210, 100, 15, 6, 3828, 3108, 1974, 490, 175, 21, 7, 25584, 30624, 12432, 5264, 980, 280, 28, 8, 270576, 230256, 137808, 37296, 11844, 1764, 420, 36, 9, 2342880, 2705760, 1151280, 459360, 93240
Offset: 1

Views

Author

Vladeta Jovovic, Aug 26 2002

Keywords

Crossrefs

Cf. A073107.

Programs

  • Maple
    G:=exp(x*y)*ln(1+x)/(1-x): Gser:=series(G,x=0,12): for n from 1 to 10 do P[n]:=n!*coeff(Gser,x^n) od: for n from 1 to 10 do seq(coeff(y*P[n],y^k),k=1..n) od; # Emeric Deutsch, Dec 14 2004

Formula

E.g.f. for k-th column is x^k/k!*log(1+x)/(1-x).
O.g.f. for n-th row is Sum_{i=0..n} binomial(n, i)*A024167(n-i)*y^i.

Extensions

More terms from Emeric Deutsch, Dec 14 2004
Showing 1-9 of 9 results.