A000165 -id:A000165 - cp's OEIS frontend

A104035 Triangle T(n,k), 0 <= k <= n, read by rows, defined by T(0,0) = 1; T(0,k) = 0 if k>0 or if k<0; T(n,k) = kT(n-1,k-1) + (k+1)T(n-1,k+1).

1, 0, 1, 1, 0, 2, 0, 5, 0, 6, 5, 0, 28, 0, 24, 0, 61, 0, 180, 0, 120, 61, 0, 662, 0, 1320, 0, 720, 0, 1385, 0, 7266, 0, 10920, 0, 5040, 1385, 0, 24568, 0, 83664, 0, 100800, 0, 40320, 0, 50521, 0, 408360, 0, 1023120, 0, 1028160, 0, 362880, 50521, 0, 1326122, 0, 6749040

Offset: 0

Philippe Deléham, Apr 06 2005

Or, triangle of coefficients (with exponents in increasing order) in polynomials Q_n(u) defined by d^n sec x / dx^n = Q_n(tan x)*sec x.
Interpolates between factorials and Euler (or secant) numbers. Related to Springer numbers.
Companion triangles are A155100 (derivative polynomials of tangent function) and A185896 (derivative polynomials of squared secant function).
A combinatorial interpretation for the polynomial Q_n(u) as the generating function for a sign change statistic on certain types of signed permutation can be found in [Verges]. A signed permutation is a sequence (x_1,x_2,...,x_n) of integers such that {|x_1|,|x_2|,...,|x_n|} = {1,2,...,n}. They form a group, the hyperoctahedral group of order 2^n*n! = A000165(n), isomorphic to the group of symmetries of the n dimensional cube.
Let x_1,...,x_n be a signed permutation. Adjoin x_0 = 0 to the front of the permutation and x_(n+1) = (-1)^n*(n+1) to the end to form x_0,x_1,...,x_n,x_(n+1). Then x_0,x_1,...,x_n,x_(n+1) is a snake of type S(n;0) when x_0 < x_1 > x_2 < ... x_(n+1). For example, 0 3 -1  2 -4 is a snake of type S(3;0).
Let sc be the number of sign changes through a snake ... sc = #{i, 0 <= i <= n, x_i*x_(i+1) < 0}. For example, the snake 0 3 -1 2 -4 has sc = 3. The polynomial Q_n(u) is the generating function for the sign change statistic on snakes of type S(n;0): ... Q_n(u) = sum {snakes in S(n;0)} u^sc. See the example section below for the cases n = 2 and n = 3.
PRODUCTION MATRIX
Let D = subdiag(1,2,3,...) be the array with the indicated sequence on the first subdiagonal and zeros elsewhere and let C = transpose(D). The production matrix for this triangle is C+D: the first row of (C+D)^n is the n-th row of this triangle. D represents the derivative operator d/dx and C represents the operator p(x) -> x*d/dx(x*p(x)) acting on the basis monomials {x^n}n>=0. See Formula (1) below.

			The polynomials Q_0(u) through Q_6(u) (with exponents in decreasing order) are:
  1
  u
  2*u^2 + 1
  6*u^3 + 5*u
  24*u^4 + 28*u^2 + 5
  120*u^5 + 180*u^3 + 61*u
  720*u^6 + 1320*u^4 + 662*u^2 + 61
Triangle begins:
  1
  0 1
  1 0 2
  0 5 0 6
  5 0 28 0 24
  0 61 0 180 0 120
  61 0 662 0 1320 0 720
  0 1385 0 7266 0 10920 0 5040
  1385 0 24568 0 83664 0 100800 0 40320
  0 50521 0 408360 0 1023120 0 1028160 0 362880
  50521 0 1326122 0 6749040 0 13335840 0 11491200 0 3628800
  0 2702765 0 30974526 0 113760240 0 185280480 0 139708800 0 39916800
  2702765 0 98329108 0 692699304 0 1979524800 0 2739623040 0 1836172800 0 479001600
Examples of sign change statistic sc on snakes of type S(n;0)
= = = = = = = = = = = = = = = = = = = = = =
.....Snakes....# sign changes sc.......u^sc
= = = = = = = = = = = = = = = = = = = = = =
n=2
...0 1 -2 3...........2.................u^2
...0 2  1 3...........0.................1
...0 2 -1 3...........2.................u^2
yields Q_2(u) = 2*u^2 + 1.
n=3
...0 1 -2  3 -4.......3.................u^3
...0 1 -3  2 -4.......3.................u^3
...0 1 -3 -2 -4.......1.................u
...0 2  1  3 -4.......1.................u
...0 2 -1  3 -4.......3.................u^3
...0 2 -3  1 -4.......3.................u^3
...0 2 -3 -2 -4.......1.................u
...0 3  1  2 -4.......1.................u
...0 3 -1  2 -4.......3.................u^3
...0 3 -2  1 -4.......3.................u^3
...0 3 -2 -1 -4.......1.................u
yields Q_3(u) = 6*u^3 + 5*u.

R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics, Addison-Wesley, Reading, MA, 2nd ed. 1998, p. 287.
S. Mukai, An Introduction to Invariants and Moduli, Cambridge, 2003; see pp. 445 and 469.

Reinhard Zumkeller, Rows n = 0..125 of triangle, flattened
K. Boyadzhiev, Derivative Polynomials for tanh, tan, sech and sec in Explicit Form, arXiv:0903.0117 [math.CA], 2009-2010.
M.-P. Grosset and A. P. Veselov, Bernoulli numbers and solitons, arXiv:math/0503175 [math.GM], 2005.
Gordon Haigh, A "natural" approach to Pick's theorem, Math. Gaz. 64 (1980), no. 429, 173-180.
Michael E. Hoffman, Derivative polynomials for tangent and secant, Amer. Math. Monthly, 102 (1995), 23-30.
Michael E. Hoffman, Derivative Polynomials, Euler Polynomials, and Associated Integer Sequences, The Electronic Journal of Combinatorics, Volume 6.1 (1999): Research paper R21, 13 p.
M. Josuat-Vergès, Enumeration of snakes and cycle-alternating permutations, arXiv:1011.0929 [math.CO], 2010.
Donald E. Knuth and Thomas J. Buckholtz, Computation of tangent, Euler and Bernoulli numbers, Math. Comp. 21 1967 663-688.

See A008294 for another version of this triangle.
Setting u=0,1,2,3,4 gives A000364, A001586, A156129, A156131, A156132.
Setting u=sqrt(2) gives A156134 and A156138; u=sqrt(3) gives A002437 and A002439.
Cf. A060187, A155100, A185896, A186492.

a104035 n k = a104035_tabl !! n !! k
a104035_row n = a104035_tabl !! n
a104035_tabl = iterate f [1] where
   f xs = zipWith (+)
     (zipWith (*) [1..] (tail xs) ++ [0,0]) ([0] ++ zipWith (*) [1..] xs)
-- Reinhard Zumkeller, Apr 27 2014

