A008277 -id:A008277 - cp's OEIS frontend

A035342 The convolution matrix of the double factorial of odd numbers (A001147).

1, 3, 1, 15, 9, 1, 105, 87, 18, 1, 945, 975, 285, 30, 1, 10395, 12645, 4680, 705, 45, 1, 135135, 187425, 82845, 15960, 1470, 63, 1, 2027025, 3133935, 1595790, 370125, 43890, 2730, 84, 1, 34459425, 58437855, 33453945, 8998290

Offset: 1

Wolfdieter Lang

Previous name was: A triangle of numbers related to the triangle A035324; generalization of Stirling numbers of second kind A008277 and Lah numbers A008297.
If one replaces in the recurrence the '2' by '0', resp. '1', one obtains the Lah-number, resp. Stirling-number of 2nd kind, triangle A008297, resp. A008277.
The product of two lower triangular Jabotinsky matrices (see A039692 for the Knuth 1992 reference) is again such a Jabotinsky matrix: J(n,m) = Sum_{j=m..n} J1(n,j)*J2(j,m). The e.g.f.s of the first columns of these triangular matrices are composed in the reversed order: f(x)=f2(f1(x)). With f1(x)=-(log(1-2*x))/2 for J1(n,m)=|A039683(n,m)| and f2(x)=exp(x)-1 for J2(n,m)=A008277(n,m) one has therefore f2(f1(x))=1/sqrt(1-2*x) - 1 = f(x) for J(n,m)=a(n,m). This proves the matrix product given below. The m-th column of a Jabotinsky matrix J(n,m) has e.g.f. (f(x)^m)/m!, m>=1.
a(n,m) gives the number of forests with m rooted ordered trees with n non-root vertices labeled in an organic way. Organic labeling means that the vertex labels along the (unique) path from the root with label 0 to any leaf (non-root vertex of degree 1) is increasing. Proof: first for m=1 then for m>=2 using the recurrence relation for a(n,m) given below. - Wolfdieter Lang, Aug 07 2007
Also the Bell transform of A001147(n+1) (adding 1,0,0,... as column 0). For the definition of the Bell transform see A264428. - Peter Luschny, Jan 19 2016

			Matrix begins:
    1;
    3,   1;
   15,   9,   1;
  105,  87,  18,   1;
  945, 975, 285,  30,   1;
  ...
Combinatoric meaning of a(3,2)=9: The nine increasing path sequences for the three rooted ordered trees with leaves labeled with 1,2,3 and the root labels 0 are: {(0,3),[(0,1),(0,2)]}; {(0,3),[(0,2),(0,1)]}; {(0,3),(0,1,2)}; {(0,1),[(0,3),(0,2)]}; [(0,1),[(0,2),(0,3)]]; [(0,2),[(0,1),(0,3)]]; {(0,2),[(0,3),(0,1)]}; {(0,1),(0,2,3)}; {(0,2),(0,1,3)}.

Reinhard Zumkeller, Rows n = 1..125 of triangle, flattened
Peter Bala, Generalized Dobinski formulas
J. Fernando Barbero G., Jesús Salas, and Eduardo J. S. Villaseñor, Bivariate Generating Functions for a Class of Linear Recurrences. I. General Structure, arXiv:1307.2010 [math.CO], 2013-2014.
Pawel Blasiak, Karol A. Penson, and Allan I. Solomon, The Boson Normal Ordering Problem and Generalized Bell Numbers, arXiv:quant-ph/0212072, 2002.
Pawel Blasiak, Karol A. Penson, and Allan I. Solomon, The general boson normal ordering problem, arXiv:quant-ph/0402027, 2004.
Richell O. Celeste, Roberto B. Corcino and Ken Joffaniel M. Gonzales. Two Approaches to Normal Order Coefficients. Journal of Integer Sequences, Vol. 20 (2017), Article 17.3.5.
Tom Copeland, A Class of Differential Operators and the Stirling Numbers, 2015.
Tom Copeland, Addendum to Mathemagical Forests, 2010.
Tom Copeland, Mathemagical Forests, 2008.
Askar Dzhumadildaev and Damir Yeliussizov, Path decompositions of digraphs and their applications to Weyl algebra, arXiv preprint arXiv:1408.6764v1 [math.CO], 2014. [Version 1 contained many references to the OEIS, which were removed in Version 2. - _N. J. A. Sloane_, Mar 28 2015]
Askar Dzhumadil'daev and Damir Yeliussizov, Walks, partitions, and normal ordering, Electronic Journal of Combinatorics, 22(4) (2015), #P4.10.
Milan Janjic, Some classes of numbers and derivatives, JIS 12 (2009) #09.8.3.
Donald E. Knuth, Convolution polynomials, arXiv:math/9207221 [math.CA], 1992; Mathematica J. 2.1 (1992), no. 4, 67-78.
Wolfdieter Lang, On generalizations of Stirling number triangles, J. Integer Seqs., Vol. 3 (2000), #00.2.4.
Wolfdieter Lang, First 10 rows.
Shi-Mei Ma, Some combinatorial sequences associated with context-free grammars, arXiv:1208.3104v2 [math.CO], 2012.
Toufik Mansour, Matthias Schork and Mark Shattuck, The Generalized Stirling and Bell Numbers Revisited, Journal of Integer Sequences, Vol. 15 (2012), #12.8.3.
Erich Neuwirth, Recursively defined combinatorial functions: Extending Galton's board, Discrete Math. (2001) Vol. 239, 33-51.
Mathias Pétréolle and Alan D. Sokal, Lattice paths and branched continued fractions. II. Multivariate Lah polynomials and Lah symmetric functions, arXiv:1907.02645 [math.CO], 2019.
Umesh Shankar, Log-concavity of rows of triangular arrays satisfying a certain super-recurrence, arXiv:2508.12467 [math.CO], 2025. See p. 4.

The column sequences are A001147, A035101, A035119, ...
Row sums: A049118(n), n >= 1.
Cf. A000108, A035324, A008277, A008297, A094638.

a035342 n k = a035342_tabl !! (n-1) !! (k-1)
a035342_row n = a035342_tabl !! (n-1)
a035342_tabl = map fst $ iterate (\(xs, i) -> (zipWith (+)
   ([0] ++ xs) $ zipWith (*) [i..] (xs ++ [0]), i + 2)) ([1], 3)
-- Reinhard Zumkeller, Mar 12 2014

T := (n,k) -> 2^(k-n)*hypergeom([k-n,k+1],[k-2*n+1],2)*GAMMA(2*n-k)/
(GAMMA(k)*GAMMA(n-k+1)); for n from 1 to 9 do seq(simplify(T(n,k)),k=1..n) od; # Peter Luschny, Mar 31 2015
T := (n, k) -> local j; 2^n*add((-1)^(k-j)*binomial(k, j)*pochhammer(j/2, n), j = 1..k)/k!: for n from 1 to 6 do seq(T(n, k), k=1..n) od;  # Peter Luschny, Mar 04 2024

a[n_, k_] := 2^(n+k)*n!/(4^n*n*k!)*Sum[(j+k)*2^(j)*Binomial[j + k - 1, k-1]*Binomial[2*n - j - k - 1, n-1], {j, 0, n-k}]; Flatten[Table[a[n,k], {n, 1, 9}, {k, 1, n}] ] [[1 ;; 40]] (* Jean-François Alcover, Jun 01 2011, after Vladimir Kruchinin *)

a(n,k):=2^(n+k)*n!/(4^n*n*k!)*sum((j+k)*2^(j)*binomial(j+k-1,k-1)*binomial(2*n-j-k-1,n-1),j,0,n-k); /* Vladimir Kruchinin, Mar 30 2011 */

# uses[bell_matrix from A264428]
# Adds a column 1,0,0,0, ... at the left side of the triangle.
print(bell_matrix(lambda n: A001147(n+1), 9)) # Peter Luschny, Jan 19 2016

a(n, m) = Sum_{j=m..n} |A039683(n, j)|*S2(j, m) (matrix product), with S2(j, m) := A008277(j, m) (Stirling2 triangle). Priv. comm. to Wolfdieter Lang by E. Neuwirth, Feb 15 2001; see also the 2001 Neuwirth reference. See the comment on products of Jabotinsky matrices.

a(n, m) = n!*A035324(n, m)/(m!*2^(n-m)), n >= m >= 1; a(n+1, m)= (2*n+m)*a(n, m)+a(n, m-1); a(n, m) := 0, n

E.g.f. of m-th column: ((x*c(x/2)/sqrt(1-2*x))^m)/m!, where c(x) = g.f. for Catalan numbers A000108.

From Vladimir Kruchinin, Mar 30 2011: (Start)

G.f. (1/sqrt(1-2*x) - 1)^k = Sum_{n>=k} (k!/n!)*a(n,k)*x^n.

a(n,k) = 2^(n+k) * n! / (4^n*n*k!) * Sum_{j=0..n-k} (j+k) * 2^(j) * binomial(j+k-1,k-1) * binomial(2*n-j-k-1,n-1). (End)

From Peter Bala, Nov 25 2011: (Start)

E.g.f.: G(x,t) = exp(t*A(x)) = 1 + t*x + (3*t + t^2)*x^2/2! + (15*t + 9*t^2 + t^3)*x^3/3! + ..., where A(x) = -1 + 1/sqrt(1-2*x) satisfies the autonomous differential equation A'(x) = (1+A(x))^3.

The generating function G(x,t) satisfies the partial differential equation t*(dG/dt+G) = (1-2*x)*dG/dx, from which follows the recurrence given above.

The row polynomials are given by D^n(exp(x*t)) evaluated at x = 0, where D is the operator (1+x)^3*d/dx. Cf. A008277 (D = (1+x)*d/dx), A105278 (D = (1+x)^2*d/dx), A035469 (D = (1+x)^4*d/dx) and A049029 (D = (1+x)^5*d/dx). (End)

The n-th row polynomial R(n,x) is given by the Dobinski-type formula R(n,x) = exp(-x)*Sum_{k>=1} k*(k+2)*...*(k+2*n-2)*x^k/k!. - Peter Bala, Jun 22 2014

T(n,k) = 2^(k-n)*hypergeom([k-n,k+1],[k-2*n+1],2)*Gamma(2*n-k)/(Gamma(k)*Gamma(n-k+1)). - Peter Luschny, Mar 31 2015

T(n,k) = 2^n*Sum_{j=1..k} ((-1)^(k-j)*binomial(k, j)*Pochhammer(j/2, n)) / k!. - Peter Luschny, Mar 04 2024

Simpler name from Peter Luschny, Mar 31 2015

A028246 Triangular array a(n,k) = (1/k)Sum_{i=0..k} (-1)^(k-i)binomial(k,i)*i^n; n >= 1, 1 <= k <= n, read by rows.

Original entry on oeis.org

1, 1, 1, 1, 3, 2, 1, 7, 12, 6, 1, 15, 50, 60, 24, 1, 31, 180, 390, 360, 120, 1, 63, 602, 2100, 3360, 2520, 720, 1, 127, 1932, 10206, 25200, 31920, 20160, 5040, 1, 255, 6050, 46620, 166824, 317520, 332640, 181440, 40320, 1, 511, 18660, 204630, 1020600, 2739240, 4233600, 3780000, 1814400, 362880

Offset: 1

N. J. A. Sloane, Doug McKenzie (mckfam4(AT)aol.com)

Let M = n X n matrix with (i,j)-th entry a(n+1-j, n+1-i), e.g., if n = 3, M = [1 1 1; 3 1 0; 2 0 0]. Given a sequence s = [s(0)..s(n-1)], let b = [b(0)..b(n-1)] be its inverse binomial transform and let c = [c(0)..c(n-1)] = M^(-1)*transpose(b). Then s(k) = Sum_{i=0..n-1} b(i)*binomial(k,i) = Sum_{i=0..n-1} c(i)*k^i, k=0..n-1. - Gary W. Adamson, Nov 11 2001
From Gary W. Adamson, Aug 09 2008: (Start)
Julius Worpitzky's 1883 algorithm generates Bernoulli numbers.
By way of example [Wikipedia]:
  B0 =  1;
  B1 = 1/1 -  1/2;
  B2 = 1/1 -  3/2 +   2/3;
  B3 = 1/1 -  7/2 +  12/3 -    6/4;
  B4 = 1/1 - 15/2 +  50/3 -   60/4 +   24/5;
  B5 = 1/1 - 31/2 + 180/3 -  390/4 +  360/5 -  120/6;
  B6 = 1/1 - 63/2 + 602/3 - 2100/4 + 3360/5 - 2520/6 + 720/7;
  ...
