A000166 -id:A000166 - cp's OEIS frontend

A001910 a(n) = na(n-1) + (n-5)a(n-2).

0, 1, 5, 31, 227, 1909, 18089, 190435, 2203319, 27772873, 378673901, 5551390471, 87057596075, 1453986832381, 25762467303377, 482626240281739, 9530573107600319, 197850855756232465, 4307357140602486869, 98125321641110663023, 2334414826276390013171

Offset: 3

N. J. A. Sloane

nonn

With offset 1, permanent of (0,1)-matrix of size n X (n+d) with d=5 and n zeros not on a line. This is a special case of Theorem 2.3 of Seok-Zun Song et al. Extremes of permanents of (0,1)-matrices, pp. 201-202. - Jaap Spies, Dec 12 2003
a(n+4)=:b(n), n>=1, enumerates the ways to distribute n beads labeled differently from 1 to n, over a set of (unordered) necklaces, excluding necklaces with exactly one bead, and k=5 indistinguishable, ordered, fixed cords, each allowed to have any number of beads. Beadless necklaces as well as a beadless cords contribute each a factor 1 in the counting, e.g., b(0):= 1*1 =1. See A000255 for the description of a fixed cord with beads.
This produces for b(n) the exponential (aka binomial) convolution of the subfactorial sequence {A000166(n)} and the sequence {A001720(n+4) = (n+4)!/4!}. See the necklaces and cords problem comment in A000153. Therefore also the recurrence b(n) = (n+4)*b(n-1) + (n-1)*b(n-2) with b(-1)=0 and b(0)=1 holds. This comment derives from a family of recurrences found by Malin Sjodahl for a combinatorial problem for certain quark and gluon diagrams (Feb 27 2010). - Wolfdieter Lang, Jun 02 2010

			Necklaces and 5 cords problem. For n=4 one considers the following weak 2 part compositions of 4: (4,0), (3,1), (2,2), and (0,4), where (1,3) does not appear because there are no necklaces with 1 bead. These compositions contribute respectively sf(4)*1, binomial(4,3)*sf(3)*c5(1), (binomial(4,2)*sf(2))*c5(2), and 1*c5(4) with the subfactorials sf(n):=A000166(n) (see the necklace comment there) and the c5(n):=A001720(n+4) numbers for the pure 5 cord problem (see the remark on the e.g.f. for the k cords problem in A000153; here for k=5: 1/(1-x)^5). This adds up as 9 + 4*2*5 + (6*1)*30 + 1680 = 1909 = b(4) = A001910(8). - _Wolfdieter Lang_, Jun 02 2010

Brualdi, Richard A. and Ryser, Herbert J., Combinatorial Matrix Theory, Cambridge NY (1991), Chapter 7.
J. Riordan, An Introduction to Combinatorial Analysis, Wiley, 1958, p. 188.
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 = 3..100
Seok-Zun Song et al., Extremes of permanents of (0,1)-matrices, Special issue on the Combinatorial Matrix Theory Conference (Pohang, 2002). Linear Algebra Appl. 373 (2003), pp. 197-210.

Cf. A000255, A000153, A000261, A001909, A001910, A055790, A090012-A090016, A086764.

A001909 (necklaces and four cords).

a := n -> `if`(n=3,0, hypergeom([6,-n+4],[],1))*(-1)^n;
seq(round(evalf(a(n),100)), n=3..20); # Peter Luschny, Sep 20 2014

t = {0, 1}; Do[AppendTo[t, n*t[[-1]] + (n - 5) t[[-2]]], {n, 5, 20}]; t (* T. D. Noe, Aug 17 2012 *)

a(n) = A086764(n+1,5), n>=3.
E.g.f. with offset -1: (exp(-x)/(1-x))*(1-x)^5 = exp(-x)/(1-x)^6. - Wolfdieter Lang, Jun 02 2010
G.f.: x*hypergeom([1,6],[],x/(x+1))/(x+1). - Mark van Hoeij, Nov 07 2011
a(n) = hypergeometric([6,-n+4],[],1)*(-1)^n for n >=4. - Peter Luschny, Sep 20 2014

A008278 Reflected triangle of Stirling numbers of 2nd kind, S(n,n-k+1), n >= 1, 1 <= k <= n.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 6, 7, 1, 1, 10, 25, 15, 1, 1, 15, 65, 90, 31, 1, 1, 21, 140, 350, 301, 63, 1, 1, 28, 266, 1050, 1701, 966, 127, 1, 1, 36, 462, 2646, 6951, 7770, 3025, 255, 1, 1, 45, 750, 5880, 22827, 42525, 34105, 9330, 511, 1

Offset: 1

N. J. A. Sloane

The n-th row also gives the coefficients of the sigma polynomial of the empty graph \bar K_n. - Eric W. Weisstein, Apr 07 2017
The n-th row also gives the coefficients of the independence polynomial of the (n-1)-triangular honeycomb bishop graph. - Eric W. Weisstein, Apr 03 2018
From Gus Wiseman, Aug 11 2020: (Start)
Conjecture: also the number of divisors of the superprimorial A006939(n - 1) that have 0 <= k <= n distinct prime factors, all appearing with distinct multiplicities. For example, row n = 4 counts the following divisors of 360:
  1  2  12  360
     3  18
     4  20
     5  24
     8  40
     9  45
        72
Equivalently, T(n,k) is the number of length-n vectors 0 <= v_i <= i with k nonzero values, all of which are distinct.
Crossrefs:
A006939 lists superprimorials or Chernoff numbers.
A022915 counts permutations of prime indices of superprimorials.
A076954 can be used instead of A006939.
A130091 lists numbers with distinct prime multiplicities.
A181796 counts divisors with distinct prime multiplicities.
A336420 is the version counting all prime factors, not just distinct ones.
Cf. A000005, A027423, A095149, A124010, A317829, A327498, A336419, A336421.
(End)
From Leonidas Liponis, Aug 26 2024: (Start)
It appears that this sequence is related to the combinatorial form of Faà di Bruno's formula. Specifically, the number of terms for the n-th derivative of a composite function y = f(g(x)) matches the number of partitions of n.
For example, consider the case where g(x) = e^x, in which all derivatives of g(x) are equal. The first 5 rows of A008278 appear as the factors of derivatives of f(x), highlighted here in brackets:
          dy/dx = [ 1 ] * f'(e^x) * e^x
          d^2y/dx^2 = [ 1 ] * f''(e^x) * e^{2x} + [ 1 ] * f'(e^x) * e^x
          d^3y/dx^3 = [ 1 ] * f'''(e^x) * e^{3x} + [ 3 ] * f''(e^x) * e^{2x} + [ 1 ] * f'(e^x) * e^x
          d^4y/dx^4 = [ 1 ] * f''''(e^x) * e^{4x} + [ 6 ] * f'''(e^x) * e^{3x} + [ 7 ] * f''(e^x) * e^{2x} + [ 1 ] * f'(e^x) * e^x
          d^5y/dx^5 = [ 1 ] * f'''''(e^x) * e^{5x} + [ 10 ] * f''''(e^x) * e^{4x} + [ 25 ] * f'''(e^x) * e^{3x} + [ 15 ] * f''(e^x) * e^{2x} + [ 1 ] * f'(e^x) * e^x