nmax = 10; t[n_, k_] := t[n, k] = k*t[n-1, k-1] + (k+1)*t[n-1, k+1]; t[0, 0] = 1; t[0, ] = 0; Flatten[ Table[t[n, k], {n, 0, nmax}, {k, 0, n}]] (* _Jean-François Alcover, Nov 14 2011 *)

T(n, n) = n!; T(n, 0) = 0 if n = 2m + 1; T(n, 0) = A000364(m) if n = 2m.
Sum_{k>=0} T(m, k)*T(n, k) = T(m+n, 0).
Sum_{k>=0} T(n, k) = A001586(n): Springer numbers.
G.f.: Sum_{n >= 0} Q_n(u)*t^n/n! = 1/(cos t - u sin t).
From Peter Bala: (Start)
RECURRENCE RELATION
For n>=0,
(1)... Q_(n+1)(u) = d/du Q_n(u) + u*d/du(u*Q_n(u))
... = (1+u^2)*d/du Q_n(u) + u*Q_n(u),
with starting condition Q_0(u) = 1. Compare with Formula (4) of A186492.
RELATION WITH TYPE B EULERIAN NUMBERS
(2)... Q_n(u) = ((u+i)/2)^n*B(n,(u-i)/(u+i)), where i = sqrt(-1) and
[B(n,u)]n>=0 = [1,1+u,1+6*u+u^2,1+23*u+23*u^2+u^3,...] is the sequence of type B Eulerian polynomials (with a factor of u removed) - see A060187.
(End)
T(n,0) = abs(A122045(n)). - Reinhard Zumkeller, Apr 27 2014

Entry revised by N. J. A. Sloane, Nov 06 2009

A132382 Lower triangular array T(n,k) generator for group of arrays related to A001147 and A102625.

Original entry on oeis.org

1, -1, 1, -1, -2, 1, -3, -3, -3, 1, -15, -12, -6, -4, 1, -105, -75, -30, -10, -5, 1, -945, -630, -225, -60, -15, -6, 1, -10395, -6615, -2205, -525, -105, -21, -7, 1, -135135, -83160, -26460, -5880, -1050, -168, -28, -8, 1, -2027025, -1216215, -374220, -79380, -13230, -1890, -252, -36, -9, 1

Offset: 0

Tom Copeland, Nov 11 2007, Nov 12 2007, Nov 19 2007, Dec 04 2007, Dec 06 2007

Let b(n) = LPT[ A001147 ] = -A001147(n-1) for n > 0 and 1 for n=0, where LPT represents the action of the list partition transform described in A133314.
Then T(n,k) = binomial(n,k) * b(n-k) .
Form the matrix of polynomials TB(n,k,t) = T(n,k) * t^(n-k) = binomial(n,k) * b(n-k) * t^(n-k) = binomial(n,k) * Pb(n-k,t),
beginning as
     1;
    -1,       1;
    -1*t,    -2,       1;
    -3*t^2,  -3*t,    -3,       1;
   -15*t^3, -12*t^2,  -6*t,    -4,    1;
  -105*t^4, -75*t^3, -30*t^2, -10*t, -5, 1;
Let Pc(n,t) = LPT(Pb(.,t)).
Then [TB(t)]^(-1) = TC(t) = [ binomial(n,k) * Pc(n-k,t) ] = LPT(TB),
whose first column is
Pc(0,t) = 1
Pc(1,t) = 1
Pc(2,t) = 2 + t
Pc(3,t) = 6 + 6*t + 3*t^2
Pc(4,t) = 24 + 36*t + 30*t^2 + 15*t^3
Pc(5,t) = 120 + 240*t + 270*t^2 + 210*t^3 + 105*t^4.
The coefficients of these polynomials are given by the reverse of A102625 with the highest order coefficients given by A001147 with an additional leading 1.
Note this is not the complete matrix TC. The complete matrix is formed by multiplying along the diagonal of the lower triangular Pascal matrix by these polynomials, embedding trees of coefficients in the matrix.
exp[Pb(.,t)*x] = 1 + [(1-2t*x)^(1/2) - 1] / (t-0) = [1 + a finite diff. of [(1-2t*x)^(1/2)] with step t] = e.g.f. of the first column of TB.
exp[Pc(.,t)*x] = 1 / { 1 + [(1-2t*x)^(1/2) - 1] / t } = 1 / exp[Pb(.,t)*x) = e.g.f. of the first column of TC.
TB(t) and TC(t), being inverse to each other, are the generators of an Abelian group.
TB(0) and TC(0) are generators for a subgroup representing the iterated Laguerre operator described in A132013 and A132014.
Let sb(t,m) and sc(t,m) be the associated sequences under the LPT to TB(t)^m = B(t,m) and TC(t)^m = C(t,m).
Let Esb(t,m) and Esc(t,m) be e.g.f.'s for sb(t,m) and sc(t,m), rB(t,m) and rC(t,m) be the row sums of B(t,m) and C(t,m) and aB(t,m) and aC(t,m) be the alternating row sums.
Then B(t,m) is the inverse of C(t,m), Esb(t,m) is the reciprocal of Esc(t,m) and sb(t,m) and sc(t,m) form a reciprocal pair under the LPT. Similar relations hold among the row sums and the alternating sign row sums and associated quantities.
All the group members have the form B(t,m) * C(u,p) = TB(t)^m * TC(u)^p = [ binomial(n,k) * s(n-k) ]
with associated e.g.f. Es(x) = exp[m * Pb(.,t) * x] * exp[p * Pc(.,u) * x] for the first column of the matrix, with terms s(n), so group multiplication is isomorphic to matrix multiplication and to multiplication of the e.g.f.'s for the associated sequences (see examples).
These results can be extended to other groups of integer-valued arrays by replacing the 2 by any natural number in the expression for exp[Pb(.,t)*x].
More generally,
[ G.f. for M = Product_{i=0..j} B[s(i),m(i)] * C[t(i),n(i)] ]
= exp(u*x) * Product_{i=0..j} { exp[m(i) * Pb(.,s(i)) * x] * exp[n(i) * Pc(.,t(i)) * x] }
= exp(u*x) * Product_{i=0..j} { 1 + [ (1 - 2*s(i)*x)^(1/2) - 1 ] / s(i) }^m(i) / { 1 + [ (1 - 2*t(i)*x)^(1/2) - 1 ] / t(i) }^n(i)
= exp(u*x) * H(x)
[ E.g.f. for M ] = I_o[2*(u*x)^(1/2)] * H(x).
M is an integer-valued matrix for m(i) and n(i) positive integers and s(i) and t(i) integers. To invert M, change B to C in Product for M.
H(x) is the e.g.f. for the first column of M and diagonally multiplying the Pascal matrix by the terms of this column generates M. See examples.
The G.f. for M, i.e., the e.g.f. for the row polynomials of M, implies that the row polynomials form an Appell sequence (see Wikipedia and Mathworld). - Tom Copeland, Dec 03 2013

			Some group members and associated arrays are