Note that in this algorithm, odd n's for the Bernoulli numbers sum to 0, not 1, and the sum for B1 = 1/2 = (1/1 - 1/2). B3 = 0 = (1 - 7/2 + 13/3 - 6/4) = 0. The summation for B4 = -1/30. (End)
Pursuant to Worpitzky's algorithm and given M = A028246 as an infinite lower triangular matrix, M * [1/1, -1/2, 1/3, ...] (i.e., the Harmonic series with alternate signs) = the Bernoulli numbers starting [1/1, 1/2, 1/6, ...]. - Gary W. Adamson, Mar 22 2012
From Tom Copeland, Oct 23 2008: (Start)
G(x,t) = 1/(1 + (1-exp(x*t))/t) = 1 + 1 x + (2 + t)*x^2/2! + (6 + 6t + t^2)*x^3/3! + ... gives row polynomials for A090582, the f-polynomials for the permutohedra (see A019538).
G(x,t-1) = 1 + 1*x + (1 + t)*x^2 / 2! + (1 + 4t + t^2)*x^3 / 3! + ... gives row polynomials for A008292, the h-polynomials for permutohedra.
G[(t+1)x,-1/(t+1)] = 1 + (1+ t) x + (1 + 3t + 2 t^2) x^2 / 2! + ... gives row polynomials for the present triangle. (End)
The Worpitzky triangle seems to be an apt name for this triangle. - Johannes W. Meijer, Jun 18 2009
If Pascal's triangle is written as a lower triangular matrix and multiplied by A028246 written as an upper triangular matrix, the product is a matrix where the (i,j)-th term is (i+1)^j. For example,
  1,0,0,0     1,1,1, 1      1,1, 1, 1
  1,1,0,0  *  0,1,3, 7   =  1,2, 4, 8
  1,2,1,0     0,0,2,12      1,3, 9,27
  1,3,3,1     0,0,0, 6      1,4,16,64
  So, numbering all three matrices' rows and columns starting at 0, the (i,j) term of the product is (i+1)^j. - Jack A. Cohen (ProfCohen(AT)comcast.net), Aug 03 2010
The Fi1 and Fi2 triangle sums are both given by sequence A000670. For the definition of these triangle sums see A180662. The mirror image of the Worpitzky triangle is A130850. - Johannes W. Meijer, Apr 20 2011
Let S_n(m) = 1^m + 2^m + ... + n^m. Then, for n >= 0, we have the following representation of S_n(m) as a linear combination of the binomial coefficients:
  S_n(m) = Sum_{i=1..n+1} a(i+n*(n+1)/2)*C(m,i). E.g., S_2(m) = a(4)*C(m,1) + a(5)*C(m,2) + a(6)*C(m,3) = C(m,1) + 3*C(m,2) + 2*C(m,3). - Vladimir Shevelev, Dec 21 2011
Given the set X = [1..n] and 1 <= k <= n, then a(n,k) is the number of sets T of size k of subset S of X such that S is either empty or else contains 1 and another element of X and such that any two elemements of T are either comparable or disjoint. - Michael Somos, Apr 20 2013
Working with the row and column indexing starting at -1, a(n,k) gives the number of k-dimensional faces in the first barycentric subdivision of the standard n-dimensional simplex (apply Brenti and Welker, Lemma 2.1). For example, the barycentric subdivision of the 2-simplex (a triangle) has 1 empty face, 7 vertices, 12 edges and 6 triangular faces giving row 4 of this triangle as (1,7,12,6). Cf. A053440. - Peter Bala, Jul 14 2014
See A074909 and above g.f.s for associations among this array and the Bernoulli polynomials and their umbral compositional inverses. - Tom Copeland, Nov 14 2014
An e.g.f. G(x,t) = exp[P(.,t)x] = 1/t - 1/[t+(1-t)(1-e^(-xt^2))] = (1-t) * x + (-2t + 3t^2 - t^3) * x^2/2! + (6t^2 - 12t^3 + 7t^4 - t^5) * x^3/3! + ... for the shifted, reverse, signed polynomials with the first element nulled, is generated by the infinitesimal generator g(u,t)d/du = [(1-u*t)(1-(1+u)t)]d/du, i.e., exp[x * g(u,t)d/du] u eval. at u=0 generates the polynomials. See A019538 and the G. Rzadkowski link below for connections to the Bernoulli and Eulerian numbers, a Ricatti differential equation, and a soliton solution to the KdV equation. The inverse in x of this e.g.f. is Ginv(x,t) = (-1/t^2)*log{[1-t(1+x)]/[(1-t)(1-tx)]} = [1/(1-t)]x + [(2t-t^2)/(1-t)^2]x^2/2 + [(3t^2-3t^3+t^4)/(1-t)^3]x^3/3  + [(4t^3-6t^4+4t^5-t^6)/(1-t)^4]x^4/4 + ... . The numerators are signed, shifted A135278 (reversed A074909), and the rational functions are the columns of A074909. Also, dG(x,t)/dx = g(G(x,t),t) (cf. A145271). (Analytic G(x,t) added, and Ginv corrected and expanded on Dec 28 2015.) - Tom Copeland, Nov 21 2014
The operator R = x + (1 + t) + t e^{-D} / [1 + t(1-e^(-D))] = x + (1+t) + t - (t+t^2) D + (t+3t^2+2t^3) D^2/2! - ... contains an e.g.f. of the reverse row polynomials of the present triangle, i.e., A123125 * A007318 (with row and column offset 1 and 1). Umbrally, R^n 1 = q_n(x;t) = (q.(0;t)+x)^n, with q_m(0;t) = (t+1)^(m+1) - t^(m+1), the row polynomials of A074909, and D = d/dx. In other words, R generates the Appell polynomials associated with the base sequence A074909. For example, R 1 = q_1(x;t) = (q.(0;t)+x) = q_1(0;t) + q__0(0;t)x = (1+2t) + x, and R^2 1 = q_2(x;t) = (q.(0;t)+x)^2 = q_2(0:t) + 2q_1(0;t)x + q_0(0;t)x^2 = 1+3t+3t^2 + 2(1+2t)x + x^2. Evaluating the polynomials at x=0 regenerates the base sequence. With a simple sign change in R, R generates the Appell polynomials associated with A248727. - Tom Copeland, Jan 23 2015
For a natural refinement of this array, see A263634. - Tom Copeland, Nov 06 2015
From Wolfdieter Lang, Mar 13 2017: (Start)
The e.g.f. E(n, x) for {S(n, m)}{m>=0} with S(n, m) = Sum{k=1..m} k^n, n >= 0, (with undefined sum put to 0) is exp(x)*R(n+1, x) with the exponential row polynomials R(n, x) = Sum_{k=1..n} a(n, k)*x^k/k!. E.g., e.g.f. for n = 2, A000330: exp(x)*(1*x/1!+3*x^2/2!+2*x^3/3!).
The o.g.f. G(n, x) for {S(n, m)}{m >=0} is then found by Laplace transform to be G(n, 1/p) = p*Sum{k=1..n} a(n+1, k)/(p-1)^(2+k).
  Hence G(n, x) = x/(1 - x)^(n+2)*Sum_{k=1..n} A008292(n,k)*x^(k-1).
  E.g., n=2: G(2, 1/p) = p*(1/(p-1)^2 + 3/(p-1)^3 + 2/(p-1)^4) = p^2*(1+p)/(p-1)^4; hence G(2, x) = x*(1+x)/(1-x)^4.
This works also backwards: from the o.g.f. to the e.g.f. of {S(n, m)}_{m>=0}. (End)
a(n,k) is the number of k-tuples of pairwise disjoint and nonempty subsets of a set of size n. - Dorian Guyot, May 21 2019
From Rajesh Kumar Mohapatra, Mar 16 2020: (Start)
a(n-1,k) is the number of chains of length k in a partially ordered set formed from subsets of an n-element set ordered by inclusion such that the first term of the chains is either the empty set or an n-element set.
Also, a(n-1,k) is the number of distinct k-level rooted fuzzy subsets of an n-set ordered by set inclusion. (End)
The relations on p. 34 of Hasan (also p. 17 of Franco and Hasan) agree with the relation between A019538 and this entry given in the formula section. - Tom Copeland, May 14 2020
T(n,k) is the size of the Green's L-classes in the D-classes of rank (k-1) in the semigroup of partial transformations on an (n-1)-set. - Geoffrey Critzer, Jan 09 2023
T(n,k) is the number of strongly connected binary relations on [n] that have period k (A367948) and index 1.  See Theorem 5.4.25(6) in Ki Hang Kim reference. - Geoffrey Critzer, Dec 07 2023

			The triangle a(n, k) starts:
n\k 1   2    3     4      5      6      7      8     9
1:  1
2:  1   1
3:  1   3    2
4:  1   7   12     6
5:  1  15   50    60     24
6:  1  31  180   390    360    120
7:  1  63  602  2100   3360   2520    720
8:  1 127 1932 10206  25200  31920  20160   5040
9:  1 255 6050 46620 166824 317520 332640 181440 40320
... [Reformatted by _Wolfdieter Lang_, Mar 26 2015]
-----------------------------------------------------
Row 5 of triangle is {1,15,50,60,24}, which is {1,15,25,10,1} times {0!,1!,2!,3!,4!}.
From _Vladimir Shevelev_, Dec 22 2011: (Start)
Also, for power sums, we have
S_0(n) = C(n,1);
S_1(n) = C(n,1) +    C(n,2);
S_2(n) = C(n,1) +  3*C(n,2) +  2*C(n,3);
S_3(n) = C(n,1) +  7*C(n,2) + 12*C(n,3) +  6*C(n,4);
S_4(n) = C(n,1) + 15*C(n,2) + 50*C(n,3) + 60*C(n,4) + 24*C(n,5); etc.
(End)
For X = [1,2,3], the sets T are {{}}, {{},{1,2}}, {{},{1,3}}, {{},{1,2,3}}, {{},{1,2},{1,2,3}}, {{},{1,3},{1,2,3}} and so a(3,1)=1, a(3,2)=3, a(3,3)=2. - _Michael Somos_, Apr 20 2013

Ki Hang Kim, Boolean Matrix Theory and Applications, Marcel Dekker, New York and Basel (1982).

Seiichi Manyama, Table of n, a(n) for n = 1..10000
V. S. Abramovich, Power sums of natural numbers, Kvant, no. 5 (1973), 22-25. (in Russian)
Peter Bala, Deformations of the Hadamard product of power series
Paul Barry, Three Études on a sequence transformation pipeline, arXiv:1803.06408 [math.CO], 2018.
Melike Baykal-Gürsoy, Marcelo Figueroa-Candia, and Zhe Duan, Completion times of jobs on two-state service processes and their asymptotic behavior, Probability in the Engineering and Information Sciences (2024), 1-29. See p. 23.
H. Belbachir, M. Rahmani, and B. Sury, Sums Involving Moments of Reciprocals of Binomial Coefficients, J. Int. Seq. 14 (2011) #11.6.6.
Hacene Belbachir and Mourad Rahmani, Alternating Sums of the Reciprocals of Binomial Coefficients, Journal of Integer Sequences, Vol. 15 (2012), #12.2.8.
F. Brenti and V. Welker, f-vectors of barycentric subdivisions, arXiv:math/0606356v1 [math.CO], Math. Z., 259(4), 849-865, 2008.
Patibandla Chanakya and Putla Harsha, Generalized Nested Summation of Powers of Natural Numbers, arXiv:1808.08699 [math.NT], 2018. See Table 1.
Tom Copeland, Generators, Inversion, and Matrix, Binomial, and Integral Transforms
Colin Defant, Troupes, Cumulants, and Stack-Sorting, arXiv:2004.11367 [math.CO], 2020.
E. Delucchi, A. Pixton and L. Sabalka. Face vectors of subdivided simplicial complexes arXiv:1002.3201v3 [math.CO], Discrete Mathematics, Volume 312, Issue 2, January 2012, Pages 248-257.
G. H. E Duchamp, N. Hoang-Nghia, and A. Tanasa, A word Hopf algebra based on the selection/quotient principle, arXiv:1207.6522 [math.CO], 2012-2013; Séminaire Lotharingien de Combinatoire 68 (2013), Article B68c.
M. Dukes and C. D. White, Web Matrices: Structural Properties and Generating Combinatorial Identities, arXiv:1603.01589 [math.CO], 2016.
Nick Early, Honeycomb tessellations and canonical bases for permutohedral blades, arXiv:1810.03246 [math.CO], 2018.
S. Franco and A. Hasan, Graded Quivers, Generalized Dimer Models and Toric Geometry , arXiv preprint arXiv:1904.07954 [hep-th], 2019
A. Hasan, Physics and Mathematics of Graded Quivers, dissertation, Graduate Center, City University of New York, 2019.
H. Hasse, Ein Summierungsverfahren für die Riemannsche Zeta-Reihe, Math. Z. 32, 458-464 (1930).
Guy Louchard, Werner Schachinger, and Mark Daniel Ward, The number of distinct adjacent pairs in geometrically distributed words: a probabilistic and combinatorial analysis, arXiv:2203.14773 [math.PR], 2022. See p. 5.
Shi-Mei Ma, A family of two-variable derivative polynomials for tangent and secant, El J. Combinat. 20(1) (2013), P11.
Richard J. Mathar, Integrals Associated with the Digamma Integral Representation, arXiv:2308.14154 [math.GM], 2023. See p. 5.
Rajesh Kumar Mohapatra and Tzung-Pei Hong, On the Number of Finite Fuzzy Subsets with Analysis of Integer Sequences, Mathematics (2022) Vol. 10, No. 7, 1161.
A. Riskin and D. Beckwith, Problem 10231, Amer. Math. Monthly, 102 (1995), 175-176.
G. Rzadkowski, Bernoulli numbers and solitons revisited, Journal of Nonlinear Mathematical Physics, Volume 17, Issue 1, 2010.
John K. Sikora, On Calculating the Coefficients of a Polynomial Generated Sequence Using the Worpitzky Number Triangles, arXiv:1806.00887 [math.NT], 2018.
G. J. Simmons, A combinatorial problem associated with a family of combination locks, Math. Mag., 37 (1964), 127-132 (but there are errors). The triangle is on page 129.
N. J. A. Sloane, Transforms
Sam Vandervelde, The Worpitzky Numbers Revisited, Amer. Math. Monthly, 125:3 (2018), 198-206.
Wikipedia, Bernoulli number.
Wikipedia, Barycentric subdivision
David C. Wood, The computation of polylogarithms (2014).

