A045379 -id:A045379 - cp's OEIS frontend

A000178 Superfactorials: product of first n factorials.

1, 1, 2, 12, 288, 34560, 24883200, 125411328000, 5056584744960000, 1834933472251084800000, 6658606584104736522240000000, 265790267296391946810949632000000000, 127313963299399416749559771247411200000000000, 792786697595796795607377086400871488552960000000000000

Offset: 0

N. J. A. Sloane

a(n) is also the Vandermonde determinant of the numbers 1,2,...,(n+1), i.e., the determinant of the (n+1) X (n+1) matrix A with A[i,j] = i^j, 1 <= i <= n+1, 0 <= j <= n. - Ahmed Fares (ahmedfares(AT)my-deja.com), May 06 2001
a(n) = (1/n!) * D(n) where D(n) is the determinant of order n in which the (i,j)-th element is i^j. - Amarnath Murthy, Jan 02 2002
Determinant of S_n where S_n is the n X n matrix S_n(i,j) = Sum_{d|i} d^j. - Benoit Cloitre, May 19 2002
Appears to be det(M_n) where M_n is the n X n matrix with m(i,j) = J_j(i) where J_k(n) denote the Jordan function of row k, column n (cf. A059380(m)). - Benoit Cloitre, May 19 2002
a(2n+1) = 1, 12, 34560, 125411328000, ... is the Hankel transform of A000182 (tangent numbers) = 1, 2, 16, 272, 7936, ...; example: det([1, 2, 16, 272; 2, 16, 272, 7936; 16, 272, 7936, 353792; 272, 7936, 353792, 22368256]) = 125411328000. - Philippe Deléham, Mar 07 2004
Determinant of the (n+1) X (n+1) matrix whose i-th row consists of terms 1 to n+1 of the Lucas sequence U(i,Q), for any Q. When Q=0, the Vandermonde matrix is obtained. - T. D. Noe, Aug 21 2004
Determinant of the (n+1) X (n+1) matrix A whose elements are A(i,j) = B(i+j) for 0 <= i,j <= n, where B(k) is the k-th Bell number, A000110(k) [I. Mezo, JIS 14 (2011) # 11.1.1]. - T. D. Noe, Dec 04 2004
The Hankel transform of the sequence A090365 is A000178(n+1); example: det([1,1,3; 1,3,11; 3,11,47]) = 12. - Philippe Deléham, Mar 02 2005
Theorem 1.3, page 2, of Polynomial points, Journal of Integer Sequences, Vol. 10 (2007), Article 07.3.6, provides an example of an Abelian quotient group of order (n-1) superfactorial, for each positive integer n. The quotient is obtained from sequences of polynomial values. - E. F. Cornelius, Jr. (efcornelius(AT)comcast.net), Apr 09 2007
Starting with offset 1 this is a 'Matryoshka doll' sequence with alpha=1, the multiplicative counterpart to the additive A000292. seq(mul(mul(i,i=alpha..k), k=alpha..n),n=alpha..12). - Peter Luschny, Jul 14 2009
For n>0, a(n) is also the determinant of S_n where S_n is the n X n matrix, indexed from 1, S_n(i,j)=sigma_i(j), where sigma_k(n) is the generalized divisor sigma function: A000203 is sigma_1(n). - Enrique Pérez Herrero, Jun 21 2010
a(n) is the multiplicative Wiener index of the (n+1)-vertex path. Example: a(4)=288 because in the path on 5 vertices there are 3 distances equal to 2, 2 distances equal to 3, and 1 distance equal to 4 (2*2*2*3*3*4=288). See p. 115 of the Gutman et al. reference. - Emeric Deutsch, Sep 21 2011
a(n-1) = Product_{j=1..n-1} j! = V(n) = Product_{1 <= i < j <= n} (j - i) (a Vandermondian V(n), see the Ahmed Fares May 06 2001 comment above), n >= 1, is in fact the determinant of any n X n matrix M(n) with entries M(n;i,j) = p(j-1,x = i), 1 <= i, j <= n, where p(m,x), m >= 0, are monic polynomials of exact degree m with p(0,x) = 1. This is a special x[i] = i choice in a general theorem given in Vein-Dale, p. 59 (written for the transposed matrix M(n;j,x_i) = p(i-1,x_j) = P_i(x_j) in Vein-Dale, and there a_{k,k} = 1, for k=1..n). See the Aug 26 2013 comment under A049310, where p(n,x) = S(n,x) (Chebyshev S). - Wolfdieter Lang, Aug 27 2013
a(n) is the number of monotonic magmas on n elements labeled 1..n with a symmetric multiplication table. I.e., Product(i,j) >= max(i,j); Product(i,j) = Product(j,i). - Chad Brewbaker, Nov 03 2013
The product of the pairwise differences of n+1 integers is a multiple of a(n) [and this does not hold for any k > a(n)]. - Charles R Greathouse IV, Aug 15 2014
a(n) is the determinant of the (n+1) X (n+1) matrix M with M(i,j) = (n+j-1)!/(n+j-i)!, 1 <= i <= n+1, 1 <= j <= n+1. - Stoyan Apostolov, Aug 26 2014
All terms are in A064807 and all terms after a(2) are in A005101. - Ivan N. Ianakiev, Sep 02 2016
Empirical: a(n-1) is the determinant of order n in which the (i,j)-th entry is the (j-1)-th derivative of x^(x+i-1) evaluated at x=1. - John M. Campbell, Dec 13 2016
Empirical: If f(x) is a smooth, real-valued function on an open neighborhood of 0 such that f(0)=1, then a(n) is the determinant of order n+1 in which the (i,j)-th entry is the (j-1)-th derivative of f(x)/((1-x)^(i-1)) evaluated at x=0. - John M. Campbell, Dec 27 2016
Also the automorphism group order of the n-triangular honeycomb rook graph. - Eric W. Weisstein, Jul 14 2017
Is the zigzag Hankel transform of A000182. That is, a(2*n+1) is the Hankel transform of A000182 and a(2*n+2) is the Hankel transform of A000182(n+1). - Michael Somos, Mar 11 2020
Except for n = 0, 1, superfactorial a(n) is never a square (proof in link Mabry and Cormick, FFF 4 p. 349); however, when k belongs to A349079 (see for further information), there exists m, 1 <= m <= k such that a(k) / m! is a square. - Bernard Schott, Nov 29 2021

			a(3) = (1/6)* | 1 1 1 | 2 4 8 | 3 9 27 |
a(7) = n! * a(n-1) = 7! * 24883200 = 125411328000.
a(12) = 1! * 2! * 3! * 4! * 5! * 6! * 7! * 8! * 9! * 10! * 11! * 12!
= 1^12 * 2^11 * 3^10 * 4^9 * 5^8 * 6^7 * 7^6 * 8^5 * 9^4 * 10^3 * 11^2 * 12^1
= 2^56 * 3^26 * 5^11 * 7^6 * 11^2.
G.f. = 1 + x + 2*x^2 + 12*x^3 + 288*x^4 + 34560*x^5 + 24883200*x^6 + ...

Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 545.
Steven R. Finch, Mathematical Constants, Cambridge, 2003, pp. 135-145.
A. Fletcher, J. C. P. Miller, L. Rosenhead and L. J. Comrie, An Index of Mathematical Tables. Vols. 1 and 2, 2nd ed., Blackwell, Oxford and Addison-Wesley, Reading, MA, 1962, Vol. 1, p. 50.
R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, p. 231.
H. J. Ryser, Combinatorial Mathematics. Mathematical Association of America, Carus Mathematical Monograph 14, 1963, p. 53.
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).
R. Vein and P. Dale, Determinants and Their Applications in Mathematical Physics, Springer, 1999.