(t,m) :: Array :: Asc. Matrix :: Asc. Sequence :: E.g.f. for sequence
..............................................................................
(0,1).::.B..::..A132013.::.(1,-1,0,0,0,0,...).....::.s(x).=.1-x
(0,1).::.C..::..A094587.::.(0!,1!,2!,3!,...)......::.1./.s(x)
(0,1).::.rB.::.~A055137.::.(1,0,-1,-2,-3,-4,...)..::.exp(x).*.s(x)
(0,1).::.rC.::....-.....::..A000522...............::.exp(x)./.s(x)
(0,1).::.aB.::....-.....::.(1,-2,3,-4,5,-6,...)...::.exp(-x).*.s(x)
(0,1).::.aC.::..A008290.::..A000166...............::.exp(-x)./.s(x)
..............................................................................
(0,2).::.B..::..A132014.::.(1,-2,2,0,0,0,0...)....::.s(x).=.(1-x)^2
(0,2).::.C..::..A132159.::.(1!,2!,3!,4!,...)......::..1./.s(x).
(0,2).::.rB.::...-......::.(1,-1,-1,1,5,11,19,29,)::.exp(x).*.s(x).
(0,2).::.rC.::...-......::..A001339...............::.exp(x)./.s(x).
(0,2).::.aB.::...-......::.(-1)^n.A002061(n+1)....::.exp(-x).*.s(x).
(0,2).::.aC.::...-......::..A000255...............::.exp(-x)./.s(x).
..............................................................................
(1,1).::.B..::..T.......::.(1,-A001147(n-1))......::.s(x).=.(1-2x)^(1/2)
(1,1).::.C..::.~A113278.::..A001147...............::.1./.s(x)...
(1,1).::.rB.::...-......::..A055142...............::.exp(x).*.s(x).
(1,1).::.rC.::...-......::..A084262...............::.exp(x)./.s(x).
(1,1).::.aB.::...-......::.(1,-2,2,-4,-4,-56,...).::.exp(-x).*.s(x).
(1,1).::.aC.::...-......::..A053871...............::.exp(-x)./.s(x).
..............................................................................
(2,1).::.B..::...-......::.(1,-A001813)...........::.s=[1+(1-4x)^(1/2)]/2....
(2,1).::.C..::...-......::..A001761...............::.1./.s(x)..
(2,1).::.rB.::...-......::.(1,0,-3,-20,-183,...)..::.exp(x).*.s(x)..
(2,1).::.rC.::...-......::.(1,2,7,46,485,...).....::.exp(x)./.s(x).
(2,1).::.aB.::...-......::.(1,-2,1,-10,-79,...)...::.exp(-x).*.s(x).
(2,1).::.aC.::...-......::.(1,0,3,20,237,...).....::.exp(-x)./.s(x)
..............................................................................
(1,2).::.B..::.~A134082.::.(1,-2,0,0,0,0,...).....::.s(x).=.1.-.2x
(1,2).::.C..::....-.....::..A000165...............::.1./.s(x)..
(1,2).::.rB.::....-.....::.(1,-1,-3,-5,-7,-9,...).::.exp(x).*.s(x).
(1,2).::.rC.::....-.....::..A010844...............::.exp(x)./.s(x)..
(1,2).::.aB.::....-.....::.(1,-3,5,-7,9,-11,...)..::.exp(-x).*.s(x).
(1,2).::.aC.::....-.....::..A000354...............::.exp(-x)./.s(x).
..............................................................................
(The tilde indicates the match is not exact--specifically, there are differences in signs from the true matrices.)
Note the row sums correspond to binomial transforms of s(x) and the alternating row sums, to inverse binomial transforms, or, finite differences.
Some additional examples:
C(1,2)*B(0,1) = B(1,-2)*C(0,-1) = [ binomial(n,k)*A002866(n-k) ] with asc. e.g.f. (1-x) / (1-2x).
B(1,2)*C(0,1) = C(1,-2)*B(0,-1) = 2I - A094587 with asc. e.g.f. (1-2x) / (1-x).

[G.f. for TB(n,k,t)] = GTB(u,x,t) = exp(u*x) * { 1 + [ (1 - 2t*x)^(1/2) - 1 ] / t } = exp[(u+Pb(.,t))*x] where TB(n,k,t) = (D_x)^n (D_u)^k /k! GTB(u,x,t) eval. at u=x=0.
[G.f. for TC(n,k,t)] = GTC(u,x,t) = exp(u*x) / { 1 + [ (1 - 2t*x)^(1/2) - 1 ] / t } = exp[(u+Pc(.,t))*x] where TC(n,k,t) = (D_x)^n (D_u)^k /k! GTC(u,x,t) eval. at u=x=0.
[E.g.f. for TB(n,k,t)] = I_o[2*(u*x)^(1/2)] * { 1 + [ (1 - 2t*x)^(1/2) - 1 ] / t } and
[E.g.f. for TC(n,k,t)] = I_o[2*(u*x)^(1/2)] / { 1 + [ (1 - 2t*x)^(1/2) - 1 ] / t }
where I_o is the zeroth modified Bessel function of the first kind, i.e.,
I_o[2*(u*x)^(1/2)] = Sum_{j>=0} (u^j/j!) * (x^j/j!).
So [e.g.f. for TB(n,k)] = I_o[2*(u*x)^(1/2)] * (1 - 2x)^(1/2).

More terms from Tom Copeland, Dec 05 2007

A051188 Sept-factorial numbers.

Original entry on oeis.org

1, 7, 98, 2058, 57624, 2016840, 84707280, 4150656720, 232436776320, 14643516908160, 1025046183571200, 78928556134982400, 6629998715338521600, 603329883095805465600, 59126328543388935628800

Offset: 0

Wolfdieter Lang

For n >= 1, a(n) is the order of the wreath product of the symmetric group S_n and the Abelian group (C_7)^n. - Ahmed Fares (ahmedfares(AT)my-deja.com), May 07 2001

Vincenzo Librandi, Table of n, a(n) for n = 0..300
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 614.
Index entries for sequences related to factorial numbers

Cf. A000142, A000165, A032031, A045754, A047053, A047058, A049209, A051186, A052562, A053106, A084947, A092516, A092750, A144827, A144739, A147585.

[7^n*Factorial(n): n in [0..20]]; // Vincenzo Librandi, Oct 05 2011

a(n) = n!*7^n; \\ Michel Marcus, Jun 08 2020