Dropping the column of 1's gives A053440.
Without the k in the denominator (in the definition), we get A019538. See also the Stirling number triangle A008277.
Cf. A087127, A087107, A087108, A087109, A087110, A087111, A084938 A075263.
Row sums give A000629(n-1) for n >= 1.
Cf. A027642, A002445. - Gary W. Adamson, Aug 09 2008
Appears in A161739 (RSEG2 triangle), A161742 and A161743. - Johannes W. Meijer, Jun 18 2009
Binomial transform is A038719. Cf. A131689.
Cf. A007318, A008292, A046802, A074909, A090582, A123125, A130850, A135278, A141618, A145271, A163626, A248727, A263634.
Cf. A119879.
From Rajesh Kumar Mohapatra, Mar 29 2020: (Start)
A000007(n-1) (column k=1), A000225(n-1) (column k=2), A028243(n-1) (column k=3), A028244(n-1) (column k=4), A028245(n-1) (column k=5), for n > 0.
Diagonal gives A000142(n-1), for n >=1.
Next-to-last diagonal gives A001710,
Third, fourth, fifth, sixth, seventh external diagonal respectively give A005460, A005461, A005462, A005463, A005464. (End)

Flat(List([1..10], n-> List([1..n], k-> Stirling2(n,k)* Factorial(k-1) ))); # G. C. Greubel, May 30 2019

[[StirlingSecond(n,k)*Factorial(k-1): k in [1..n]]: n in [1..10]]; // G. C. Greubel, May 30 2019

a := (n,k) -> add((-1)^(k-i)*binomial(k,i)*i^n, i=0..k)/k;
seq(print(seq(a(n,k),k=1..n)),n=1..10);
T := (n,k) -> add(eulerian1(n,j)*binomial(n-j,n-k), j=0..n):
seq(print(seq(T(n,k),k=0..n)),n=0..9); # Peter Luschny, Jul 12 2013

a[n_, k_] = Sum[(-1)^(k-i) Binomial[k,i]*i^n, {i,0,k}]/k; Flatten[Table[a[n, k], {n, 10}, {k, n}]] (* Jean-François Alcover, May 02 2011 *)

{T(n, k) = if( k<0 || k>n, 0, n! * polcoeff( (x / log(1 + x + x^2 * O(x^n) ))^(n+1), n-k))}; /* Michael Somos, Oct 02 2002 */

{T(n,k) = stirling(n,k,2)*(k-1)!}; \\ G. C. Greubel, May 31 2019

# Assuming offset (n, k) = (0, 0).
def T(n, k):
    if k >  n: return 0
    if k == 0: return 1
    return k*T(n - 1, k - 1) + (k + 1)*T(n - 1, k)
for n in range(9):
    print([T(n, k) for k in range(n + 1)])  # Peter Luschny, Apr 26 2022

def A163626_row(n) :
    x = polygen(ZZ,'x')
    A = []
    for m in range(0, n, 1) :
        A.append((-x)^m)
        for j in range(m, 0, -1):
            A[j - 1] = j * (A[j - 1] - A[j])
    return list(A[0])
for i in (1..7) : print(A163626_row(i))  # Peter Luschny, Jan 25 2012

[[stirling_number2(n,k)*factorial(k-1) for k in (1..n)] for n in (1..10)] # G. C. Greubel, May 30 2019

