cp's OEIS Frontend

Number of ways n competitors can rank in a competition, allowing for the possibility of ties.

Also number of asymmetric generalized weak orders on n points.

Also called the ordered Bell numbers.

A weak order is a relation that is transitive and complete.

Called Fubini numbers by Comtet: counts formulas in Fubini theorem when switching the order of summation in multiple sums. - Olivier Gérard, Sep 30 2002 [Named after the Italian mathematician Guido Fubini (1879-1943). - Amiram Eldar, Jun 17 2021]

If the points are unlabeled then the answer is a(0) = 1, a(n) = 2^(n-1) (cf. A011782).

For n>0, a(n) is the number of elements in the Coxeter complex of type A_{n-1}. The corresponding sequence for type B is A080253 and there one can find a worked example as well as a geometric interpretation. - Tim Honeywill and Paul Boddington, Feb 10 2003

Also number of labeled (1+2)-free posets. - Detlef Pauly, May 25 2003

Also the number of chains of subsets starting with the empty set and ending with a set of n distinct objects. - Andrew Niedermaier, Feb 20 2004

From Michael Somos, Mar 04 2004: (Start)

Stirling transform of A007680(n) = [3,10,42,216,...] gives [3,13,75,541,...].

Stirling transform of a(n) = [1,3,13,75,...] is A083355(n) = [1,4,23,175,...].

Stirling transform of A000142(n) = [1,2,6,24,120,...] is a(n) = [1,3,13,75,...].

Stirling transform of A005359(n-1) = [1,0,2,0,24,0,...] is a(n-1) = [1,1,3,13,75,...].

Stirling transform of A005212(n-1) = [0,1,0,6,0,120,0,...] is a(n-1) = [0,1,3,13,75,...].

(End)

Unreduced denominators in convergent to log(2) = lim_{n->infinity} n*a(n-1)/a(n).

a(n) is congruent to a(n+(p-1)p^(h-1)) (mod p^h) for n >= h (see Barsky).

Stirling-Bernoulli transform of 1/(1-x^2). - Paul Barry, Apr 20 2005

This is the sequence of moments of the probability distribution of the number of tails before the first head in a sequence of fair coin tosses. The sequence of cumulants of the same probability distribution is A000629. That sequence is twice the result of deletion of the first term of this sequence. - Michael Hardy (hardy(AT)math.umn.edu), May 01 2005

With p(n) = the number of integer partitions of n, p(i) = the number of parts of the i-th partition of n, d(i) = the number of different parts of the i-th partition of n, p(j,i) = the j-th part of the i-th partition of n, m(i,j) = multiplicity of the j-th part of the i-th partition of n, one has: a(n) = Sum_{i=1..p(n)} (n!/(Product_{j=1..p(i)} p(i,j)!)) * (p(i)!/(Product_{j=1..d(i)} m(i,j)!)). - Thomas Wieder, May 18 2005

The number of chains among subsets of [n]. The summed term in the new formula is the number of such chains of length k. - Micha Hofri (hofri(AT)wpi.edu), Jul 01 2006

Occurs also as first column of a matrix-inversion occurring in a sum-of-like-powers problem. Consider the problem for any fixed natural number m>2 of finding solutions to the equation Sum_{k=1..n} k^m = (k+1)^m. Erdős conjectured that there are no solutions for n, m > 2. Let D be the matrix of differences of D[m,n] := Sum_{k=1..n} k^m - (k+1)^m. Then the generating functions for the rows of this matrix D constitute a set of polynomials in n (for varying n along columns) and the m-th polynomial defining the m-th row. Let GF_D be the matrix of the coefficients of this set of polynomials. Then the present sequence is the (unsigned) first column of GF_D^-1. - Gottfried Helms, Apr 01 2007

Assuming A = log(2), D is d/dx and f(x) = x/(exp(x)-1), we have a(n) = (n!/2*A^(n+1)) Sum_{k=0..n} (A^k/k!) D^n f(-A) which gives Wilf's asymptotic value when n tends to infinity. Equivalently, D^n f(-a) = 2*( A*a(n) - 2*a(n-1) ). - Martin Kochanski (mjk(AT)cardbox.com), May 10 2007

List partition transform (see A133314) of (1,-1,-1,-1,...). - Tom Copeland, Oct 24 2007

First column of A154921. - Mats Granvik, Jan 17 2009

