|
| |
|
|
A055462
|
|
Superduperfactorials: product of first n superfactorials.
|
|
22
|
|
|
|
1, 1, 2, 24, 6912, 238878720, 5944066965504000, 745453331864786829312000000, 3769447945987085350501386572267520000000000, 6916686207999802072984424331678589933649915805696000000000000000
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
|
OFFSET
|
0,3
|
|
|
COMMENTS
|
Next term has 92 digits and is too large to display.
Starting with offset 1, a(n) is a 'Matryoshka doll' sequence with alpha=1, the multiplicative counterpart to the additive A000332. The sequence for m with alpha<=m<=L is then computed as Prod_{n=alpha..m}(Prod_{k=alpha..n}(Prod_{i=alpha..k}(i))). - Peter Luschny, Jul 14 2009
|
|
|
LINKS
|
|
|
|
FORMULA
|
a(n) = a(n-1)*A000178(n) = Product_{i=1..n} (i!)^(n-i+1) = Product_{i=1..n} i^((n-i+1)*(n-i+2)/2).
a(n) = exp((6 + 13 n + 9 n^2 + 2 n^3 - 8*(n + 2)*log(A)-2*(n + 2)^2*log(2*Pi) + 4*(2 n + 1)*logG(n + 2) - 4*(n + 1)^2*logGamma(n + 2) + 8*psi(-3, n + 2))/8) where A is the Glaisher-Kinkelin constant, logG(z) is the logarithm of the Barnes G function (A000178), and psi(-3, z) is a polygamma function of negative order (i.e., the second iterated integral of logGamma(z)). - Jan Mangaldan, Mar 21 2013
a(n) ~ exp(Zeta(3)/(8*Pi^2) - (2*n+3)*(11*n^2 + 24*n - 3)/72) * n^((2*n+3)*(2*n^2 + 6*n + 3)/24) * (2*Pi)^((n+1)*(n+2)/4) / A^(n+3/2), where A = A074962 = 1.28242712910062263687... is the Glaisher-Kinkelin constant and Zeta(3) = A002117 = 1.2020569031595942853997... . - Vaclav Kotesovec, Feb 20 2015
|
|
|
EXAMPLE
|
a(4) = 1!2!3!4!*1!2!3!*1!2!*1! = 288*12*2*1 = 6912.
|
|
|
MAPLE
|
seq(mul(mul(mul(i, i=alpha..k), k=alpha..n), n=alpha..m), m=alpha..10); # Peter Luschny, Jul 14 2009
|
|
|
MATHEMATICA
|
Table[Product[BarnesG[j], {j, k + 1}], {k, 10}] (* Jan Mangaldan, Mar 21 2013 *)
Table[Round[Exp[(n+2)*(n+3)*(2*n+5)/8] * Exp[PolyGamma[-3, n+3]] * BarnesG[n+3]^(n+3/2) / (Glaisher^(n+3) * (2*Pi)^((n+3)^2/4) * Gamma[n+3]^((n+2)^2/2))], {n, 0, 10}] (* Vaclav Kotesovec, Feb 20 2015 after Jan Mangaldan *)
Nest[FoldList[Times, #]&, Range[0, 15]!, 2] (* Harvey P. Dale, Jul 14 2023 *)
|
|
|
PROG
|
(Magma) [n eq 0 select 1 else (&*[j^Binomial(n-j+2, 2): j in [1..n]]): n in [0..10]]; // G. C. Greubel, Jan 31 2024
(SageMath) [product(j^binomial(n-j+2, 2) for j in range(1, n+1)) for n in range(11)] # G. C. Greubel, Jan 31 2024
|
|
|
CROSSREFS
|
|
|
|
KEYWORD
|
nonn,easy
|
|
|
AUTHOR
|
|
|
|
EXTENSIONS
|
|
|
|
STATUS
|
approved
|
| |
|
|