E.g.f.: -log(1-y*(exp(x)-1)). - Vladeta Jovovic, Sep 28 2003
a(n, k) = S2(n, k)*(k-1)! where S2(n, k) is a Stirling number of the second kind (cf. A008277). Also a(n,k) = T(n,k)/k, where T(n, k) = A019538.
Essentially same triangle as triangle [1, 0, 2, 0, 3, 0, 4, 0, 5, 0, 6, 0, 7, ...] DELTA [1, 1, 2, 2, 3, 3, 4, 4, 5, 5, ...] where DELTA is Deléham's operator defined in A084938, but the notation is different.
Sum of terms in n-th row = A000629(n) - Gary W. Adamson, May 30 2005
The row generating polynomials P(n, t) are given by P(1, t)=t, P(n+1, t) = t(t+1)(d/dt)P(n, t) for n >= 1 (see the Riskin and Beckwith reference). - Emeric Deutsch, Aug 09 2005
From Gottfried Helms, Jul 12 2006: (Start)
Delta-matrix as can be read from H. Hasse's proof of a connection between the zeta-function and Bernoulli numbers (see link below).
Let P = lower triangular matrix with entries P[row,col] = binomial(row,col).
Let J = unit matrix with alternating signs J[r,r]=(-1)^r.
Let N(m) = column matrix with N(m)(r) = (r+1)^m, N(1)--> natural numbers.
Let V = Vandermonde matrix with V[r,c] = (r+1)^c.
V is then also N(0)||N(1)||N(2)||N(3)... (indices r,c always beginning at 0).
Then Delta = P*J * V and B' = N(-1)' * Delta, where B is the column matrix of Bernoulli numbers and ' means transpose, or for the single k-th Bernoulli number B_k with the appropriate column of Delta,
B_k = N(-1)' * Delta[ *,k ] = N(-1)' * P*J * N(k).
Using a single column instead of V and assuming infinite dimension, H. Hasse showed that in x = N(-1) * P*J * N(s), where s can be any complex number and s*zeta(1-s) = x.
His theorem reads: s*zeta(1-s) = Sum_{n>=0..inf} (n+1)^-1*delta(n,s), where delta(n,s) = Sum_{j=0..n} (-1)^j * binomial(n,j) * (j+1)^s.
(End)
a(n,k) = k*a(n-1,k) + (k-1)*a(n-1,k-1) with a(n,1) = 1 and a(n,n) = (n-1)!. - Johannes W. Meijer, Jun 18 2009
Rephrasing the Meijer recurrence above: Let M be the (n+1)X(n+1) bidiagonal matrix with M(r,r) = M(r,r+1) = r, r >= 1, in the two diagonals and the rest zeros. The row a(n+1,.) of the triangle is row 1 of M^n. - Gary W. Adamson, Jun 24 2011
From Tom Copeland, Oct 11 2011: (Start)
With e.g.f.. A(x,t) = G[(t+1)x,-1/(t+1)]-1 (from 2008 comment) = -1 + 1/[1-(1+t)(1-e^(-x))] = (1+t)x + (1+3t+2t^2)x^2/2! + ..., the comp. inverse in x is
B(x,t)= -log(t/(1+t)+1/((1+t)(1+x))) = (1/(1+t))x - ((1+2t)/(1+t)^2)x^2/2 + ((1+3t+3t^2)/(1+t)^3)x^3/3 + .... The numerators are the row polynomials of A074909, and the rational functions are (omitting the initial constants) signed columns of the re-indexed Pascal triangle A007318.
Let h(x,t)= 1/(dB/dx) = (1+x)(1+t(1+x)), then the row polynomial P(n,t) = (1/n!)(h(x,t)*d/dx)^n x, evaluated at x=0, A=exp(x*h(y,t)*d/dy) y, eval. at y=0, and dA/dx = h(A(x,t),t), with P(1,t)=1+t. (Series added Dec 29 2015.)(End)
Let  denote the Eulerian numbers A173018(n,k), then T(n,k) = Sum_{j=0..n} *binomial(n-j,n-k). - Peter Luschny, Jul 12 2013
Matrix product A007318 * A131689. The n-th row polynomial R(n,x) = Sum_{k >= 1} k^(n-1)*(x/(1 + x))^k, valid for x in the open interval (-1/2, inf). Cf A038719. R(n,-1/2) = (-1)^(n-1)*(2^n - 1)*Bernoulli(n)/n. - Peter Bala, Jul 14 2014
a(n,k) = A141618(n,k) / C(n,k-1). - Tom Copeland, Oct 25 2014
For the row polynomials, A028246(n,x) = A019538(n-1,x) * (1+x). - Tom Copeland, Dec 28 2015
A248727 = A007318*(reversed A028246) = A007318*A130850 = A007318*A123125*A007318 = A046802*A007318. - Tom Copeland, Nov 14 2016
n-th row polynomial R(n,x) = (1+x) o (1+x) o ... o (1+x) (n factors), where o denotes the black diamond multiplication operator of Dukes and White. See example E11 in the Bala link. - Peter Bala, Jan 12 2018
From Dorian Guyot, May 21 2019: (Start)
Sum_{i=0..k} binomial(k,i) * a(n,i) = (k+1)^n.
Sum_{k=0..n} a(n,k) = 2*A000670(n).
(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 A123125. Then the row polynomials of this entry, A028246, are given by x^n * A_n(1 + 1/x;0). Other specializations of A_n(x;y) give A046802, A090582, A119879, A130850, and A248727. - Tom Copeland, Jan 24 2020
The row generating polynomials R(n,x) = Sum_{i=1..n} a(n,i) * x^i satisfy the recurrence equation R(n+1,x) = R(n,x) + Sum_{k=0..n-1} binomial(n-1,k) * R(k+1,x) * R(n-k,x) for n >= 1 with initial value R(1,x) = x. - Werner Schulte, Jun 17 2021

Definition corrected by Li Guo, Dec 16 2006
Typo in link corrected by Johannes W. Meijer, Oct 17 2009
Error in title corrected by Johannes W. Meijer, Sep 24 2010
Edited by M. F. Hasler, Oct 29 2014

A032443 a(n) = Sum_{i=0..n} binomial(2*n, i).

Original entry on oeis.org

1, 3, 11, 42, 163, 638, 2510, 9908, 39203, 155382, 616666, 2449868, 9740686, 38754732, 154276028, 614429672, 2448023843, 9756737702, 38897306018, 155111585372, 618679078298, 2468152192772, 9848142504068, 39301087452632, 156861290196878, 626155256640188

Offset: 0

J. H. Conway, Dec 11 1999

Array interpretation: first row is filled with 1's, first column with powers of 2, b(i,j) = b(i-1,j) + b(i,j-1); then a(n) = b(n,n). - Benoit Cloitre, Apr 01 2002
   1   1   1   1   1   1   1 ...
   2   3   4   5   6   7   8 ...
   4   7  11  16  22 ...
   8  15  26  42  64 ...
  16  31  ..  99 163 ...
From Gary W. Adamson, Dec 27 2008: (Start)
This is an analog of the Catalan sequence: Let M denote an infinite Cartan matrix (-1's in the super and subdiagonals and (2,2,2,...) in the main diagonal which we modify to (1,2,2,2,...)). Then A000108 can be generated by accessing the leftmost term in M^n * [1,0,0,0,...]. Change the operation to M^n * [1,2,3,...] to get this sequence. Or, take iterates M * [1,2,3,...] -> M * ANS, -> M * ANS,...; accessing the leftmost term. (End)
Convolved with the Catalan sequence, A000108: (1, 1, 2, 5, 14, ...) = powers of 4, A000302: (1, 4, 16, 64, ...). - Gary W. Adamson, May 15 2009
Row sums of A094527. - Paul Barry, Sep 07 2009
Hankel transform of the aeration of this sequence is C(n+2, 2). - Paul Barry, Sep 26 2009
Number of 4-ary words of length n in which the number of 1's does not exceed the number of 0's. - David Scambler, Aug 14 2012
Number of options available to a voter who has up to n (0-n) votes to allot among 2n candidates. - Lorraine Lee, Apr 27 2019
2*a(n-1) is the number of all triangulations that can be obtained from a Möbius strip with n chosen points on its edge. See Bazier-Matte et al. - Michel Marcus, Sep 15 2020

			G.f. = 1 + 3*x + 11*x^2 + 42*x^3 + 163*x^4 + 638*x^5 + 2510*x^6 + 9908*x^7 + ...
According to the second formula, we see the fourth row of Pascal's triangle has the terms 1,4,6,4,1 and the partial sums are 1,5,11,15,16. Using these we get 1*1 + 4*5 + 6*11 + 4*15*1*16 = 1 + 20 + 66 + 60 + 16 = 163 = a(4). - _J. M. Bergot_, Apr 29 2014

D. Phulara and L. W. Shapiro, Descendants in ordered trees with a marked vertex, Congressus Numerantium, 205 (2011), 121-128.

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Véronique Bazier-Matte, Ruiyan Huang, and Hanyi Luo, Number of Triangulations of a Möbius Strip, arXiv:2009.05785 [math.CO], 2020.
A. Bernini, F. Disanto, R. Pinzani and S. Rinaldi, Permutations defining convex permutominoes, J. Int. Seq. 10 (2007) # 07.9.7.
Celeste Damiani, Paul Martin, and Eric C. Rowell, Generalisations of Hecke algebras from Loop Braid Groups, arXiv:2008.04840 [math.GT], 2020.
M. Klazar, Twelve countings with rooted plane trees, European Journal of Combinatorics 18 (1997), 195-210; Addendum, 18 (1997), 739-740.
Mircea Merca, A Note on Cosine Power Sums J. Integer Sequences, Vol. 15 (2012), Article 12.5.3.

Binomial transform of A027914. Hankel transform is {1, 2, 3, 4, ..., n, ...}.

Cf. A000108, A000302, A110166.

[(4^n+Binomial(2*n, n))/2:n in [0.. 30]]; // Vincenzo Librandi, Oct 18 2014

A032443:=n->(4^n + binomial(2*n, n))/2; seq(A032443(n), n=0..30); # Wesley Ivan Hurt, Apr 15 2014

a[ n_] := If[ n<-3, 0, (4^n + Binomial[2 n, n]) / 2]; Table[a[n], {n, 0, 30}] (* Michael Somos, Jan 25 2014 *)

A032443 = n->sum(i=0,n,binomial(2*n,i)) \\ M. F. Hasler, Jan 02 2014

A032443 = n->binomial(2*n-1,n)+4^n\2 \\ M. F. Hasler, Jan 02 2014

{a(n) = if( n<-3, 0, (4^n + binomial(2*n, n)) / 2)} /* Michael Somos, Jan 25 2014 */

a(n) = (4^n+binomial(2*n, n))/2. - David W. Wilson
a(n) = Sum_{0 <= i_1 <= i_2 <= n} binomial(n, i_2) * binomial(n, i_1+i_2). - Benoit Cloitre, Oct 14 2004
Sequence with interpolated zeros has a(n) = Sum_{k=0..floor(n/2)} (if((n-2k) mod 2)=0, C(n, k), 0). - Paul Barry, Jan 14 2005
a(n) = Sum_{k=0..n} C(n+k-1, k)*2^(n-k). - Paul Barry, Sep 28 2007
E.g.f.: exp(2*x)*(exp(2*x) + BesselI(0,2*x))/2. For BesselI see Abramowitz-Stegun (reference and link under A008277), p. 375, eq. 9.6.10. See also A000984 for its e.g.f. - Wolfdieter Lang, Jan 16 2012
From Sergei N. Gladkovskii, Aug 13 2012: (Start)
G.f.: (1/sqrt(1-4*x) + 1/(1-4*x))/2 = G(0)/2 where G(k) = 1 + ((2*k)!)/(k!)^2/(4^k - 4*x*(16^k)/( 4*x*(4^k) + ((2*k)!)/(k!)^2/G(k+1))); (continued fraction).
E.g.f.: G(0)/2 where G(k) = 1 + ((2*k)!)/(k!)^2/(4^k - 4*x*(16^k)/( 4*x*(4^k) + (k+1)*((2*k)!)/(k!)^2/G(k+1))); (continued fraction).
(End)
O.g.f.: (1 - x*(2 + c(x)))/(1 - 4*x)^(3/2), with c the o.g.f. of A000108 (Catalan). - Wolfdieter Lang, Nov 22 2012
D-finite with recurrence: n*a(n) + 2*(-4*n+3)*a(n-1) + 8*(2*n-3)*a(n-2) = 0. - R. J. Mathar, Dec 04 2012
a(n) = binomial(2*n-1,n) + floor(4^n/2), or a(n+1) = A001700(n) + A004171(n), for all n >= 0. See A000346 for the difference. - M. F. Hasler, Jan 02 2014
0 = a(n) * (256*a(n+1) - 224*a(n+2) + 40*a(n+3)) + a(n+1) * (-32*a(n+1) + 56*a(n+2) - 14*a(n+3)) + a(n+2) * (-2*a(n+2) + a(n+3)) if n > -4. - Michael Somos, Jan 25 2014
a(n) = coefficient of x^n in (4*x + 1 / (1 + x))^n. - Michael Somos, Jan 25 2014
Binomial transform is A110166. - Michael Somos, Jan 25 2014
Asymptotics: a(n) ~ 2^(2*n-1)*(1+1/sqrt(Pi*n)). - Fung Lam, Apr 13 2014
Self-convolution is A240879. - Fung Lam, Apr 13 2014
a(0) = 1, a(n+1) = A001700(n) + 2^(2n+1). - Philippe Deléham, Oct 11 2014
exp( Sum_{n >= 1} a(n)*x^n/n ) = 1 + 3*x + 10*x^2 + 35*x^3 + 126*x^4 + ... is the o.g.f. for A001700. - Peter Bala, Jul 21 2015
a(n) = 4*a(n-1) - A000108(n-1). - Bob Selcoe, Apr 02 2017
a(n) = [x^n] 1/((1 - x)^n*(1 - 2*x)). - Ilya Gutkovskiy, Oct 12 2017
The Gauss congruences a(n*p^k) == a(n*p^(k-1)) (mod p^k) hold for prime p and positive integers n and k. - Peter Bala, Jan 08 2022
O.g.f.: ( hypergeom([1/4, 3/4], [1/2], 4*x) )^2. - Peter Bala, Mar 04 2022
a(n) = binomial(2*n, n) * hypergeom([1, -n], [n + 1], -1). - Peter Luschny, Oct 06 2023

A094638 Triangle read by rows: T(n,k) = |s(n,n+1-k)|, where s(n,k) are the signed Stirling numbers of the first kind A008276 (1 <= k <= n; in other words, the unsigned Stirling numbers of the first kind in reverse order).

Original entry on oeis.org

1, 1, 1, 1, 3, 2, 1, 6, 11, 6, 1, 10, 35, 50, 24, 1, 15, 85, 225, 274, 120, 1, 21, 175, 735, 1624, 1764, 720, 1, 28, 322, 1960, 6769, 13132, 13068, 5040, 1, 36, 546, 4536, 22449, 67284, 118124, 109584, 40320, 1, 45, 870, 9450, 63273, 269325, 723680, 1172700, 1026576, 362880

Offset: 1

André F. Labossière, May 17 2004

Triangle of coefficients of the polynomial (x+1)(x+2)...(x+n), expanded in decreasing powers of x. - T. D. Noe, Feb 22 2008
Row n also gives the number of permutation of 1..n with complexity 0,1,...,n-1. See the comments in A008275. - N. J. A. Sloane, Feb 08 2019
T(n,k) is the number of deco polyominoes of height n and having k columns. A deco polyomino is a directed column-convex polyomino in which the height, measured along the diagonal, is attained only in the last column. Example: T(2,1)=1 and T(2,2)=1 because the deco polyominoes of height 2 are the vertical and horizontal dominoes, having, respectively, 1 and 2 columns. - Emeric Deutsch, Aug 14 2006
Sum_{k=1..n} k*T(n,k) = A121586. - Emeric Deutsch, Aug 14 2006
Let the triangle U(n,k), 0 <= k <= n, read by rows, be given by [1,0,1,0,1,0,1,0,1,0,1,...] DELTA [1,1,2,2,3,3,4,4,5,5,6,...] where DELTA is the operator defined in A084938; then T(n,k) = U(n-1,k-1). - Philippe Deléham, Jan 06 2007
From Tom Copeland, Dec 15 2007: (Start)
Consider c(t) = column vector(1, t, t^2, t^3, t^4, t^5, ...).
Starting at 1 and sampling every integer to the right, we obtain (1,2,3,4,5,...). And T * c(1) = (1, 1*2, 1*2*3, 1*2*3*4,...), giving n! for n > 0. Call this sequence the right factorial (n+)!.
Starting at 1 and sampling every integer to the left, we obtain (1,0,-1,-2,-3,-4,-5,...). And T * c(-1) = (1, 1*0, 1*0*-1, 1*0*-1*-2,...) = (1, 0, 0, 0, ...), the left factorial (n-)!.
Sampling every other integer to the right, we obtain (1,3,5,7,9,...). T * c(2) = (1, 1*3, 1*3*5, ...) = (1,3,15,105,945,...), giving A001147 for n > 0, the right double factorial, (n+)!!.
Sampling every other integer to the left, we obtain (1,-1,-3,-5,-7,...). T * c(-2) = (1, 1*-1, 1*-1*-3, 1*-1*-3*-5,...) = (1,-1,3,-15,105,-945,...) = signed A001147, the left double factorial, (n-)!!.
Sampling every 3 steps to the right, we obtain (1,4,7,10,...). T * c(3) = (1, 1*4, 1*4*7,...) = (1,4,28,280,...), giving A007559 for n > 0, the right triple factorial, (n+)!!!.
Sampling every 3 steps to the left, we obtain (1,-2,-5,-8,-11,...), giving T * c(-3) = (1, 1*-2, 1*-2*-5, 1*-2*-5*-8,...) = (1,-2,10,-80,880,...) = signed A008544, the left triple factorial, (n-)!!!.
The list partition transform A133314 of [1,T * c(t)] gives [1,T * c(-t)] with all odd terms negated; e.g., LPT[1,T*c(2)] = (1,-1,-1,-3,-15,-105,-945,...) = (1,-A001147). And e.g.f. for [1,T * c(t)] = (1-xt)^(-1/t).
The above results hold for t any real or complex number. (End)
Let R_n(x) be the real and I_n(x) the imaginary part of Product_{k=0..n} (x + I*k). Then, for n=1,2,..., we have R_n(x) = Sum_{k=0..floor((n+1)/2)}(-1)^k*Stirling1(n+1,n+1-2*k)*x^(n+1-2*k), I_n(x) = Sum_{k=0..floor(n/2)}(-1)^(k+1)*Stirling1(n+1,n-2*k)*x^(n-2*k). - Milan Janjic, May 11 2008
T(n,k) is also the number of permutations of n with "reflection length" k (i.e., obtained from 12..n by k not necessarily adjacent transpositions). For example, when n=3, 132, 213, 321 are obtained by one transposition, while 231 and 312 require two transpositions. - Kyle Petersen, Oct 15 2008
From Tom Copeland, Nov 02 2010: (Start)
[x^(y+1) D]^n = x^(n*y) [T(n,1)(xD)^n + T(n,2)y (xD)^(n-1) + ... + T(n,n)y^(n-1)(xD)], with D the derivative w.r.t. x.
E.g., [x^(y+1) D]^4 = x^(4*y) [(xD)^4 + 6 y(xD)^3 + 11 y^2(xD)^2 + 6 y^3(xD)].
(xD)^m can be further expanded in terms of the Stirling numbers of the second kind and operators of the form x^j D^j. (End)
With offset 0, 0 <= k <= n: T(n,k) is the sum of products of each size k subset of {1,2,...,n}. For example, T(3,2) = 11 because there are three subsets of size two: {1,2},{1,3},{2,3}. 1*2 + 1*3 + 2*3 = 11. - Geoffrey Critzer, Feb 04 2011
The Kn11, Fi1 and Fi2 triangle sums link this triangle with two sequences, see the crossrefs. For the definitions of these triangle sums see A180662. The mirror image of this triangle is A130534. - Johannes W. Meijer, Apr 20 2011
T(n+1,k+1) is the elementary symmetric function a_k(1,2,...,n), n >= 0, k >= 0, (a_0(0):=1). See the T. D. Noe and Geoffrey Critzer comments given above. For a proof see the Stanley reference, p. 19, Second Proof. - Wolfdieter Lang, Oct 24 2011
Let g(t) = 1/d(log(P(j+1,-t)))/dt (see Tom Copeland's 2007 formulas). The Mellin transform (t to s) of t*Dirac[g(t)] gives Sum_{n=1..j} n^(-s), which as j tends to infinity gives the Riemann zeta function for Re(s) > 1. Dirac(x) is the Dirac delta function. The complex contour integral along a circle of radius 1 centered at z=1 of z^s/g(z) gives the same result. - Tom Copeland, Dec 02 2011
Rows are coefficients of the polynomial expansions of the Pochhammer symbol, or rising factorial, Pch(n,x) = (x+n-1)!/(x-1)!. Expansion of Pch(n,xD) = Pch(n,Bell(.,:xD:)) in a polynomial with terms :xD:^k=x^k*D^k gives the Lah numbers A008297. Bell(n,x) are the unsigned Bell polynomials or Stirling polynomials of the second kind A008277. - Tom Copeland, Mar 01 2014
From Tom Copeland, Dec 09 2016: (Start)
The Betti numbers, or dimension, of the pure braid group cohomology. See pp. 12 and 13 of the Hyde and Lagarias link.
Row polynomials and their products appear in presentation of the Jack symmetric functions of R. Stanley. See Copeland link on the Witt differential generator.
(End)
From Tom Copeland, Dec 16 2019: (Start)
The e.g.f. given by Copeland in the formula section appears in a combinatorial Dyson-Schwinger equation of quantum field theory in Yeats in Thm. 2 on p. 62 related to a Hopf algebra of rooted trees. See also the Green function on p. 70.
Per comments above, this array contains the coefficients in the expansion in polynomials of the Euler, or state number, operator xD of the rising factorials Pch(n,xD) = (xD+n-1)!/(xD-1)! = x [:Dx:^n/n!]x^{-1} = L_n^{-1}(-:xD:), where :Dx:^n = D^n x^n and :xD:^n = x^n D^n. The polynomials L_n^{-1} are the Laguerre polynomials of order -1, i.e., normalized Lah polynomials.
The Witt differential operators L_n = x^(n+1) D and the row e.g.f.s appear in Hopf and dual Hopf algebra relations presented by Foissy. The Witt operators satisfy L_n L_k - L_k L_n = (k-n) L_(n+k), as for the dual Hopf algebra. (End)

			Triangle starts:
  1;
  1,  1;
  1,  3,  2;
  1,  6, 11,  6;
  1, 10, 35, 50, 24;
...

M. Miyata and J. W. Son, On the complexity of permutations and the metric space of bijections, Tensor, 60 (1998), No. 1, 109-116 (MR1768839).
Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 18, equations 18:4:2 - 18:4:8 at page 151.
R. P. Stanley, Enumerative Combinatorics, Vol. 1, Cambridge University Press, 1997.

T. D. Noe, Rows n = 1..51 of triangle, flattened
E. Barcucci, A. Del Lungo and R. Pinzani, "Deco" polyominoes, permutations and random generation, Theoretical Computer Science, 159, 1996, 29-42.
F. Bergeron, Ph. Flajolet and B. Salvy, Varieties of Increasing Trees, Lecture Notes in Computer Science vol. 581, ed. J.-C. Raoult, Springer 1992, pp. 24-48. Added Mar 01 2014
F. Bergeron, Philippe Flajolet and Bruno Salvy, Varieties of increasing trees, HAL, Rapport De Recherche Inria. Added Mar 01 2014
Tom Copeland, Addendum to Mathemagical Forests
Tom Copeland, A Class of Differential Operators and the Stirling Numbers
Tom Copeland, Mellin Interpolation of Differential Ops and Associated Infinigens and Appell Polynomials: The Ordered, Laguerre, and Scherk-Witt-Lie Diff Ops
Tom Copeland, Witt differential generator for special Jack symmetric functions / polynomials
FindStat - Combinatorial Statistic Finder, The absolute length of a permutation
L. Foissy, Faa di Bruno subalgebras of the Hopf algebra of planar trees from combinatorial Dyson-Schwinger equations, arXiv:0707.1204 [math.RA], (2007).
O. Furdui, T. Trif, On the Summation of Certain Iterated Series, J. Int. Seq. 14 (2011) #11.6.1
F. Hivert, J.-C. Novelli and J.-Y. Thibon, The Algebra of Binary Search Trees, Theoretical Computer Science, 339 (2005), 129-165.
T. Hyde and J. Lagarias Polynomial splitting measures and cohomology of the pure braid group, arXiv preprint arXiv:1604.05359 [math.RT], 2016.
MathOverflow, Motivation of Virasoro algebra, an answer by Tom Copeland to an MO question posed in 2012.
Romeo Mestrovic, Lucas' theorem: its generalizations, extensions and applications (1878--2014), arXiv preprint arXiv:1409.3820 [math.NT], 2014.
Robert E. Moritz, On the sum of products of n consecutive integers, Univ. Washington Publications in Math., 1 (No. 3, 1926), 44-49 [Annotated scanned copy]
M. D. Schmidt, Generalized j-Factorial Functions, Polynomials, and Applications , J. Int. Seq. 13 (2010), 10.6.7.
M. Z. Spivey, On Solutions to a General Combinatorial Recurrence, J. Int. Seq. 14 (2011) # 11.9.7.
K. Yeats, A Combinatorial Perspective on Quantum Field Theory, SpringerBriefs in Mathematical Physics, Vol. 15, 2017.

A008276 gives the (signed) Stirling numbers of the first kind.
Cf. A000108, A014137, A001246, A033536, A000984, A094639, A006134, A082894, A002897, A079727, A000217 (2nd column), A000914 (3rd column), A001303 (4th column), A000915 (5th column), A053567 (6th column), A000142 (row sums).
Triangle sums (see the comments): A124380 (Kn11), A001710 (Fi1, Fi2). - Johannes W. Meijer, Apr 20 2011
Cf. A121586, A130534, A143491.

Flat(List([1..10], n-> List([1..n], k-> Stirling1(n,n-k+1) ))); # G. C. Greubel, Dec 29 2019

a094638 n k = a094638_tabl !! (n-1) !! (k-1)
a094638_row n = a094638_tabl !! (n-1)
a094638_tabl = map reverse a130534_tabl
-- Reinhard Zumkeller, Aug 01 2014

[(-1)^(k+1)*StirlingFirst(n,n-k+1): k in [1..n], n in [1..10]]; // G. C. Greubel, Dec 29 2019

T:=(n,k)->abs(Stirling1(n,n+1-k)): for n from 1 to 10 do seq(T(n,k),k=1..n) od; # yields sequence in triangular form. # Emeric Deutsch, Aug 14 2006

Table[CoefficientList[Series[Product[1 + i x, {i,n}], {x,0,20}], x], {n,0,6}] (* Geoffrey Critzer, Feb 04 2011 *)
Table[Abs@StirlingS1[n, n-k+1], {n, 10}, {k, n}]//Flatten (* Michael De Vlieger, Aug 29 2015 *)

create_list(abs(stirling1(n+1,n-k+1)),n,0,10,k,0,n); /* Emanuele Munarini, Jun 01 2012 */

{T(n,k)=if(n<1 || k>n,0,(n-1)!*polcoeff(polcoeff(x*y/(1 - x*y+x*O(x^n))^(1 + 1/y),n,x),k,y))} /* Paul D. Hanna, Jul 21 2011 */

[[stirling_number1(n, n-k+1) for k in (1..n)] for n in (1..10)] # G. C. Greubel, Dec 29 2019

With P(n,t) = Sum_{k=0..n-1} T(n,k+1) * t^k = 1*(1+t)*(1+2t)...(1+(n-1)*t) and P(0,t)=1, exp[P(.,t)*x] = (1-tx)^(-1/t). T(n,k+1) = (1/k!) (D_t)^k (D_x)^n [ (1-tx)^(-1/t) - 1 ] evaluated at t=x=0. (1-tx)^(-1/t) - 1 is the e.g.f. for a plane m-ary tree when t=m-1. See Bergeron et al. in "Varieties of Increasing Trees". - Tom Copeland, Dec 09 2007
First comment and formula above rephrased as o.g.f. for row n: Product_{i=0...n} (1+i*x). - Geoffrey Critzer, Feb 04 2011
n-th row polynomials with alternate signs are the characteristic polynomials of the (n-1)x(n-1) matrices with 1's in the superdiagonal, (1,2,3,...) in the main diagonal, and the rest zeros. For example, the characteristic polynomial of [1,1,0; 0,2,1; 0,0,3] is x^3 - 6*x^2 + 11*x - 6. - Gary W. Adamson, Jun 28 2011
E.g.f.: A(x,y) = x*y/(1 - x*y)^(1 + 1/y) = Sum_{n>=1, k=1..n} T(n,k)*x^n*y^k/(n-1)!. - Paul D. Hanna, Jul 21 2011
With F(x,t) = (1-t*x)^(-1/t) - 1 an e.g.f. for the row polynomials P(n,t) of A094638 with P(0,t)=0, G(x,t)= [1-(1+x)^(-t)]/t is the comp. inverse in x. Consequently, with H(x,t) = 1/(dG(x,t)/dx) = (1+x)^(t+1),
  P(n,t) = [(H(x,t)*d/dx)^n] x evaluated at x=0; i.e.,
  F(x,t) = exp[x*P(.,t)] = exp[x*H(u,t)*d/du] u, evaluated at u = 0.
  Also, dF(x,t)/dx = H(F(x,t),t). - Tom Copeland, Sep 20 2011
T(n,k) = |A008276(n,k)|. - R. J. Mathar, May 19 2016
The row polynomials of this entry are the reversed row polynomials of A143491 multiplied by (1+x). E.g., (1+x)(1 + 5x + 6x^2) = (1 + 6x + 11x^2 + 6x^3). - Tom Copeland, Dec 11 2016
Regarding the row e.g.f.s in Copeland's 2007 formulas, e.g.f.s for A001710, A001715, and A001720 give the compositional inverses of the e.g.f. here for t = 2, 3, and 4 respectively. - Tom Copeland, Dec 28 2019

Edited by Emeric Deutsch, Aug 14 2006

A217900 O.g.f.: Sum_{n>=0} n^n * (n+1)^(n-1) * exp(-n(n+1)x) * x^n / n!.

Original entry on oeis.org

1, 1, 4, 38, 576, 12052, 322848, 10564304, 408903680, 18288706544, 928575662400, 52780935007968, 3321208845997056, 229232635832433664, 17221699990084108288, 1399139700462119135232, 122235936429355565580288, 11428226675376971405577984, 1138551595285580854471388160

Offset: 0

Paul D. Hanna, Oct 14 2012

nonn

Compare the g.f. to the LambertW identity:
1 = Sum_{n>=0} (n+1)^(n-1) * exp(-(n+1)*x) * x^n/n!.
More generally, if we define a(n) for fixed integers m, t, and s>=0, by:
(0) Sum_{n>=0} m * n^(s*n) * (n*t+m)^(n-1) * exp(-n^s*(n*t+m)*x) * x^n/n! = Sum_{n>=0} a(n)*x^n
then the coefficients a(n) are integral and may be expressed by:
(1) a(n) = 1/n! * Sum_{k=0..n} m*(-1)^(n-k)*binomial(n,k) * k^(s*n) * (k*t+m)^(n-1).
(2) a(n) = 1/n! * [x^n] Sum_{k>=0} m*k^(s*k)*(k*t+m)^(k-1)*x^k / (1 + k^s*(k*t+m)*x)^(k+1).
(3) a(n) = 1/t^((s-1)*n) * [x^(s*n)] 1 + m*x*(1+m*x)^(n-1) / Product_{k=1..n} (1-k*t*x).
(4) a(n) = 1/t^((s-1)*n) * [x^(s*n)] 1 + m*x*(1-m*x)^(s*n) / Product_{k=1..n} (1-(k*t+m)*x).

			O.g.f.: A(x) = 1 + x + 4*x^2 + 38*x^3 + 576*x^4 + 12052*x^5 + 322848*x^6 +...
where
A(x) = 1 + 1^1*2^0*x*exp(-1*2*x) + 2^2*3^1*exp(-2*3*x)*x^2/2! + 3^3*4^2*exp(-3*4*x)*x^3/3! + 4^4*5^3*exp(-4*5*x)*x^4/4! + 5^5*6^4*exp(-5*6*x)*x^5/5! +...
simplifies to a power series in x with integer coefficients.

G. C. Greubel, Table of n, a(n) for n = 0..345

Cf. A078739, A217899, A217901, A217902, A217903, A217904, A217905, A217910, A217913, A218300.

a[n_] := 1/n!*Sum[(-1)^(n-k)*Binomial[n, k]*k^n*(k+1)^(n-1), {k, 0, n}]; a[0]=1; Table[a[n], {n, 0, 18}] (* Jean-François Alcover, Mar 06 2013 *)

{a(n)=polcoeff(sum(m=0,n,m^m*(m+1)^(m-1)*x^m*exp(-m*(m+1)*x+x*O(x^n))/m!),n)}

{a(n)=(1/n!)*polcoeff(sum(k=0, n, k^k*(k+1)^(k-1)*x^k/(1+k*(k+1)*x +x*O(x^n))^(k+1)), n)}

{a(n)=1/n!*sum(k=0,n, (-1)^(n-k)*binomial(n,k)*k^n*(k+1)^(n-1))}

{a(n)=polcoeff(1+x*(1+x)^(n-1)/prod(k=0, n, 1-k*x +x*O(x^n)), n)}

{a(n)=polcoeff(1+x*(1-x)^n/prod(k=0, n, 1-(k+1)*x +x*O(x^n)), n)}
for(n=0,30,print1(a(n),", "))

{Stirling2(n, k)=n!*polcoeff(((exp(x+x*O(x^n))-1)^k)/k!, n)}
{a(n)=if(n==0,1,sum(k=0,n-1, binomial(n-1,k) * Stirling2(2*n-k-1,n)))} \\ Paul D. Hanna, Nov 13 2012
/* PARI Programs for the General Case (START) ...................... */

{a(n,m=1,t=1,s=1)=polcoeff(sum(k=0, n, m*k^(s*k)*(t*k+m)^(k-1)*exp(-k^s*(t*k+m)*x+x*O(x^n))*x^k/k!), n)}

{a(n,m=1,t=1,s=1)=(1/n!)*polcoeff(sum(k=0, n, m*k^(s*k)*(t*k+m)^(k-1)*x^k/(1+k^s*(t*k+m)*x +x*O(x^n))^(k+1)), n)}

{a(n,m=1,t=1,s=1)=1/n!*sum(k=0, n, m*(-1)^(n-k)*binomial(n, k)*k^(s*n)*(t*k+m)^(n-1))}

{a(n,m=1,t=1,s=1)=(1/t^((s-1)*n))*polcoeff(1+m*x*(1+m*x)^(n-1)/prod(k=0, n, 1-t*k*x +x*O(x^(s*n))), s*n)}

{a(n,m=1,t=1,s=1)=(1/t^((s-1)*n))*polcoeff(1+m*x*(1-m*x)^(s*n)/prod(k=0, n, 1-(t*k+m)*x +x*O(x^(s*n))), s*n)}
/* (END) ........................................................... */

a(n) = 1/n! * Sum_{k=0..n} (-1)^(n-k)*binomial(n,k) * k^n * (k+1)^(n-1).
a(n) = 1/n! * [x^n] Sum_{k>=0} k^k*(k+1)^(k-1)*x^k / (1 + k*(k+1)*x)^(k+1).
a(n) = [x^n] 1 + x*(1+x)^(n-1) / Product_{k=1..n} (1-k*x).
a(n) = [x^n] 1 + x*(1-x)^(n-1) / Product_{k=1..n} (1-(k+1)*x).
a(n) = A078739(n,n) for n>=1.
a(n) = Sum_{k=0..n-1} binomial(n-1,k) * Stirling2(2*n-k-1,n) for n>0, where Stirling2(n,k) = A008277(n,k). - Paul D. Hanna, Nov 13 2012
a(n) ~ 2^(2*n-1) * n^(n-3/2) / (sqrt(Pi*(1-c)) * exp(n) * (2-c)^(n-1) * c^(n+1/2)), where c = -LambertW(-2*exp(-2)) = 0.4063757399599599... = 2*A106533. - Vaclav Kotesovec, May 09 2014

A001296 4-dimensional pyramidal numbers: a(n) = (3n+1)binomial(n+2, 3)/4. Also Stirling2(n+2, n).

Original entry on oeis.org

0, 1, 7, 25, 65, 140, 266, 462, 750, 1155, 1705, 2431, 3367, 4550, 6020, 7820, 9996, 12597, 15675, 19285, 23485, 28336, 33902, 40250, 47450, 55575, 64701, 74907, 86275, 98890, 112840, 128216, 145112, 163625, 183855, 205905, 229881, 255892, 284050, 314470

Offset: 0

N. J. A. Sloane

Permutations avoiding 12-3 that contain the pattern 31-2 exactly once.
Kekulé numbers for certain benzenoids. - Emeric Deutsch, Nov 18 2005
Partial sums of A002411. - Jonathan Vos Post, Mar 16 2006
If Y is a 3-subset of an n-set X then, for n>=6, a(n-5) is the number of 6-subsets of X having at least two elements in common with Y. - Milan Janjic, Nov 23 2007
Starting with 1 = binomial transform of [1, 6, 12, 10, 3, 0, 0, 0, ...]. Equals row sums of triangle A143037. - Gary W. Adamson, Jul 18 2008
Rephrasing the Perry formula of 2003: a(n) is the sum of all products of all two  numbers less than or equal to n, including the squares. Example: for n=3 the sum of these products is 1*1 + 1*2 + 1*3 + 2*2 + 2*3 + 3*3 = 25. - J. M. Bergot, Jul 16 2011
Half of the partial sums of A011379. [Jolley, Summation of Series, Dover (1961), page 12 eq (66).] - R. J. Mathar, Oct 03 2011
Also the number of (w,x,y,z) with all terms in {1,...,n+1} and w < x >= y > z (see A211795). - Clark Kimberling, May 19 2012
Convolution of A000027 with A000326. - Bruno Berselli, Dec 06 2012
This sequence is related to A000292 by a(n) = n*A000292(n) - Sum_{i=0..n-1} A000292(i) for n>0. - Bruno Berselli, Nov 23 2017
a(n-2) is the maximum number of intersections made from the perpendicular bisectors of all pair combinations of n points. - Ian Tam, Dec 22 2020

			G.f. = x + 7*x^2 + 25*x^3 + 65*x^4 + 140*x^5 + 266*x^6 + 462*x^7 + 750*x^8 + 1155*x^9 + ...

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 835.
A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 195.
L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 227, #16.
S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988 (see p. 166, Table 10.4/I/3).
F. N. David, M. G. Kendall and D. E. Barton, Symmetric Function and Allied Tables, Cambridge, 1966, p. 223.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe, Table of n, a(n) for n = 0..1000
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
P. Aluffi, Degrees of projections of rank loci, arXiv:1408.1702 [math.AG], 2014. ["After compiling the results of many explicit computations, we noticed that many of the numbers d_{n,r,S} appear in the existing literature in contexts far removed from the enumerative geometry of rank conditions; we owe this surprising (to us) observation to perusal of [Slo14]."]
S. Butler and P. Karasik, A note on nested sums, J. Int. Seq. 13 (2010), 10.4.4, page 5.
M. Griffiths, Remodified Bessel Functions via Coincidences and Near Coincidences, Journal of Integer Sequences, Vol. 14 (2011), Article 11.7.1.
L. Hogben, Choice and Chance by Cardpack and Chessboard, Vol. 1, Chanticleer Press, NY, 1950, p. 36.
C. Krishnamachaki, The operator (xD)^n, J. Indian Math. Soc., 15 (1923),3-4. [Annotated scanned copy]
Feihu Liu, Guoce Xin, and Chen Zhang, Ehrhart Polynomials of Order Polytopes: Interpreting Combinatorial Sequences on the OEIS, arXiv:2412.18744 [math.CO], 2024. See p. 29.
T. Mansour, Restricted permutations by patterns of type 2-1, arXiv:math/0202219 [math.CO], 2002.
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992
Eric Weisstein's World of Mathematics, Stirling numbers of the 2nd kind.
Index to sequences related to pyramidal numbers
Index entries for linear recurrences with constant coefficients, signature (5,-10,10,-5,1).

a(n)=f(n, 2) where f is given in A034261.
Cf. A000217, A000326, A001297, A001298, A002411, A008277, A008517, A094262.
a(n)= A093560(n+3, 4), (3, 1)-Pascal column.
Cf. A220212 for a list of sequences produced by the convolution of the natural numbers with the k-gonal numbers.
Cf. similar sequences listed in A241765 and A254142.
Cf. A000914.

/* A000027 convolved with A000326: */ A000326:=func; [&+[(n-i+1)*A000326(i): i in [0..n]]: n in [0..40]]; // Bruno Berselli, Dec 06 2012

[(3*n+1)*Binomial(n+2,3)/4: n in [0..40]]; // Vincenzo Librandi, Jul 30 2014

A001296:=-(1+2*z)/(z-1)**5; # Simon Plouffe in his 1992 dissertation for sequence without the leading zero

Table[n*(1+n)*(2+n)*(1+3*n)/24, {n, 0, 100}]
CoefficientList[Series[x (1 + 2 x)/(1 - x)^5, {x, 0, 40}], x] (* Vincenzo Librandi, Jul 30 2014 *)
Table[StirlingS2[n+2, n], {n, 0, 40}] (* Jean-François Alcover, Jun 24 2015 *)
Table[ListCorrelate[Accumulate[Range[n]],Range[n]],{n,0,40}]//Flatten (* or *) LinearRecurrence[{5,-10,10,-5,1},{0,1,7,25,65},40] (* Harvey P. Dale, Aug 14 2017 *)

t(n)=n*(n+1)/2
for(i=1,30,print1(","sum(j=1,i,j*t(j))))

{a(n) = n * (n+1) * (n+2) * (3*n+1) / 24}; /* Michael Somos, Sep 04 2017 */

[stirling_number2(n+2,n) for n in range(0,38)] # Zerinvary Lajos, Mar 14 2009

a(n) = n*(1+n)*(2+n)*(1+3*n)/24. - T. D. Noe, Jan 21 2008
G.f.: x*(1+2*x)/(1-x)^5. - Paul Barry, Jul 23 2003
a(n) = Sum_{j=0..n} j*A000217(j). - Jon Perry, Jul 28 2003
E.g.f. with offset -1: exp(x)*(1*(x^2)/2! + 4*(x^3)/3! + 3*(x^4)/4!). For the coefficients [1, 4, 3] see triangle A112493.
E.g.f. x*exp(x)*(24 + 60*x + 28*x^2 + 3*x^3)/24 (above e.g.f. differentiated).
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4) + 3. - Kieren MacMillan, Sep 29 2008
a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5). - Jaume Oliver Lafont, Nov 23 2008
O.g.f. is D^2(x/(1-x)) = D^3(x), where D is the operator x/(1-x)*d/dx. - Peter Bala, Jul 02 2012
a(n) = A153978(n)/2. - J. M. Bergot, Aug 09 2013
a(n) = A002817(n) + A000292(n-1). - J. M. Bergot, Aug 29 2013; [corrected by Cyril Damamme, Feb 26 2018]
a(n) = A000914(n+1) - 2 * A000330(n+1). - Antal Pinter, Dec 31 2015
a(n) = A080852(3,n-1). - R. J. Mathar, Jul 28 2016
a(n) = 1*(1+2+...+n) + 2*(2+3+...+n) + ... + n*n. For example, a(6) = 266 = 1(1+2+3+4+5+6) + 2*(2+3+4+5+6) + 3*(3+4+5+6) + 4*(4+5+6) + 5*(5+6) + 6*(6).- J. M. Bergot, Apr 20 2017
a(n) = A000914(-2-n) for all n in Z. - Michael Somos, Sep 04 2017
a(n) = A000292(n) + A050534(n+1). - Cyril Damamme, Feb 26 2018
From Amiram Eldar, Jul 02 2020: (Start)
Sum_{n>=1} 1/a(n) = (6/5) * (47 - 3*sqrt(3)*Pi - 27*log(3)).
Sum_{n>=1} (-1)^(n+1)/a(n) = (6/5) * (16*log(2) + 6*sqrt(3)*Pi - 43). (End)

A074909 Running sum of Pascal's triangle (A007318), or beheaded Pascal's triangle read by beheaded rows.

Original entry on oeis.org

1, 1, 2, 1, 3, 3, 1, 4, 6, 4, 1, 5, 10, 10, 5, 1, 6, 15, 20, 15, 6, 1, 7, 21, 35, 35, 21, 7, 1, 8, 28, 56, 70, 56, 28, 8, 1, 9, 36, 84, 126, 126, 84, 36, 9, 1, 10, 45, 120, 210, 252, 210, 120, 45, 10, 1, 11, 55, 165, 330, 462, 462, 330, 165, 55, 11

Offset: 0

Wouter Meeussen, Oct 01 2002

This sequence counts the "almost triangular" partitions of n. A partition is triangular if it is of the form 0+1+2+...+k. Examples: 3=0+1+2, 6=0+1+2+3. An "almost triangular" partition is a triangular partition with at most 1 added to each of the parts. Examples: 7 = 1+1+2+3 = 0+2+2+3 = 0+1+3+3 = 0+1+2+4. Thus a(7)=4. 8 = 1+2+2+3 = 1+1+3+3 = 1+1+2+4 = 0+2+3+3 = 0+2+2+4 = 0+1+3+4 so a(8)=6. - Moshe Shmuel Newman, Dec 19 2002
The "almost triangular" partitions are the ones cycled by the operation of "Bulgarian solitaire", as defined by Martin Gardner.
Start with A007318 - I (I = Identity matrix), then delete right border of zeros. - Gary W. Adamson, Jun 15 2007
Also the number of increasing acyclic functions from {1..n-k+1} to {1..n+2}. A function f is acyclic if for every subset B of the domain the image of B under f does not equal B. For example, T(3,1)=4 since there are exactly 4 increasing acyclic functions from {1,2,3} to {1,2,3,4,5}: f1={(1,2),(2,3),(3,4)}, f2={(1,2),(2,3),(3,5)}, f3={(1,2),(2,4),(3,5)} and f4={(1,3),(2,4),(4,5)}. - Dennis P. Walsh, Mar 14 2008
Second Bernoulli polynomials are (from A164555 instead of A027641) B2(n,x) = 1; 1/2, 1; 1/6, 1, 1; 0, 1/2, 3/2, 1; -1/30, 0, 1, 2, 1; 0, -1/6, 0, 5/3, 5/2, 1; ... . Then (B2(n,x)/A002260) = 1; 1/2, 1/2; 1/6, 1/2, 1/3; 0, 1/4, 1/2, 1/4; -1/30, 0, 1/3, 1/2, 1/5; 0, -1/12, 0, 5/12, 1/2, 1/6; ... . See (from Faulhaber 1631) Jacob Bernoulli Summae Potestatum (sum of powers) in A159688. Inverse polynomials are 1; -1, 2; 1, -3, 3; -1, 4, -6, 4; ... = A074909 with negative even diagonals. Reflected A053382/A053383 = reflected B(n,x) = RB(n,x) = 1; -1/2, 1; 1/6, -1, 1; 0, 1/2, -3/2, 1; ... . A074909 is inverse of RB(n,x)/A002260 = 1; -1/2, 1/2; 1/6, -1/2, 1/3; 0, 1/4, -1/2, 1/4; ... . - Paul Curtz, Jun 21 2010
A054143 is the fission of the polynomial sequence (p(n,x)) given by p(n,x) = x^n + x^(n-1) + ... + x + 1 by the polynomial sequence ((x+1)^n). See A193842 for the definition of fission. - Clark Kimberling, Aug 07 2011
Reversal of A135278. - Philippe Deléham, Feb 11 2012
For a closed-form formula for arbitrary left and right borders of Pascal-like triangles see A228196. - Boris Putievskiy, Aug 19 2013
For a closed-form formula for generalized Pascal's triangle see A228576. - Boris Putievskiy, Sep 09 2013
From A238363, the operator equation d/d(:xD:)f(xD)={exp[d/d(xD)]-1}f(xD) = f(xD+1)-f(xD) follows. Choosing f(x) = x^n and using :xD:^n/n! = binomial(xD,n) and (xD)^n = Bell(n,:xD:), the Bell polynomials of A008277, it follows that the lower triangular matrix [padded A074909]
  A) = [St2]*[dP]*[St1] = A048993*A132440*[padded A008275]
  B) = [St2]*[dP]*[St2]^(-1)
  C) = [St1]^(-1)*[dP]*[St1],
  where [St1]=padded A008275 just as [St2]=A048993=padded A008277 whereas [padded A074909]=A007318-I with I=identity matrix. - Tom Copeland, Apr 25 2014
