A103194 -id:A103194 - cp's OEIS frontend

A105278 Triangle read by rows: T(n,k) = binomial(n,k)*(n-1)!/(k-1)!.

1, 2, 1, 6, 6, 1, 24, 36, 12, 1, 120, 240, 120, 20, 1, 720, 1800, 1200, 300, 30, 1, 5040, 15120, 12600, 4200, 630, 42, 1, 40320, 141120, 141120, 58800, 11760, 1176, 56, 1, 362880, 1451520, 1693440, 846720, 211680, 28224, 2016, 72, 1, 3628800, 16329600

Offset: 1

Miklos Kristof, Apr 25 2005

T(n,k) is the number of partially ordered sets (posets) on n elements that consist entirely of k chains. For example, T(4, 3)=12 since there are exactly 12 posets on {a,b,c,d} that consist entirely of 3 chains. Letting ab denote a<=b and using a slash "/" to separate chains, the 12 posets can be given by a/b/cd, a/b/dc, a/c/bd, a/c/db, a/d/bc, a/d/cb, b/c/ad, b/c/da, b/d/ac, b/d/ca, c/d/ab and c/d/ba, where the listing of the chains is arbitrary (e.g., a/b/cd = a/cd/b =...cd/b/a). - Dennis P. Walsh, Feb 22 2007
Also the matrix product |S1|.S2 of Stirling numbers of both kinds.
This Lah triangle is a lower triangular matrix of the Jabotinsky type. See the column e.g.f. and the D. E. Knuth reference given in A008297. - Wolfdieter Lang, Jun 29 2007
The infinitesimal matrix generator of this matrix is given in A132710. See A111596 for an interpretation in terms of circular binary words and generalized factorials. - Tom Copeland, Nov 22 2007
Three combinatorial interpretations: T(n,k) is (1) the number of ways to split [n] = {1,...,n} into a collection of k nonempty lists ("partitions into sets of lists"), (2) the number of ways to split [n] into an ordered collection of n+1-k nonempty sets that are noncrossing ("partitions into lists of noncrossing sets"), (3) the number of Dyck n-paths with n+1-k peaks labeled 1,2,...,n+1-k in some order. - David Callan, Jul 25 2008
Given matrices A and B with A(n,k) = T(n,k)*a(n-k) and B(n,k) = T(n,k)*b(n-k), then A*B = D where D(n,k) = T(n,k)*[a(.)+b(.)]^(n-k), umbrally. - Tom Copeland, Aug 21 2008
An e.g.f. for the row polynomials of A(n,k) = T(n,k)*a(n-k) is exp[a(.)* D_x * x^2] exp(x*t) = exp(x*t) exp[(.)!*Lag(.,-x*t,1)*a(.)*x], umbrally, where [(.)! Lag(.,x,1)]^n = n! Lag(n,x,1) is a normalized Laguerre polynomial of order 1. - Tom Copeland, Aug 29 2008
Triangle of coefficients from the Bell polynomial of the second kind for f = 1/(1-x). B(n,k){x1,x2,x3,...} = B(n,k){1/(1-x)^2,...,(j-1)!/(1-x)^j,...} = T(n,k)/(1-x)^(n+k). - Vladimir Kruchinin, Mar 04 2011
The triangle, with the row and column offset taken as 0, is the generalized Riordan array (exp(x), x) with respect to the sequence n!*(n+1)! as defined by Wang and Wang (the generalized Riordan array (exp(x), x) with respect to the sequence n! is Pascal's triangle A007318, and with respect to the sequence n!^2 is A021009 unsigned). - Peter Bala, Aug 15 2013
For a relation to loop integrals in QCD, see p. 33 of Gopakumar and Gross and Blaizot and Nowak. - Tom Copeland, Jan 18 2016
Also the Bell transform of (n+1)!. For the definition of the Bell transform see A264428. - Peter Luschny, Jan 27 2016
Also the number of k-dimensional flats of the n-dimensional Shi arrangement. - Shuhei Tsujie, Apr 26 2019
The numbers T(n,k) appear as coefficients when expanding the rising factorials (x)^k = x(x+1)...(x+k-1) in the basis of falling factorials (x)k = x(x-1)...(x-k+1). Specifically, (x)^n = Sum{k=1..n} T(n,k) (x)k. - _Jeremy L. Martin, Apr 21 2021

			T(1,1) = C(1,1)*0!/0! = 1,
T(2,1) = C(2,1)*1!/0! = 2,
T(2,2) = C(2,2)*1!/1! = 1,
T(3,1) = C(3,1)*2!/0! = 6,
T(3,2) = C(3,2)*2!/1! = 6,
T(3,3) = C(3,3)*2!/2! = 1,
Sheffer a-sequence recurrence: T(6,2)= 1800 = (6/3)*120 + 6*240.
B(n,k) =
   1/(1-x)^2;
   2/(1-x)^3,  1/(1-x)^4;
   6/(1-x)^4,  6/(1-x)^5,  1/(1-x)^6;
  24/(1-x)^5, 36/(1-x)^6, 12/(1-x)^7, 1/(1-x)^8;
The triangle T(n,k) begins:
  n\k      1       2       3      4      5     6    7  8  9 ...
  1:       1
  2:       2       1
  3:       6       6       1
  4:      24      36      12      1
  5:     120     240     120     20      1
  6:     720    1800    1200    300     30     1
  7:    5040   15120   12600   4200    630    42    1
  8:   40320  141120  141120  58800  11760  1176   56  1
  9:  362880 1451520 1693440 846720 211680 28224 2016 72  1
  ...
Row n=10: [3628800, 16329600, 21772800, 12700800, 3810240, 635040, 60480, 3240, 90, 1]. - _Wolfdieter Lang_, Feb 01 2013
From _Peter Bala_, Feb 24 2025: (Start)
The array factorizes as an infinite product (read from right to left):
  /  1                \        /1             \^m /1           \^m /1           \^m
  |  2    1            |      | 0   1          |  |0  1         |  |1  1         |
  |  6    6   1        | = ...| 0   0   1      |  |0  1  1      |  |0  2  1      |
  | 24   36  12   1    |      | 0   0   1  1   |  |0  0  2  1   |  |0  0  3  1   |
  |120  240 120  20   1|      | 0   0   0  2  1|  |0  0  0  3  1|  |0  0  0  4  1|
  |...                 |      |...             |  |...          |  |...          |
where m = 2. Cf. A008277 (m = 1), A035342 (m = 3), A035469 (m = 4), A049029 (m = 5) A049385 (m = 6), A092082 (m = 7), A132056 (m = 8), A223511 - A223522 (m = 9 through 20), A001497 (m = -1), A004747 (m = -2), A000369 (m = -3), A011801 (m = -4), A013988 (m = -5). (End)

Reinhard Zumkeller, Rows n = 1..100 of triangle, flattened
Peter Bala, Factorising (r,b)-Stirling arrays
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.
Paul Barry, Eulerian polynomials as moments, via exponential Riordan arrays, arXiv preprint arXiv:1105.3043 [math.CO], 2011, J. Int. Seq. 14 (2011) # 11.9.5
Jean-Paul Blaizot and Maciej A. Nowak, Large N_c confinement and turbulence, arXiv:0801.1859 [hep-th], 2008.
David Callan, Sets, Lists and Noncrossing Partitions, arXiv:0711.4841 [math.CO], 2007-2008.
Pietro Codara, Ottavio M. D'Antona, and Pavol Hell, A simple combinatorial interpretation of certain generalized Bell and Stirling numbers, arXiv preprint arXiv:1308.1700 [cs.DM], 2013.
Tom Copeland, Mathemagical Forests, Addendum to Mathemagical Forests, The Inverse Mellin Transform, Bell Polynomials, a Generalized Dobinski Relation, and the Confluent Hypergeometric Functions, A Class of Differential Operators and the Stirling Numbers
Siad Daboul, Jan Mangaldan, Michael Z. Spivey and Peter Taylor, The Lah Numbers and the n-th Derivative of exp(1/x), Math. Mag., 86 (2013), 39-47.
Bérénice Delcroix-Oger and Clément Dupont, Lie-operads and operadic modules from poset cohomology, arXiv:2505.06094 [math.CO], 2025. See p. 33.
G. H. E. Duchamp et al., Feynman graphs and related Hopf algebras, J. Phys. (Conf Ser) 30 (2006) 107-118.
Rajesh Gopakumar and David J. Gross, Mastering the master field, arXiv:hep-th/9411021, 1994.
Gábor Hetyei, Meixner polynomials of the second kind and quantum algebras representing su(1,1), arXiv preprint arXiv:0909.4352 [math.QA], 2009, p. 4. - From _Tom Copeland_, Oct 01 2015
Milan Janjic, Some classes of numbers and derivatives, JIS 12 (2009) 09.8.3.
Donald E. Knuth, Convolution polynomials, The Mathematica J., 2 (1992), 67-78.
Shi-Mei Ma, Some combinatorial sequences associated with context-free grammars, arXiv:1208.3104v2 [math.CO], 2012. - From _N. J. A. Sloane_, Aug 21 2012
MacTutor History of Mathematics archive: Biography of Ivo Lah.
Robert S. Maier, Boson Operator Ordering Identities from Generalized Stirling and Eulerian Numbers, arXiv:2308.10332 [math.CO], 2023. See p. 19.
Norihiro Nakashima and Shuhei Tsujie, Enumeration of Flats of the Extended Catalan and Shi Arrangements with Species, arXiv:1904.09748 [math.CO], 2019. See p. 18.
Michael Penn, Lah Numbers and an appearance of exponential generating functions, YouTube video, 2025.
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. See p. 4.
Tilman Piesk, Illustration of the first four rows
Kornelia Ufniarz and Grzegorz Siudem, Combinatorial origins of the canonical ensemble, arXiv:2008.00244 [math-ph], 2020. See p. 5.
Weiping Wang and Tianming Wang, Generalized Riordan arrays, Discrete Mathematics, Vol. 308, No. 24, 6466-6500.
Wikipedia, Lah number

Cf. A000262, A103194, A105220.
Triangle of Lah numbers (A008297) unsigned.
Cf. A111596 (signed triangle with extra n=0 row and m=0 column).
Cf. A130561 (for a natural refinement).
Cf. A094638 (for differential operator representation).
Cf. A248045 (central terms), A002868 (row maxima).
Cf, A059110.
Cf. A001263, A008277.
Cf. A089231 (triangle with mirrored rows).
Cf. A271703 (triangle with extra n=0 row and m=0 column).

Flat(List([1..10],n->List([1..n],k->Binomial(n,k)*Factorial(n-1)/Factorial(k-1)))); # Muniru A Asiru, Jul 25 2018

a105278 n k = a105278_tabl !! (n-1) !! (k-1)
a105278_row n = a105278_tabl !! (n-1)
a105278_tabl = [1] : f [1] 2 where
   f xs i = ys : f ys (i + 1) where
     ys = zipWith (+) ([0] ++ xs) (zipWith (*) [i, i + 1 ..] (xs ++ [0]))
-- Reinhard Zumkeller, Sep 30 2014, Mar 18 2013

/* As triangle */ [[Binomial(n,k)*Factorial(n-1)/Factorial(k-1): k in [1..n]]: n in [1.. 15]]; // Vincenzo Librandi, Oct 31 2014

The triangle: for n from 1 to 13 do seq(binomial(n,k)*(n-1)!/(k-1)!,k=1..n) od;
the sequence: seq(seq(binomial(n,k)*(n-1)!/(k-1)!,k=1..n),n=1..13);
# The function BellMatrix is defined in A264428.
# Adds (1, 0, 0, 0, ...) as column 0.
BellMatrix(n -> (n+1)!, 9); # Peter Luschny, Jan 27 2016

nn = 9; a = x/(1 - x); f[list_] := Select[list, # > 0 &]; Flatten[Map[f, Drop[Range[0, nn]! CoefficientList[Series[Exp[y a], {x, 0, nn}], {x, y}], 1]]] (* Geoffrey Critzer, Dec 11 2011 *)
nn = 9; Flatten[Table[(j - k)! Binomial[j, k] Binomial[j - 1, k - 1], {j, nn}, {k, j}]] (* Jan Mangaldan, Mar 15 2013 *)
rows = 10;
t = Range[rows]!;
T[n_, k_] := BellY[n, k, t];
Table[T[n, k], {n, 1, rows}, {k, 1, n}] // Flatten (* Jean-François Alcover, Jun 23 2018, after Peter Luschny *)
T[n_, n_] := 1; T[n_, k_] /;0Oliver Seipel, Dec 06 2024 *)

use ntheory ":all"; say join ", ", map { my $n=$; map { stirling($n,$,3) } 1..$n; } 1..9; # Dana Jacobsen, Mar 16 2017

T(n,k) = Sum_{m=n..k} |S1(n,m)|*S2(m,k), k>=n>=1, with Stirling triangles S2(n,m):=A048993 and S1(n,m):=A048994.
T(n,k) = C(n,k)*(n-1)!/(k-1)!.
Sum_{k=1..n} T(n,k) = A000262(n).
n*Sum_{k=1..n} T(n,k) = A103194(n) = Sum_{k=1..n} T(n,k)*k^2.
E.g.f. column k: (x^(k-1)/(1-x)^(k+1))/(k-1)!, k>=1.
Recurrence from Sheffer (here Jabotinsky) a-sequence [1,1,0,...] (see the W. Lang link under A006232): T(n,k)=(n/k)*T(n-1,m-1) + n*T(n-1,m). - Wolfdieter Lang, Jun 29 2007
The e.g.f. is, umbrally, exp[(.)!* L(.,-t,1)*x] = exp[t*x/(1-x)]/(1-x)^2 where L(n,t,1) = Sum_{k=0..n} T(n+1,k+1)*(-t)^k = Sum_{k=0..n} binomial(n+1,k+1)* (-t)^k / k! is the associated Laguerre polynomial of order 1. - Tom Copeland, Nov 17 2007
For this Lah triangle, the n-th row polynomial is given umbrally by
  n! C(B.(x)+1+n,n) = (-1)^n C(-B.(x)-2,n), where C(x,n)=x!/(n!(x-n)!),
  the binomial coefficient, and B_n(x)= exp(-x)(xd/dx)^n exp(x), the n-th Bell / Touchard / exponential polynomial (cf. A008277). E.g.,
  2! C(-B.(-x)-2,2) = (-B.(x)-2)(-B.(x)-3) = B_2(x) + 5*B_1(x) + 6 = 6 + 6x + x^2.
  n! C(B.(x)+1+n,n) = n! e^(-x) Sum_{j>=0} C(j+1+n,n)x^j/j! is a corresponding Dobinski relation. See the Copeland link for the relation to inverse Mellin transform. - Tom Copeland, Nov 21 2011
The row polynomials are given by D^n(exp(x*t)) evaluated at x = 0, where D is the operator (1+x)^2*d/dx. Cf. A008277 (D = (1+x)*d/dx), A035342 (D = (1+x)^3*d/dx), A035469 (D = (1+x)^4*d/dx) and A049029 (D = (1+x)^5*d/dx). - Peter Bala, Nov 25 2011
T(n,k) = Sum_{i=k..n} A130534(n-1,i-1)*A008277(i,k). - Reinhard Zumkeller, Mar 18 2013
Let E(x) = Sum_{n >= 0} x^n/(n!*(n+1)!). Then a generating function is exp(t)*E(x*t) = 1 + (2 + x)*t + (6 + 6*x + x^2)*t^2/(2!*3!) + (24 + 36*x + 12*x^2 + x^3)*t^3/(3!*4!) + ... . - Peter Bala, Aug 15 2013
P_n(x) = L_n(1+x) = n!*Lag_n(-(1+x);1), where P_n(x) are the row polynomials of A059110; L_n(x), the Lah polynomials of A105278; and Lag_n(x;1), the Laguerre polynomials of order 1. These relations follow from the relation between the iterated operator (x^2 D)^n and ((1+x)^2 D)^n with D = d/dx. - Tom Copeland, Jul 23 2018
Dividing each n-th diagonal by n!, where the main diagonal is n=1, generates the Narayana matrix A001263. - Tom Copeland, Sep 23 2020
T(n,k) = A089231(n,n-k). - Ron L.J. van den Burg, Dec 12 2021
T(n,k) = T(n-1,k-1) + (n+k-1)*T(n-1,k). - Bérénice Delcroix-Oger, Jun 25 2025

Stirling comments and e.g.f.s from Wolfdieter Lang, Apr 11 2007

A206703 Triangular array read by rows. T(n,k) is the number of partial permutations (injective partial functions) of {1,2,...,n} that have exactly k elements in a cycle. The k elements are not necessarily in the same cycle. A fixed point is considered to be in a cycle.

Original entry on oeis.org

1, 1, 1, 3, 2, 2, 13, 9, 6, 6, 73, 52, 36, 24, 24, 501, 365, 260, 180, 120, 120, 4051, 3006, 2190, 1560, 1080, 720, 720, 37633, 28357, 21042, 15330, 10920, 7560, 5040, 5040, 394353, 301064, 226856, 168336, 122640, 87360, 60480, 40320, 40320

Offset: 0

Geoffrey Critzer, Feb 11 2012

			     1;
     1,     1;
     3,     2,     2;
    13,     9,     6,     6;
    73,    52,    36,    24,    24;
   501,   365,   260,   180,   120,  120;
  4051,  3006,  2190,  1560,  1080,  720,   720;
  ...

Mohammad K. Azarian, On the Fixed Points of a Function and the Fixed Points of its Composite Functions, International Journal of Pure and Applied Mathematics, Vol. 46, No. 1, 2008, pp. 37-44. Mathematical Reviews, MR2433713 (2009c:65129), March 2009. Zentralblatt MATH, Zbl 1160.65015.
Mohammad K. Azarian, Fixed Points of a Quadratic Polynomial, Problem 841, College Mathematics Journal, Vol. 38, No. 1, January 2007, p. 60. Solution published in Vol. 39, No. 1, January 2008, pp. 66-67.

Alois P. Heinz, Rows n = 0..140, flattened
Philippe Flajolet and Robert Sedgewick, Analytic Combinatorics, Cambridge Univ. Press, 2009, page 132.

Columns k = 0..1 give: A000262, A006152.
Main diagonal gives A000142.
Row sums give A002720.
T(2n,n) gives A088026.
Cf. A002720, A052852, A103194, A105219, A331725.

b:= proc(n) option remember; `if`(n=0, 1, add((p-> p+x^j*
      coeff(p, x, 0))(b(n-j)*binomial(n-1, j-1)*j!), j=1..n))
    end:
T:= n-> (p-> seq(coeff(p, x, i), i=0..n))(b(n)):
seq(T(n), n=0..10);  # Alois P. Heinz, Feb 19 2022

nn = 7; a = 1/(1 - x); ay = 1/(1 - y x); f[list_] := Select[list, # > 0 &]; Map[f, Range[0, nn]! CoefficientList[Series[Exp[a x] ay, {x, 0, nn}], {x, y}]] // Flatten

E.g.f.: exp(x/(1-x))/(1-y*x).
From Alois P. Heinz, Feb 19 2022: (Start)
Sum_{k=1..n} T(n,k) = A052852.
Sum_{k=0..n} k * T(n,k) = A103194(n).
Sum_{k=0..n} (n-k) * T(n,k) =  A105219(n).
Sum_{k=0..n} (-1)^k * T(n,k) = A331725(n). (End)

A317277 a(n) = Sum_{k=0..n} binomial(n-1,k-1)k^nn!/k!; a(0) = 1.

Original entry on oeis.org

1, 1, 6, 81, 1828, 60565, 2734926, 160109005, 11724156648, 1045312448841, 111114793839610, 13845807451708441, 1994597720747571468, 328351264019737949341, 61162428777982281583302, 12782305566531823350524805, 2975150384583838798131401296, 766253903501365584725344992529

Offset: 0

Ilya Gutkovskiy, Jul 25 2018

nonn

a(n) is the n-th term of the Lah transform of the n-th powers.

G. C. Greubel, Table of n, a(n) for n = 0..250
N. J. A. Sloane, Transforms

Cf. A000262, A052852, A103194, A293145, A317278, A317279.

[1]cat[(&+[Binomial(n-1,j-1)*Binomial(n,j)*Factorial(n-j)*j^n: j in [0..n]]): n in [1..30]]; // G. C. Greubel, Mar 09 2021

A317277:= n-> `if`(n=0,1, add(binomial(n-1,j-1)*binomial(n,j)*(n-j)!*j^n, j=0..n)); seq(A317277(n), n=0..30); # G. C. Greubel, Mar 09 2021

Join[{1}, Table[Sum[Binomial[n - 1, k - 1] k^n n!/k!, {k, n}], {n, 17}]]
Join[{1}, Table[n! SeriesCoefficient[Sum[k^n (x/(1 - x))^k/k!, {k, n}], {x, 0, n}], {n, 17}]]

a(n) = if (n==0, 1, sum(k=0, n, binomial(n-1, k-1)*k^n*n!/k!)); \\ Michel Marcus, Mar 10 2021; corrected Jun 15 2022

[1]+[sum(binomial(n-1,j-1)*binomial(n,j)*factorial(n-j)*j^n for j in (0..n)) for n in (1..30)] # G. C. Greubel, Mar 09 2021

a(n) = n! * [x^n] Sum_{k>=0} k^n*(x/(1 - x))^k/k!.

Name edited by Michel Marcus, Jun 15 2022

A341200 Square array T(n,k), n >= 0, k >= 0, read by antidiagonals, where T(n,k) = Sum_{j=0..n} j^k * (n-j)! * binomial(n,j)^2.

Original entry on oeis.org

1, 0, 2, 0, 1, 7, 0, 1, 6, 34, 0, 1, 8, 39, 209, 0, 1, 12, 63, 292, 1546, 0, 1, 20, 117, 544, 2505, 13327, 0, 1, 36, 243, 1168, 5225, 24306, 130922, 0, 1, 68, 549, 2800, 12525, 55656, 263431, 1441729, 0, 1, 132, 1323, 7312, 33425, 145836, 653023, 3154824, 17572114

Offset: 0

Seiichi Manyama, Feb 06 2021

			Square array begins:
     1,    0,    0,     0,     0,     0, ...
     2,    1,    1,     1,     1,     1, ...
     7,    6,    8,    12,    20,    36, ...
    34,   39,   63,   117,   243,   549, ...
   209,  292,  544,  1168,  2800,  7312, ...
  1546, 2505, 5225, 12525, 33425, 97125, ...

Columns k=0..4 gives A002720, A103194, A105219, A105218, A341196.
Main diagonal gives A341197.
Cf. A289192.

T[n_, k_] := Sum[If[j == k == 0, 1, j^k] * (n - j)! * Binomial[n, j]^2, {j, 0, n}]; Table[T[k, n - k], {n, 0, 9}, {k, 0, n}] // Flatten (* Amiram Eldar, Feb 06 2021 *)

T(n, k) = sum(j=0, n, j^k*(n-j)!*binomial(n, j)^2);

About e.g.f. of column k, see A105218 or A105219 comment.

A256467 Inverse Lah transform of the squares.

Original entry on oeis.org

0, 1, 2, -9, 28, -55, -234, 5047, -59464, 620433, -6210710, 60312791, -552386988, 4291343641, -14786103682, -469083221865, 17904311480176, -458594711604703, 10473023418660306, -228670491372982217, 4899169866194557580, -104056906653521654679, 2196053393686810460902

Offset: 0

Peter Luschny, Mar 30 2015

sign

Cf. A103194.

a := n -> `if`(n=0,0, -(-1)^n*n!*hypergeom([2, 1-n], [1, 1], 1)):
seq(simplify(a(n)), n=0..22);

a(n) = Sum_{k=0..n}(-1)^(n-k)*(n-k)!*C(n,n-k)*C(n-1,n-k)*k^2.
a(n) = (-1)^(n+1)*n!*hypergeom([2, 1-n], [1, 1], 1) for n>=1.
D-finite with recurrence +(-n+1)*a(n) +(-2*n^2+3*n+4)*a(n-1) -(n-1)*(n-2)*(n+1)*a(n-2)=0. - R. J. Mathar, Jul 27 2022

A300159 Number of ways of converting one set of lists containing n elements to another set of lists containing n elements by removing the last element from one of the lists and either appending it to an existing list or treating it as a new list.

Original entry on oeis.org

0, 0, 4, 30, 240, 2140, 21300, 235074, 2853760, 37819800, 543445380, 8416452550, 139753069104, 2476581106740, 46648575724660, 930581784937770, 19597766647728000, 434455097953799344, 10112163333554834820, 246539064280189932270, 6282671083849941925360

Offset: 0

Mitchell Keith Bloch, Feb 26 2018

nonn

All terms are even.

			a(0) = 0 since for 0 lists, 0 conversions are possible.
a(1) = 0 since for the 1 set of 1 list of length 1, there exist no possible conversions.
a(2) = 4 since for the 2 sets of 1 list of length 2, there exists only 1 conversion, and for the 1 set of 2 lists of length 1, there exist 2 conversions.
a(3) = 30 since for the 6 sets of 1 list of length 3, there exists 1 conversion, for the 6 sets of 1 list of length 2 and 1 list of length 1, there exist 3 conversions, and for the 1 set of 3 lists of length 1, there exists 6 conversions.

Alois P. Heinz, Table of n, a(n) for n = 0..443 (first 71 terms from Mitchell Keith Bloch)
Mitchell Keith Bloch, Program C++
Mitchell Keith Bloch, Program C++ with Boost
Index entries for sequences related to Laguerre polynomials

Extends A000262 to count conversions in addition to sets of lists.

Cf. A000041, A006152, A052852, A103194.

l:= func< n,b | Evaluate(LaguerrePolynomial(n,1), b) >;
[0,0,4]cat[Factorial(n)*( 2*l(n-2,-1) - l(n-3,-1) ): n in [3..30]]; // G. C. Greubel, Mar 09 2021

b:= proc(n, t, c) option remember; `if`(n=0, t^2-c, add(j!*
      binomial(n-1, j-1)*b(n-j, t+1, c+`if`(j=1, 1, 0)), j=1..n))
    end:
a:= n-> b(n, 0$2):
seq(a(n), n=0..25);  # Alois P. Heinz, Mar 05 2018
# second Maple program:
a:= proc(n) option remember; `if`(n<6, [0$2, 4, 30, 240, 2140][n+1],
     (n*(2*n^2-13*n+16)*a(n-1)-n*(n-1)*(n-3)*(n-4)*a(n-2))/((n-2)*(n-5)))
    end:
seq(a(n), n=0..25);  # Alois P. Heinz, Mar 05 2018

(* First program *)
b[n_, t_, c_]:= b[n,t,c]= If[n==0, t^2 -c, Sum[j! Binomial[n-1, j-1]b[n-j,t+1,c + If[j==1, 1, 0]], {j,n}]];
a[n_]:= b[n, 0, 0];
a/@ Range[0, 25] (* Jean-François Alcover, Nov 24 2020, after Alois P. Heinz *)
(* Second program *)
Table[If[n<2, 0, n!*(2*LaguerreL[n-2,1,-1] -LaguerreL[n-3,1,-1])], {n,0,30}] (* G. C. Greubel, Mar 09 2021 *)

[0,0,4]+[factorial(n)*(2*gen_laguerre(n-2,1,-1) - gen_laguerre(n-3,0,-1)) for n in (3..30)] # G. C. Greubel, Mar 09 2021

a(n) = Sum_{i=1..p(n)} (n!/(Product_{j=1..n} k(i,j)!) * ((Sum_{j=1..n} k(i,j))^2 - k(i,1))) (where p(n) is the number of partitions A000041 and k(i,j) is the number of partitions of size j in partitioning i).
From Alois P. Heinz, Mar 05 2018: (Start)
E.g.f.: x^2*(2-x)*exp(x/(1-x))/(x-1)^2.
a(n) = (n*(2*n^2-13*n+16)*a(n-1) - n*(n-1)*(n-3)*(n-4)*a(n-2)) / ((n-2)*(n-5)) for n>5. (End)
a(n) ~ n^(n + 3/4) * exp(2*sqrt(n) - n - 1/2) / sqrt(2). - Vaclav Kotesovec, Jun 02 2018
a(n) = n!*( 2*LaguerreL(n-2,1,-1) - LaguerreL(n-3,1,-1) ) for n > 1, with a(0) = a(1) = 0. - G. C. Greubel, Mar 09 2021

More terms from Mitchell Keith Bloch, Mar 05 2018

cp's OEIS Frontend

A105278 Triangle read by rows: T(n,k) = binomial(n,k)*(n-1)!/(k-1)!.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

Perl

Formula

Extensions

A206703 Triangular array read by rows. T(n,k) is the number of partial permutations (injective partial functions) of {1,2,...,n} that have exactly k elements in a cycle. The k elements are not necessarily in the same cycle. A fixed point is considered to be in a cycle.

Views

Author

Keywords

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A317277 a(n) = Sum_{k=0..n} binomial(n-1,k-1)*k^n*n!/k!; a(0) = 1.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Sage

Formula

Extensions

A341200 Square array T(n,k), n >= 0, k >= 0, read by antidiagonals, where T(n,k) = Sum_{j=0..n} j^k * (n-j)! * binomial(n,j)^2.

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica

PARI

Formula

A256467 Inverse Lah transform of the squares.

Views

Author

Keywords

Crossrefs

Programs

Maple

Formula

A300159 Number of ways of converting one set of lists containing n elements to another set of lists containing n elements by removing the last element from one of the lists and either appending it to an existing list or treating it as a new list.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Sage

Formula

Extensions

A323770 Expansion of e.g.f. x*(2 - x)*exp(x/(1 - x))/(2*(1 - x)^2).

Views

A317277 a(n) = Sum_{k=0..n} binomial(n-1,k-1)k^nn!/k!; a(0) = 1.

A323770 Expansion of e.g.f. x(2 - x)exp(x/(1 - x))/(2*(1 - x)^2).