Boris Hostnik, Table of n, a(n) for n = 0..46
Christian Aebi and Grant Cairns, Generalizations of Wilson's Theorem for Double-, Hyper-, Sub-and Superfactorials, The American Mathematical Monthly 122.5 (2015): 433-443.
Andreas B. G. Blobel, On convolution powers of 1/x, arXiv:2203.09519 [math.CO], 2022.
E. F. Cornelius, Jr. and Phill Schultz, Polynomial points , Journal of Integer Sequences, Vol. 10 (2007), Article 07.3.6.
Selden Crary, Factorization of the Determinant of the Gaussian-Covariance Matrix of Evenly Spaced Points Using an Inter-dimensional Multiset Duality, arXiv preprint arXiv:1406.6326 [math.ST], 2014-2019.
N. Destainville, R. Mosseri and F. Bailly, Configurational Entropy of Codimension-One Tilings and Directed Membranes, J. Stat. Phys. 87, Nos 3/4, 697 (1997).
J. East and R. D. Gray, Idempotent generators in finite partition monoids and related semigroups, arXiv preprint arXiv:1404.2359 [math.GR], 2014.
Richard Ehrenborg, The Hankel determinant of exponential polynomials, Amer. Math. Monthly, 107 (2000), 557-560.
William Q. Erickson and Jan Kretschmann, The structure and normalized volume of Monge polytopes, arXiv:2311.07522 [math.CO], 2023. See p. 7.
Steven R. Finch, Glaisher-Kinkelin Constant (gives asymptotic expressions for A002109, A000178) [Broken link]
Steven R. Finch, Glaisher-Kinkelin Constant (gives asymptotic expressions for A002109, A000178) [From the Wayback machine]
Shyam Sunder Gupta, Fascinating Factorials, Exploring the Beauty of Fascinating Numbers, Springer (2025) Ch. 16, 411-442.
Ivan Gutman, Wolfgang Linert, István Lukovits and Željko Tomović, The multiplicative version of the Wiener index, J. Chem. Inf. Comput. Sci., Vol. 40, No. 1 (2000), pp. 113-116.
Brady Haran and Sophie Maclean, What's special about 288?, Numberphile video (2023).
Aoife Hennessy, A Study of Riordan Arrays with Applications to Continued Fractions, Orthogonal Polynomials and Lattice Paths, Ph. D. Thesis, Waterford Institute of Technology, Oct. 2011.
Nick Hobson, Python program for this sequence.
A. M. Ibrahim, Extension of factorial concept to negative numbers, Notes on Number Theory and Discrete Mathematics, Vol. 19, 2013, 2, 30-42.
Pavel L. Krapivsky, Jean-Marc Luck and Kirone Mallick, Quantum return probability of a system of N non-interacting lattice fermions, Journal of Statistical Mechanics: Theory and Experiment, Vol. 2018, No. 2 (2018), 023104; arXiv preprint, arXiv:1710.08178 [cond-mat.mes-hall], 2017-2018.
Jeffrey C. Lagarias and Harsh Mehta, Products of binomial coefficients and unreduced Farey fractions, International Journal of Number Theory, Vol. 12, No. 1 (2016), pp. 57-91; arXiv preprint, arXiv:1409.4145 [math.NT], 2014-2015.
Mogens Esrom Larsen, Wronskian Harmony, Mathematics Magazine, vol. 63, no. 1, 1990, pp. 33-37.
John W. Layman, The Hankel Transform and Some of its Properties, J. Integer Sequences, 4 (2001), #01.1.5.
Rick Mabry and Laura McCormick, Square products of punctured sequences of factorials, Gaz. Aust. Math. Soc., 2009.
Rémy Mosseri and Francis Bailly, Configurational Entropy in Octagonal Tiling Models, Int. J. Mod. Phys. B, Vol. 7, No. 6-7 (1993), pp. 1427-1436.
Rémy Mosseri, F. Bailly and C. Sire, Configurational Entropy in Random Tiling Models, J. Non-Cryst. Solids, Vol. 153-154 (1993), pp. 201-204.
Amarnath Murthy, Miscellaneous Results and Theorems on Smarandache terms and factor partitions, Smarandache Notions Journal, Vol. 11, No. 1-2-3, Spring 2000.
Amarnath Murthy and Charles Ashbacher, Generalized Partitions and Some New Ideas on Number Theory and Smarandache Sequences, Hexis, Phoenix; USA 2005. See Section 3.14.
Christian Radoux, Query 145, Notices Amer. Math. Soc., 25-3 (1978), p. 197.
Christian Radoux, Déterminants de Hankel et théorème de Sylvester, Séminaire Lotharingien de Combinatoire, B28b (1992), 9 pp.
Vignesh Raman, The Generalized Superfactorial, Hyperfactorial and Primorial functions, arXiv:2012.00882 [math.NT], 2020.
Michel Waldschmidt, Schanuel Property for Elliptic and Quasi-Elliptic Functions, arXiv:2504.14041 [math.NT], 2025. Mentions this sequence, see p. 10.
Eric Weisstein's World of Mathematics, Barnes G-Function.
Eric Weisstein's World of Mathematics, Bell Number.
Eric Weisstein's World of Mathematics, Factorial Products.
Eric Weisstein's World of Mathematics, Graph Automorphism.
Eric Weisstein's World of Mathematics, Lucas Sequence.
Eric Weisstein's World of Mathematics, Superfactorial.
Eric Weisstein's World of Mathematics, Vandermonde Determinant.
Index to divisibility sequences.
Index entries for sequences related to factorial numbers.
Index to sequences related to Olympiads and other Mathematical competitions.

Partial products of A000142.
Cf. A002109, A036561, A000292, A098694, A098695, A113271, A087316, A113208, A113231, A113257, A113258, A113320, A113336, A113498, A113173, A113170, A113475, A113492, A113497, A113533, A113534, A113535, A113153, A113154, A113122, A114045, A055462, A137986, A137987.
A002109(n)*A000178(n-1) = (n!)^n = A036740(n) for n >= 1.
A000178 is the Hankel transform (see A001906 for definition) of A000085, A000110, A000296, A005425, A005493, A005494 and A045379. - John W. Layman, Jul 28 2000
Cf. A255322, A255358, A255359, A255360.
Cf. A051675, A255321, A255323, A255344, A287013.
Cf. A348692, A349079.

[&*[Factorial(k): k in [0..n]]: n in [0..20]]; // Bruno Berselli, Mar 11 2015

A000178 := proc(n)
    mul(i!,i=1..n) ;
end proc:
seq(A000178(n),n=0..10) ; # R. J. Mathar, Oct 30 2015

a[0] := 1; a[1] := 1; a[n_] := n!*a[n - 1]; Table[a[n], {n, 1, 12}] (* Stefan Steinerberger, Mar 10 2006 *)
Table[BarnesG[n], {n, 2, 14}] (* Zerinvary Lajos, Jul 16 2009 *)
FoldList[Times,1,Range[20]!] (* Harvey P. Dale, Mar 25 2011 *)
RecurrenceTable[{a[n] == n! a[n - 1], a[0] == 1}, a, {n, 0, 12}] (* Ray Chandler, Jul 30 2015 *)
BarnesG[Range[2, 20]] (* Eric W. Weisstein, Jul 14 2017 *)