T(n,k) generated by m-gon expansions in the case of odd m with "vertex to side" version or even m with "vertex to vertes" version. Refer to triangle expansions in A061777 and A101946 (and their companions for m-gons) which are "vertex to vertex" and "vertex to side" versions respectively. The label values at each iteration can be arranged as a triangle. Any m-gon can also be arranged as the same triangle with conditions: (i) m is odd and expansion is "vertex to side" version or (ii) m is even and expansion is "vertex to vertex" version. m*Sum_{i=1..k} T(n,k) gives the total label value at the n-th iteration. See also A247976. Vertex to vertex: A061777, A247618, A247619, A247620. Vertex to side: A101946, A247903, A247904, A247905. - Kival Ngaokrajang Sep 28 2014
From Tom Copeland, Nov 12 2014: (Start)
With P(n,x) = [(x+1)^(n+1)-x^(n+1)], the row polynomials of this entry, Up(n,x) = P(n,x)/(n+1) form an Appell sequence of polynomials that are the umbral compositional inverses of the Bernoulli polynomials B(n,x), i.e., B[n,Up(.,x)] = x^n = Up[n,B(.,x)] under umbral substitution, e.g., B(.,x)^n = B(n,x).
The e.g.f. for the Bernoulli polynomials is [t/(e^t - 1)] e^(x*t), and for Up(n,x) it's exp[Up(.,x)t] = [(e^t - 1)/t] e^(x*t).
Another g.f. is G(t,x) = log[(1-x*t)/(1-(1+x)*t)] = log[1 + t /(1 + -(1+x)t)] = t/(1-t*Up(.,x)) = Up(0,x)*t + Up(1,x)*t^2 + Up(2,x)*t^3 + ... = t + (1+2x)/2 t^2 + (1+3x+3x^2)/3 t^3 + (1+4x+6x^2+4x^3)/4 t^4 + ... = -log(1-t*P(.,x)), expressed umbrally.
The inverse, Ginv(t,x), in t of the g.f. may be found in A008292 from Copeland's list of formulas (Sep 2014) with a=(1+x) and b=x. This relates these two sets of polynomials to algebraic geometry, e.g., elliptic curves, trigonometric expansions, Chebyshev polynomials, and the combinatorics of permutahedra and their duals.
Ginv(t,x) = [e^((1+x)t) - e^(xt)] / [(1+x) * e^((1+x)t) - x * e^(xt)] = [e^(t/2) - e^(-t/2)] / [(1+x)e^(t/2) - x*e^(-t/2)] = (e^t - 1) / [1 + (1+x) (e^t - 1)] = t - (1 + 2 x) t^2/2! + (1 + 6 x + 6 x^2) t^3/3! - (1 + 14 x + 36 x^2 + 24 x^3) t^4/4! + ... = -exp[-Perm(.,x)t], where Perm(n,x) are the reverse face polynomials, or reverse f-vectors, for the permutahedra, i.e., the face polynomials for the duals of the permutahedra. Cf. A090582, A019538, A049019, A133314, A135278.
With L(t,x) = t/(1+t*x) with inverse L(t,-x) in t, and Cinv(t) = e^t - 1 with inverse C(t) = log(1 + t). Then Ginv(t,x) = L[Cinv(t),(1+x)] and G(t,x) = C[L[t,-(1+x)]]. Note L is the special linear fractional (Mobius) transformation.
Connections among the combinatorics of the permutahedra, simplices (cf. A135278), and the associahedra can be made through the Lagrange inversion formula (LIF) of A133437 applied to G(t,x) (cf. A111785 and the Schroeder paths A126216 also), and similarly for the LIF A134685 applied to Ginv(t,x) involving the simplicial Whitehouse complex, phylogenetic trees, and other structures. (See also the LIFs A145271 and A133932). (End)
R = x - exp[-[B(n+1)/(n+1)]D] = x - exp[zeta(-n)D] is the raising operator for this normalized sequence UP(n,x) = P(n,x) / (n+1), that is, R UP(n,x) = UP(n+1,x), where D = d/dx, zeta(-n) is the value of the Riemann zeta function evaluated at -n, and B(n) is the n-th Bernoulli number, or constant B(n,0) of the Bernoulli polynomials. The raising operator for the Bernoulli polynomials is then x + exp[-[B(n+1)/(n+1)]D]. [Note added Nov 25 2014: exp[zeta(-n)D] is abbreviation of exp(a.D) with (a.)^n = a_n = zeta(-n)]. - Tom Copeland, Nov 17 2014
The diagonals T(n, n-m), for n >= m, give the m-th iterated partial sum of the positive integers; that is A000027(n+1), A000217(n), A000292(n-1), A000332(n+1), A000389(n+1), A000579(n+1), A000580(n+1), A000581(n+1), A000582(n+1), ... . - Wolfdieter Lang, May 21 2015
The transpose gives the numerical coefficients of the Maurer-Cartan form matrix for the general linear group GL(n,1) (cf. Olver, but note that the formula at the bottom of p. 6 has an error--the 12 should be a 15). - Tom Copeland, Nov 05 2015
The left invariant Maurer-Cartan form polynomial on p. 7 of the Olver paper for the group GL^n(1) is essentially a binomial convolution of the row polynomials of this entry with those of A133314, or equivalently the row polynomials generated by the product of the e.g.f. of this entry with that of A133314, with some reindexing. - Tom Copeland, Jul 03 2018
From Tom Copeland, Jul 10 2018: (Start)
The first column of the inverse matrix is the sequence of Bernoulli numbers, which follows from the umbral definition of the Bernoulli polynomials (B.(0) + x)^n = B_n(x) evaluated at x = 1 and the relation B_n(0) = B_n(1) for n > 1 and -B_1(0) = 1/2 = B_1(1), so the Bernoulli numbers can be calculated using Cramer's rule acting on this entry's matrix and, therefore, from the ratios of volumes of parallelepipeds determined by the columns of this entry's square submatrices. - Tom Copeland, Jul 10 2018
Umbrally composing the row polynomials with B_n(x), the Bernoulli polynomials, gives (B.(x)+1)^(n+1) - (B.(x))^(n+1) = d[x^(n+1)]/dx = (n+1)*x^n, so multiplying this entry as a lower triangular matrix (LTM) by the LTM of the coefficients of the Bernoulli polynomials gives the diagonal matrix of the natural numbers. Then the inverse matrix of this entry has the elements B_(n,k)/(k+1), where B_(n,k) is the coefficient of x^k for B_n(x), and the e.g.f. (1/x) (e^(xt)-1)/(e^t-1). (End)

			T(4,2) = 0+0+1+3+6 = 10 = binomial(5, 2).