This pattern is observed in Mathematica for the first 10 cases, using the code below.
(End)

			The e.g.f. of [0,0,1,7,25,65,...], the k=3 column of A008278, but with offset n=0, is exp(x)*(1*(x^2)/2! + 4*(x^3)/3! + 3*(x^4)/4!).
Triangle starts:
  1;
  1,  1;
  1,  3,   1;
  1,  6,   7,    1;
  1, 10,  25,   15,    1;
  1, 15,  65,   90,   31,    1;
  1, 21, 140,  350,  301,   63,    1;
  1, 28, 266, 1050, 1701,  966,  127,   1;
  1, 36, 462, 2646, 6951, 7770, 3025, 255, 1;
  ...

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.
F. N. David, M. G. Kendall and D. E. Barton, Symmetric Function and Allied Tables, Cambridge, 1966, p. 223.
R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics, Addison-Wesley, 2nd ed., 1994.

T. D. Noe, Rows n = 0..100 of triangle, flattened
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].
Noureddine Chair, Exact two-point resistance, and the simple random walk on the complete graph minus N edges, Ann. Phys. 327, No. 12, 3116-3129 (2012), eq. (27).
Xi Chen, Bishal Deb, Alexander Dyachenko, Tomack Gilmore, and Alan D. Sokal, Coefficientwise total positivity of some matrices defined by linear recurrences, arXiv:2012.03629 [math.CO], 2020.
T. Copeland, Generators, Inversion, and Matrix, Binomial, and Integral Transforms, 2015
U. N. Katugampola, A new Fractional Derivative and its Mellin Transform, arXiv preprint arXiv:1106.0965 [math.CA], 2011.
Eric Weisstein's World of Mathematics, Bell Polynomial
Eric Weisstein's World of Mathematics, Empty Graph
Eric Weisstein's World of Mathematics, Independence Polynomial
Eric Weisstein's World of Mathematics, Sigma Polynomial
Eric Weisstein's World of Mathematics, Stirling Number of the Second Kind

See A008277 and A048993, which are the main entries for this triangle of numbers.

Cf. A094262, A008277, A008276, A003422, A000166, A000110, A000204, A000045, A000108, A213735.