A slightly more transparent interpretation of a(n) is as the number of 'factor sequences' of N for the case in which N is a product of n distinct primes. A factor sequence of N of length k is of the form 1 = x(1), x(2), ..., x(k) = N, where {x(i)} is an increasing sequence such that x(i) divides x(i+1), i=1,2,...,k-1. For example, N=70 has the 13 factor sequences {1,70}, {1,2,70}, {1,5,70}, {1,7,70}, {1,10,70}, {1,14,70}, {1,35,70}, {1,2,10,70}, {1,2,14,70}, {1,5,10,70}, {1,5,35,70}, {1,7,14,70}, {1,7,35,70}. - Martin Griffiths, Mar 25 2009

Starting (1, 3, 13, 75, ...) = row sums of triangle A163204. - Gary W. Adamson, Jul 23 2009

Equals double inverse binomial transform of A007047: (1, 3, 11, 51, ...). - Gary W. Adamson, Aug 04 2009

If f(x) = Sum_{n>=0} c(n)*x^n converges for every x, then Sum_{n>=0} f(n*x)/2^(n+1) = Sum_{n>=0} c(n)*a(n)*x^n. Example: Sum_{n>=0} exp(n*x)/2^(n+1) = Sum_{n>=0} a(n)*x^n/n! = 1/(2-exp(x)) = e.g.f. - Miklos Kristof, Nov 02 2009

Hankel transform is A091804. - Paul Barry, Mar 30 2010

It appears that the prime numbers greater than 3 in this sequence (13, 541, 47293, ...) are of the form 4n+1. - Paul Muljadi, Jan 28 2011

The Fi1 and Fi2 triangle sums of A028246 are given by the terms of this sequence. For the definitions of these triangle sums, see A180662. - Johannes W. Meijer, Apr 20 2011

The modified generating function A(x) = 1/(2-exp(x))-1 = x + 3*x^2/2! + 13*x^3/3! + ... satisfies the autonomous differential equation A' = 1 + 3*A + 2*A^2 with initial condition A(0) = 0. Applying [Bergeron et al., Theorem 1] leads to two combinatorial interpretations for this sequence: (A) a(n) gives the number of plane-increasing 0-1-2 trees on n vertices, where vertices of outdegree 1 come in 3 colors and vertices of outdegree 2 come in 2 colors. (B) a(n) gives the number of non-plane-increasing 0-1-2 trees on n vertices, where vertices of outdegree 1 come in 3 colors and vertices of outdegree 2 come in 4 colors. Examples are given below. - Peter Bala, Aug 31 2011

