A377023 -id:A377023 - cp's OEIS frontend

A213080 Decimal expansion of Product_{n>=1} n! /(sqrt(2Pin) * (n/e)^n * (1+1/n)^(1/12)).

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

Offset: 1

Michael David Hirschhorn, Jun 04 2012

Just as Stirling's formula for the asymptotic expansion of n! involves the constant sqrt{2 Pi}, the asymptotic expansion of the product of all binomial coefficients in a row of Pascal's triangle involves a constant, the reciprocal of the constant C defined and evaluated here.
From Bernd C. Kellner, Oct 13 2024: (Start)
It turns out that 1/C is not the complete asymptotic constant for the product of the binomial coefficients in a row of Pascal's triangle. A constant factor of (2*Pi)^(-1/4) was overlooked in the asymptotic expansion of that product given by Hirschhorn in 2013. The correct asymptotic constant is A377023.
However, the constant C equals the constant F(1) as introduced before in Kellner 2009. The constants F(1), F(2), ... occur in the same context of asymptotic constants related to asymptotic products of factorials as well as of binomial and multinomial coefficients. Moreover, the sequence (F(k)){k >= 1} is strictly decreasing with limit 1. For example, for k >= 1 the asymptotic product Prod{v >= 1} (k*v)! has the asymptotic constant F(k)*A^k*(2*Pi)^(1/4), where A = A074962 denotes the Glaisher-Kinkelin constant. Let gamma = A001620 be Euler's constant and Gamma(x) be the gamma function.
For k >= 1, the constants F(k) can be computed by an explicit formula and a divergent series expansion, as follows. We have log(F(k)) = (1/(12*k))*(1-log(k)) + (k/4)*log(2*Pi) - ((k^2+1)/k)*log(A) - Sum_{v=1..k-1} (v/k)*log(Gamma(v/k)) = gamma/(12*k) - t*zeta(3)/(360*k^3) with some t in (0,1), respectively.
It follows that log(F(1)) = 1/12 + log(2*Pi)/4 - 2*log(A) = gamma/12 - t*zeta(3)/360 with some t in (0,1), and so F(1) lies in the interval (1.0457...,1.0492...) (see Kellner 2009 and 2024). (End)

			1.04633506677050318098095065697776037101974218113264442441587534042035751563744...

G. C. Greubel, Table of n, a(n) for n = 1..10000
Michael D. Hirschhorn, On the asymptotic behavior of Product_{k=0..n} C(n,k), Fib. Q., 51 (2013), 163-173.
Bernd C. Kellner, On asymptotic constants related to products of Bernoulli numbers and factorials, Integers 9 (2009), Article #A08, 83-106; alternative link; arXiv:0604505 [math.NT], 2006.
Bernd C. Kellner, Asymptotic products of binomial and multinomial coefficients revisited, Integers 24 (2024), Article #A59, 10 pp.; arXiv:2312.11369 [math.CO], 2023.

Cf. A074962, A000178, A084448, A241140, A272097.

Cf. A001620, A002117, A377023, A377024.

exp(2*Zeta(1,-1)-1/12)*(2*Pi)^(1/4); evalf(%,100); # Peter Luschny, Jun 22 2012

RealDigits[(Exp[1]^(1/12) (2 Pi)^(1/4))/Glaisher^2, 10, 100][[1]] (*Peter Luschny, Jun 22 2012 *)

exp(2*zeta'(-1)-1/12)*(2*Pi)^(1/4) \\ Charles R Greathouse IV, Dec 12 2013

import mpmath
mpmath.mp.pretty=True; mpmath.mp.dps = 200 #precision
mpmath.exp(2*mpmath.zeta(-1,1,1)-1/12)*(2*pi)^(1/4) # Peter Luschny, Jun 22 2012

Equals (exp(1)^(1/12)*(2*Pi)^(1/4))/A^2 where A denotes the Glaisher-Kinkelin constant.
Equals exp(2*zeta'(-1)-1/12)*(2*Pi)^(1/4).
A closely related constant is K = Product_{n>=1} (n!*(e/n)^(n+1/2))/ ((1+1/(n+1/2))^(1/12)*sqrt(2*Pi*e)) = (2^(1/6)*(3*e)^(1/12)*Pi^(1/4))/A^2 = exp(2*zeta'(-1)-1/12)*2^(1/6)*3^(1/12)*Pi^(1/4) = 1.082293504658977773529439... - Peter Luschny, Jun 22 2012
The sqrt of the constant equals Limit_{n>=1} (Product_{k=1..n-1} k!) / f(n) where f(n) = (2*Pi)^(n/2-1/8)*exp(1/24-3/4*n^2)*n^(1/2*n^2-1/12). - Peter Luschny, Jun 23 2012

A377024 Decimal expansion of the constant F(2) related to asymptotic products of factorials.

Original entry on oeis.org

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

Offset: 1

Bernd C. Kellner, Oct 13 2024

The constants F(1) = A213080, F(2), ... occur in the context of asymptotic constants related to asymptotic products of factorials as well as of binomial and multinomial coefficients. Moreover, the sequence (F(k)){k >= 1} is strictly decreasing with limit 1. For example, for k >= 1 the asymptotic product Prod{v >= 1} (k*v)! has the asymptotic constant F(k)*A^k*(2*Pi)^(1/4), where A = A074962 denotes the Glaisher-Kinkelin constant. Let gamma = A001620 be Euler's constant and Gamma(x) be the gamma function.
For k >= 1, the constants F(k) can be computed by an explicit formula and a divergent series expansion, as follows. We have log(F(k)) = (1/(12*k))*(1-log(k)) + (k/4)*log(2*Pi) - ((k^2+1)/k)*log(A) - Sum_{v=1..k-1} (v/k)*log(Gamma(v/k)) = gamma/(12*k) - t*zeta(3)/(360*k^3) with some t in (0,1), respectively.
It follows that log(F(2)) = 1/24 + log(2*Pi)/4 + (5/24)*log(2) - (5/2)*log(A) = gamma/24 - t*zeta(3)/2880 with some t in (0,1), and so F(2) lies in the interval (1.023914..., 1.024342...) (see Kellner 2009 and 2024).

			1.02393741163711840157795078258621780080376098043644005129469909513476924124007...

Bernd C. Kellner, Table of n, a(n) for n = 1..10000
Bernd C. Kellner, On asymptotic constants related to products of Bernoulli numbers and factorials, Integers 9 (2009), Article #A08, 83-106; alternative link; arXiv:0604505 [math.NT], 2006.
Bernd C. Kellner, Asymptotic products of binomial and multinomial coefficients revisited, Integers 24 (2024), Article #A59, 10 pp.; arXiv:2312.11369 [math.CO], 2023.

Cf. A001620, A002117, A074962, A213080, A377023.

exp(-1/6+5/2*Zeta(1, -1))*(2*Pi)^(1/4)*2^(5/24); evalf(%, 100);

RealDigits[Exp[1/24] (2 Pi)^(1/4) 2^(5/24) / Glaisher^(5/2), 10, 100][[1]]

default(realprecision, 100);
exp(-1/6+5/2*zeta'(-1))*(2*Pi)^(1/4)*2^(5/24)

import mpmath
mpmath.mp.pretty = True; mpmath.mp.dps = 100
mpmath.exp(-1/6+5/2*mpmath.zeta(-1, 1, 1))*(2*pi)^(1/4)*2^(5/24)

Equals exp(1/24)*(2*Pi)^(1/4)*2^(5/24)/A^(5/2) where A = A074962.

Equals exp(-1/6+(5/2)*zeta'(-1))*(2*Pi)^(1/4)*2^(5/24).

cp's OEIS Frontend

A213080 Decimal expansion of Product_{n>=1} n! /(sqrt(2Pin) * (n/e)^n * (1+1/n)^(1/12)).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Sage

Formula

A377024 Decimal expansion of the constant F(2) related to asymptotic products of factorials.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Sage

Formula

cp's OEIS Frontend

A213080 Decimal expansion of Product_{n>=1} n! /(sqrt(2*Pi*n) * (n/e)^n * (1+1/n)^(1/12)).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Sage

Formula

A377024 Decimal expansion of the constant F(2) related to asymptotic products of factorials.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Sage

Formula

A213080 Decimal expansion of Product_{n>=1} n! /(sqrt(2Pin) * (n/e)^n * (1+1/n)^(1/12)).