Triangle T(n,k) begins:
n\k 0  1  2   3   4   5   6   7   8   9 10 11
0:  1
1:  1  2
2:  1  3  3
3:  1  4  6   4
4:  1  5 10  10   5
5:  1  6 15  20  15   6
6:  1  7 21  35  35  21   7
7:  1  8 28  56  70  56  28   8
8:  1  9 36  84 126 126  84  36  9
9:  1 10 45 120 210 252 210 120 45   10
10: 1 11 55 165 330 462 462 330 165  55 11
11: 1 12 66 220 495 792 924 792 495 220 66 12
... Reformatted. - _Wolfdieter Lang_, Nov 04 2014
.
Can be seen as the square array A(n, k) = binomial(n + k + 1, n) read by descending antidiagonals. A(n, k) is the number of monotone nondecreasing functions f: {1,2,..,k} -> {1,2,..,n}. - _Peter Luschny_, Aug 25 2019
[0]  1,  1,   1,   1,    1,    1,     1,     1,     1, ... A000012
[1]  2,  3,   4,   5,    6,    7,     8,     9,    10, ... A000027
[2]  3,  6,  10,  15,   21,   28,    36,    45,    55, ... A000217
[3]  4, 10,  20,  35,   56,   84,   120,   165,   220, ... A000292
[4]  5, 15,  35,  70,  126,  210,   330,   495,   715, ... A000332
[5]  6, 21,  56, 126,  252,  462,   792,  1287,  2002, ... A000389
[6]  7, 28,  84, 210,  462,  924,  1716,  3003,  5005, ... A000579
[7]  8, 36, 120, 330,  792, 1716,  3432,  6435, 11440, ... A000580
[8]  9, 45, 165, 495, 1287, 3003,  6435, 12870, 24310, ... A000581
[9] 10, 55, 220, 715, 2002, 5005, 11440, 24310, 48620, ... A000582