Starting with offset 1 = the eigensequence of A074909 (the beheaded Pascal's triangle), and row sums of triangle A208744. - Gary W. Adamson, Mar 05 2012

a(n) = number of words of length n on the alphabet of positive integers for which the letters appearing in the word form an initial segment of the positive integers. Example: a(2) = 3 counts 11, 12, 21. The map "record position of block containing i, 1<=i<=n" is a bijection from lists of sets on [n] to these words. (The lists of sets on [2] are 12, 1/2, 2/1.) - David Callan, Jun 24 2013

This sequence was the subject of one of the earliest uses of the database. Don Knuth, who had a computer printout of the database prior to the publication of the 1973 Handbook, wrote to N. J. A. Sloane on May 18, 1970, saying: "I have just had my first real 'success' using your index of sequences, finding a sequence treated by Cayley that turns out to be identical to another (a priori quite different) sequence that came up in connection with computer sorting." A000670 is discussed in Exercise 3 of Section 5.3.1 of The Art of Computer Programming, Vol. 3, 1973. - N. J. A. Sloane, Aug 21 2014

Ramanujan gives a method of finding a continued fraction of the solution x of an equation 1 = x + a2*x^2 + ... and uses log(2) as the solution of 1 = x + x^2/2 + x^3/6 + ... as an example giving the sequence of simplified convergents as 0/1, 1/1, 2/3, 9/13, 52/75, 375/541, ... of which the sequence of denominators is this sequence, while A052882 is the numerators. - Michael Somos, Jun 19 2015

For n>=1, a(n) is the number of Dyck paths (A000108) with (i) n+1 peaks (UD's), (ii) no UUDD's, and (iii) at least one valley vertex at every nonnegative height less than the height of the path. For example, a(2)=3 counts UDUDUD (of height 1 with 2 valley vertices at height 0), UDUUDUDD, UUDUDDUD. These paths correspond, under the "glove" or "accordion" bijection, to the ordered trees counted by Cayley in the 1859 reference, after a harmless pruning of the "long branches to a leaf" in Cayley's trees. (Cayley left the reader to infer the trees he was talking about from examples for small n and perhaps from his proof.) - David Callan, Jun 23 2015

From David L. Harden, Apr 09 2017: (Start)

Fix a set X and define two distance functions d,D on X to be metrically equivalent when d(x_1,y_1) <= d(x_2,y_2) iff D(x_1,y_1) <= D(x_2,y_2) for all x_1, y_1, x_2, y_2 in X.

Now suppose that we fix a function f from unordered pairs of distinct elements of X to {1,...,n}. Then choose positive real numbers d_1 <= ... <= d_n such that d(x,y) = d_{f(x,y)}; the set of all possible choices of the d_i's makes this an n-parameter family of distance functions on X. (The simplest example of such a family occurs when n is a triangular number: When that happens, write n = (k 2). Then the set of all distance functions on X, when |X| = k, is such a family.) The number of such distance functions, up to metric equivalence, is a(n).

It is easy to see that an equivalence class of distance functions gives rise to a well-defined weak order on {d_1, ..., d_n}. To see that any weak order is realizable, choose distances from the set of integers {n-1, ..., 2n-2} so that the triangle inequality is automatically satisfied. (End)

a(n) is the number of rooted labeled forests on n nodes that avoid the patterns 213, 312, and 321. - Kassie Archer, Aug 30 2018

From A.H.M. Smeets, Nov 17 2018: (Start)

Also the number of semantic different assignments to n variables (x_1, ..., x_n) including simultaneous assignments. From the example given by Joerg Arndt (Mar 18 2014), this is easily seen by replacing

"{i}" by "x_i := expression_i(x_1, ..., x_n)",

"{i, j}" by "x_i, x_j := expression_i(x_1, .., x_n), expression_j(x_1, ..., x_n)", i.e., simultaneous assignment to two different variables (i <> j),

similar for simultaneous assignments to more variables, and

"<" by ";", i.e., the sequential constructor. These examples are directly related to "Number of ways n competitors can rank in a competition, allowing for the possibility of ties." in the first comment.

From this also the number of different mean definitions as obtained by iteration of n different mean functions on n initial values. Examples:

the AGM(x1,x2) = AGM(x2,x1) is represented by {arithmetic mean, geometric mean}, i.e., simultaneous assignment in any iteration step;

Archimedes's scheme (for Pi) is represented by {geometric mean} < {harmonic mean}, i.e., sequential assignment in any iteration step;

the geometric mean of two values can also be observed by {arithmetic mean, harmonic mean};

the AGHM (as defined in A319215) is represented by {arithmetic mean, geometric mean, harmonic mean}, i.e., simultaneous assignment, but there are 12 other semantic different ways to assign the values in an AGHM scheme.

By applying power means (also called Holder means) this can be extended to any value of n. (End)

Total number of faces of all dimensions in the permutohedron of order n. For example, the permutohedron of order 3 (a hexagon) has 6 vertices + 6 edges + 1 2-face = 13 faces, and the permutohedron of order 4 (a truncated octahedron) has 24 vertices + 36 edges + 14 2-faces + 1 3-face = 75 faces. A001003 is the analogous sequence for the associahedron. - Noam Zeilberger, Dec 08 2019

Number of odd multinomial coefficients N!/(a_1!*a_2!*...*a_k!). Here each a_i is positive, and Sum_{i} a_i = N (so 2^{N-1} multinomial coefficients in all), where N is any positive integer whose binary expansion has n 1's. - Richard Stanley, Apr 05 2022 (edited Oct 19 2022)

From Peter Bala, Jul 08 2022: (Start)

Conjecture: Let k be a positive integer. The sequence obtained by reducing a(n) modulo k is eventually periodic with the period dividing phi(k) = A000010(k). For example, modulo 16 we obtain the sequence [1, 1, 3, 13, 11, 13, 11, 13, 11, 13, ...], with an apparent period of 2 beginning at a(4). Cf. A354242.

More generally, we conjecture that the same property holds for integer sequences having an e.g.f. of the form G(exp(x) - 1), where G(x) is an integral power series. (End)

a(n) is the number of ways to form a permutation of [n] and then choose a subset of its descent set. - Geoffrey Critzer, Apr 29 2023

This is the Akiyama-Tanigawa transform of A000079, the powers of two. - Shel Kaphan, May 02 2024

Counts binary rooted trees (with out-degree <= 2), embedded in plane, with n labeled end nodes of degree 1. Unlabeled version gives Catalan numbers A000108.

Define a "downgrade" to be the permutation which places the items of a permutation in descending order. We are concerned with permutations that are identical to their downgrades. Only permutations of order 4n and 4n+1 can have this property; the number of permutations of length 4n having this property are equinumerous with those of length 4n+1. If a permutation p has this property then the reversal of this permutation also has it. a(n) = number of permutations of length 4n and 4n+1 that are identical to their downgrades. - Eugene McDonnell (eemcd(AT)mac.com), Oct 26 2003

Number of broadcast schemes in the complete graph on n+1 vertices, K_{n+1}. - Calin D. Morosan (cd_moros(AT)alumni.concordia.ca), Nov 28 2008

Hankel transform is A137565. - Paul Barry, Nov 25 2009

The e.g.f. of 1/a(n) = n!/(2*n)! is (exp(sqrt(x)) + exp(-sqrt(x)) )/2. - Wolfdieter Lang, Jan 09 2012

From Tom Copeland, Nov 15 2014: (Start)

Aerated with intervening zeros (1,0,2,0,12,0,120,...) = a(n) (cf. A123023 and A001147), the e.g.f. is e^(t^2), so this is the base for the Appell sequence with e.g.f. e^(t^2) e^(x*t) = exp(P(.,x),t) (reverse A059344, cf. A099174, A066325 also). P(n,x) = (a. + x)^n with (a.)^n = a_n and comprise the umbral compositional inverses for e^(-t^2)e^(x*t) = exp(UP(.,x),t), i.e., UP(n,P(.,t)) = x^n = P(n,UP(.,t)), e.g., (P(.,t))^n = P(n,t).

Equals A000407*2 with leading 1 added. (End)

a(n) is also the number of square roots of any permutation in S_{4*n} whose disjoint cycle decomposition consists of 2*n transpositions. - Luis Manuel Rivera Martínez, Mar 04 2015

Self-convolution gives A076729. - Vladimir Reshetnikov, Oct 11 2016

For n > 1, it follows from the formula dated Aug 07 2013 that a(n) is a Zumkeller number (A083207). - Ivan N. Ianakiev, Feb 28 2017

For n divisible by 4, a(n/4) is the number of ways to place n points on an n X n grid with pairwise distinct abscissae, pairwise distinct ordinates, and 90-degree rotational symmetry. For n == 1 (mod 4), the number of ways is a((n-1)/4) because the center point can be considered "fixed". For 180-degree rotational symmetry see A006882, for mirror symmetry see A000085, A135401, and A297708. - Manfred Scheucher, Dec 29 2017

C(n,k) is the number of lattice paths from (0,0) to (n,k) using steps (1,0) and (0,1). - Joerg Arndt, Jul 01 2011

Row sums are A001700.

T(n, k) is also the number of order-preserving full transformations (of an n-chain) of waist k (waist(alpha) = max(Im(alpha))). - Abdullahi Umar, Oct 02 2008

If T(r,c), r=0,1,2,..., c=1,2,...,(r+1), are the triangle elements, then for r > 0, T(r,c) = binomial(r+c-1,c-1) = M(r,c) is the number of monotonic mappings from an ordered set of r elements into an ordered set of c elements. For example, there are 15 monotonic mappings from an ordered set of 4 elements into an ordered set of 3 elements. For c > r+1, use the identity M(r,c) = M(c-1,r+1) = T(c-1,r+1). For example, there are 210 monotonic mappings from an ordered set of 4 elements into an ordered set of 7 elements, because M(4,7) = T(6,5) = 210. Number of monotonic endomorphisms in a set of r elements, M(r,r), therefore appear on the second diagonal of the triangle which coincides with A001700. - Stanislav Sykora, May 26 2012

Start at the origin. Flip a fair coin to determine steps of (1,0) or (0,1). Stop when you are a (perpendicular) distance of n steps from the x axis or the y axis. For k = 0,1,...,n-1, C(n-1,k)/2^(n+k) is the probability that you will stop on the point (n,k). This is equal to the probability that you will stop on the point (k,n). Hence, Sum_{k=0..n} C(n,k)/2^(n+k) = 1. - Geoffrey Critzer, May 13 2017

The g.f. for the sequence {b(N,p)}, with b(N,p) the number of cyclically inequivalent two-color, N bead necklaces with p beads of one color and N-p beads of the other color is, for prime numbers p, G(p(n),x):=P(p(n)-1,x)/((1-x)^(p(n)-1)*(1-x^p(n))), with the numerator polynomial P(p(n)-1,x):= sum(r(n,k)*x^k,k=0..p(n)-1) and the row polynomials of this triangle r(n,k):=sum(a(k,m)*p(n)^m,m=0..k). p(n)=A000040(n) (prime numbers).

Row sums (signed) give A000142(k)=k!. Row sums (unsigned) coincide with A007680(k)=(2*k+1)*k!, k>=0.

The (unsigned) column sequences are, for m=0..10: A000774, A081052, A103719-A103727.

The coefficients in the upper right triangle of the ED1 array (m > n) were found with the a(n,m) formula while the coefficients in the lower left triangle of the ED1 array (m <= n) were found with the recurrence relation, see below. We use for the array rows the letter n (>= 1) and for the array columns the letter m (>= 1).

Our procedure for finding the coefficients in the lower left triangle can be compared with the procedure that De Smit and Lenstra used to fill in the hole in the middle of 'The Print Gallery' by M. C. Escher, see the links. In this lithograph Escher made use of the so-called Droste effect, hence we propose to call this square array of numbers the ED1 array.

For the ED2, ED3 and ED4 arrays see A167560, A167572 and A167584.

Row sums are the factorials. T(n,n) = (n-2)!*(2n-3) = A007680(n-2) (for n>=2). T(n,1) = Fibonacci(n+1)-1 = A000071(n+1). Sum_{k=0..n-1} k*T(n,k) = A130153(n). For the statistic max(p(i)-i) see A056151.

Denominators of coefficients of the Taylor series of sinh(sqrt(2*x))/(sqrt(2*x)). - J. Zurita (jrzurita(AT)inaoep.mx), Dec 01 2007

The row polynomials t(n,x):=Sum_{k=0..n} T(n,k)*x^k satisfy the recurrence relation t(n,x) = x*t(n-1,x) + ((n-1)^2)*t(n-2,x); t(-1,x)=0, t(0,x)=1. - Wolfdieter Lang, see above.

This is an example of a Sheffer triangle (coefficient triangle for Sheffer polynomials). In the umbral calculus (see the Roman reference given under A048854) s(n,x) := Sum_{k=0..n} T(n,k)*x^k would be called Sheffer polynomials for (1/cosh(t),tanh(t)), which translates to the e.g.f. for column number k>=0 given by (1/sqrt(1-x^2))*((arctanh(x))^k)/k!. The e.g.f. given below is rewritten in this Sheffer context as (1/sqrt(1-x^2))*exp(y*log(sqrt((1+x)/(1-x))))= (1/sqrt(1-x^2))*exp(y*arctanh(x)). The rows of the Jabotinsky type triangle |A049218| provide the coefficients of the associated polynomials. - Wolfdieter Lang, Feb 24 2005

The solution of the differential-difference relation f(n+1,x)= (d/dx)f(n,x) + (n^2)*f(n-1,x), n >= 1, with inputs f(0,x) and f(1,x) = (d/dx)f(0,x) is f(n,x) = t(n,d_x)*f(0,x), with the differential operator d_x:=d/dx and the row polynomials t(n,x) defined above. This problem appears in a computation of thermo field dynamics where f(0,x)=1/cosh(x). See the triangle A060081. - Wolfdieter Lang, Feb 24 2005

The inverse of the Sheffer matrix T with elements T(n,k) is the Sheffer matrix A060081. - Wolfdieter Lang, Jul 22 2005

T(n,k)=0 if n-k= 1(mod 2), else T(n,k) = sum of M2(n,p), p from {1,...,A000041(n)} restricted to partitions with exactly k odd parts and any nonnegative number of even parts. For the M2-multinomial numbers in A-St order see A036039(n,p). - Wolfdieter Lang, Aug 07 2007

Previous name was: A simple grammar.

Stirling transform of A005359(n-1)=[0,0,2,0,24,0,...] is a(n-1)=[0,0,2,12,74,...]. - Michael Somos, Mar 04 2004

Stirling transform of -(-1)^n*A052566(n-1)=[1,-1,4,-6,48,...] is a(n-1)=[1,0,2,12,74,...]. - Michael Somos, Mar 04 2004

Stirling transform of A000142(n)=[0,2,6,24,120,...] is a(n)=[0,2,12,74,...]. - Michael Somos, Mar 04 2004

Stirling transform of A007680(n)=[2,10,42,216,...] is a(n+1)=[2,12,74,...]. - Michael Somos, Mar 04 2004

a(n) is the number of chains in the power set of {1,2,...,n} that do not contain the empty set and do not contain {1,2,...,n}. Equivalently, a(n) is the number of ordered set partitions of {1,2,...,n} into at least 2 classes. - Geoffrey Critzer, Sep 01 2014