A000178(n):=prod(k!,k,0,n)$ makelist(A000178(n),n,0,30); /* Martin Ettl, Oct 23 2012 */

A000178(n)=prod(k=2,n,k!) \\ M. F. Hasler, Sep 02 2007

a(n)=polcoeff(1-sum(k=0, n-1, a(k)*x^k/prod(j=1, k+1, (1+j!*x+x*O(x^n)) )), n) \\ Paul D. Hanna, Oct 02 2013

for(j=1,13, print1(prod(k=1,j,k^(j-k)),", ")) \\ Hugo Pfoertner, Apr 09 2020

A000178_list, n, m = [1], 1,1
for i in range(1,100):
    m *= i
    n *= m
    A000178_list.append(n) # Chai Wah Wu, Aug 21 2015

from math import prod
def A000178(n): return prod(i**(n-i+1) for i in range(2,n+1)) # Chai Wah Wu, Nov 26 2023

def mono_choices(a,b,n)
    n - [a,b].max
end
def comm_mono_choices(n)
    accum =1
    0.upto(n-1) do |i|
        i.upto(n-1) do |j|
            accum = accum * mono_choices(i,j,n)
        end
    end
    accum
end
1.upto(12) do |k|
    puts comm_mono_choices(k)
end # Chad Brewbaker, Nov 03 2013

