OFFSET
0,1
FORMULA
a(n) ~ c * n^(n^2 + 2*n + 5/6) * (2*Pi)^(n+1) / (A^2 * exp(3*n^2/2 + 2*n - 1/6)), where c = Product_{k>=0} (1 + 1/k!^2) = 5.1481781945902396880952694880498895... and A is the Glaisher-Kinkelin constant A074962.
MATHEMATICA
Table[Product[k!^2 + 1, {k, 0, n}], {n, 0, 12}]
Table[BarnesG[n+2]^2 * Product[1 + 1/k!^2, {k, 0, n}], {n, 0, 12}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, Mar 26 2019
STATUS
approved