a(n) = n!*7^n =: (7*n)(!^7).
a(n) = 7*A034834(n) = Product_{k=1..n} 7*k, n >= 1.
E.g.f.: 1/(1 - 7*x).
G.f.: 1/(1 - 7*x/(1 - 7*x/(1 - 14*x/(1 - 14*x/(1 - 21*x/(1 - 21*x/(1 - 28*x/(1 - 28*x/(1 - ... (continued fraction). - Philippe Deléham, Jan 08 2012
From Amiram Eldar, Jun 25 2020: (Start)
Sum_{n>=0} 1/a(n) = e^(1/7) (A092516).
Sum_{n>=0} (-1)^n/a(n) = e^(-1/7) (A092750). (End)

A155100 Triangle read by rows: coefficients in polynomials P_n(u) arising from the expansion of D^(n-1) (tan x) in increasing powers of tan x for n>=1 and 1 for n=0.

Original entry on oeis.org

1, 0, 1, 1, 0, 1, 0, 2, 0, 2, 2, 0, 8, 0, 6, 0, 16, 0, 40, 0, 24, 16, 0, 136, 0, 240, 0, 120, 0, 272, 0, 1232, 0, 1680, 0, 720, 272, 0, 3968, 0, 12096, 0, 13440, 0, 5040, 0, 7936, 0, 56320, 0, 129024, 0, 120960, 0, 40320, 7936, 0, 176896, 0, 814080, 0, 1491840

Offset: 0

N. J. A. Sloane, Nov 05 2009

The definition is d^(n-1) tan x / dx^n = P_n(tan x) for n>=1 and 1 for n=0.
Interpolates between factorials and tangent numbers.
From Peter Bala, Mar 02 2011: (Start)
Companion triangles are A104035 and A185896.
A combinatorial interpretation for the polynomial P_n(t) as the generating function for a sign change statistic on certain types of signed permutation can be found in [Verges].
A signed permutation is a sequence (x_1,x_2,...,x_n) of integers such that {|x_1|,|x_2|,...,|x_n|} = {1,2,...,n}.
They form a group, the hyperoctahedral group of order 2^n*n! = A000165(n), isomorphic to the group of symmetries of the n dimensional cube.
Let x_1,...,x_n be a signed permutation and put x_0 = -(n+1) and x_(n+1) = (-1)^n*(n+1). Then x_0,x_1,...,x_n,x_(n+1) is a snake of type S(n) when x_0 < x_1 > x_2 < ... x_(n+1). For example, -5 4 -3 -1 -2 5 is a snake of type S(4).
Let sc be the number of sign changes through a snake sc = #{i, 0 <= i <= n, x_i*x_(i+1) < 0}. For example, the snake -5 4 -3 -1 -2 5 has sc = 3.
The polynomial P_(n+1)(t) is the generating function for the sign change statistic on snakes of type S(n): P_(n+1)(t) = sum {snakes in S(n)} t^sc.
See the example section below for the cases n=1 and n=2.
(End)
Equals A107729 when the first column is removed. - Georg Fischer, Jul 26 2023

			The polynomials P_{-1}(u) through P_6(u) with exponents in decreasing order:
      1
      u
      u^2 +    1
    2*u^3 +    2*u
    6*u^4 +    8*u^2 +    2
   24*u^5 +   40*u^3 +   16*u
  120*u^6 +  240*u^4 +  136*u^2 +  16
  720*u^7 + 1680*u^5 + 1232*u^3 + 272*u
  ...
Triangle begins:
  1
  0, 1
  1, 0, 1
  0, 2, 0, 2
  2, 0, 8, 0, 6
  0, 16, 0, 40, 0, 24
  16, 0, 136, 0, 240, 0, 120
  0, 272, 0, 1232, 0, 1680, 0, 720
  272, 0, 3968, 0, 12096, 0, 13440, 0, 5040
  0, 7936, 0, 56320, 0, 129024, 0, 120960, 0, 40320
  7936, 0, 176896, 0, 814080, 0, 1491840, 0, 1209600, 0, 362880
  0, 353792, 0, 3610112, 0, 12207360, 0, 18627840, 0, 13305600, 0, 3628800
  ...
From _Peter Bala_, Feb 07 2011: (Start)
Examples of sign change statistic sc on snakes of type S(n):
    Snakes     # sign changes sc  t^sc
  ===========  =================  ====
n=1:
  -2  1 -2 ........... 2 ........ t^2
  -2 -1 -2 ........... 0 ........ 1
                  yields P_2(t) = 1 + t^2;
n=2:
  -3  1 -2  3 ........ 3 ........ t^3
  -3  2  1  3 ........ 1 ........ t
  -3  2 -1  3 ........ 3 ........ t^3
  -3 -1 -2  3 ........ 1 ........ t
                  yields P_3(t) = 2*t + 2*t^3. (End)

R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics, Addison-Wesley, Reading, MA, 2nd ed. 1998, p. 287.

G. C. Greubel, Rows n = 0..101 of triangle, flattened
K. Boyadzhiev, Derivative Polynomials for tanh, tan, sech and sec in Explicit Form, arXiv:0903.0117 [math.CA], 2009-2010.
M.-P. Grosset and A. P. Veselov, Bernoulli numbers and solitons, arXiv:math/0503175 [math.GM], 2005.
Gordon Haigh, A "natural" approach to Pick's theorem, Math. Gaz. 64 (1980), no. 429, 173-180.
Michael E. Hoffman, Derivative polynomials for tangent and secant, Amer. Math. Monthly, 102 (1995), 23-30.
Michael E. Hoffman, Derivative Polynomials, Euler Polynomials, and Associated Integer Sequences, The Electronic Journal of Combinatorics, Volume 6.1 (1999): Research paper R21, 13 p.
Donald E. Knuth and Thomas J. Buckholtz, Computation of tangent, Euler and Bernoulli numbers, Math. Comp. 21 1967 663-688.
Shi-Mei Ma, Qi Fang, Toufik Mansour, and Yeong-Nan Yeh, Alternating Eulerian polynomials and left peak polynomials, arXiv:2104.09374 [math.CO], 2021.

For other versions of this triangle see A008293, A101343.
A104035 is a companion triangle.
Highest order coefficients give factorials A000142. Constant terms give tangent numbers A000182. Other coefficients: A002301.
Setting u=1 in P_n gives A000831, u=2 gives A156073, u=3 gives A156075, u=4 gives A156076, u=1/2 gives A156102.
Setting u=sqrt(2) in P_n gives A156108 and A156122; setting u=sqrt(3) gives A156103 and A000436.
Cf. A008292, A019538, A107729, A185896.

P:=proc(n) option remember;
if n=-1 then RETURN(1); elif n=0 then RETURN(u); else RETURN(expand((u^2+1)*diff(P(n-1),u))); fi;
end;
for n from -1 to 12 do t1:=series(P(n),u,20); lprint(seriestolist(t1)); od:
# Alternatively:
with(PolynomialTools): seq(print(CoefficientList(`if`(i=0,1,D@@(i-1))(tan),tan)), i=0..7); # Peter Luschny, May 19 2015

p[n_, u_] := D[Tan[x], {x, n}] /. Tan[x] -> u /. Sec[x] -> Sqrt[1 + u^2] // Expand; p[-1, u_] = 1; Flatten[ Table[ CoefficientList[ p[n, u], u], {n, -1, 9}]] (* Jean-François Alcover, Jun 28 2012 *)
T[ n_, k_] := Which[n<0, Boole[n==-1 && k==0], n==0, Boole[k==1], True, (k-1)*T[n-1, k-1] + (k+1)*T[n-1, k+1]]; (* Michael Somos, Jul 09 2024 *)

{T(n, k) = if(n<0, n==-1 && k==0, n==0, k==1, (k-1)*T(n-1, k-1) + (k+1)*T(n-1, k+1))}; /* Michael Somos, Jul 09 2024 */

If the polynomials are denoted by P_n(u), we have the recurrence P_{-1}=1, P_0 = u, P_n = (u^2+1)*dP_{n-1}/du.
G.f.: Sum_{n >= 0} P_n(u) t^n/n! = (sin t + u*cos t)/(cos t - u sin t). [Hoffman]
From Peter Bala, Feb 07 2011: (Start)
RELATION WITH BERNOULLI NUMBERS A000367 AND A002445
Put T(n,t) = P_n(i*t), where i = sqrt(-1). We have the definite integral evaluation, valid when both m and n are >=1 and m+n >= 4:
int( T(m,t)*T(n,t)/(1-t^2), t = -1..1) = (-1)^((m-n)/2)*2^(m+n-1)*Bernoulli(m+n-2).
The case m = n is equivalent to the result of [Grosset and Veselov]. The methods used there extend to the general case.
RELATION WITH OTHER ROW POLYNOMIALS
The following three identities hold for n >= 1:
P_(n+1)(t) = (1+t^2)*R(n-1,t) where R(n,t) is the n-th row polynomial of A185896.
P_(n+1)(t) = (-2*i)^n*(t-i)*R(n,-1/2+1/2*i*t), where i = sqrt(-1) and R(n,x) is an ordered Bell polynomial, that is, the n-th row polynomial of A019538.
P_(n+1)(t) = (t-i)*(t+i)^n*A(n,(t-i)/(t+i)), where {A(n,t)}n>=1 = [1,1+t,1+4*t+t^2,1+11*t+11*t^2+t^3,...] is the sequence of Eulerian polynomials - see A008292. (End)
T(n,k) = cos((n+k)*Pi/2) * Sum_{p=0..n-1} A008292(n-1,p+1) Sum_{j=0..k}(-1)^(p+j+1) * binomial(p+1,k-j) *binomial(n-p-1,j) for n>1. - Ammar Khatab, Aug 15 2024

Name clarified by Peter Luschny, May 25 2015

A037223 Number of solutions to non-attacking rooks problem on n X n board that are invariant under 180-degree rotation.

Original entry on oeis.org

1, 1, 2, 2, 8, 8, 48, 48, 384, 384, 3840, 3840, 46080, 46080, 645120, 645120, 10321920, 10321920, 185794560, 185794560, 3715891200, 3715891200, 81749606400, 81749606400, 1961990553600, 1961990553600, 51011754393600, 51011754393600, 1428329123020800, 1428329123020800

Offset: 0

Miklos SZABO (mike(AT)ludens.elte.hu)

This is just A000165 doubled up. Normally such sequences do not get their own entry in the OEIS. This is an exception. - N. J. A. Sloane, Sep 23 2006
Also the number of permutations of (1,2,3,...,n) for which the reverse of the inverse is the same as the inverse of the reverse. - Ian Duff, Mar 09 2007
Conjecture: a(n) = Product_{1<=i<=n and phi(i)<=floor(i/2)}i. - Enrique Pérez Herrero, May 31 2012. This conjecture is WRONG, counterexample is n=105. [Vaclav Kotesovec, Sep 07 2012]

E. Lucas, Theorie des nombres, Gauthiers-Villars, Paris, 1891, Vol 1, p. 221.

Alois P. Heinz, Table of n, a(n) for n = 0..500
E. Lucas, Théorie des nombres, Gauthiers-Villars, Paris, 1891, Vol 1, p. 221.
R. W. Robinson, Counting arrangements of bishops, pp. 198-214 of Combinatorial Mathematics IV (Adelaide 1975), Lect. Notes Math., 560 (1976).
M. Szabo, Non-attacking Queens Problem Page

Cf. A000165, A033148, A037224, A032522, A037223.

[Factorial((n div 2) -1)*2^((n div 2)-1): n   in [2..35]]; // Vincenzo Librandi, Nov 17 2018

For Maple program see A000903.
# second Maple program:
a:= n-> (r-> r!*2^r)(iquo(n, 2)):
seq(a(n), n=0..30);  # Alois P. Heinz, Dec 23 2013

f[n_]:=Times@@Select[Range[n],EulerPhi[#]<=Floor[#/2]&]; Table[f[n],{n,1,30}] (* Conjectured: Enrique Pérez Herrero, May 31 2012 *)(* This conjecture and also program is WRONG for n=105, Vaclav Kotesovec, Sep 07 2012 *)
a[n_] := (2*Floor[n/2])!!; Table[a[n], {n, 0, 27}] (* Jean-François Alcover, Dec 23 2013, after N. J. A. Sloane's comment *)

a(2n) = a(2n+1) = n!*2^n.
E.g.f.: 1 + x + (1 + x + x^2)*exp(x^2/2)*sqrt(Pi/2)*erf(x/sqrt(2)), where erf denotes the error function. - Antonio G. Astudillo (afg_astudillo(AT)hotmail.com), Nov 01 2002
For asymptotics see the Robinson paper.
E.g.f.: Q(0)  where Q(k)= 1 + x/(2*k + 1 - x*(2*k+1)/(x+1/Q(k+1))); (continued fraction, 3-step). - Sergei N. Gladkovskii, Sep 21 2012
E.g.f.: 1/(W(0)-x)  where W(k)= x + 1/(1 + x/(2*k + 1 - x*(2*k+1)/W(k+1))); (continued fraction, 3-step). - Sergei N. Gladkovskii, Sep 22 2012
a(n) = Product_{i=1..floor(n/2)} 2*i. - Wesley Ivan Hurt, Oct 19 2014
D-finite with recurrence: a(n) +a(n-1) -n*a(n-2) +(-n+2)*a(n-3)=0. - R. J. Mathar, Feb 20 2020

More terms from Antonio G. Astudillo (afg_astudillo(AT)hotmail.com), Nov 01 2002

Edited by N. J. A. Sloane, Sep 23 2006

A075497 Stirling2 triangle with scaled diagonals (powers of 2).

Original entry on oeis.org

1, 2, 1, 4, 6, 1, 8, 28, 12, 1, 16, 120, 100, 20, 1, 32, 496, 720, 260, 30, 1, 64, 2016, 4816, 2800, 560, 42, 1, 128, 8128, 30912, 27216, 8400, 1064, 56, 1, 256, 32640, 193600, 248640, 111216, 21168, 1848, 72, 1

Offset: 1

Wolfdieter Lang, Oct 02 2002

This is a lower triangular infinite matrix of the Jabotinsky type. See the D. E. Knuth reference given in A039692 for exponential convolution arrays.
The row polynomials p(n,x) := Sum_{m=1..n} a(n,m)x^m, n >= 1, have e.g.f. J(x; z)= exp((exp(2*z) - 1)*x/2) - 1.
Subtriangle of (0, 2, 0, 4, 0, 6, 0, 8, 0, 10, 0, 12, ...) DELTA (1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Feb 13 2013
Also the inverse Bell transform of the double factorial of even numbers Product_ {k=0..n-1} (2*k+2) (A000165). For the definition of the Bell transform see A264428 and for cross-references A265604. - Peter Luschny, Dec 31 2015
This is the exponential Riordan array [exp(2*x), (exp(2*x) - 1)/2] belonging to the derivative subgroup of the exponential Riordan group. In the notation of Corcino, this is the triangle of (2, 2)-Stirling numbers of the second kind. A factorization of the array as an infinite product is given in the example section. - Peter Bala, Feb 20 2025

			Triangle begins:
  [1];
  [2,1];
  [4,6,1]; p(3,x) = x*(4 + 6*x + x^2).
  ...;
Triangle (0, 2, 0, 4, 0, 6, 0, 8, ...) DELTA (1, 0, 1, 0, 1, 0, 1, 0, ...) begins:
  1
  0,  1
  0,  2,   1
  0,  4,   6,   1
  0,  8,  28,  12,  1
  0, 16, 120, 100, 20, 1. - _Philippe Deléham_, Feb 13 2013
From _Peter Bala_, Feb 23 2025: (Start)
The array factorizes as
/ 1               \       /1             \ /1             \ /1            \
| 2    1           |     | 2   1          ||0  1           ||0  1          |
| 4    6   1       |  =  | 4   4   1      ||0  2   1       ||0  0  1       | ...
| 8   28  12   1   |     | 8  12   6  1   ||0  4   4  1    ||0  0  2  1    |
|16  120 100  20  1|     |16  32  24  8  1||0  8  12  6  1 ||0  0  4  4  1 |
|...               |     |...             ||...            ||...           |
where, in the infinite product on the right-hand side, the first array is the Riordan array (1/(1 - 2*x), x/(1 - 2*x)) = P^2, where P denotes Pascal's triangle. See A038207. Cf. A143494. (End)

Alois P. Heinz, Rows n = 1..141, flattened
Peter Bala, The white diamond product of power series
Peter Bala, Factorising (r,b)-Stirling arrays
Paul Barry, Three Études on a sequence transformation pipeline, arXiv:1803.06408 [math.CO], 2018.
John R. Britnell and Mark Wildon, Bell numbers, partition moves and the eigenvalues of the random-to-top shuffle in types A, B and D, arXiv 1507.04803 [math.CO], 2015.
Roberto B. Corcino, The (r, β)-Stirling Numbers, The Mindanao Forum, Vol. XIV, No.2, pp. 91-99, 1999.
Roberto B. Corcino and Maribeth B. Montero, The (r, β)-Stirling Numbers in the Context of 0-1 Tableau, Jour. Math. Soc. of the Philippines, ISSN 0115-6926, Vol. 32, No. 1 (2009), pp. 45-52
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.
Wolfdieter Lang, First 10 rows.
Toufik Mansour, Generalization of some identities involving the Fibonacci numbers, arXiv:math/0301157 [math.CO], 2003.
Emanuele Munarini, Characteristic, admittance and matching polynomials of an antiregular graph, Appl. Anal. Discrete Math 3 (1) (2009) 157-176.

Columns 1-7 are A000079, A006516, A016283, A025966, A075510-A075512.
Row sums are A004211.
Cf. A008277, A075498-A075505.

with(combinat):
b:= proc(n, i) option remember; expand(`if`(n=0, 1,
       `if`(i<1, 0, add(x^j*multinomial(n, n-i*j, i$j)/j!*add(
        binomial(i, 2*k), k=0..i/2)^j*b(n-i*j, i-1), j=0..n/i))))
    end:
T:= n-> (p-> seq(coeff(p, x, i), i=1..n))(b(n$2)):
seq(T(n), n=1..12);  # Alois P. Heinz, Aug 13 2015
# Alternatively, giving the triangle in the form displayed in the Example section:
gf := exp(x*exp(z)*sinh(z)):
X := n -> series(gf, z, n+2):
Z := n -> n!*expand(simplify(coeff(X(n), z, n))):
A075497_row := n -> op(PolynomialTools:-CoefficientList(Z(n), x)):
seq(A075497_row(n), n=0..9); # Peter Luschny, Jan 14 2018

Table[(2^(n - m)) StirlingS2[n, m], {n, 9}, {m, n}] // Flatten (* Michael De Vlieger, Dec 31 2015 *)

for(n=1, 11, for(m=1, n, print1(2^(n - m) * stirling(n, m, 2),", ");); print();) \\ Indranil Ghosh, Mar 25 2017

# uses[inverse_bell_transform from A265605]
multifact_2_2 = lambda n: prod(2*k + 2 for k in (0..n-1))
inverse_bell_matrix(multifact_2_2, 9) # Peter Luschny, Dec 31 2015

a(n, m) = (2^(n-m)) * Stirling2(n, m).
a(n, m) = (Sum_{p=0..m-1} A075513(m, p)*((p+1)*2)^(n-m))/(m-1)! for n >= m >= 1, else 0.
a(n, m) = 2*m*a(n-1, m) + a(n-1, m-1), n >= m >= 1, else 0, with a(n, 0) := 0 and a(1, 1)=1.
G.f. for m-th column: (x^m)/Product_{k=1..m}(1-2*k*x), m >= 1.
E.g.f. for m-th column: (((exp(2*x)-1)/2)^m)/m!, m >= 1.
The row polynomials in t are given by D^n(exp(x*t)) evaluated at x = 0, where D is the operator (1+2*x)*d/dx. Cf. A008277. - Peter Bala, Nov 25 2011
From Peter Bala, Jan 13 2018: (Start)
n-th row polynomial R(n,x)= x o x o ... o x (n factors), where o is the deformed Hadamard product of power series defined in Bala, section 3.1.
R(n+1,x)/x = (x + 2) o (x + 2) o...o (x + 2) (n factors).
R(n+1,x) = x*Sum_{k = 0..n} binomial(n,k)*2^(n-k)*R(k,x).
Dobinski-type formulas: R(n,x) = exp(-x/2)*Sum_{i >= 0} (2*i)^n* (x/2)^i/i!; 1/x*R(n+1,x) = exp(-x/2)*Sum_{i >= 0} (2 + 2*i)^n* (x/2)^i/i!. (End)

A051189 Octo-factorial numbers.

Original entry on oeis.org

1, 8, 128, 3072, 98304, 3932160, 188743680, 10569646080, 676457349120, 48704929136640, 3896394330931200, 342882701121945600, 32916739307706777600, 3423340888001504870400, 383414179456168545484800

Offset: 0

Wolfdieter Lang

For n >= 1, a(n) is the order of the wreath product of the symmetric group S_n and the Abelian group (C_8)^n. - Ahmed Fares (ahmedfares(AT)my-deja.com), May 07 2001
Number of n X n monomial matrices whose nonzero entries are unit quaternions.
Number of ways of reassembling n slices of toast or of binding n square pages. - Donald S. McDonald, Sep 24 2005

Vincenzo Librandi, Table of n, a(n) for n = 0..300
Index entries for sequences related to factorial numbers

Cf. A000165, A034976, A047058, A051188, A053115.

Shifted absolute values are column 1 of A051187.

[8^n*Factorial(n): n in [0..20]]; // Vincenzo Librandi, Oct 05 2011

Table[n! 8^n,{n,0,20}] (* Harvey P. Dale, Aug 14 2021 *)

[8^n*factorial(n) for n in range(40)] # G. C. Greubel, Oct 21 2022

a(n) = 8*A034976(n) = Product_{k=1..n} 8*k, n >= 1; a(0) = 1.
a(n) = n!*8^n.
E.g.f.: 1/(1-8*x).
G.f.: 1/(1 - 8*x/(1 - 8*x/(1 - 16*x/(1 - 16*x/(1 - 24*x/(1 - 24*x/(1 - 32*x/(1 - 32*x/(1 - ... (continued fraction). - Philippe Deléham, Jan 07 2012
From Amiram Eldar, Jun 25 2020: (Start)
Sum_{n>=0} 1/a(n) = e^(1/8).
Sum_{n>=0} (-1)^n/a(n) = e^(-1/8). (End)

A092605 Decimal expansion of e^(-1/2) or 1/sqrt(e).

Original entry on oeis.org

6, 0, 6, 5, 3, 0, 6, 5, 9, 7, 1, 2, 6, 3, 3, 4, 2, 3, 6, 0, 3, 7, 9, 9, 5, 3, 4, 9, 9, 1, 1, 8, 0, 4, 5, 3, 4, 4, 1, 9, 1, 8, 1, 3, 5, 4, 8, 7, 1, 8, 6, 9, 5, 5, 6, 8, 2, 8, 9, 2, 1, 5, 8, 7, 3, 5, 0, 5, 6, 5, 1, 9, 4, 1, 3, 7, 4, 8, 4, 2, 3, 9, 9, 8, 6, 4, 7, 6, 1, 1, 5, 0, 7, 9, 8, 9, 4, 5, 6, 0, 2, 6, 4, 2, 3

Offset: 0

Mohammad K. Azarian, Apr 22 2004

For x = e^(-1/2), the largest prime factor of a random integer n is equally likely to be above or below n^x. - Charles R Greathouse IV, May 25 2009

Siegel's conjecture: this constant gives the density of regular primes among all the primes (see Ribenboim and Siegel). - Stefano Spezia, Apr 22 2025

			0.6065306597126334...

Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See p. 225.
C. L. Siegel, Zu zwei Bemerkungen Kummers. Nachr. Akad. d. Wiss. Göttingen, Math. Phys. Kl., II, 1964, 51-62. Reprinted in Gesammelte Abhandlungen (edited by K. Chandrasekharan and H. Maas), Vol. III, 436-442. Springer-Verlag, Berlin, 1966.

Michael I. Shamos, A catalog of the real numbers, (2007). See p. 502.
Index entries for transcendental numbers.

Cf. A000165, A001113, A019774, A068985, A290506.

RealDigits[E^-(1/2),10,120][[1]] (* Harvey P. Dale, Jul 23 2012 *)

exp(-.5) \\ Charles R Greathouse IV, Oct 02 2022

Equals Sum_{k>=0} (-1)^k/(2^k * k!) = Sum_{k>=0} (-1)^k/A000165(k). - Amiram Eldar, Aug 15 2020
From Peter Bala, Jan 16 2022; (Start)
Equals 16*Sum_{n >= 0} (-1)^n*n^2/((4*n^2 - 1)*(4*n^2 - 9)*(2^n)*n!).
Equals 8*Sum_{n >= 0} (-1)^n/(p(n)*p(n+1)*(2^n)*n!), where p(n) = 4*n^2 + 8*n + 1.
Equals 48*Sum_{n >= 0} (-1)^n/(q(n)*q(n+1)*(2^n)*n!), where q(n) = 8*n^3 + 36*n^2 + 34*n + 1. (End)
Equals i^(i/Pi), where i denotes the imaginary unit. - Stefano Spezia, Mar 01 2025
Equals 1 - A290506. - Amiram Eldar, Apr 22 2025

A002454 Central factorial numbers: a(n) = 4^n * (n!)^2.

Original entry on oeis.org

1, 4, 64, 2304, 147456, 14745600, 2123366400, 416179814400, 106542032486400, 34519618525593600, 13807847410237440000, 6682998146554920960000, 3849406932415634472960000, 2602199086312968903720960000, 2040124083669367620517232640000, 1836111675302430858465509376000000

Offset: 0

N. J. A. Sloane

Denominators in the series for Bessel's J0(x) = 1 - x^2/4 + x^4/64 - x^6/2304 + ...
a(n) is the unreduced numerator in Product_{k=1..n} (4*k^2)/(4*k^2-1), therefore a(n)/A079484(n) = Pi/2 as n -> oo. - Daniel Suteu, Dec 02 2016
From Zhi-Wei Sun, Jun 26 2022: (Start)
Conjecture: Let zeta be a primitive 2n+1-th root of unity. Then the permanent of the 2n X 2n matrix [m(j,k)]_{j,k=1..2n} is a(n)/(2n+1) = ((2n)!!)^2/(2n+1), where m(j,k) is 1 or (1+zeta^(j-k))/(1-zeta^(j-k)) according as j = k or not.
The determinant of the matrix [m(j,k)]_{j,k=1..2n} was shown to be (-1)^(n-1)*((2n)!!)^2/(2n(2n+1)) by Han Wang and Zhi-Wei Sun in 2022. (End)

Richard Bellman, A Brief Introduction to Theta Functions, Dover, 2013 (20.1).
Bronstein-Semendjajew, Taschenbuch der Mathematik, 7th german ed. 1965, ch. 4.4.7
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. 110.
E. L. Ince, Ordinary Differential Equations, Dover, NY, 1956; see p. 173.
J. Riordan, Combinatorial Identities, Wiley, 1968, p. 217.
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).
Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapters 49 and 52, equations 49:6:1 and 52:6:2 at pages 483, 513.

T. D. Noe, Table of n, a(n) for n = 0..50
T. R. Van Oppolzer, Lehrbuch zur Bahnbestimmung der Kometen und Planeten, Vol. 2, Engelmann, Leipzig, 1880, p. 7.
Han Wang and Zhi-Wei Sun, Proof of a conjecture involving derangements and roots of unity, arXiv:2206.02589 [math.CO], 2022.
Index to divisibility sequences.
Index entries for sequences related to factorial numbers.

Cf. A000165, A001818, A079484, A197036, A334380.

J1: A002474, J2: A002506, J3: A014401.

[4^n*Factorial(n)^2: n in [0..15]]; // Vincenzo Librandi, Mar 15 2019

Array[4^# (#!)^2 &, 14, 0] (* Michael De Vlieger, Nov 01 2017 *)

a(n) = 4^n*(n!)^2; \\ Michel Marcus, Mar 13 2019

(-1)^n*a(n) is the coefficient of x^1 in Product_{k=0..2*n} (x+2*k-2*n). - Benoit Cloitre and Michael Somos, Nov 22 2002
E.g.f.: A(x) = arcsin(x)*sec(arcsin(x)). - Vladimir Kruchinin, Sep 12 2010
E.g.f.: arcsin(x)*sec(arcsin(x)) = arcsin(x)/sqrt(1-x^2) = x/G(0); G(k) = 2k*(x^2+1)+1-x^2*(2k+1)*(2k+2)/G(k+1); (continued fraction). - Sergei N. Gladkovskii, Nov 20 2011
G.f.: 1 + x*(G(0) - 1)/(x-1) where G(k) = 1 - (2*k+2)^2/(1-x/(x - 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jan 15 2013
From Ilya Gutkovskiy, Dec 02 2016: (Start)
a(n) ~ Pi*2^(2*n+1)*n^(2*n+1)/exp(2*n).
Sum_{n>=0} 1/a(n) = BesselI(0,1) = A197036. (End)
From Daniel Suteu, Dec 02 2016: (Start)
a(n) ~ 2^(2*n) * gamma(n+1/2) * gamma(n+3/2).
a(n) ~ Pi*(2*n+1)*(4*n^2-1)^n/exp(2*n). (End)
2*a(n)/(2*n+1)! = A101926(n) / A001803(n). - Daniel Suteu, Feb 03 2017
Limit_{n->oo} n*a(n)/((2n+1)!!)^2 = Pi/4. - Daniel Suteu, Nov 01 2017
Sum_{n>=0} (-1)^n/a(n) = BesselJ(0, 1) (A334380). - Amiram Eldar, Apr 09 2022
Limit_{n->oo} a(n) / (n * A001818(n)) = Pi. - Daniel Suteu, Apr 09 2022

A114799 Septuple factorial, 7-factorial, n!7, n!!!!!!!, a(n) = n*a(n-7) if n > 1, else 1.

Original entry on oeis.org

1, 1, 2, 3, 4, 5, 6, 7, 8, 18, 30, 44, 60, 78, 98, 120, 288, 510, 792, 1140, 1560, 2058, 2640, 6624, 12240, 19800, 29640, 42120, 57624, 76560, 198720, 379440, 633600, 978120, 1432080, 2016840, 2756160, 7352640, 14418720, 24710400, 39124800

Offset: 0

Jonathan Vos Post, Feb 18 2006

Many of the terms yield multifactorial primes a(n) + 1, e.g.: a(2) + 1 = 3, a(4) + 1 = 5, a(6) + 1 = 7, a(9) + 1 = 19, a(10) + 1 = 31, a(12) + 1 = 61, a(13) + 1 = 79, a(24) + 1 = 12241, a(25) + 1 = 19801, a(26) + 1 = 29641, a(29) + 1 = 76561, a(31) + 1 = 379441, a(35) + 1 = 2016841, a(36) + 1 = 2756161, ...

Equivalently, product of all positive integers <= n congruent to n (mod 7). - M. F. Hasler, Feb 23 2018

			a(40) = 40 * a(40-7) = 40 * a(33) = 40 * (33*a(26)) = 40 * 33 * (26*a(19)) = 40 * 33 * 26 * (19*a(12)) = 40 * 33 * 26 * 19 * (12*a(5)) = 40 * 33 * 26 * 19 * 12 5 = 39124800.

G. C. Greubel, Table of n, a(n) for n = 0..1000
Eric Weisstein's World of Mathematics, Multifactorial.
Index entries for sequences related to factorial numbers

Cf. A045754, A084947, A288094.

Cf. k-fold factorials: A000142, A001147 (and A000165, A006882), A007559 (and A032031, A008544, A007661), A007696 (and A001813, A008545, A047053, A007662), A008548 (and A052562, A047055, A085157), A085158 (and A008542, A047058, A047657), A045755.

a:= function(n)
    if n<1 then return 1;
    else return n*a(n-7);
    fi;
  end;
List([0..40], n-> a(n) ); # G. C. Greubel, Aug 20 2019

b:= func< n | (n lt 8) select n else n*Self(n-7) >;
[1] cat [b(n): n in [1..40]]; // G. C. Greubel, Aug 20 2019

A114799 := proc(n)
    option remember;
    if n < 1 then
        1;
    else
        n*procname(n-7) ;
    end if;
end proc:
seq(A114799(n),n=0..40) ; # R. J. Mathar, Jun 23 2014
A114799 := n -> product(n-7*k,k=0..(n-1)/7); # M. F. Hasler, Feb 23 2018

a[n_]:= If[n<1, 1, n*a[n-7]]; Table[a[n], {n,0,40}] (* G. C. Greubel, Aug 20 2019 *)

A114799(n,k=7)=prod(j=0,(n-1)\k,n-j*k) \\ M. F. Hasler, Feb 23 2018

def a(n):
    if (n<1): return 1
    else: return n*a(n-7)
[a(n) for n in (0..40)] # G. C. Greubel, Aug 20 2019

a(n) = 1 for n <= 1, else a(n) = n*a(n-7).

Sum_{n>=0} 1/a(n) = A288094. - Amiram Eldar, Nov 10 2020

Edited by M. F. Hasler, Feb 23 2018

cp's OEIS Frontend

A104035 Triangle T(n,k), 0 <= k <= n, read by rows, defined by T(0,0) = 1; T(0,k) = 0 if k>0 or if k<0; T(n,k) = k*T(n-1,k-1) + (k+1)*T(n-1,k+1).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Mathematica

Formula

Extensions

A132382 Lower triangular array T(n,k) generator for group of arrays related to A001147 and A102625.

Views

Author

Keywords

Comments

Examples

Formula

Extensions

A051188 Sept-factorial numbers.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

PARI

Formula

A155100 Triangle read by rows: coefficients in polynomials P_n(u) arising from the expansion of D^(n-1) (tan x) in increasing powers of tan x for n>=1 and 1 for n=0.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A037223 Number of solutions to non-attacking rooks problem on n X n board that are invariant under 180-degree rotation.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Formula

Extensions

A075497 Stirling2 triangle with scaled diagonals (powers of 2).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Sage

Formula

A051189 Octo-factorial numbers.

Views

Author

A104035 Triangle T(n,k), 0 <= k <= n, read by rows, defined by T(0,0) = 1; T(0,k) = 0 if k>0 or if k<0; T(n,k) = kT(n-1,k-1) + (k+1)T(n-1,k+1).