a008278 n k = a008278_tabl !! (n-1) !! (k-1)
a008278_row n = a008278_tabl !! (n-1)
a008278_tabl = iterate st2 [1] where
  st2 row = zipWith (+) ([0] ++ row') (row ++ [0])
            where row' = reverse $ zipWith (*) [1..] $ reverse row
-- Reinhard Zumkeller, Jun 22 2013

rows = 10; Flatten[Table[StirlingS2[n, k], {n, 1, rows}, {k, n, 1, -1}]] (* Jean-François Alcover, Nov 17 2011, *)
Table[CoefficientList[x^n BellB[n, 1/x], x], {n, 10}] // Flatten (* Eric W. Weisstein, Apr 05 2017 *)
n = 5; Grid[Prepend[Transpose[{Range[1, n], Table[D[f[Exp[x]], {x, i}], {i, 1, n}]}], {"Order","Derivative"}], Frame -> All, Spacings -> {2, 1}] (* Leonidas Liponis, Aug 27 2024 *)

for(n=1,10,for(k=1,n,print1(stirling(n,n-k+1,2),", "))) \\ Hugo Pfoertner, Aug 30 2020

T(n, k)=0 if n < k, T(n, 0)=0, T(1, 1)=1, T(n, k) = (n-k+1)*T(n-1, k-1) + T(n-1, k) otherwise.
O.g.f. for the k-th column: 1/(1-x) if k=1 and A(k,x):=((x^k)/(1-x)^(2*k+1))*Sum_{m=0..k-1} A008517(k,m+1)*x^m if k >= 2. A008517 is the second-order Eulerian triangle. Cf. p. 257, eq. (6.43) of the R. L. Graham et al. book. - Wolfdieter Lang, Oct 14 2005
E.g.f. for the k-th column (with offset n=0): E(k,x):=exp(x)*Sum_{m=0..k-1} A112493(k-1,m)*(x^(k-1+m))/(k-1+m)! if k >= 1. - Wolfdieter Lang, Oct 14 2005
a(n) = abs(A213735(n-1)). - Hugo Pfoertner, Sep 07 2020

Name edited by Gus Wiseman, Aug 11 2020

A008306 Triangle T(n,k) read by rows: associated Stirling numbers of first kind (n >= 2, 1 <= k <= floor(n/2)).

Original entry on oeis.org

1, 2, 6, 3, 24, 20, 120, 130, 15, 720, 924, 210, 5040, 7308, 2380, 105, 40320, 64224, 26432, 2520, 362880, 623376, 303660, 44100, 945, 3628800, 6636960, 3678840, 705320, 34650, 39916800, 76998240, 47324376, 11098780, 866250, 10395

Offset: 2

N. J. A. Sloane

Also, T(n,k) is the number of derangements (permutations with no fixed points) of {1..n} with k cycles.

The sum of the n-th row is the n-th subfactorial: A000166(n). - Gary Detlefs, Jul 14 2010

			Rows 2 through 7 are:
    1;
    2;
    6,   3;
   24,  20;
  120, 130,  15;
  720, 924, 210;

L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 256.
J. Riordan, An Introduction to Combinatorial Analysis, Wiley, 1958, p. 75.

Reinhard Zumkeller, Rows n = 2..125 of table, flattened
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.
W. Carlitz, On some polynomials of Tricomi, Bollettino dell'Unione Matematica Italiana, Serie 3, Vol. 13, (1958), n. 1, p. 58-64
Tom Copeland, Generators, Inversion, and Matrix, Binomial, and Integral Transforms
W. Gautschi, The incomplete gamma functions since Tricomi (Cf. p. 206-207.)
P. Gniewek and B. Jeziorski, Convergence properties of the multipole expansion of the exchange contribution to the interaction energy, arXiv preprint arXiv:1601.03923 [physics.chem-ph], 2016.
S. Karlin and J. McGregor, Many server queuing processes with Poisson input and exponential service times, Pacific Journal of Mathematics, Vol. 8, No. 1, p. 87-118, March (1958). Cf. p. 117.
R. Paris, A uniform asymptotic expansion for the incomplete gamma function, Journal of Computational and Applied Mathematics, 148 (2002), p. 223-239. (See 333. From Tom Copeland, Jan 03 2016)
M. Z. Spivey, On Solutions to a General Combinatorial Recurrence, J. Int. Seq. 14 (2011) # 11.9.7.
N. Temme, A class of polynomials related to those of Laguerre
N. Temme, Traces to Tricomi in recent work on special functions and asymptotics of integrals
A. Topuzoglu, The Carlitz rank of permutations of finite fields: A survey, Journal of Symbolic Computation, Online, Dec 07, 2013.
Eric Weisstein's World of Mathematics, Permutation Cycle
Eric Weisstein's World of Mathematics, Stirling Number of the First Kind
Shawn L. Witte, Link Nomenclature, Random Grid Diagrams, and Markov Chain Methods in Knot Theory, Ph. D. Dissertation, University of California-Davis (2020).

Cf. A000166, A106828 (another version), A079510 (rearranged triangle), A235706 (specializations).
Diagonals: A000142, A000276, A000483.
Diagonals give reversed rows of A111999.

a008306 n k = a008306_tabf !! (n-2) !! (k-1)
a008306_row n = a008306_tabf !! (n-2)
a008306_tabf = map (fst . fst) $ iterate f (([1], [2]), 3) where
   f ((us, vs), x) =
     ((vs, map (* x) $ zipWith (+) ([0] ++ us) (vs ++ [0])), x + 1)
-- Reinhard Zumkeller, Aug 05 2013

A008306 := proc(n,k) local j;
add(binomial(j,n-2*k)*A008517(n-k,j),j=0..n-k) end;
seq(print(seq(A008306(n,k),k=1..iquo(n,2))),n=2..12):
# Peter Luschny, Apr 20 2011

t[0, 0] = 1; t[n_, 0] = 0; t[n_, k_] /; k > n/2 = 0; t[n_, k_] := t[n, k] = (n - 1)*(t[n - 1, k] + t[n - 2, k - 1]); A008306 = Flatten[ Table[ t[n, k], {n, 2, 12}, {k, 1, Quotient[n, 2]}]] (* Jean-François Alcover, Jan 25 2012, after David Callan *)

{ A008306(n,k) = (-1)^(n+k) * sum(i=0,k, (-1)^i * binomial(n,i) * stirling(n-i,k-i,1) ); } \\ Max Alekseyev, Sep 08 2018

T(n,k) = Sum_{i=0..k} (-1)^i * binomial(n,i) * |stirling1(n-i,k-i)| = (-1)^(n+k) * Sum_{i=0..k} (-1)^i * binomial(n,i) * A008275(n-i,k-i). - Max Alekseyev, Sep 08 2018
E.g.f.: 1 + Sum_{1 <= 2*k <= n} T(n, k)*t^n*u^k/n! = exp(-t*u)*(1-t)^(-u).
Recurrence: T(n, k) = (n-1)*(T(n-1, k) + T(n-2, k-1)) for 1 <= k <= n/2 with boundary conditions T(0,0) = 1, T(n,0) = 0 for n >= 1, and T(n,k) = 0 for k > n/2. - David Callan, May 16 2005
E.g.f. for column k: B(A(x)) where A(x) = log(1/1-x)-x and B(x) = x^k/k!.
From Tom Copeland, Jan 05 2016: (Start)
The row polynomials of this signed array are the orthogonal NL(n,x;x-n) = n! Sum_{k=0..n} binomial(x,n-k)*(-x)^k/k!, the normalized Laguerre polynomials of order (x-n) as discussed in Gautschi (the Temme, Carlitz, and Karlin and McGregor references come from this paper) in regard to asymptotic expansions of the upper incomplete gamma function--Tricomi's Cinderella of special functions.
e^(x*t)*(1-t)^x = Sum_{n>=0} NL(n,x;x-n)*x^n/n!.
The first few are
NL(0,x) = 1
NL(1,x) = 0
NL(2,x) = -x
NL(3,x) = 2*x
NL(4,x) = -6*x + 3*x^2.
With D=d/dx, :xD:^n = x^n D^n, :Dx:^n = D^n x^n, and K(a,b,c), the Kummer confluent hypergeometric function, NL(n,x;y-n) = n!*e^x binomial(xD+y,n)*e^(-x) = n!*e^x Sum_{k=0..n} binomial(k+y,n) (-x)^k/k! = e^x x^(-y+n) D^n (x^y e^(-x)) = e^x x^(-y+n) :Dx:^n x^(y-n)*e^(-x) = e^x*x^(-y+n)*n!*L(n,:xD:,0)*x^(y-n)*e^(-x) = n! binomial(y,n)*K(-n,y-n+1,x) = n!*e^x*(-1)^n*binomial(-xD-y+n-1,n)*e^(-x). Evaluate these expressions at y=x after the derivative operations to obtain NL(n,x;x-n). (End)

More terms from Larry Reeves (larryr(AT)acm.org), Feb 16 2001

A068106 Euler's difference table: triangle read by rows, formed by starting with factorial numbers (A000142) and repeatedly taking differences. T(n,n) = n!, T(n,k) = T(n,k+1) - T(n-1,k).

Original entry on oeis.org

1, 0, 1, 1, 1, 2, 2, 3, 4, 6, 9, 11, 14, 18, 24, 44, 53, 64, 78, 96, 120, 265, 309, 362, 426, 504, 600, 720, 1854, 2119, 2428, 2790, 3216, 3720, 4320, 5040, 14833, 16687, 18806, 21234, 24024, 27240, 30960, 35280, 40320, 133496, 148329, 165016, 183822, 205056, 229080, 256320, 287280, 322560, 362880

Offset: 0

N. J. A. Sloane, Apr 12 2002

Triangle T(n,k) (n >= 1, 1 <= k <= n) giving number of ways of winning with (n-k+1)st card in the generalized "Game of Thirteen" with n cards.
From Emeric Deutsch, Apr 21 2009: (Start)
T(n-1,k-1) is the number of non-derangements of {1,2,...,n} having largest fixed point equal to k. Example: T(3,1)=3 because we have 1243, 4213, and 3241.
Mirror image of A047920.
(End)

			Triangle begins:
[0]    1;
[1]    0,    1;
[2]    1,    1,    2;
[3]    2,    3,    4,    6;
[4]    9,   11,   14,   18,   24;
[5]   44,   53,   64,   78,   96,  120;
[6]  265,  309,  362,  426,  504,  600,  720;
[7] 1854, 2119, 2428, 2790, 3216, 3720, 4320, 5040.

Reinhard Zumkeller, Rows n = 0..150 of triangle, flattened
W. Y. C. Chen et al., Higher-order log-concavity in Euler's difference table, Discrete Math., 311 (2011), 2128-2134.
P. R. de Montmort, On the Game of Thirteen (1713), reprinted in Annotated Readings in the History of Statistics, ed. H. A. David and A. W. F. Edwards, Springer-Verlag, 2001, pp. 25-29.
Emeric Deutsch and S. Elizalde, The largest and the smallest fixed points of permutations, arXiv:0904.2792 [math.CO], 2009.
D. Dumont, Matrices d'Euler-Seidel, Sem. Loth. Comb. B05c (1981) 59-78.
Philip Feinsilver and John McSorley, Zeons, Permanents, the Johnson scheme, and Generalized Derangements, arXiv:1710.00788 [math.CO], (2017); see page 29.
P. Feinsilver and J. McSorley, Zeons, Permanents, the Johnson scheme, and Generalized Derangements, International Journal of Combinatorics, 2011 (2011).
Fanja Rakotondrajao, k-Fixed-Points-Permutations, Integers: Electronic journal of combinatorial number theory 7 (2007) A36.
Index entries for sequences related to factorial numbers

Row sums give A002467.
Diagonals give A000142, A001563, A001564, A001565, A001688, A001689, A023043, A023044, A023045, A023046, A023047 (factorials and k-th differences, k=1..10).
See A047920 and A086764 for other versions.
T(2*n, n) is A033815.
Columns k=0..10 give A000166, A000255, A055790, A277609, A277563, A280425, A280920, A284204, A284205, A284206, A284207.

a068106 n k = a068106_tabl !! n !! k
a068106_row n = a068106_tabl !! n
a068106_tabl = map reverse a047920_tabl
-- Reinhard Zumkeller, Mar 05 2012

d[0] := 1: for n to 15 do d[n] := n*d[n-1]+(-1)^n end do: T := proc (n, k) if k <= n then sum(binomial(k, j)*d[n-j], j = 0 .. k) else 0 end if end proc: for n from 0 to 9 do seq(T(n, k), k = 0 .. n) end do; # yields sequence in triangular form; Emeric Deutsch, Jul 18 2009

t[n_, k_] := Sum[(-1)^j*Binomial[n-k, j]*(n-j)!, {j, 0, n}]; Flatten[ Table[ t[n, k], {n, 0, 9}, {k, 0, n}]] (* Jean-François Alcover, Feb 21 2012, after Philippe Deléham *)
T[n_, k_] := n! HypergeometricPFQ[{k-n}, {-n}, -1];
Table[T[n, k], {n,0,9}, {k,0,n}] // Flatten (* Peter Luschny, Oct 05 2017 *)

T(n, k) = Sum_{j>= 0} (-1)^j*binomial(n-k, j)*(n-j)!. - Philippe Deléham, May 29 2005
From Emeric Deutsch, Jul 18 2009: (Start)
T(n,k) = Sum_{j=0..k} d(n-j)*binomial(k, j), where d(i) = A000166(i) are the derangement numbers.
Sum_{k=0..n} (k+1)*T(n,k) = A000166(n+2) (the derangement numbers). (End)
T(n, k) = n!*hypergeom([k-n], [-n], -1). - Peter Luschny, Oct 05 2017
D-finite recurrence for columns: T(n,k) = n*T(n-1,k) + (n-k)*T(n-2,k). - Georg Fischer, Aug 13 2022

More terms from Antonio G. Astudillo (afg_astudillo(AT)lycos.com), Apr 01 2003

Edited by N. J. A. Sloane, Sep 24 2011

A188674 Stack polyominoes with square core.

Original entry on oeis.org

1, 1, 0, 0, 1, 2, 3, 4, 5, 7, 9, 13, 17, 24, 31, 42, 54, 71, 90, 117, 147, 188, 236, 298, 371, 466, 576, 716, 882, 1088, 1331, 1633, 1987, 2422, 2935, 3557, 4290, 5177, 6216, 7465, 8932, 10682, 12731, 15169, 18016, 21387, 25321, 29955, 35353, 41696, 49063, 57689, 67698, 79375, 92896, 108633, 126817, 147922, 172272

Offset: 0

Emanuele Munarini, Apr 08 2011

nonn

a(n) is the number of stack polyominoes of area n with square core.
The core of stack is the set of all maximal columns.
The core is a square when the number of columns is equal to their height.
Equivalently, a(n) is the number of unimodal compositions of n, where the number of the parts of maximum value equal the maximum value itself. For instance, for n = 10, we have the following stacks:
(1,3,3,3), (3,3,3,1), (1,1,1,1,1,1,2,2), (1,1,1,1,1,2,2,1), (1,1,1,1,2,2,1,1), (1,1,1,2,2,1,1,1), (1,1,2,2,1,1,1,1), (1,2,2,1,1,1,1,1), (2,2,1,1,1,1,1,1).
From Gus Wiseman, Apr 06 2019 and May 21 2022: (Start)
Also the number of integer partitions of n with final part in their inner lining partition equal to 1, where the k-th part of the inner lining partition of a partition is the number of squares in its Young diagram that are k diagonal steps from the lower-right boundary. For example, the a(4) = 1 through a(10) = 9 partitions are:
  (22)  (32)   (42)    (52)     (62)      (72)       (82)
        (221)  (321)   (421)    (521)     (333)      (433)
               (2211)  (3211)   (4211)    (621)      (721)
                       (22111)  (32111)   (5211)     (3331)
                                (221111)  (42111)    (6211)
                                          (321111)   (52111)
                                          (2211111)  (421111)
                                                     (3211111)
                                                     (22111111)
Also partitions that have a fixed point and a conjugate fixed point, ranked by A353317. The strict case is A352829. For example, the a(0) = 0 through a(9) = 7 partitions are:
  ()  .  .  (21)  (31)   (41)    (51)     (61)      (71)
                  (211)  (311)   (411)    (511)     (332)
                         (2111)  (3111)   (4111)    (611)
                                 (21111)  (31111)   (5111)
                                          (211111)  (41111)
                                                    (311111)
                                                    (2111111)
Also partitions of n + 1 without a fixed point or conjugate fixed point.
(End)

Alois P. Heinz, Table of n, a(n) for n = 0..1000
Brian Hopkins, James A. Sellers, and Dennis Stanton, Dyson's Crank and the Mex of Integer Partitions, arXiv:2009.10873 [math.CO], 2020. Mentions this sequence.

Cf. A001523 (stacks).
Cf. A006918, A252464, A325135, A325163, A325165.
Positive crank: A001522, ranked by A352874.
Zero crank: A064410, ranked by A342192.
Nonnegative crank: A064428, ranked by A352873.
Fixed point but no conjugate fixed point: A118199, ranked by A353316.
A000041 counts partitions, strict A000009.
A002467 counts permutations with a fixed point, complement A000166.
A115720/A115994 count partitions by Durfee square, rank statistic A257990.
A238352 counts reversed partitions by fixed points, rank statistic A352822.
A238394 counts reversed partitions without a fixed point, ranked by A352830.
A238395 counts reversed partitions with a fixed point, ranked by A352872.
A352833 counts partitions by fixed points.
Cf. A000700, A088902, A114088, A300788, A330644, A352828, A352829, A352832.

a[n_]:=CoefficientList[Series[1+Sum[x^((k+1)^2)/Product[(1-x^i)^2,{i,1,k}],{k,0,n}],{x,0,n}],x]
(* second program *)
pml[ptn_]:=If[ptn=={},{},FixedPointList[If[#=={},{},DeleteCases[Rest[#]-1,0]]&,ptn][[-3]]];
Table[Length[Select[IntegerPartitions[n],pml[#]=={1}&]],{n,0,30}] (* Gus Wiseman, Apr 06 2019 *)

G.f.: 1 + sum(k>=0, x^((k+1)^2)/((1-x)^2*(1-x^2)^2*...*(1-x^k)^2)).

A055209 a(n) = Product_{i=0..n} i!^2.

Original entry on oeis.org

1, 1, 4, 144, 82944, 1194393600, 619173642240000, 15728001190723584000000, 25569049282962188245401600000000, 3366980847587422591723894776791040000000000, 44337041641882947649156022595410930014617600000000000000

Offset: 0

N. J. A. Sloane, Jul 18 2000

a(n) is the discriminant of the polynomial x(x+1)(x+2)...(x+n). - Yuval Dekel (dekelyuval(AT)hotmail.com), Nov 13 2003
This is the Hankel transform (see A001906 for definition) of the sequence: 1, 0, 1, 0, 5, 0, 61, 0, 1385, 0, 50521, ... (see A000364: Euler numbers). - Philippe Deléham, Apr 06 2005
Also, for n>0, the quotient of (-1)^(n-1)S(u)^(n^2)/S(un) and the determinant of the (n-1) X (n-1) square matrix [P'(u), P''(u), ..., P^(n-1)(u); P''(u), P'''(u), ..., P^(n)(u); P'''(u), P^(4)(u), ..., P^(n+1)(u); ...; P^(n-1)(u), P^(n)(u), ..., P^(2n-3)(u)] where S and P are the Weierstrass Sigma and The Weierstrass P-function, respectively and f^(n) is the n-th derivative of f. See the King and Schwarz & Weierstrass references. - Balarka Sen, Jul 31 2013
a(n) is the number of idempotent monotonic labeled magmas. That is, prod(i,j) >= max(i,j) and prod(i,i) = i. - Chad Brewbaker, Nov 03 2013
Ramanujan's infinite nested radical sqrt(1+2*sqrt(1+3*sqrt(1+...))) = 3 can be written sqrt(1+sqrt(4+sqrt(144+...))) = sqrt(a(1)+sqrt(a(2)+sqrt(a(3)+...))). Vijayaraghavan used that to prove convergence of Ramanujan's formula. - Petros Hadjicostas and Jonathan Sondow, Mar 22 2014
a(n) is the determinant of the (n+1)-th order Hankel matrix whose (i,j)-entry is equal to A000142(i+j), i,j = 0,1,...,n. - Michael Shmoish, Sep 02 2020

R. Bruce King, Beyond The Quartic Equation, Birkhauser Boston, Berlin, 1996, p. 72.
Srinivasa Ramanujan, J. Indian Math. Soc., III (1911), 90 and IV (1912), 226.
T. Vijayaraghavan, in Collected Papers of Srinivasa Ramanujan, G.H. Hardy, P.V. Seshu Aiyar and B.M. Wilson, eds., Cambridge Univ. Press, 1927, p. 348; reprinted by Chelsea, 1962.

G. C. Greubel, Table of n, a(n) for n = 0..32
Paul Barry, A Note on Three Families of Orthogonal Polynomials defined by Circular Functions, and Their Moment Sequences, Journal of Integer Sequences, Vol. 15 (2012), #12.7.2.
Richard Ehrenborg, The Hankel determinant of exponential polynomials, Amer. Math. Monthly, Vol. 107, No. 6 (2000), pp. 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.
John W. Layman, The Hankel Transform and Some of its Properties, J. Integer Sequences, 4 (2001), #01.1.5.
Christian Radoux, Déterminants de Hankel et théorème de Sylvester, Séminaire Lotharingien de Combinatoire, B28b (1992), 9 pp.
H. A. Schwarz and K. Weierstrass, Formeln und Lehrsätze zum Gebrauche der elliptischen Functionen, Springer, Berlin, 1893, p. 19.
Jonathan Sondow and Petros Hadjicostas, The generalized-Euler-constant function gamma(z) and a generalization of Somos's quadratic recurrence constant, J. Math. Anal. Appl., Vol. 332, No. 1 (2007), pp. 292-314; see pp. 305-306.
Michael Z. Spivey and Laura L. Steil, The k-Binomial Transforms and the Hankel Transform, Journal of Integer Sequences, Vol. 9 (2006), Article 06.1.1.
Index entries for sequences related to factorial numbers

Cf. A055209 is the Hankel transform (see A001906 for definition) of A000023, A000142, A000166, A000522, A003701, A010842, A010843, A051295, A052186, A053486, A053487.

Cf. A112302, A258619.

[1] cat [(&*[(Factorial(k))^2: k in [1..n]]): n in [1..10]]; // G. C. Greubel, Oct 14 2018

seq(mul(mul(j^2,j=1..k), k=0..n), n=0..10); # Zerinvary Lajos, Sep 21 2007

Table[Product[(i!)^2,{i,n}],{n,0,11}] (* Harvey P. Dale, Jul 06 2011 *)
Table[BarnesG[n + 2]^2, {n, 0, 11}] (* Jan Mangaldan, May 07 2014 *)

a(n)=prod(i=1,n,i!)^2 \\ Charles R Greathouse IV, Jan 12 2012

def A055209(n) :
   return prod(factorial(i)^(2) for i in (0..n))
[A055209(n) for n in (0..11)] # Jani Melik, Jun 06 2015

a(n) = A000178(n)^2. - Philippe Deléham, Mar 06 2004
a(n) = Product_{i=0..n} i^(2*n - 2*i + 2). - Charles R Greathouse IV, Jan 12 2012
Asymptotic: a(n) ~ exp(2*zeta'(-1)-3/2*(1+n^2)-3*n)*(2*Pi)^(n+1)*(n+1)^ (n^2+2*n+5/6). - Peter Luschny, Jun 23 2012
lim_{n->infinity} a(n)^(2^(-(n+1))) = 1. - Vaclav Kotesovec, Jun 06 2015
Sum_{n>=0} 1/a(n) = A258619. - Amiram Eldar, Nov 17 2020

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

Original entry on oeis.org

1, 0, 2, 4, 24, 128, 880, 6816, 60032, 589312, 6384384, 75630080, 972387328, 13483769856, 200571078656, 3185540657152, 53800242216960, 962741176500224, 18195808235880448, 362183230599856128, 7572922094360723456, 165945771111208714240, 3802923921298533384192, 90965940197460917878784, 2267151124921333646884864

Offset: 0

Gordon F. Royle, Oct 28 2003

Permanent of an (n+1) X (n+1) (+1, -1)-matrix with exactly n -1's on the diagonal and 1's everywhere else.
It is conjectured by Kräuter and Seifter that for n >= 5 a(n-1) is the maximal possible value for the permanent of a nonsingular n X n (+1, -1)-matrix. I do not know for which values of n this has been confirmed - compare A087982. - N. J. A. Sloane
The Kräuter conjecture on permanents is true (see Budrevich and Guterman). - Sergei Shteiner, Jan 17 2020
The maximal possible value for the permanent of a singular n X n (+1, -1)-matrix is obviously n!.
Degree of the "hyperdeterminant" of a multilinear polynomial on (\P^1(\C))^n, or equivalently of an element of (\C^2)^{⊗ n}: see Gelfand, Kapranov and Zelevinsky. - Eric Rains, Mar 15 2004
(-1)^n * a(n) = Polynomials in A010027 evaluated at -1. - Ralf Stephan, Dec 15 2004
a(n) is the number of n X n (-1, 0, 1)-matrices containing in every row and every column exactly one -1 and one 1 such that the main diagonal does not contain 0's. - Vladimir Shevelev, Apr 01 2010
a(n) is the number of colored permutations with no fixed points of n elements where each cycle is one of two colors. - Michael Somos, Jan 19 2011
Binomial transform is A000255. Hankel transform is A059332. - Paul Barry, Apr 11 2011
Exponential self-convolution of subfactorials (A000166). - Vladimir Reshetnikov, Oct 07 2016

			G.f. = 1 + 2*x^2 + 4*x^3 + 24*x^4 + 128*x^5 + 880*x^6 + 6816*x^7 + ...
Since a(1) = 0, then, for n = 2, we have a(2) = -(-2)^3/4 = 2; further, for n = 3, we find a(3) = (3*6/5)*2 - (-2)^4/5 = 36/5 - 16/5 = 4. - _Vladimir Shevelev_, Apr 01 2010
a(4) = 24 because there are 6 derangements with one 4-cycle with 2^1 ways to color each derangement and 3 derangements with two 2-cycles with 2^2 ways to color each derangement. - _Michael Somos_, Jan 19 2011

I. M. Gelfand, M. M. Kapranov, and A. V. Zelevinsky, Discriminants, Resultants and Multidimensional Determinants, Birkhauser, 1994; see Corollary 2.10 in Chapter 14 (p. 457).

T. D. Noe, Table of n, a(n) for n = 0..100
Mikhail V. Budrevich and Alexander E. Guterman, Kräuter conjecture on permanents is true, arXiv:1810.04439 [math.CO], 2018.
Steven Finch, Rounds, Color, Parity, Squares, arXiv:2111.14487 [math.CO], 2021.
I. M. Gelfand, M. M. Kapranov, and A. V. Zelevinsky, Hyperdeterminants, Advances in Mathematics 96(2) (1992), 226-263; see Corollary 3.10 (p. 246).
Arnold R. Kräuter, Permanenten - Ein kurzer Überblick, Séminaire Lotharingien de Combinatoire, B09b (1983), 34 pp.
Arnold R. Kräuter and Norbert Seifter, Some properties of the permanent of (1,-1)-matrices, Linear and Multilinear Algebra 15 (1984), 207-223.
Norbert Seifter, Upper bounds for permanents of (1,-1)-matrices, Israel J. Math. 48 (1984), 69-78.
Edward Tzu-Hsia Wang, On permanents of (1,-1)-matrices, Israel J. Math. 18 (1974), 353-361.
Wikipedia, Hyperdeterminant
Index entries for sequences related to binary matrices

Cf. A087982, A087983, A176097.

Cf. A000007, A000023, A000166, A024000, A137775. - Michael Somos, Jan 19 2011

seq(simplify(KummerU(-n, -n-1, -2)), n = 0..24); # Peter Luschny, May 10 2022

Range[0, 20]! CoefficientList[Series[Exp[-2 x]/(1 - x)^2, {x, 0, 20}], x]
Table[(-2)^n HypergeometricPFQ[{2, -n}, {}, 1/2], {n, 0, 20}] (* Vladimir Reshetnikov, Oct 07 2016 *)

{a(n) = if( n<0, 0, n! * polcoeff( exp( -2 * x + x * O(x^n) ) / ( 1 - x )^2, n ) )} /* Michael Somos, Jan 19 2011 */

Krauter and Seifter prove that the permanent of an n X n {-1, 1} matrix is divisible by 2^{n - [log_2(n)] - 1}.
Let c(n) be the permanent of the {-1, +1}-matrix of order n X n with n diagonal -1's only. Let a(n) be the permanent of the {-1, +1}-matrix of order (n+1) X (n+1) with n diagonal -1's only. Then by expanding along the first row (like determinant, but with no sign) we get c(n+1) = -c(n) + n a(n-1), a(n) = c(n) + n a(n-1), with c(2) = 2, a(2) = 2. {c(n)} has e.g.f. exp(-2x)/(1-x), see A000023. Also a(n) = c(n+1) + 2*c(n).
The following 4 formulas hold: a(n) = Sum_{k = 0..n} C(n, k)*D_k*D_{n-k}, where D_n = A000166(n); a(n) = n!*Sum_{j = 0..n} (n+1-j)*(-2)^j/j!; a(0) = 1, a(1) = 0 and, for n > 0, a(n+1) = n*(a(n) + 2*a(n-1)); a(0) = 1 and, for n > 0, a(n) = (n*(n+3)/(n+2))*a(n-1) - (-2)^(n+1)/(n+2). - Vladimir Shevelev, Apr 01 2010 [edited by Michael Somos, Jan 19 2011]
G.f.: 1/(1-2x^2/(1-2x-6x^2/(1-4x-12x^2/(1-6x-20x^2/(1-.../(1-2n*x-(n+1)(n+2)x^2/(1-... (continued fraction). - Paul Barry, Apr 11 2011
E.g.f.: 1/U(0) where U(k)=  1 - 2*x/( 1 + x/(2 - x - 4/( 2 + x*(k+1)/U(k+1)))) ; (continued fraction, 3rd kind, 4-step). - Sergei N. Gladkovskii, Oct 28 2012
G.f.: 1/Q(0), where Q(k) = 1 + 2*x - x*(k+2)/(1 - x*(k+1)/Q(k+1)); (continued fraction). - Sergei N. Gladkovskii, Apr 22 2013
G.f.: 1/Q(0) where Q(k) = 1 - 2*k*x - x^2*(k + 1)*(k+2)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, May 10 2013
G.f.: S(x)/x - 1/x = G(0)/x - 1/x, where S(x) = sum(k >= 0, k!*(x/(1+2*x))^k  ), G(k) = 1 + (2*k + 1)*x/( 1+2*x - 2*x*(1+2*x)*(k+1)/(2*x*(k+1) + (1+2*x)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 26 2013
a(n) = (-2)^n*hypergeom([2, -n], [], 1/2) = 4*(-2)^n*(1 - 2*hypergeom([1, -n-3], [], 1/2))/(n^2+3*n+2) = (4*(-2)^n + Gamma(n+4, -2)*exp(-2))/(n^2+3*n+2). - Vladimir Reshetnikov, Oct 07 2016
a(n) ~ sqrt(2*Pi) * n^(n+3/2) / exp(n+2). - Vaclav Kotesovec, Oct 08 2016
a(n) = KummerU(-n, -n - 1, -2). - Peter Luschny, May 10 2022

More terms from Jaap Spies, Oct 28 2003
Further terms from Gordon F. Royle, Oct 29 2003
Definition via e.g.f. from Eric Rains, Mar 15 2004
Changed the offset and terms to correspond to e.g.f, Michael Somos, Jan 19 2011

A180191 Number of permutations of [n] having at least one succession. A succession of a permutation p is a position i such that p(i+1)-p(i) = 1.

Original entry on oeis.org

0, 1, 3, 13, 67, 411, 2921, 23633, 214551, 2160343, 23897269, 288102189, 3760013027, 52816397219, 794536751217, 12744659120521, 217140271564591, 3916221952414383, 74539067188152941, 1493136645424092773, 31400620285465593339, 691708660911435955579

Offset: 1

Emeric Deutsch, Sep 07 2010

nonn

a(n) = A180190(n,1).
a(n+2) = p(n+2) where p(x) is the unique degree-n polynomial such that p(k) = k! for k = 1, ..., n+1. - Michael Somos, Jan 05 2012
From Jon Perry, Jan 04 2013: (Start)
Number of permutations of {1,...,n-1,n+1} with at least one indexed point p(k)=k with 1<=k<=n. Note that this means p(k)=n+1 is never an indexed point as kFor n>1, a(n) is the number of permutations of [n+1] that have a fixed point and contain 12; for example the a(3)=3 such permutations of {1,2,3,4} are 1234, 1243, and 3124.
(End)
For n > 0: row sums of triangle A116853. - Reinhard Zumkeller, Aug 31 2014

			x^2 + 3*x^3 + 13*x^4 + 67*x^5 + 411*x^6 + 2921*x^7 + 23633*x^8 + ...
a(3) = 3 because we have 123, 312, and 231; the permutations 132, 213, and 321 have no successions.
a(4) = 13 since p(x) = (3*x^2 - 7*x + 6) / 2 interpolates p(1) = 1, p(2) = 2, p(3) = 6, and p(4) = 13. - _Michael Somos_, Jan 05 2012

Alois P. Heinz, Table of n, a(n) for n = 1..200
Uriel Feige, Tighter bounds for online bipartite matching, 2018.

Cf. A000142, A000166, A000255, A002467, A028417, A180190.
Cf. A116853, A276975.
Column k=1 of A306234, A306461, and of A324362(n-1).

a180191 n = if n == 1 then 0 else sum $ a116853_row (n - 1)
-- Reinhard Zumkeller, Aug 31 2014

d[0] := 1: for n to 50 do d[n] := n*d[n-1]+(-1)^n end do: seq(factorial(n)-d[n]-d[n-1], n = 1 .. 22);

f[n_] := Sum[ -(-1)^k (n - k)! Binomial[n - 1, k], {k, 1, n}]; Array[f, 20] (* Robert G. Wilson v, Oct 16 2010 *)

{a(n) = if( n<2, 0, n--; subst( polinterpolate( vector( n, k, k!)), x, n+1))} /* Michael Somos, Jan 05 2012 */

a(n) = n! - d(n) - d(n-1), where d(j) = A000166(j) are the derangement numbers.
a(n) = n! - A000255(n-1) = A002467(n) - A000166(n-1). - Jon Perry, Jan 05 2013
a(n) = (n-1)! [x^(n-1)] (1-exp(-x))/(1-x)^2. - Alois P. Heinz, Feb 23 2019

A109043 a(n) = lcm(n,2).

Original entry on oeis.org

0, 2, 2, 6, 4, 10, 6, 14, 8, 18, 10, 22, 12, 26, 14, 30, 16, 34, 18, 38, 20, 42, 22, 46, 24, 50, 26, 54, 28, 58, 30, 62, 32, 66, 34, 70, 36, 74, 38, 78, 40, 82, 42, 86, 44, 90, 46, 94, 48, 98, 50, 102, 52, 106, 54, 110, 56, 114, 58, 118, 60, 122, 62, 126, 64, 130, 66, 134

Offset: 0

Mitch Harris, Jun 18 2005

Exponent of the dihedral group D(2n) = . - Arkadiusz Wesolowski, Sep 10 2013
Second column of table A210530. - Boris Putievskiy, Jan 29 2013
For n > 1, the basic period of A000166(k) (mod n) (Miska, 2016). - Amiram Eldar, Mar 03 2021

Vincenzo Librandi, Table of n, a(n) for n = 0..10000
Dorin Andrica, Sorin Rădulescu, and George Cătălin Ţurcaş, The Exponent of a Group: Properties, Computations and Applications, Disc. Math. and Applications, Springer, Cham (2020), 57-108.
Piotr Miska, Arithmetic properties of the sequence of derangements, Journal of Number Theory, Vol. 163 (2016), pp. 114-145; arXiv preprint, arXiv:1508.01987 [math.NT], 2015. See p. 124 (p. 14 in the preprint).
Index entries for linear recurrences with constant coefficients, signature (0,2,0,-1).
Index entries for sequences related to lcm's.

Cf. A000166, A109042, A152749 (partial sums).
Cf. A066043 (essentially the same), A000034 (=a(n)/n), A026741 (=a(n)/2).
Cf. A186421, A210530.

a109043 = (lcm 2)
a109043_list = zipWith (*) [0..] a000034_list
-- Reinhard Zumkeller, Mar 31 2012

[0, 2, 2] cat [Exponent(DihedralGroup(n)) : n in [3..65]]; // Arkadiusz Wesolowski, Sep 10 2013

LCM[Range[0,70],2] (* Harvey P. Dale, Aug 19 2012 *)

a(n)=lcm(n,2) \\ Charles R Greathouse IV, Sep 24 2015

def A109043(n): return n<<1 if n&1 else n # Chai Wah Wu, Aug 05 2024

[lcm(n,2) for n in range(0, 68)] # Zerinvary Lajos, Jun 07 2009

For n > 0: a(n) = A186421(n) + A186421(n+2).
a(n) = n*2 / gcd(n, 2).
a(n) = -(n*((-1)^n-3))/2. - Stephen Crowley, Feb 11 2007
From R. J. Mathar, Aug 20 2008: (Start)
a(n) = A066043(n), n > 1.
a(n) = 2*A026741(n).
G.f.: 2*x(1+x+x^2)/((1-x)^2*(1+x)^2). (End)
a(n) = n*A000034(n). - Paul Curtz, Mar 25 2011
E.g.f.: x*(2*cosh(x) + sinh(x)). - Stefano Spezia, May 09 2021
Sum_{k=1..n} a(k) ~ (3/4) * n^2. - Amiram Eldar, Nov 26 2022

A086764 Triangle T(n, k), read by row, related to Euler's difference table A068106 (divide column k of A068106 by k!).

Original entry on oeis.org

1, 0, 1, 1, 1, 1, 2, 3, 2, 1, 9, 11, 7, 3, 1, 44, 53, 32, 13, 4, 1, 265, 309, 181, 71, 21, 5, 1, 1854, 2119, 1214, 465, 134, 31, 6, 1, 14833, 16687, 9403, 3539, 1001, 227, 43, 7, 1, 133496, 148329, 82508, 30637, 8544, 1909, 356, 57, 8, 1

Offset: 0

Philippe Deléham, Aug 02 2003

The k-th column sequence, k >= 0, without leading zeros, enumerates the ways to distribute n beads, n >= 1, labeled differently from 1 to n, over a set of (unordered) necklaces, excluding necklaces with exactly one bead, and k+1 indistinguishable, ordered, fixed cords, each allowed to have any number of beads. Beadless necklaces as well as beadless cords each contribute a factor 1, hence for n=0 one has 1. See A000255 for the description of a fixed cord with beads. This comment derives from a family of recurrences found by Malin Sjodahl for a combinatorial problem for certain quark and gluon diagrams (Feb 27 2010). - Wolfdieter Lang, Jun 02 2010

			Formatted as a square array:
      1      3     7    13   21   31  43 57 ... A002061;
      2     11    32    71  134  227 356    ... A094792;
      9     53   181   465 1001 1909        ... A094793;
     44    309  1214  3539 8544             ... A094794;
    265   2119  9403 30637                  ... A023043;
   1854  16687 82508                        ... A023044;
  14833 148329                              ... A023045;
Formatted as a triangular array (mirror of A076731):
       1;
       0      1;
       1      1     1;
       2      3     2     1;
       9     11     7     3    1;
      44     53    32    13    4    1;
     265    309   181    71   21    5    1;
    1854   2119  1214   465  134   31    6   1;
   14833  16687  9403  3539 1001  227   43   7   1;
  133496 148329 82508 30637 8544 1909  356  57   8   1;

Indranil Ghosh, Rows 0..50, flattened
W. Y. C. Chen et al., Higher-order log-concavity in Euler's difference table, Discrete Math., 311 (2011), 2128-2134. (These are the numbers d^k_n.)
Fanja Rakotondrajao, k-Fixed-Points-Permutations, Integers: Electronic journal of combinatorial number theory 7 (2007) A36.

Columns: A000166, A000155, A000153, A000261, A001909, A001910, A176732, A176733, A176734, A176735, A176736.
Mirror image of A076731.
Cf. A002061, A003470, A008288, A023043, A023044, A023045, A068106.
Cf. A076731, A094792, A094793, A094794, A095000, A108625, A143007.
Cf. A143409, A143410, A143411, A143413.

A086764:= func< n,k | (&+[(-1)^j*Binomial(n-k,j)*Factorial(n-j): j in [0..n]])/Factorial(k) >;
[A086764(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Oct 05 2023

T[n_,k_]:=(1/k!)*Sum[(-1)^j*Binomial[n-k,j]*(n-j)!,{j,0,n}];Flatten[Table[T[n,k],{n,0,11},{k,0,n}]] (* Indranil Ghosh, Feb 20 2017 *)
T[n_, k_] := (n!/k!) HypergeometricPFQ[{k-n},{-n},-1];
Table[T[n,k], {n,0,9}, {k,0,n}] // Flatten (* Peter Luschny, Oct 05 2017 *)

def A086764(n,k): return sum((-1)^j*binomial(n-k,j)*factorial(n-j) for j in range(n+1))//factorial(k)
flatten([[A086764(n,k) for k in range(n+1)] for n in range(13)]) # G. C. Greubel, Oct 05 2023

T(n, n) = 1; T(n+1, n) = n.
T(n+2, n) = A002061(n+1) = n^2 + n + 1; T(n+3, n) = n^3 + 3*n^2 + 5*n + 2.
T(n, k) = (k + 1)*T(n, k + 1) - T(n-1, k); T(n, n) = 1; T(n, k) = 0, if k > n.
T(n, k) = (n-1)*T(n-1, k) + (n-k-1)*T(n-2, k).
k!*T(n, k) = A068106(n, k). [corrected by Georg Fischer, Aug 13 2022]
Sum_{k>=0} T(n, k) = A003470(n+1).
T(n, k) = (1/k!) * Sum_{j>=0} (-1)^j*binomial(n-k, j)*(n-j)!. - Philippe Deléham, Jun 13 2005
From Peter Bala, Aug 14 2008: (Start)
The following remarks all relate to the array read as a square array: e.g.f for column k: exp(-y)/(1-y)^(k+1); e.g.f. for array: exp(-y)/(1-x-y) = (1 + x + x^2 + x^3 + ...) + (x + 2*x^2 + 3*x^3 + 4*x^4 + ...)*y + (1 + 3*x + 7*x^2 + 13*x^3 + ...)*y^2/2! + ... .
This table is closely connected to the constant e. The row, column and diagonal entries of this table occur in series formulas for e.
Row n for n >= 2: e = n!*(1/T(n,0) + (-1)^n*[1/(1!*T(n,0)*T(n,1)) + 1/(2!*T(n,1)*T(n,2)) + 1/(3!*T(n,2)*T(n,3)) + ...]). For example, row 3 gives e = 6*(1/2 - 1/(1!*2*11) - 1/(2!*11*32) - 1/(3!*32*71) - ...). See A095000.
Column 0: e = 2 + Sum_{n>=2} (-1)^n*n!/(T(n,0)*T(n+1,0)) = 2 + 2!/(1*2) - 3 !/(2*9) + 4!/(9*44) - ... .
Column k, k >= 1: e = (1 + 1/1! + 1/2! + ... + 1/k!) + 1/k!*Sum_{n >= 0} (-1)^n*n!/(T(n,k)*T(n+1,k)). For example, column 3 gives e = 8/3 + 1/6*(1/(1*3) - 1/(3*13) + 2/(13*71) - 6/(71*465) + ...).
Main diagonal: e = 1 + 2*(1/(1*1) - 1/(1*7) + 1/(7*71) - 1/(71*1001) + ...).
First subdiagonal: e = 8/3 + 5/(3*32) - 7/(32*465) + 9/(465*8544) - ... .
Second subdiagonal: e = 2*(1 + 2^2/(1*11) - 3^2/(11*181) + 4^2/(181*3539) - ...). See A143413.
Third subdiagonal: e = 3 - (2*3*5)/(2*53) + (3*4*7)/(53*1214) - (4*5*9)/(1214*30637) + ... .
For the corresponding results for the constants 1/e, sqrt(e) and 1/sqrt(e) see A143409, A143410 and A143411 respectively. For other arrays similarly related to constants see A008288 (for log(2)), A108625 (for zeta(2)) and A143007 (for zeta(3)). (End)
G.f. for column k is hypergeom([1,k+1],[],x/(x+1))/(x+1). - Mark van Hoeij, Nov 07 2011
T(n, k) = (n!/k!)*hypergeom([k-n], [-n], -1). - Peter Luschny, Oct 05 2017

More terms from David Wasserman, Mar 28 2005
Additional comments from Zerinvary Lajos, Mar 30 2006
Edited by N. J. A. Sloane, Sep 24 2011

Previous Showing 81-90 of 544 results. Next

cp's OEIS Frontend

A001910 a(n) = n*a(n-1) + (n-5)*a(n-2).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A008278 Reflected triangle of Stirling numbers of 2nd kind, S(n,n-k+1), n >= 1, 1 <= k <= n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Formula

Extensions

A008306 Triangle T(n,k) read by rows: associated Stirling numbers of first kind (n >= 2, 1 <= k <= floor(n/2)).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

Extensions

A068106 Euler's difference table: triangle read by rows, formed by starting with factorial numbers (A000142) and repeatedly taking differences. T(n,n) = n!, T(n,k) = T(n,k+1) - T(n-1,k).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Formula

Extensions

A188674 Stack polyominoes with square core.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

A055209 a(n) = Product_{i=0..n} i!^2.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Magma

Maple

A001910 a(n) = na(n-1) + (n-5)a(n-2).