Reinhard Zumkeller, Rows n = 0..150 of triangle, flattened
Feryal Alayont and Evan Henning, Edge Covers of Caterpillars, Cycles with Pendants, and Spider Graphs, J. Int. Seq. (2023) Vol. 26, Art. 23.9.4.
Paul Barry, On the f-Matrices of Pascal-like Triangles Defined by Riordan Arrays, arXiv:1805.02274 [math.CO], 2018.
Tom Copeland, Appell polynomials, cumulants, noncrossing partitions, Dyck lattice paths, and inversion, 2014.
Tom Copeland, Generators, Inversion, and Matrix, Binomial, and Integral Transforms, 2015.
Sela Fried, The expected degree of noninvertibility of compositions of functions and a related combinatorial identity, arXiv:2202.13061 [math.CO], 2022.
J. R. Griggs, The Cycling of Partitions and Compositions under Repeated Shifts, Advances in Applied Mathematics, Volume 21, Issue 2, August 1998, Pages 205-227.
P. Olver, The canonical contact form p. 7.
D. P. Walsh, A short note on increasing acyclic functions

Cf. A007318, A181971, A228196, A228576.
Row sums are A000225, diagonal sums are A052952.
The number of acyclic functions is A058127.
Cf. A008292, A090582, A019538, A049019, A133314, A135278, A133437, A111785, A126216, A134685, A133932, A248727, A033282, A134264.
Cf. A000027, A000217, A000292, A000332, A000389, A000579, A000580, A000581, A000582.
Cf. A325191.