a(0) = 1, a(n) = n!*a(n-1). - Lee Hae-hwang, May 13 2003, corrected by Ilya Gutkovskiy, Jul 30 2016
a(0) = 1, a(n) = 1^n * 2^(n-1) * 3^(n-2) * ... * n = Product_{r=1..n} r^(n-r+1). - Amarnath Murthy, Dec 12 2003 [Formula corrected by Derek Orr, Jul 27 2014]
a(n) = sqrt(A055209(n)). - Philippe Deléham, Mar 07 2004
a(n) = Product_{i=1..n} Product_{j=0..i-1} (i-j). - Paul Barry, Aug 02 2008
log a(n) = 0.5*n^2*log n - 0.75*n^2 + O(n*log n). - Charles R Greathouse IV, Jan 13 2012
Asymptotic: a(n) ~ exp(zeta'(-1) - 3/4 - (3/4)*n^2 - (3/2)*n)*(2*Pi)^(1/2 + (1/2)*n)*(n+1)^((1/2)*n^2 + n + 5/12). For example, a(100) is approx. 0.270317...*10^6941. (See A213080.) - Peter Luschny, Jun 23 2012
G.f.: 1 + x/(U(0) - x) where U(k) = 1 + x*(k+1)! - x*(k+2)!/U(k+1); (continued fraction). - Sergei N. Gladkovskii, Oct 02 2012
G.f.: G(0)/2, where G(k) = 1 + 1/(1 - 1/(1 + 1/((k+1)!*x*G(k+1)))); (continued fraction). - Sergei N. Gladkovskii, Jun 14 2013
G.f.: 1 = Sum_{n>=0} a(n)*x^n / Product_{k=1..n+1} (1 + k!*x). - Paul D. Hanna, Oct 02 2013
A203227(n+1)/a(n) -> e, as n -> oo. - Daniel Suteu, Jul 30 2016
From Ilya Gutkovskiy, Jul 30 2016: (Start)
a(n) = G(n+2), where G(n) is the Barnes G-function.
a(n) ~ exp(1/12 - n*(3*n+4)/4)*n^(n*(n+2)/2 + 5/12)*(2*Pi)^((n+1)/2)/A, where A is the Glaisher-Kinkelin constant (A074962).
Sum_{n>=0} (-1)^n/a(n) = A137986. (End)
0 = a(n)*a(n+2)^3 + a(n+1)^2*a(n+2)^2 - a(n+1)^3*a(n+3) for all n in Z (if a(-1)=1). - Michael Somos, Mar 11 2020
Sum_{n>=0} 1/a(n) = A287013 = 1/A137987. - Amiram Eldar, Nov 19 2020
a(n) = Wronskian(1, x, x^2, ..., x^n). - Mohammed Yaseen, Aug 01 2023
From Andrea Pinos, Apr 04 2024: (Start)
a(n) = e^(Sum_{k=1..n} (Integral_{x=1..k+1} Psi(x) dx)).
a(n) = e^(Integral_{x=1..n+1} (log(sqrt(2*Pi)) - (x-1/2) + x*Psi(x)) dx).
a(n) = e^(Integral_{x=1..n+1} (log(sqrt(2*Pi)) - (x-1/2) + (n+1)*Psi(x) - log(Gamma(x))) dx).
Psi(x) is the digamma function. (End)

A056857 Triangle read by rows: T(n,c) = number of successive equalities in set partitions of n.

Original entry on oeis.org

1, 1, 1, 2, 2, 1, 5, 6, 3, 1, 15, 20, 12, 4, 1, 52, 75, 50, 20, 5, 1, 203, 312, 225, 100, 30, 6, 1, 877, 1421, 1092, 525, 175, 42, 7, 1, 4140, 7016, 5684, 2912, 1050, 280, 56, 8, 1, 21147, 37260, 31572, 17052, 6552, 1890, 420, 72, 9, 1, 115975, 211470, 186300, 105240, 42630, 13104, 3150, 600, 90, 10, 1

Offset: 1

Winston C. Yang (winston(AT)cs.wisc.edu), Aug 31 2000

Number of successive equalities s_i = s_{i+1} in a set partition {s_1, ..., s_n} of {1, ..., n}, where s_i is the subset containing i, s(1) = 1 and s(i) <= 1 + max of previous s(j)'s.
T(n,c) = number of set partitions of the set {1,2,...,n} in which the size of the block containing the element 1 is k+1. Example: T(4,2)=3 because we have 123|4, 124|3 and 134|2. - Emeric Deutsch, Nov 10 2006
Let P be the lower-triangular Pascal-matrix (A007318), Then this is exp(P) / exp(1). - Gottfried Helms, Mar 30 2007. [This comment was erroneously attached to A011971, but really belongs here. - N. J. A. Sloane, May 02 2015]
From David Pasino (davepasino(AT)yahoo.com), Apr 15 2009: (Start)
As an infinite lower-triangular matrix (with offset 0 rather than 1, so the entries would be B(n - c)*binomial(n, c), B() a Bell number, rather than B(n - 1 - c)*binomial(n - 1, c) as below), this array is S P S^-1 where P is the Pascal matrix A007318, S is the Stirling2 matrix A048993, and S^-1 is the Stirling1 matrix A048994.
Also, S P S^-1 = (1/e)*exp(P). (End)
Exponential Riordan array [exp(exp(x)-1), x]. Equal to A007318*A124323. - Paul Barry, Apr 23 2009
Equal to A049020*A048994 as infinite lower triangular matrices. - Philippe Deléham, Nov 19 2011
Build a superset Q[n] of set partitions of {1,2,...,n} by distinguishing "some" (possibly none or all) of the singletons. Indexed from n >= 0, 0 <= k <= n, T(n,k) is the number of elements in Q[n] that have exactly k distinguished singletons. A singleton is a subset containing one element. T(3,1) = 6 because we have {{1}'{2,3}}, {{1,2}{3}'}, {{1,3}{2}'}, {{1}'{2}{3}}, {{1}{2}'{3}}, {{1}{2}{3}'}. - Geoffrey Critzer, Nov 10 2012
Let Bell(n,x) denote the n-th Bell polynomial, the n-th row polynomial of A048993. Then this is the triangle of connection constants when expressing the basis polynomials Bell(n,x + 1) in terms of the basis polynomials Bell(n,x). For example, row 3 is (5, 6, 3, 1) and 5 + 6*Bell(1,x) + 3*Bell(2,x) + Bell(3,x) = 5 + 6*x + 3*(x + x^2) + (x + 3*x^2 + x^3) = 5 + 10*x + 6*x^2 + x^3 = (x + 1) + 3*(x + 1)^2 + (x + 1)^3 = Bell(3,x + 1). - Peter Bala, Sep 17 2013

			For example {1, 2, 1, 2, 2, 3} is a set partition of {1, 2, 3, 4, 5, 6} and has 1 successive equality, at i = 4.
Triangle begins:
    1;
    1,   1;
    2,   2,   1;
    5,   6,   3,   1;
   15,  20,  12,   4,   1;
   52,  75,  50,  20,   5,   1;
  203, 312, 225, 100,  30,   6,   1;
  ...
From _Paul Barry_, Apr 23 2009: (Start)
Production matrix is
  1,  1;
  1,  1,  1;
  1,  2,  1,  1;
  1,  3,  3,  1,  1;
  1,  4,  6,  4,  1,  1;
  1,  5, 10, 10,  5,  1,  1;
  1,  6, 15, 20, 15,  6,  1,  1;
  1,  7, 21, 35, 35, 21,  7,  1,  1;
  1,  8, 28, 56, 70, 56, 28,  8,  1,  1; ... (End)

W. C. Yang, Conjectures on some sequences involving set partitions and Bell numbers, preprint, 2000. [Apparently unpublished]

Alois P. Heinz, Rows n = 1..141, flattened
H. W. Becker, Rooks and rhymes, Math. Mag., 22 (1948/49), 23-26. See Table III.
H. W. Becker, Rooks and rhymes, Math. Mag., 22 (1948/49), 23-26. [Annotated scanned copy]
Fufa Beyene, Jörgen Backelin, Roberto Mantaci, and Samuel A. Fufa, Set Partitions and Other Bell Number Enumerated Objects, J. Int. Seq., Vol. 26 (2023), Article 23.1.8.
A. Hennessy and Paul Barry, Generalized Stirling Numbers, Exponential Riordan Arrays, and Orthogonal Polynomials, J. Int. Seq. 14 (2011) # 11.8.2, Corollary 17.
G. Hurst and A. Schultz, An elementary (number theory) proof of Touchard's congruence, arXiv:0906.0696v2 [math.CO], 2009.
A. O. Munagi, Set partitions with successions and separations, Intl. J. Math. Math. Sci. 2005 (2005) 451-463.
M. Spivey, A generalized recurrence for Bell numbers, J. Int. Seq., 11 (2008), no. 2, Article 08.2.5
W. Yang, Bell numbers and k-trees, Disc. Math. 156 (1996) 247-252.

Cf. Bell numbers A000110 (column c=0), A052889 (c=1), A105479 (c=2), A105480 (c=3).
Cf. A056858-A056863. Essentially same as A056860, where the rows are read from right to left.
Cf. also A007318, A005493, A270953.
See A259691 for another version.
T(2n+1,n+1) gives A124102.
T(2n,n) gives A297926.

with(combinat): A056857:=(n,c)->binomial(n-1,c)*bell(n-1-c): for n from 1 to 11 do seq(A056857(n,c),c=0..n-1) od; # yields sequence in triangular form; Emeric Deutsch, Nov 10 2006
with(linalg): # Yields sequence in matrix form:
A056857_matrix := n -> subs(exp(1)=1, exponential(exponential(
matrix(n,n,[seq(seq(`if`(j=k+1,j,0),k=0..n-1),j=0..n-1)])))):
A056857_matrix(8); # Peter Luschny, Apr 18 2011

t[n_, k_] := BellB[n-1-k]*Binomial[n-1, k]; Flatten[ Table[t[n, k], {n, 1, 11}, {k, 0, n-1}]](* Jean-François Alcover, Apr 25 2012, after Emeric Deutsch *)

B(n) = sum(k=0, n, stirling(n, k, 2));
tabl(nn)={for(n=1, nn, for(k=0, n - 1, print1(B(n - 1 - k) * binomial(n - 1, k),", ");); print(););};
tabl(12); \\ Indranil Ghosh, Mar 19 2017

from sympy import bell, binomial
for n in range(1,12):
    print([bell(n - 1 - k) * binomial(n - 1, k) for k in range(n)]) # Indranil Ghosh, Mar 19 2017

def a(n): return (-1)^n / factorial(n)
@cached_function
def p(n, m):
    R = PolynomialRing(QQ, "x")
    if n == 0: return R(a(m))
    return R((m + x)*p(n - 1, m) - (m + 1)*p(n - 1, m + 1))
for n in range(11): print(p(n, 0).list())  # Peter Luschny, Jun 18 2023

T(n,c) = B(n-1-c)*binomial(n-1, c), where T(n,c) is the number of set partitions of {1, ..., n} that have c successive equalities and B() is a Bell number.
E.g.f.: exp(exp(x)+x*y-1). - Vladeta Jovovic, Feb 13 2003
G.f.: 1/(1-xy-x-x^2/(1-xy-2x-2x^2/(1-xy-3x-3x^2/(1-xy-4x-4x^2/(1-... (continued fraction). - Paul Barry, Apr 23 2009
Considered as triangle T(n,k), 0 <= k <= n: T(n,k) = A007318(n,k)*A000110(n-k) and Sum_{k=0..n} T(n,k)*x^k = A000296(n), A000110(n), A000110(n+1), A005493(n), A005494(n), A045379(n) for x = -1, 0, 1, 2, 3, 4 respectively. - Philippe Deléham, Dec 13 2009
Let R(n,x) denote the n-th row polynomial of the triangle. Then A000110(n+j) = Bell(n+j,1) = Sum_{k = 1..n} R(j,k)*Stirling2(n,k) (Spivey). - Peter Bala, Sep 17 2013

More terms from David Wasserman, Apr 22 2002

A005494 3-Bell numbers: E.g.f.: exp(3*z + exp(z) - 1).

Original entry on oeis.org

1, 4, 17, 77, 372, 1915, 10481, 60814, 372939, 2409837, 16360786, 116393205, 865549453, 6713065156, 54190360453, 454442481041, 3952241526188, 35590085232519, 331362825860749, 3185554606447814, 31581598272055879, 322516283206446897, 3389017736055752178

Offset: 0

N. J. A. Sloane, Simon Plouffe

nonn

For further information, references, programs, etc. for r-Bell numbers see A005493. - N. J. A. Sloane, Nov 27 2013
From expansion of falling factorials (binomial transform of A005493).
Row sums of Sheffer triangle (exp(3*x), exp(x)-1). - Wolfdieter Lang, Sep 29 2011

			G.f. = 1 + 4*x + 17*x^2 + 77*x^3 + 372*x^4 + 1915*x^5 + 10481*x^6 + 60814*x^7 + ...

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Eldar Fischer, Johann A. Makowsky, and Vsevolod Rakita, MC-finiteness of restricted set partition functions, arXiv:2302.08265 [math.CO], 2023.
R. Jakimczuk, Successive Derivatives and Integer Sequences, J. Int. Seq. 14 (2011) # 11.7.3.
J. W. Layman, The Hankel Transform and Some of its Properties, J. Integer Sequences, 4 (2001), #01.1.5.
Toufik Mansour and Mark Shattuck, A recurrence related to the Bell numbers, INTEGERS 11 (2011), #A67.
I. Mezo, The r-Bell numbers, J. Int. Seq. 14 (2011) # 11.1.1.
J. Riordan, Letter, Oct 31 1977
N. J. A. Sloane, Transforms
Earl Glen Whitehead Jr., Stirling number identities from chromatic polynomials, J. Combin. Theory, A 24 (1978), 314-317.

Cf. A000110, A005493, A045379, A108087, A124323, A196834.

A row or column of the array A108087.

A005494:= func< n | (&+[Binomial(n,j)*3^(n-j)*Bell(j): j in [0..n]]) >;
[A005494(n): n in [0..30]]; // G. C. Greubel, Dec 01 2022

seq(add(3^(n-i)*combinat:-bell(i)*binomial(n,i),i=0..n), n=0..50); # Robert Israel, Dec 16 2014
# second Maple program:
b:= proc(n, m) option remember; `if`(n=0,
      m^2, m*b(n-1, m)+b(n-1, m+1))
    end:
a:= n-> b(n+1, 0)-b(n, 0):
seq(a(n), n=0..23);  # Alois P. Heinz, Aug 03 2025

Range[0, 40]! CoefficientList[Series[Exp[3 x + Exp[x] - 1], {x, 0, 40}], x] (* Vincenzo Librandi, Mar 04 2014 *)

def A005494(n): return sum( 3^(n-j)*bell_number(j)*binomial(n,j) for j in range(n+1))
[A005494(n) for n in range(31)] # G. C. Greubel, Dec 01 2022

a(n) = Sum_{i=0..n} 3^(n-i)*B(i)*binomial(n,i) where B(n) = Bell numbers A000110(n). - Fred Lunnon, Aug 04 2007
a(n) = exp(-1)*Sum_{k>=0} ((k+3)^n)/k!. - Gerald McGarvey, Jun 03 2004. May be rewritten as a(n) = Sum_{k>=3} (k^n*(k-1)*(k-2)/k!)/exp(1), which is a Dobinski-type relation for this sequence. - Karol A. Penson, Aug 18 2006
Define f_1(x), f_2(x), ... such that f_1(x) = x^2*e^x, f_{n+1}(x) = (d/dx)(x*f_n(x)), for n=2,3,.... Then a(n-1) = e^(-1)*f_n(1). - Milan Janjic, May 30 2008
Let A be the upper Hessenberg matrix of order n defined by: A[i,i-1]=-1, A[i,j]=binomial(j-1,i-1), (i <= j), and A[i,j]=0 otherwise. Then, for n >= 1, a(n) = (-1)^(n)charpoly(A,-3). - Milan Janjic, Jul 08 2010
a(n) = Sum_{k=3..n+3} A143495(n+3,k), n >= 0. - Wolfdieter Lang, Sep 29 2011
G.f.: 1/U(0) where U(k)= 1 - x*(k+4) - x^2*(k+1)/U(k+1); (continued fraction, 1-step). - Sergei N. Gladkovskii, Oct 11 2012
G.f.: Sum_{k>0} x^(k-1) / ((1 - 3*x) * (1 - 4*x) * ... * (1 - (k+2)*x)). - Michael Somos, Feb 26 2014
G.f.: Sum_{k>0} k * x^(k-1) / ((1 - 2*x) * (1 - 3*x) * ... * (1 - (k+1)*x)). - Michael Somos, Feb 26 2014
a(n) ~ exp(n/LambertW(n) - n - 1) * n^(n + 3) / LambertW(n)^(n + 7/2). - Vaclav Kotesovec, Jun 10 2020
a(0) = 1; a(n) = 3 * a(n-1) + Sum_{k=0..n-1} binomial(n-1,k) * a(k). - Ilya Gutkovskiy, Jul 02 2020
a(n) = Sum_{k=0..n} 4^k*A124323(n, k). - Mélika Tebni, Jun 10 2022

A143496 4-Stirling numbers of the second kind.

Original entry on oeis.org

1, 4, 1, 16, 9, 1, 64, 61, 15, 1, 256, 369, 151, 22, 1, 1024, 2101, 1275, 305, 30, 1, 4096, 11529, 9751, 3410, 545, 39, 1, 16384, 61741, 70035, 33621, 7770, 896, 49, 1, 65536, 325089, 481951, 305382, 95781, 15834, 1386, 60, 1, 262144, 1690981, 3216795

Offset: 4

Peter Bala, Aug 20 2008

This is the case r = 4 of the r-Stirling numbers of the second kind. The 4-Stirling numbers of the second kind count the ways of partitioning the set {1,2,...,n} into k nonempty disjoint subsets with the restriction that the elements 1, 2, 3 and 4 belong to distinct subsets. For remarks on the general case see A143494 (r = 2). The corresponding array of 4-Stirling numbers of the first kind is A143493. The theory of r-Stirling numbers of both kinds is developed in [Broder]. For 4-Lah numbers refer to A143499.
From Wolfdieter Lang, Sep 29 2011: (Start)
T(n,k) = S(n,k,4), n >= k >= 4, in Mikhailov's first paper, eq.(28) or (A3). E.g.f. column k from (A20) with k->4, r->k. Therefore, with offset [0,0], this triangle is the Sheffer triangle (exp(4*x),exp(x)-1) with e.g.f. of column no. m >= 0: exp(4*x)*((exp(x)-1)^m)/m!. See one of the formulas given below. For Sheffer matrices see the W. Lang link under A006232 with the S. Roman reference, also found in A132393.
(End)

			Triangle begins
n\k|.....4.....5.....6.....7.....8.....9
========================================
4..|.....1
5..|.....4.....1
6..|....16.....9.....1
7..|....64....61....15.....1
8..|...256...369...151....22.....1
9..|..1024..2101..1275...305....30.....1
...
T(6,5) = 9. The set {1,2,3,4,5,6} can be partitioned into five subsets such that 1, 2, 3 and 4 belong to different subsets in 9 ways: {{1,5}{2}{3}{4}{6}}, {{1,6}{2}{3}{4}{5}}, {{2,5}{1}{3}{4}{6}}, {{2,6}{1}{3}{4}{5}}, {{3,5}{1}{2}{4}{6}}, {{3,6}{1}{2}{4}{5}}, {{4,5}{1}{2}{3}{6}}, {{4,6}{1}{2}{3}{5}} and {{5,6}{1}{2}{3}{4}}.

Andrei Z. Broder, The r-Stirling numbers, Report Number: CS-TR-82-949, Stanford University, Department of Computer Science; see also, Discrete Math. 49, 241-259 (1984).
Askar Dzhumadil'daev 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.
Paweł Hitczenko, A class of polynomial recurrences resulting in (n/log n, n/log^2 n)-asymptotic normality, arXiv:2403.03422 [math.CO], 2024. See p. 8.
V. V. Mikhailov, Ordering of some boson operator functions, J. Phys A: Math. Gen. 16 (1983) 3817-3827.
V. V. Mikhailov, Normal ordering and generalised Stirling numbers, J. Phys A: Math. Gen. 18 (1985) 231-235.
Erich Neuwirth, Recursively defined combinatorial functions: Extending Galton's board, Discrete Math. 239 No. 1-3, 33-51 (2001)
Lily L. Liu and Yi Wang, A unified approach to polynomial sequences with only real zeros, arXiv:math/0509207 [math.CO], 2005-2006.
Michael J. Schlosser and Meesue Yoo, Elliptic Rook and File Numbers, Electronic Journal of Combinatorics, 24(1) (2017), #P1.31.

Cf. A003468 (column 7), A005060 (column 5), A008277, A016103 (column 6), A045379 (row sums), A049459 (matrix inverse), A143493, A143494, A143495, A143499.

with combinat: T := (n, k) -> 1/(k-4)!*add ((-1)^(k-i)*binomial(k-4,i)*(i+4)^(n-4),i = 0..k-4): for n from 4 to 13 do seq(T(n, k), k = 4..n) end do;

t[n_, k_] := StirlingS2[n, k] - 6*StirlingS2[n-1, k] + 11*StirlingS2[n-2, k] - 6*StirlingS2[n-3, k]; Flatten[ Table[ t[n, k], {n, 4, 13}, {k, 4, n}]] (* Jean-François Alcover, Dec 02 2011 *)

T(n+4,k+4) = (1/k!)*Sum_{i = 0..k} (-1)^(k-i)*C(k,i)*(i+4)^n, n,k >= 0.
T(n,k) = Stirling2(n,k) - 6*Stirling2(n-1,k) + 11*Stirling2(n-2,k) - 6*Stirling2(n-3,k) for n,k >= 4.
Recurrence relation: T(n,k) = T(n-1,k-1) + k*T(n-1,k) for n > 4 with boundary conditions: T(n,3) = T(3,n) = 0 for all n; T(4,4) = 1; T(4,k) = 0 for k > 4. Special cases: T(n,4) = 4^(n-4); T(n,5) = 5^(n-4) - 4^(n-4).
E.g.f. (k+4)-th column (with offset 4): (1/k!)*exp(4*x)*(exp(x)-1)^k.
O.g.f. k-th column: Sum_{n>=k} T(n,k)*x^n = x^k/((1-4*x)*(1-5*x)*...*(1-k*x)).
E.g.f.: exp(4*t + x*(exp(t)-1)) = Sum_{n = 0..infinity} Sum_(k = 0..n) T(n+4,k+4)*x^k*t^n/n! = Sum_{n = 0..infinity} B_n(4;x)*t^n/n! = 1 + (4+x)*t/1! + (16+9*x+x^2)*t^2/2! + ..., where the row polynomials, B_n(4;x) := Sum_{k = 0..n} T(n+4,k+4)*x^k, may be called the 4-Bell polynomials.
Dobinski-type identities: Row polynomial B_n(4;x) = exp(-x)*Sum_{i = 0..infinity} (i+4)^n*x^i/i!; Sum_{k = 0..n} k!*T(n+4,k+4)*x^k = Sum_{i = 0..infinity} (i+4)^n*x^i/(1+x)^(i+1).
The T(n,k) are the connection coefficients between the falling factorials and the shifted monomials (x+4)^(n-4). For example, 16 + 9*x + x*(x-1) = (x+4)^2; 64 + 61*x + 15*x*(x-1) + x*(x-1)*(x-2) = (x+4)^3.
This array is the matrix product P^3 * S, where P denotes Pascal's triangle, A007318 and S denotes the lower triangular array of Stirling numbers of the second kind, A008277 (apply Theorem 10 of [Neuwirth]).
The inverse array is A049459, the signed 4-Stirling numbers of the first kind.
From Peter Bala, Sep 19 2008: (Start)
Let D be the derivative operator d/dx and E the Euler operator x*d/dx. Then x^(-4)*E^n*x^4 = Sum_{k = 0..n} T(n+4,k+4)*x^k*D^k.
The row generating polynomials R_n(x) := Sum_{k=4..n} T(n,k)*x^k satisfy the recurrence R_(n+1)(x) = x*R_n(x) + x*d/dx(R_n(x)) with R_4(x) = x^4. It follows that the polynomials R_n(x) have only real zeros (apply Corollary 1.2. of [Liu and Wang]).
Relation with the 4-Eulerian numbers E_4(n,j) := A144698(n,j): T(n,k) = 4!/k!*Sum_{j = n-k..n-4} E_4(n,j)*binomial(j,n-k) for n >= k >= 4.
(End)

A108087 Array, read by antidiagonals, where A(n,k) = exp(-1)*Sum_{i>=0} (i+k)^n/i!.

Original entry on oeis.org

1, 1, 1, 2, 2, 1, 5, 5, 3, 1, 15, 15, 10, 4, 1, 52, 52, 37, 17, 5, 1, 203, 203, 151, 77, 26, 6, 1, 877, 877, 674, 372, 141, 37, 7, 1, 4140, 4140, 3263, 1915, 799, 235, 50, 8, 1, 21147, 21147, 17007, 10481, 4736, 1540, 365, 65, 9, 1, 115975, 115975, 94828, 60814, 29371, 10427, 2727, 537, 82, 10, 1

Offset: 0

Gerald McGarvey, Jun 05 2005

The column for k=0 is A000110 (Bell or exponential numbers). The column for k=1 is A000110 starting at offset 1. The column for k=2 is A005493 (Sum_{k=0..n} k*Stirling2(n,k).). The column for k=3 is A005494 (E.g.f.: exp(3*z+exp(z)-1).). The column for k=4 is A045379 (E.g.f.: exp(4*z+exp(z)-1).). The row for n=0 is 1's sequence, the row for n=1 is the natural numbers. The row for n=2 is A002522 (n^2 + 1.). The row for n=3 is A005491 (n^3 + 3n + 1.). The row for n=4 is A005492.
Number of ways of placing n labeled balls into n+k boxes, where k of the boxes are labeled and the rest are indistinguishable. - Bradley Austin (artax(AT)cruzio.com), Apr 24 2006
The column for k = -1 (not shown) is A000296 (Number of partitions of an n-set into blocks of size >1. Also number of cyclically spaced (or feasible) partitions.). - Gerald McGarvey, Oct 08 2006
Equals antidiagonals of an array in which (n+1)-th column is the binomial transform of n-th column, with leftmost column = the Bell sequence, A000110. - Gary W. Adamson, Apr 16 2009
Number of partitions of [n+k] where at least k blocks contain their own index element.  A(2,2) = 10: 134|2, 13|24, 13|2|4, 14|23, 1|234, 1|23|4, 14|2|3, 1|24|3, 1|2|34, 1|2|3|4. - Alois P. Heinz, Jan 07 2022

			Array A(n,k) begins:
   1,   1,   1,    1,    1,     1,     1,     1,     1,      1, ... A000012;
   1,   2,   3,    4,    5,     6,     7,     8,     9,     10, ... A000027;
   2,   5,  10,   17,   26,    37,    50,    65,    82,    101, ... A002522;
   5,  15,  37,   77,  141,   235,   365,   537,   757,   1031, ... A005491;
  15,  52, 151,  372,  799,  1540,  2727,  4516,  7087,  10644, ... A005492;
  52, 203, 674, 1915, 4736, 10427, 20878, 38699, 67340, 111211, ... ;
Antidiagonal triangle, T(n, k), begins as:
     1;
     1,    1;
     2,    2,    1;
     5,    5,    3,    1;
    15,   15,   10,    4,   1;
    52,   52,   37,   17,   5,   1;
   203,  203,  151,   77,  26,   6,  1;
   877,  877,  674,  372, 141,  37,  7,  1;
  4140, 4140, 3263, 1915, 799, 235, 50,  8,  1;

F. Ruskey, Combinatorial Generation, preprint, 2001.

Alois P. Heinz, Antidiagonals n = 0..140, flattened
I. Mezo, The r-Bell numbers, J. Int. Seq. 14 (2011) # 11.1.1, Figure 1.
J. Riordan, Letter, Oct 31 1977, The array is on the second page.
F. Ruskey, Combinatorial Generation, 2003.
F. Ruskey, Lexicographic Algorithms [Broken link]

Cf. A000012, A000027, A000110, A002522, A005490, A005491.
Cf. A005492, A005493, A005494, A045379, A086659, A124824.
Main diagonal gives A134980.
Antidiagonal sums give A347420.

A108087:= func< n,k | (&+[Binomial(n-k,j)*k^j*Bell(n-k-j): j in [0..n-k]]) >;
[A108087(n,k): k in [0..n], n in [0..13]]; // G. C. Greubel, Dec 02 2022

with(combinat):
A:= (n, k)-> add(binomial(n, i) * k^i * bell(n-i), i=0..n):
seq(seq(A(d-k, k), k=0..d), d=0..12);  # Alois P. Heinz, Jul 18 2012

Unprotect[Power]; 0^0 = 1; A[n_, k_] := Sum[Binomial[n, i] * k^i * BellB[n - i], {i, 0, n}]; Table[Table[A[d - k, k], {k, 0, d}], {d, 0, 12}] // Flatten (* Jean-François Alcover, Nov 05 2015, after Alois P. Heinz *)

f(n,k)=round (suminf(i=0,(i+k)^n/i!)/exp(1));
g(n,k)=for(k=0,k,print1(f(n,k),",")) \\ prints k+1 terms of n-th row

def A108087(n,k): return sum( k^j*bell_number(n-k-j)*binomial(n-k,j) for j in range(n-k+1))
flatten([[A108087(n,k) for k in range(n+1)] for n in range(14)]) # G. C. Greubel, Dec 02 2022

For n> 1, A(n, k) = k^n + sum_{i=0..n-2} A086659(n, i)*k^i. (A086659 is set partitions of n containing k-1 blocks of length 1, with e.g.f: exp(x*y)*(exp(exp(x)-1-x)-1).)
A(n, k) = k * A(n-1, k) + A(n-1, k+1), A(0, k) = 1. - Bradley Austin (artax(AT)cruzio.com), Apr 24 2006
A(n,k) = Sum_{i=0..n} C(n,i) * k^i * Bell(n-i). - Alois P. Heinz, Jul 18 2012
Sum_{k=0..n-1} A(n-k,k) = A005490(n). - Alois P. Heinz, Jan 05 2022
From G. C. Greubel, Dec 02 2022: (Start)
T(n, n) = A000012(n).
T(n, n-1) = A000027(n).
T(n, n-2) = A002522(n-1).
T(n, n-3) = A005491(n-2).
T(n, n-4) = A005492(n+1).
T(2*n, n) = A134980(n).
T(2*n, n+1) = A124824(n), n >= 1.
Sum_{k=0..n} T(n, k) = A347420(n). (End)

A346739 Expansion of e.g.f.: exp(exp(x) - 4*x - 1).

Original entry on oeis.org

1, -3, 10, -35, 127, -472, 1787, -6855, 26572, -103765, 407695, -1608378, 6369117, -25271183, 100542930, -400114103, 1597052419, -6359524256, 25481982047, -101103395443, 409291679676, -1592903606657, 6729506287091, -23748796926026, 123501587468073, -227183793907851

Offset: 0

Ilya Gutkovskiy, Jul 31 2021

sign

Seiichi Manyama, Table of n, a(n) for n = 0..580

Cf. A000110, A000296, A045379, A126617, A193684, A290219, A346738, A346740.

R:=PowerSeriesRing(Rationals(), 30);
Coefficients(R!(Laplace( Exp(Exp(x) -4*x -1) ))) // G. C. Greubel, Jun 12 2024

nmax = 25; CoefficientList[Series[Exp[Exp[x] - 4 x - 1], {x, 0, nmax}], x] Range[0, nmax]!
Table[Sum[Binomial[n, k] (-4)^(n - k) BellB[k], {k, 0, n}], {n, 0, 25}]
a[0] = 1; a[n_] := a[n] = -4 a[n - 1] + Sum[Binomial[n - 1, k] a[k], {k, 0, n - 1}]; Table[a[n], {n, 0, 25}]

[factorial(n)*( exp(exp(x) -4*x -1) ).series(x, n+1).list()[n] for n in (0..30)] # G. C. Greubel, Jun 12 2024

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

A196834 Row sums of Sheffer triangle A193685 (5-restricted Stirling2 numbers).

Original entry on oeis.org

1, 6, 37, 235, 1540, 10427, 73013, 529032, 3967195, 30785747, 247126450, 2050937445, 17585497797, 155666739742, 1421428484337, 13377704321695, 129659127547372, 1293095848212799, 13259069937250169, 139671750579429512, 1510382932875294447, 16754464511605466311

Offset: 0

Wolfdieter Lang, Oct 07 2011

			a(2) = 25 + 11 + 1 = 37.

Seiichi Manyama, Table of n, a(n) for n = 0..500

Cf. A000110, A005493, A005494, A045379, A196835 (alternating row sums).

b:= proc(n, m) option remember;
     `if`(n=0, 1, m*b(n-1, m)+b(n-1, m+1))
    end:
a:= n-> b(n, 5):
seq(a(n), n=0..23);  # Alois P. Heinz, Aug 22 2021

nmax = 20; CoefficientList[Series[E^(E^x + 5*x - 1), {x, 0, nmax}], x] * Range[0, nmax]! (* Vaclav Kotesovec, Jun 10 2020 *)

a(n) = Sum_{m=0..n} A193685(n,m).
E.g.f.: exp(exp(x)+5*x-1).
a(n) ~ exp(n/LambertW(n) - n - 1) * n^(n + 5) / LambertW(n)^(n + 11/2). - Vaclav Kotesovec, Jun 10 2020
a(0) = 1; a(n) = 5 * a(n-1) + Sum_{k=0..n-1} binomial(n-1,k) * a(k). - Ilya Gutkovskiy, Jul 03 2020

A343674 a(0) = 1; a(n) = 4 * n * a(n-1) + Sum_{k=0..n-1} binomial(n,k) * a(k).

Original entry on oeis.org

1, 5, 51, 781, 15947, 407021, 12466251, 445452813, 18191122219, 835737327661, 42661645147403, 2395510523568845, 146739531459316587, 9737742346694258157, 695911661109898805323, 53286006304099668950413, 4352120920347139791200171, 377674509364714706139413933, 34702277449656625185428239755

Offset: 0

Ilya Gutkovskiy, Apr 25 2021

nonn

Cf. A000670, A006155, A032032, A045379, A343672, A343673.

a[0] = 1; a[n_] := a[n] = 4 n a[n - 1] + Sum[Binomial[n, k] a[k], {k, 0, n - 1}]; Table[a[n], {n, 0, 18}]
nmax = 18; CoefficientList[Series[1/(2 (1 - 2 x) - Exp[x]), {x, 0, nmax}], x] Range[0, nmax]!

E.g.f.: 1 / (2 * (1 - 2*x) - exp(x)).

a(n) ~ n! * 2^(n-1) / ((1 + LambertW(exp(1/2)/4)) * (1 - 2*LambertW(exp(1/2)/4))^(n+1)). - Vaclav Kotesovec, Jun 20 2022

A347204 a(n) = a(f(n)/2) + a(floor((n+f(n))/2)) for n > 0 with a(0) = 1 where f(n) = A129760(n).

Original entry on oeis.org

1, 2, 3, 5, 4, 7, 10, 15, 5, 9, 13, 20, 17, 27, 37, 52, 6, 11, 16, 25, 21, 34, 47, 67, 26, 43, 60, 87, 77, 114, 151, 203, 7, 13, 19, 30, 25, 41, 57, 82, 31, 52, 73, 107, 94, 141, 188, 255, 37, 63, 89, 132, 115, 175, 235, 322, 141, 218, 295, 409, 372, 523, 674

Offset: 0

Mikhail Kurkov, Aug 23 2021 [verification needed]

Modulo 2 binomial transform of A243499(n).

Antti Karttunen, Table of n, a(n) for n = 0..8191
Kevin Ryde, PARI/GP Code

Cf. A000110, A000120, A002720, A005493, A005494, A007814, A008277, A035009, A045379, A053645, A129760, A143494, A143495, A196834, A265649, A295989.

Similar recurrences: A124758, A243499, A284005, A329369, A341392.

function a = A347204(max_n)
    a(1) = 1;
    a(2) = 2;
    for nloop = 3:max_n
        n = nloop-1;
        s = 0;
        for k = 0:floor(log2(n))-1
            s = s + a(1+A053645(n)-2^k*(mod(floor(n/(2^k)),2)));
        end
        a(nloop) = 2*a(A053645(n)+1) + s;
    end
end
function a_n = A053645(n)
    a_n = n - 2^floor(log2(n));
end % Thomas Scheuerle, Oct 25 2021

f[n_] := BitAnd[n, n - 1]; a[0] = 1; a[n_] := a[n] = a[f[n]/2] + a[Floor[(n + f[n])/2]]; Array[a, 100, 0] (* Amiram Eldar, Nov 19 2021 *)

f(n) = bitand(n, n-1); \\ A129760
a(n) = if (n<=1, n+1, if (n%2, a(n\2)+a(n-1), a(f(n/2)) + a(n/2+f(n/2)))); \\ Michel Marcus, Oct 25 2021

PARI
```
\\ Also see links.
```

A129760(n) = bitand(n, n-1);
memoA347204 = Map();
A347204(n) = if (n<=1, n+1, my(v); if(mapisdefined(memoA347204,n,&v), v, v = if(n%2, A347204(n\2)+A347204(n-1), A347204(A129760(n/2)) + A347204(n/2+A129760(n/2))); mapput(memoA347204,n,v); (v))); \\ (Memoized version of Michel Marcus's program given above) - Antti Karttunen, Nov 20 2021

a(n) = a(n - 2^f(n)) + (1 + f(n))*a((n - 2^f(n))/2) for n > 0 with a(0) = 1 where f(n) = A007814(n).
a(2n+1) = a(n) + a(2n) for n >= 0.
a(2n) = a(n - 2^f(n)) + a(2n - 2^f(n)) for n > 0 with a(0) = 1 where f(n) = A007814(n).
a(n) = 2*a(f(n)) + Sum_{k=0..floor(log_2(n))-1} a(f(n) - 2^k*T(n,k)) for n > 1 with a(0) = 1, a(1) = 2, and where f(n) = A053645(n), T(n,k) = floor(n/2^k) mod 2.
Sum_{k=0..2^n - 1} a(k) = A035009(n+1) for n >= 0.
a((4^n - 1)/3) = A002720(n) for n >= 0.
a(2^n - 1) = A000110(n+1),
a(2*(2^n - 1)) = A005493(n),
a(2^2*(2^n - 1)) = A005494(n),
a(2^3*(2^n - 1)) = A045379(n),
a(2^4*(2^n - 1)) = A196834(n),
a(2^m*(2^n-1)) = T(n,m+1) is the n-th (m+1)-Bell number for n >= 0, m >= 0 where T(n,m) = m*T(n-1,m) + Sum_{k=0..n-1} binomial(n-1,k)*T(k,m) with T(0,m) = 1.
a(n) = Sum_{j=0..2^A000120(n)-1} A243499(A295989(n,j)) for n >= 0. Also A243499(n) = Sum_{j=0..2^f(n)-1} (-1)^(f(n)-f(j)) a(A295989(n,j)) for n >= 0 where f(n) = A000120(n). In other words, a(n) = Sum_{j=0..n} (binomial(n,j) mod 2)*A243499(j) and A243499(n) = Sum_{j=0..n} (-1)^(f(n)-f(j))*(binomial(n,j) mod 2)*a(j) for n >= 0 where f(n) = A000120(n).
Generalization:
b(n, x) = (1/x)*b((n - 2^f(n))/2, x) + (-1)^n*b(floor((2n - 2^f(n))/2), x) for n > 0 with b(0, x) = 1 where f(n) = A007814(n).
Sum_{k=0..2^n - 1} b(k, x) = (1/x)^n for n >= 0.
b((4^n - 1)/3, x) = (1/x)^n*n!*L_{n}(x) for n >= 0 where L_{n}(x) is the n-th Laguerre polynomial.
b((8^n - 1)/7, x) = (1/x)^n*Sum_{k=0..n} (-x)^k*A265649(n, k) for n >= 0.
b(2^n - 1, x) = (1/x)^n*Sum_{k=0..n} (-x)^k*A008277(n+1, k+1),
b(2*(2^n - 1), x) = (1/x)^n*Sum_{k=0..n} (-x)^k*A143494(n+2, k+2),
b(2^2*(2^n - 1), x) = (1/x)^n*Sum_{k=0..n} (-x)^k*A143495(n+3, k+3),
b(2^m*(2^n - 1), x) = (1/x)^n*Sum_{k=0..n} (-x)^k*T(n+m+1, k+m+1, m+1) for n >= 0, m >= 0 where T(n,k,m) is m-Stirling numbers of the second kind.

cp's OEIS Frontend

A000178 Superfactorials: product of first n factorials.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Maxima

PARI

PARI

PARI

Python

Python

Ruby

Formula

A056857 Triangle read by rows: T(n,c) = number of successive equalities in set partitions of n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

SageMath

Formula

Extensions

A005494 3-Bell numbers: E.g.f.: exp(3*z + exp(z) - 1).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

SageMath

Formula

A143496 4-Stirling numbers of the second kind.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A108087 Array, read by antidiagonals, where A(n,k) = exp(-1)*Sum_{i>=0} (i+k)^n/i!.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI