Products of factorials and binomial coefficients
DefinitionGlaisher–Kinkelin constant#✓fully proved
The Glaisher–Kinkelin constant is the real number
where \(\zeta'\) is the complex derivative of the Riemann zeta function \(\zeta\). Since \(\zeta'(-1)\) is real, \(\log A = \frac{1}{12} - \zeta'(-1)\).
public noncomputable def glaisherKinkelinConstant : ℝ :=
Real.exp (1 / 12 - (deriv riemannZeta (-1 : ℂ)).re)
Sources
- Bernd C. Kellner, Asymptotic Products of Binomial and Multinomial Coefficients Revisited, INTEGERS 24 (2024), Article A59, Section 2.
Used by 8
TheoremLogarithm of the Glaisher–Kinkelin constant#✓fully proved
The Glaisher–Kinkelin constant \(A\) satisfies \(\log A = \frac{1}{12} - \operatorname{Re} \zeta'(-1)\), where \(\zeta'\) is the complex derivative of the Riemann zeta function.
public theorem log_glaisherKinkelinConstant :
Real.log glaisherKinkelinConstant = 1 / 12 - (deriv riemannZeta (-1 : ℂ)).re
public noncomputable def glaisherKinkelinConstant : ℝ :=
Real.exp (1 / 12 - (deriv riemannZeta (-1 : ℂ)).re)
WantedExt/NumberTheory/BinomialProductAsymptoticsWanted.lean:27
Sources
- Bernd C. Kellner, Asymptotic Products of Binomial and Multinomial Coefficients Revisited, INTEGERS 24 (2024), Article A59, Section 2.
Uses 1
DefinitionFactorial-product constants#✓fully proved
For a natural number \(k \ge 1\), the constant \(F_k\) is the positive real number defined by
where \(A\) is the Glaisher–Kinkelin constant and \(\Gamma\) is the gamma function. For \(k = 0\) the definition gives \(F_0 = 1\).
The constant \(F_k\) appears in the asymptotic formula for \(\prod_{\nu=1}^{n} (k\nu)!\) as \(n \to \infty\); that formula is a separate statement.
public noncomputable def factorialProductConstant (k : ℕ) : ℝ :=
Real.exp ((k : ℝ) / 4 * Real.log (2 * Real.pi) -
((k : ℝ) ^ 2 + 1) / k * Real.log glaisherKinkelinConstant +
(1 - Real.log k) / (12 * k) -
∑ ν ∈ Finset.Ico 1 k, (ν : ℝ) / k * Real.log (Real.Gamma ((ν : ℝ) / k)))
public noncomputable def glaisherKinkelinConstant : ℝ :=
Real.exp (1 / 12 - (deriv riemannZeta (-1 : ℂ)).re)
WantedExt/NumberTheory/BinomialProductAsymptoticsWanted.lean:27
Sources
- Bernd C. Kellner, Asymptotic Products of Binomial and Multinomial Coefficients Revisited, INTEGERS 24 (2024), Article A59, equation (3) (Kellner 2009, Theorems 12 and 13).
Uses 1 · Used by 10
TheoremLogarithm of the first factorial product constant#✓fully proved
The factorial product constant \(F_1\) satisfies \(\log F_1 = \frac{1}{4} \log(2\pi) + \frac{1}{12} - 2 \log A\), where \(A\) is the Glaisher–Kinkelin constant.
Equivalently, \(F_1 = (2\pi)^{1/4} e^{1/12} A^{-2}\).
public theorem log_factorialProductConstant_one :
Real.log (factorialProductConstant 1) =
1 / 4 * Real.log (2 * Real.pi) + 1 / 12 - 2 * Real.log glaisherKinkelinConstant
public noncomputable def factorialProductConstant (k : ℕ) : ℝ :=
Real.exp ((k : ℝ) / 4 * Real.log (2 * Real.pi) -
((k : ℝ) ^ 2 + 1) / k * Real.log glaisherKinkelinConstant +
(1 - Real.log k) / (12 * k) -
∑ ν ∈ Finset.Ico 1 k, (ν : ℝ) / k * Real.log (Real.Gamma ((ν : ℝ) / k)))
WantedExt/NumberTheory/BinomialProductAsymptoticsWanted.lean:46
public noncomputable def glaisherKinkelinConstant : ℝ :=
Real.exp (1 / 12 - (deriv riemannZeta (-1 : ℂ)).re)
WantedExt/NumberTheory/BinomialProductAsymptoticsWanted.lean:27
Sources
- Bernd C. Kellner, Asymptotic Products of Binomial and Multinomial Coefficients Revisited, INTEGERS 24 (2024), Article A59, Table 1.
Uses 2
TheoremPositivity of the factorial-product constants#✓fully proved
For every natural number \(k\), the factorial-product constant \(F_k\) is positive.
public theorem factorialProductConstant_pos (k : ℕ) : 0 < factorialProductConstant k
public noncomputable def factorialProductConstant (k : ℕ) : ℝ :=
Real.exp ((k : ℝ) / 4 * Real.log (2 * Real.pi) -
((k : ℝ) ^ 2 + 1) / k * Real.log glaisherKinkelinConstant +
(1 - Real.log k) / (12 * k) -
∑ ν ∈ Finset.Ico 1 k, (ν : ℝ) / k * Real.log (Real.Gamma ((ν : ℝ) / k)))
WantedExt/NumberTheory/BinomialProductAsymptoticsWanted.lean:46
Sources
- Bernd C. Kellner, Asymptotic Products of Binomial and Multinomial Coefficients Revisited, INTEGERS 24 (2024), Article A59, Section 2.
Uses 2
DefinitionComparison function for products of factorials#✓fully proved
For a natural number \(k\), the comparison function \(P_k \colon \mathbb{R} \to \mathbb{R}\) is defined by
public noncomputable def factorialProductComparison (k : ℕ) (x : ℝ) : ℝ :=
((k : ℝ) ^ k) ^ (x * (x + 1) / 2) * (2 * Real.pi * k / Real.exp 1 * x) ^ (x / 2) *
x ^ ((3 * (k : ℝ) + 1) / (12 * k))
Sources
- Bernd C. Kellner, Asymptotic Products of Binomial and Multinomial Coefficients Revisited, INTEGERS 24 (2024), Article A59, Theorem 5.
TheoremValue of the factorial-product comparison function at \(1\)#✓fully proved
For every natural number \(k\), the factorial-product comparison function \(P_k\) satisfies \(P_k(1) = k^k (2\pi k/e)^{1/2}\).
public theorem factorialProductComparison_one (k : ℕ) :
factorialProductComparison k 1 =
(k : ℝ) ^ k * (2 * Real.pi * k / Real.exp 1) ^ ((1 : ℝ) / 2)
public noncomputable def factorialProductComparison (k : ℕ) (x : ℝ) : ℝ :=
((k : ℝ) ^ k) ^ (x * (x + 1) / 2) * (2 * Real.pi * k / Real.exp 1 * x) ^ (x / 2) *
x ^ ((3 * (k : ℝ) + 1) / (12 * k))
WantedExt/NumberTheory/BinomialProductAsymptoticsWanted.lean:74
Sources
- Bernd C. Kellner, Asymptotic Products of Binomial and Multinomial Coefficients Revisited, INTEGERS 24 (2024), Article A59, Theorem 5.
Uses 1
DefinitionComparison function for products of binomial coefficients#✓fully proved
For natural numbers \(a\) and \(b\), let \(c = a - b\) if \(b \le a\) and \(c = 0\) otherwise. The comparison function \(P_{a,b} \colon \mathbb{R} \to \mathbb{R}\) is defined by
public noncomputable def binomialProductComparison (a b : ℕ) (x : ℝ) : ℝ :=
((a : ℝ) ^ a / ((b : ℝ) ^ b * ((a - b : ℕ) : ℝ) ^ (a - b))) ^ (x * (x + 1) / 2) *
x ^ ((1 / (a : ℝ) - a / (b * ((a - b : ℕ) : ℝ))) / 12 - 1 / 4) /
(2 * Real.pi / Real.exp 1 * (b * ((a - b : ℕ) : ℝ) / a) * x) ^ (x / 2)
Sources
- Bernd C. Kellner, Asymptotic Products of Binomial and Multinomial Coefficients Revisited, INTEGERS 24 (2024), Article A59, Theorem 6.
TheoremValue of the binomial-product comparison function at \(1\)#✓fully proved
For natural numbers \(a\) and \(b\), let \(c = a - b\) if \(b \le a\) and \(c = 0\) otherwise. Then the binomial-product comparison function \(P_{a,b}\) satisfies
public theorem binomialProductComparison_one (a b : ℕ) :
binomialProductComparison a b 1 =
(a : ℝ) ^ a / ((b : ℝ) ^ b * ((a - b : ℕ) : ℝ) ^ (a - b)) /
(2 * Real.pi / Real.exp 1 * (b * ((a - b : ℕ) : ℝ) / a)) ^ ((1 : ℝ) / 2)
public noncomputable def binomialProductComparison (a b : ℕ) (x : ℝ) : ℝ :=
((a : ℝ) ^ a / ((b : ℝ) ^ b * ((a - b : ℕ) : ℝ) ^ (a - b))) ^ (x * (x + 1) / 2) *
x ^ ((1 / (a : ℝ) - a / (b * ((a - b : ℕ) : ℝ))) / 12 - 1 / 4) /
(2 * Real.pi / Real.exp 1 * (b * ((a - b : ℕ) : ℝ) / a) * x) ^ (x / 2)
WantedExt/NumberTheory/BinomialProductAsymptoticsWanted.lean:92
Sources
- Bernd C. Kellner, Asymptotic Products of Binomial and Multinomial Coefficients Revisited, INTEGERS 24 (2024), Article A59, Theorem 6.
Uses 1
TheoremAsymptotics of products of factorials of multiples#●ready to prove
Let \(k \ge 1\) be a natural number, let \(F_k\) be the \(k\)-th factorial product constant and \(A\) the Glaisher–Kinkelin constant. As \(n \to \infty\), the ratio of \(\prod_{\nu=1}^{n} (k\nu)!\) to
tends to \(F_k\).
The divergent asymptotic expansion (4) of the source is not part of this statement.
public theorem_wanted tendsto_prod_factorial_mul_div (k : ℕ) (hk : 1 ≤ k) :
Tendsto (fun n : ℕ => (∏ ν ∈ Finset.Icc 1 n, ((k * ν).factorial : ℝ)) /
(glaisherKinkelinConstant ^ k * (2 * Real.pi) ^ ((1 : ℝ) / 4) *
((k : ℝ) * Real.exp (-3 / 2) * n) ^ ((k : ℝ) * (n * (n + 1) / 2)) *
(2 * Real.pi * k * Real.exp ((k : ℝ) / 2 - 1) * n) ^ ((n : ℝ) / 2) *
(n : ℝ) ^ (((k : ℝ) ^ 2 + 3 * k + 1) / (12 * k))))
atTop (𝓝 (factorialProductConstant k))
public noncomputable def glaisherKinkelinConstant : ℝ :=
Real.exp (1 / 12 - (deriv riemannZeta (-1 : ℂ)).re)
WantedExt/NumberTheory/BinomialProductAsymptoticsWanted.lean:27
public noncomputable def factorialProductConstant (k : ℕ) : ℝ :=
Real.exp ((k : ℝ) / 4 * Real.log (2 * Real.pi) -
((k : ℝ) ^ 2 + 1) / k * Real.log glaisherKinkelinConstant +
(1 - Real.log k) / (12 * k) -
∑ ν ∈ Finset.Ico 1 k, (ν : ℝ) / k * Real.log (Real.Gamma ((ν : ℝ) / k)))
WantedExt/NumberTheory/BinomialProductAsymptoticsWanted.lean:46
Sources
- Bernd C. Kellner, Asymptotic Products of Binomial and Multinomial Coefficients Revisited, INTEGERS 24 (2024), Article A59, Theorem 1 (Kellner 2009, Theorems 12 and 13).
Uses 2
TheoremThe factorial-product constants decrease strictly to \(1\)#●ready to prove
The sequence \((F_k)_{k \ge 1}\) of factorial-product constants is strictly decreasing, and \(F_k \to 1\) as \(k \to \infty\).
public theorem_wanted factorialProductConstant_strictAnti_tendsto_one :
StrictAnti (fun k : ℕ => factorialProductConstant (k + 1)) ∧
Tendsto factorialProductConstant atTop (𝓝 1)
public noncomputable def factorialProductConstant (k : ℕ) : ℝ :=
Real.exp ((k : ℝ) / 4 * Real.log (2 * Real.pi) -
((k : ℝ) ^ 2 + 1) / k * Real.log glaisherKinkelinConstant +
(1 - Real.log k) / (12 * k) -
∑ ν ∈ Finset.Ico 1 k, (ν : ℝ) / k * Real.log (Real.Gamma ((ν : ℝ) / k)))
WantedExt/NumberTheory/BinomialProductAsymptoticsWanted.lean:46
Sources
- Bernd C. Kellner, Asymptotic Products of Binomial and Multinomial Coefficients Revisited, INTEGERS 24 (2024), Article A59, Corollary 2.
Uses 1
TheoremBounds for \(1/a + 1/b - 1/(a+b)\)#●ready to prove
Let \(a\) and \(b\) be natural numbers with \(a, b \ge 1\) and \(a + b \ge 3\), and put \(g(a,b) = 1/a + 1/b - 1/(a+b)\). Then \(0 < g(a,b) \le 7/6\).
public theorem_wanted inv_add_inv_sub_inv_add_bounds (a b : ℕ) (ha : 1 ≤ a) (hb : 1 ≤ b)
(hab : 3 ≤ a + b) :
0 < (1 / a + 1 / b - 1 / (a + b) : ℚ) ∧ (1 / a + 1 / b - 1 / (a + b) : ℚ) ≤ 7 / 6
Sources
- Bernd C. Kellner, Asymptotic Products of Binomial and Multinomial Coefficients Revisited, INTEGERS 24 (2024), Article A59, Lemma 3.
TheoremBounds for the binomial-product constants#●ready to prove
Let \(a\) and \(b\) be natural numbers with \(1 \le b < a\), and let \(F_k\) denote the factorial-product constants. Then \(c_{a,b} = F_a / (F_b F_{a-b})\) satisfies
that is, \(c_{2,1} \le c_{a,b} < 1\).
public theorem_wanted factorialProductConstant_ratio_mem_Ico (a b : ℕ) (hb : 1 ≤ b)
(hab : b < a) :
factorialProductConstant a / (factorialProductConstant b * factorialProductConstant (a - b))
∈ Set.Ico (factorialProductConstant 2 / factorialProductConstant 1 ^ 2) 1
public noncomputable def factorialProductConstant (k : ℕ) : ℝ :=
Real.exp ((k : ℝ) / 4 * Real.log (2 * Real.pi) -
((k : ℝ) ^ 2 + 1) / k * Real.log glaisherKinkelinConstant +
(1 - Real.log k) / (12 * k) -
∑ ν ∈ Finset.Ico 1 k, (ν : ℝ) / k * Real.log (Real.Gamma ((ν : ℝ) / k)))
WantedExt/NumberTheory/BinomialProductAsymptoticsWanted.lean:46
Sources
- Bernd C. Kellner, Asymptotic Products of Binomial and Multinomial Coefficients Revisited, INTEGERS 24 (2024), Article A59, Lemma 4.
Uses 1
TheoremAsymptotics of products of multinomial coefficients#●ready to prove
Let \(r \ge 2\), let \(b_1, \dots, b_r \ge 1\) be natural numbers with \(a = b_1 + \dots + b_r\), let \(F_k\) denote the factorial product constants and \(P_k\) the factorial comparison functions, and put \(c_{a,b} = F_a / (F_{b_1} \cdots F_{b_r})\). As \(m \to \infty\), the ratio of \(\prod_{n=1}^{m} \binom{an}{b_1 n, \dots, b_r n}\) to \(P_a(m) / (P_{b_1}(m) \cdots P_{b_r}(m))\) tends to \(c_{a,b} (2\pi)^{(1-r)/4}\); moreover \(0 < c_{a,b} < 1\).
public theorem_wanted tendsto_prod_multinomial_div (r : ℕ) (hr : 2 ≤ r) (b : Fin r → ℕ)
(hb : ∀ i, 1 ≤ b i) :
Tendsto (fun m : ℕ => (∏ n ∈ Finset.Icc 1 m,
(Nat.multinomial Finset.univ (fun i => b i * n) : ℝ)) /
(factorialProductComparison (∑ i, b i) m /
∏ i, factorialProductComparison (b i) m))
atTop (𝓝 (factorialProductConstant (∑ i, b i) / (∏ i, factorialProductConstant (b i)) *
(2 * Real.pi) ^ ((1 - (r : ℝ)) / 4))) ∧
factorialProductConstant (∑ i, b i) / (∏ i, factorialProductConstant (b i)) ∈
Set.Ioo 0 1
public noncomputable def factorialProductComparison (k : ℕ) (x : ℝ) : ℝ :=
((k : ℝ) ^ k) ^ (x * (x + 1) / 2) * (2 * Real.pi * k / Real.exp 1 * x) ^ (x / 2) *
x ^ ((3 * (k : ℝ) + 1) / (12 * k))
WantedExt/NumberTheory/BinomialProductAsymptoticsWanted.lean:74
public noncomputable def factorialProductConstant (k : ℕ) : ℝ :=
Real.exp ((k : ℝ) / 4 * Real.log (2 * Real.pi) -
((k : ℝ) ^ 2 + 1) / k * Real.log glaisherKinkelinConstant +
(1 - Real.log k) / (12 * k) -
∑ ν ∈ Finset.Ico 1 k, (ν : ℝ) / k * Real.log (Real.Gamma ((ν : ℝ) / k)))
WantedExt/NumberTheory/BinomialProductAsymptoticsWanted.lean:46
Sources
- Bernd C. Kellner, Asymptotic Products of Binomial and Multinomial Coefficients Revisited, INTEGERS 24 (2024), Article A59, Theorem 5.
Uses 2
TheoremAsymptotics of products of binomial coefficients \(\binom{an}{bn}\)#●ready to prove
Let \(a > b \ge 1\) be natural numbers, let \(F_k\) denote the factorial product constants and \(P_{a,b}\) the binomial comparison function, and put \(c_{a,b} = F_a / (F_b F_{a-b})\) and \(C_{a,b} = c_{a,b} (2\pi)^{-1/4}\). Then \(\prod_{n=1}^{m} \binom{an}{bn} / P_{a,b}(m) \to C_{a,b}\) as \(m \to \infty\); moreover \(C_{2,1} \le C_{a,b} < (2\pi)^{-1/4}\) and \(c_{2,1} \le c_{a,b} < 1\), where \(c_{2,1} = F_2 / F_1^2\) and \(C_{2,1} = c_{2,1} (2\pi)^{-1/4}\).
The decimal approximations of the constants are not part of this statement.
public theorem_wanted tendsto_prod_choose_div (a b : ℕ) (hb : 1 ≤ b) (hab : b < a) :
Tendsto (fun m : ℕ => (∏ n ∈ Finset.Icc 1 m, ((a * n).choose (b * n) : ℝ)) /
binomialProductComparison a b m)
atTop (𝓝 (factorialProductConstant a /
(factorialProductConstant b * factorialProductConstant (a - b)) *
(2 * Real.pi) ^ (-(1 : ℝ) / 4))) ∧
factorialProductConstant a /
(factorialProductConstant b * factorialProductConstant (a - b)) *
(2 * Real.pi) ^ (-(1 : ℝ) / 4) ∈
Set.Ico (factorialProductConstant 2 / factorialProductConstant 1 ^ 2 *
(2 * Real.pi) ^ (-(1 : ℝ) / 4)) ((2 * Real.pi) ^ (-(1 : ℝ) / 4)) ∧
factorialProductConstant a /
(factorialProductConstant b * factorialProductConstant (a - b)) ∈
Set.Ico (factorialProductConstant 2 / factorialProductConstant 1 ^ 2) 1
public noncomputable def binomialProductComparison (a b : ℕ) (x : ℝ) : ℝ :=
((a : ℝ) ^ a / ((b : ℝ) ^ b * ((a - b : ℕ) : ℝ) ^ (a - b))) ^ (x * (x + 1) / 2) *
x ^ ((1 / (a : ℝ) - a / (b * ((a - b : ℕ) : ℝ))) / 12 - 1 / 4) /
(2 * Real.pi / Real.exp 1 * (b * ((a - b : ℕ) : ℝ) / a) * x) ^ (x / 2)
WantedExt/NumberTheory/BinomialProductAsymptoticsWanted.lean:92
public noncomputable def factorialProductConstant (k : ℕ) : ℝ :=
Real.exp ((k : ℝ) / 4 * Real.log (2 * Real.pi) -
((k : ℝ) ^ 2 + 1) / k * Real.log glaisherKinkelinConstant +
(1 - Real.log k) / (12 * k) -
∑ ν ∈ Finset.Ico 1 k, (ν : ℝ) / k * Real.log (Real.Gamma ((ν : ℝ) / k)))
WantedExt/NumberTheory/BinomialProductAsymptoticsWanted.lean:46
Sources
- Bernd C. Kellner, Asymptotic Products of Binomial and Multinomial Coefficients Revisited, INTEGERS 24 (2024), Article A59, Theorem 6.
Uses 2
TheoremAsymptotics of products of central binomial coefficients and Catalan numbers#●ready to prove
Let \(A\) be the Glaisher–Kinkelin constant, let \(\mathrm{Cat}_n\) denote the \(n\)-th Catalan number, and put \(C_{2,1} = 2^{5/24} A^{3/2} / \left((2\pi)^{1/2} e^{1/8}\right)\). As \(m \to \infty\), the ratio of \(\prod_{n=1}^{m} \binom{2n}{n}\) to \(2^{m^2} (4e/\pi)^{m/2} / m^{m/2 + 3/8}\) tends to \(C_{2,1}\), and the ratio of \(\prod_{n=1}^{m} \mathrm{Cat}_n\) to \(2^{m^2} (4e^3/\pi)^{m/2} / m^{3m/2 + 15/8}\) tends to \(C_{2,1} (2\pi)^{-1/2}\).
The decimal approximations of the constants are not part of this statement.
public theorem_wanted tendsto_prod_centralBinom_catalan_div :
Tendsto (fun m : ℕ => (∏ n ∈ Finset.Icc 1 m, (n.centralBinom : ℝ)) /
(2 ^ ((m : ℝ) ^ 2) * (4 * Real.exp 1 / Real.pi) ^ ((m : ℝ) / 2) /
(m : ℝ) ^ ((m : ℝ) / 2 + 3 / 8)))
atTop (𝓝 (2 ^ ((5 : ℝ) / 24) * glaisherKinkelinConstant ^ ((3 : ℝ) / 2) /
((2 * Real.pi) ^ ((1 : ℝ) / 2) * Real.exp (1 / 8)))) ∧
Tendsto (fun m : ℕ => (∏ n ∈ Finset.Icc 1 m, (catalan n : ℝ)) /
(2 ^ ((m : ℝ) ^ 2) * (4 * Real.exp 3 / Real.pi) ^ ((m : ℝ) / 2) /
(m : ℝ) ^ (3 * (m : ℝ) / 2 + 15 / 8)))
atTop (𝓝 (2 ^ ((5 : ℝ) / 24) * glaisherKinkelinConstant ^ ((3 : ℝ) / 2) /
((2 * Real.pi) ^ ((1 : ℝ) / 2) * Real.exp (1 / 8)) * (2 * Real.pi) ^ (-(1 : ℝ) / 2)))
public noncomputable def glaisherKinkelinConstant : ℝ :=
Real.exp (1 / 12 - (deriv riemannZeta (-1 : ℂ)).re)
WantedExt/NumberTheory/BinomialProductAsymptoticsWanted.lean:27
Sources
- Bernd C. Kellner, Asymptotic Products of Binomial and Multinomial Coefficients Revisited, INTEGERS 24 (2024), Article A59, Corollary 7.
Uses 1
TheoremAsymptotics of the product of a row of Pascal's triangle#●ready to prove
Let \(F_1\) be the first factorial product constant and \(A\) the Glaisher–Kinkelin constant, and put \(C_{\mathrm{row}} = \left(F_1 (2\pi)^{1/4}\right)^{-1}\). As \(n \to \infty\), the ratio of \(\prod_{\nu=0}^{n} \binom{n}{\nu}\) to \(e^{n^2/2 + n} / \left((2\pi)^{n/2} n^{n/2 + 1/3}\right)\) tends to \(C_{\mathrm{row}}\), and \(C_{\mathrm{row}} = A^2 / \left((2\pi)^{1/2} e^{1/12}\right)\).
public theorem_wanted tendsto_prod_choose_row_div :
Tendsto (fun n : ℕ => (∏ ν ∈ Finset.range (n + 1), (n.choose ν : ℝ)) /
(Real.exp ((n : ℝ) ^ 2 / 2 + n) /
((2 * Real.pi) ^ ((n : ℝ) / 2) * (n : ℝ) ^ ((n : ℝ) / 2 + 1 / 3))))
atTop (𝓝 ((factorialProductConstant 1 * (2 * Real.pi) ^ ((1 : ℝ) / 4))⁻¹)) ∧
(factorialProductConstant 1 * (2 * Real.pi) ^ ((1 : ℝ) / 4))⁻¹ =
glaisherKinkelinConstant ^ 2 / ((2 * Real.pi) ^ ((1 : ℝ) / 2) * Real.exp (1 / 12))
public noncomputable def factorialProductConstant (k : ℕ) : ℝ :=
Real.exp ((k : ℝ) / 4 * Real.log (2 * Real.pi) -
((k : ℝ) ^ 2 + 1) / k * Real.log glaisherKinkelinConstant +
(1 - Real.log k) / (12 * k) -
∑ ν ∈ Finset.Ico 1 k, (ν : ℝ) / k * Real.log (Real.Gamma ((ν : ℝ) / k)))
WantedExt/NumberTheory/BinomialProductAsymptoticsWanted.lean:46
public noncomputable def glaisherKinkelinConstant : ℝ :=
Real.exp (1 / 12 - (deriv riemannZeta (-1 : ℂ)).re)
WantedExt/NumberTheory/BinomialProductAsymptoticsWanted.lean:27
Sources
- Bernd C. Kellner, Asymptotic Products of Binomial and Multinomial Coefficients Revisited, INTEGERS 24 (2024), Article A59, Theorem 10.
Uses 2
TheoremPascal-row asymptotics as a limit of binomial-product asymptotics#●ready to prove
Let \(F_k\) denote the factorial product constants and \(P_{a,b}\) the binomial comparison functions, and let \(n \ge 1\) be a natural number and \(x > 0\) a real number. Then, as \(a \to \infty\) through the natural numbers, \(\frac{F_a}{F_1 F_{a-1}} (2\pi)^{-1/4} \to \left(F_1 (2\pi)^{1/4}\right)^{-1}\), \(a^{-n(n+1)/2} \prod_{\nu=1}^{n} \binom{a\nu}{\nu} \to \prod_{\nu=1}^{n} \binom{n}{\nu}\), and \(a^{-x(x+1)/2} P_{a,1}(x) \to e^{x^2/2 + x} / \left((2\pi)^{x/2} x^{x/2 + 1/3}\right)\).
The source does not state the range of \(x\); the hypothesis \(x > 0\) is added here.
public theorem_wanted tendsto_binomialProduct_limit_row (n : ℕ) (hn : 1 ≤ n) (x : ℝ)
(hx : 0 < x) :
Tendsto (fun a : ℕ => factorialProductConstant a /
(factorialProductConstant 1 * factorialProductConstant (a - 1)) *
(2 * Real.pi) ^ (-(1 : ℝ) / 4))
atTop (𝓝 ((factorialProductConstant 1 * (2 * Real.pi) ^ ((1 : ℝ) / 4))⁻¹)) ∧
Tendsto (fun a : ℕ => (a : ℝ) ^ (-((n : ℝ) * (n + 1) / 2)) *
∏ ν ∈ Finset.Icc 1 n, ((a * ν).choose ν : ℝ))
atTop (𝓝 (∏ ν ∈ Finset.Icc 1 n, (n.choose ν : ℝ))) ∧
Tendsto (fun a : ℕ => (a : ℝ) ^ (-(x * (x + 1) / 2)) * binomialProductComparison a 1 x)
atTop (𝓝 (Real.exp (x ^ 2 / 2 + x) / ((2 * Real.pi) ^ (x / 2) * x ^ (x / 2 + 1 / 3))))
public noncomputable def factorialProductConstant (k : ℕ) : ℝ :=
Real.exp ((k : ℝ) / 4 * Real.log (2 * Real.pi) -
((k : ℝ) ^ 2 + 1) / k * Real.log glaisherKinkelinConstant +
(1 - Real.log k) / (12 * k) -
∑ ν ∈ Finset.Ico 1 k, (ν : ℝ) / k * Real.log (Real.Gamma ((ν : ℝ) / k)))
WantedExt/NumberTheory/BinomialProductAsymptoticsWanted.lean:46
public noncomputable def binomialProductComparison (a b : ℕ) (x : ℝ) : ℝ :=
((a : ℝ) ^ a / ((b : ℝ) ^ b * ((a - b : ℕ) : ℝ) ^ (a - b))) ^ (x * (x + 1) / 2) *
x ^ ((1 / (a : ℝ) - a / (b * ((a - b : ℕ) : ℝ))) / 12 - 1 / 4) /
(2 * Real.pi / Real.exp 1 * (b * ((a - b : ℕ) : ℝ) / a) * x) ^ (x / 2)
WantedExt/NumberTheory/BinomialProductAsymptoticsWanted.lean:92
Sources
- Bernd C. Kellner, Asymptotic Products of Binomial and Multinomial Coefficients Revisited, INTEGERS 24 (2024), Article A59, Theorem 11.
Uses 2
TheoremHirschhorn's constant#●ready to prove
Let \(A\) be the Glaisher–Kinkelin constant and \(F_1\) the first factorial product constant. As \(n \to \infty\),
and this limit equals \(F_1\).
Hirschhorn's asymptotic expansion with divergent Bernoulli sums is not part of this statement.
public theorem_wanted tendsto_hirschhorn_constant :
Tendsto (fun n : ℕ => (n : ℝ) ^ (-(1 : ℝ) / 12) *
∏ k ∈ Finset.Icc 1 n, ((k.factorial : ℝ) /
(Real.sqrt (2 * Real.pi * k) * ((k : ℝ) / Real.exp 1) ^ k)))
atTop (𝓝 ((2 * Real.pi) ^ ((1 : ℝ) / 4) * Real.exp (1 / 12) /
glaisherKinkelinConstant ^ 2)) ∧
(2 * Real.pi) ^ ((1 : ℝ) / 4) * Real.exp (1 / 12) / glaisherKinkelinConstant ^ 2 =
factorialProductConstant 1
public noncomputable def glaisherKinkelinConstant : ℝ :=
Real.exp (1 / 12 - (deriv riemannZeta (-1 : ℂ)).re)
WantedExt/NumberTheory/BinomialProductAsymptoticsWanted.lean:27
public noncomputable def factorialProductConstant (k : ℕ) : ℝ :=
Real.exp ((k : ℝ) / 4 * Real.log (2 * Real.pi) -
((k : ℝ) ^ 2 + 1) / k * Real.log glaisherKinkelinConstant +
(1 - Real.log k) / (12 * k) -
∑ ν ∈ Finset.Ico 1 k, (ν : ℝ) / k * Real.log (Real.Gamma ((ν : ℝ) / k)))
WantedExt/NumberTheory/BinomialProductAsymptoticsWanted.lean:46
Sources
- Bernd C. Kellner, Asymptotic Products of Binomial and Multinomial Coefficients Revisited, INTEGERS 24 (2024), Article A59, Theorem 12 (Hirschhorn) and the discussion after it.
Uses 2
Loading dependency explorer…
The dependency explorer needs JavaScript and its graph data.