Flat(List([0..10],n->List([0..n],k->Binomial(n+1,k)))); # Muniru A Asiru, Jul 10 2018

a074909 n k = a074909_tabl !! n !! k
a074909_row n = a074909_tabl !! n
a074909_tabl = iterate
   (\row -> zipWith (+) ([0] ++ row) (row ++ [1])) [1]
-- Reinhard Zumkeller, Feb 25 2012

/* As triangle */ [[Binomial(n+1,k): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Jul 22 2018

A074909 := proc(n,k)
    if k > n or k < 0 then
        0;
    else
        binomial(n+1,k) ;
    end if;
end proc: # Zerinvary Lajos, Nov 09 2006

Flatten[Join[{1}, Table[Sum[Binomial[k, m], {k, 0, n}], {n, 0, 12}, {m, 0, n}] ]] (* or *) Flatten[Join[{1}, Table[Binomial[n, m], {n, 12}, {m, n}]]]

print1(1);for(n=1,10,for(k=1,n,print1(", "binomial(n,k)))) \\ Charles R Greathouse IV, Mar 26 2013

from math import comb, isqrt
def A074909(n): return comb(r:=(m:=isqrt(k:=n+1<<1))+(k>m*(m+1)),n-comb(r,2)) # Chai Wah Wu, Nov 12 2024

T(n, k) = Sum_{i=0..n} C(i, n-k) = C(n+1, k).
Row n has g.f. (1+x)^(n+1)-x^(n+1).
E.g.f.: ((1+x)*e^t - x) e^(x*t). The row polynomials p_n(x) satisfy dp_n(x)/dx = (n+1)*p_(n-1)(x). - Tom Copeland, Jul 10 2018
T(n, k) = T(n-1, k-1) + T(n-1, k) for k: 0Reinhard Zumkeller, Apr 18 2005
T(n,k) = T(n-1,k) + 2*T(n-1,k-1) - T(n-2,k-1) - T(n-2,k-2), T(0,0)=1, T(1,0)=1, T(1,1)=2, T(n,k)=0 if k<0 or if k>n. - Philippe Deléham, Dec 27 2013
G.f. for column k (with leading zeros): x^(k-1)*(1/(1-x)^(k+1)-1), k >= 0. - Wolfdieter Lang, Nov 04 2014
Up(n, x+y) = (Up(.,x)+ y)^n = Sum_{k=0..n} binomial(n,k) Up(k,x)*y^(n-k), where Up(n,x) = ((x+1)^(n+1)-x^(n+1)) / (n+1) = P(n,x)/(n+1) with P(n,x) the n-th row polynomial of this entry. dUp(n,x)/dx = n * Up(n-1,x) and dP(n,x)/dx = (n+1)*P(n-1,x). - Tom Copeland, Nov 14 2014
The o.g.f. GF(x,t) = x / ((1-t*x)*(1-(1+t)x)) = x + (1+2t)*x^2 + (1+3t+3t^2)*x^3 + ... has the inverse GFinv(x,t) = (1+(1+2t)x-sqrt(1+(1+2t)*2x+x^2))/(2t(1+t)x) in x about 0, which generates the row polynomials (mod row signs) of A033282. The reciprocal of the o.g.f., i.e., x/GF(x,t), gives the free cumulants (1, -(1+2t) , t(1+t) , 0, 0, ...) associated with the moments defined by GFinv, and, in fact, these free cumulants generate these moments through the noncrossing partitions of A134264. The associated e.g.f. and relations to Grassmannians are described in A248727, whose polynomials are the basis for an Appell sequence of polynomials that are umbral compositional inverses of the Appell sequence formed from this entry's polynomials (distinct from the one described in the comments above, without the normalizing reciprocal). - Tom Copeland, Jan 07 2015
T(n, k) = (1/k!) * Sum_{i=0..k} Stirling1(k,i)*(n+1)^i, for 0<=k<=n. - Ridouane Oudra, Oct 23 2022

I added an initial 1 at the suggestion of Paul Barry, which makes the triangle a little nicer but may mean that some of the formulas will now need adjusting. - N. J. A. Sloane, Feb 11 2003

Formula section edited, checked and corrected by Wolfdieter Lang, Nov 04 2014

Previous Showing 71-80 of 644 results. Next

Number of palindromic structures using exactly three different symbols; Mobius transform: A056279. - Marks R. Nester
Number of ways of placing n labeled balls into k=3 indistinguishable boxes. - Thomas Wieder, Nov 30 2004
With two leading zeros, this is the second binomial transform of cosh(x)-1 and the binomial transform of A000225 (with extra leading zero). - Paul Barry, May 13 2003
Let [m] denote the first m positive integers. Then a(n) is the number of functions f from [n] to [n+1] that satisfy (i) f(x) > x for all x, (ii) f(x) = n+1 for exactly 3 elements and (iii) f(f(x)) = n+1 for the remaining n-3 elements of [n]. For example, a(4)=6 since there are exactly 6 functions from {1,2,3,4} to {1,2,3,4,5} such that f(x) > x, f(x) = 5 for 3 elements and f(f(x)) = 5 for the remaining element. The functions are f1 = {(1,5), (2,5), (3,4), (4,5)}, f2 = {(1,5), (2,3), (3,5), (4,5)}, f3 = {(1,5), (2,4), (3,5), (4,5)}, f4 = {(1,2), (2,5), (3,5), (4,5)}, f5 = {(1,3), (2,5), (3,5), (4,5)}, f6 = {(1,4), (2,5), (3,5), (4,5)}. - Dennis P. Walsh, Feb 20 2007
Conjecture. Let S(1)={1} and, for n > 1, let S(n) be the smallest set containing x, 2x and 3x for each element x in S(n-1). Then a(n+2) is the sum of the elements in S(n). (It is easy to prove that the number of elements in S(n) is the n-th triangular number given by A001952.) See A122554 for a sequence defined in this way. - John W. Layman, Nov 21 2007; corrected (a(n) to a(n+2) due to offset change) by Fred Daniel Kline, Oct 02 2014
Let P(A) be the power set of an n-element set A. Then a(n+1) = the number of pairs of elements {x,y} of P(A) for which x and y are disjoint and for which x is not a subset of y and y is not a subset of x. Wieder calls these "disjoint strict 2-combinations". - Ross La Haye, Jan 11 2008; corrected by Ross La Haye, Oct 29 2008
Also, let P(A) be the power set of an n-element set A. Then a(n+2) = the number of pairs of elements {x,y} of P(A) for which either 0) x and y are disjoint and for which either x is a subset of y or y is a subset of x, or 1) x and y are disjoint and for which x is not a subset of y and y is not a subset of x, or 2) x and y are intersecting and for which either x is a proper subset of y or y is a proper subset of x. - Ross La Haye, Jan 11 2008
3 * a(n+1) = p(n+1) where p(x) is the unique degree-n polynomial such that p(k) = a(k+1) for k = 0, 1, ..., n. - Michael Somos, Apr 29 2012
John W. Layman's conjecture that a(n+2) is the sum of elements in S(n) follows from the identification of S(n) with the first n rows of A036561, whose row sums are A001047. - Fred Daniel Kline, Oct 02 2014
From M. Sinan Kul, Sep 08 2016: (Start)
Let m be equal to the product of n-1 distinct primes. Then a(n) is equal to the number of distinct fractions >=1 that may be created by dividing a divisor of m by another divisor. For example for m = 2*3*5 = 30, we would have the following 6 fractions: 6/5, 3/2, 5/3, 5/2, 10/3, 15/2.
Here finding the number of fractions would be equivalent to distributing n-1 balls (distinct primes) to two bins (numerator and denominator) with no empty bins which can be found by Stirling numbers of the second kind. So another definition for a(n) is a(n) = Sum_{i=2..n-1} Stirling2(i,2)*binomial(n-1,i).
Also for n > 0, a(n) = (d(m^2)+1)/2 - d(m) where m is equal to the product of n-1 distinct primes. Example for a(5): m = 2*3*5*7 = 210 (product of four distinct primes) so a(5) = (d(210^2)+1)/2 - d(210) = 41 - 16 = 25. (End)
6*a(n) is the number of ternary strings of length n that contain at least one of each of the 3 symbols on which they are defined. For example, for n=4, the strings are the 12 permutations of 0012, the 12 permutations of 0112, and the 12 permutations of 0122. - Enrique Navarrete, Aug 23 2021
A simpler form of La Haye's first comment is: a(n+1) is the number of ways we can form disjoint unions of two nonempty subsets of [n] (see example below). Cf. A001047 for the requirement that the union contains n. - Enrique Navarrete, Aug 24 2021
As partial sums of the Nicomachus triangle's rows and the differences of the powers of 3 and 2 (A001047), each iteration corresponds to two figurate variations of the Sierpinski triangle (3^n) with cross-correlation to the Nicomachus triangle, see illustrations in links. The Sierpinski half-hexagons of (A001047) stack and conform to the footprint of 2^n - 1 triangular numbers. The 3^n Sierpinski triangle minus its 2^n bottom row, also correlates to the Nicomachus triangle according to each Sierpinski triangular sub-row. - John Elias, Oct 04 2021

Chan and Manna prove that a(n) is odd if and only if n is in A003714. - Jason Kimberley, Sep 14 2009
The number of ways to partition a set of 2*n elements into n disjoint subsets. - Vladimir Reshetnikov, Oct 10 2016
Conjecture: a(2*n+1) is divisible by (2*n + 1)^2. - Peter Bala, Mar 30 2025

Coefficients of (D^k f)(g(t))*(D g(t))^k1*(D^2 g(t))^k2*... in the Faa di Bruno formula for D^m(f(g(t))) where k = k1 + k2 + ..., m = 1*k1 + 2*k2 + ....

Number of set partitions whose block sizes are the prime indices of n (i.e., the integer partition with Heinz number n). - Gus Wiseman, Sep 12 2018

cp's OEIS Frontend

A000392 Stirling numbers of second kind S(n,3).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

Maple

Mathematica

PARI

Sage

Formula

Extensions

A007820 Stirling numbers of second kind S(2n,n).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

Maxima

PARI

PARI

PARI

Python

Sage

Formula

Extensions

A124794 Coefficients of incomplete Bell polynomials in the prime factorization order.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A035342 The convolution matrix of the double factorial of odd numbers (A001147).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Maxima

Sage

Formula

Extensions

A028246 Triangular array a(n,k) = (1/k)*Sum_{i=0..k} (-1)^(k-i)*binomial(k,i)*i^n; n >= 1, 1 <= k <= n, read by rows.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

Magma

Maple

A028246 Triangular array a(n,k) = (1/k)Sum_{i=0..k} (-1)^(k-i)binomial(k,i)*i^n; n >= 1, 1 <= k <= n, read by rows.

A217900 O.g.f.: Sum_{n>=0} n^n * (n+1)^(n-1) * exp(-n(n+1)x) * x^n / n!.

A001296 4-dimensional pyramidal numbers: a(n) = (3n+1)binomial(n+2, 3)/4. Also Stirling2(n+2, n).