Higher-order Bernoulli numbers
DefinitionHigher-order Bernoulli numbers#✓fully proved
For a complex number \(\alpha\) and a natural number \(n\), the higher-order Bernoulli number \(B_n^{(\alpha)}\) is \(n!\) times the coefficient of \(t^n\) in the formal power series \(\left(t / (e^t - 1)\right)^{\alpha}\). Here \(t / (e^t - 1)\) is the Bernoulli generating series \(\sum_{k \ge 0} B_k t^k / k!\) (with \(B_1 = -1/2\)), and its \(\alpha\)-th power is the binomial series \((1 + u)^{\alpha} = \sum_{j \ge 0} \binom{\alpha}{j} u^j\) with \(u = t / (e^t - 1) - 1\) substituted.
Equivalently, \(B_n^{(\alpha)} = B_n^{(\alpha)}(0)\), where the generalized Bernoulli polynomials \(B_n^{(\alpha)}(x)\) are given by the exponential generating function \(\sum_{n \ge 0} B_n^{(\alpha)}(x) t^n / n! = \left(t / (e^t - 1)\right)^{\alpha} e^{tx}\).
public noncomputable def higherOrderBernoulliNumber (α : ℂ) (n : ℕ) : ℂ :=
(n.factorial : ℂ) * PowerSeries.coeff n
(PowerSeries.subst (bernoulliPowerSeries ℂ - 1) (PowerSeries.binomialSeries ℂ α))
Sources
- https://cs.uwaterloo.ca/journals/JIS/VOL23/Chellal/chellal7.tex, lines 99-112.
TheoremHigher-order Bernoulli numbers of natural order#✓fully proved
Let \(B_k\) denote the Bernoulli numbers, with \(B_1 = -1/2\). For natural numbers \(d\) and \(n\), the higher-order Bernoulli number \(B_n^{(d)}\) equals \(n!\) times the coefficient of \(t^n\) in the \(d\)-th power of the Bernoulli generating series \(\sum_{k \ge 0} B_k t^k / k!\).
public theorem higherOrderBernoulliNumber_nat (d n : ℕ) :
higherOrderBernoulliNumber (d : ℂ) n =
(n.factorial : ℂ) * PowerSeries.coeff n (bernoulliPowerSeries ℂ ^ d)
public noncomputable def higherOrderBernoulliNumber (α : ℂ) (n : ℕ) : ℂ :=
(n.factorial : ℂ) * PowerSeries.coeff n
(PowerSeries.subst (bernoulliPowerSeries ℂ - 1) (PowerSeries.binomialSeries ℂ α))
MathlibExt/NumberTheory/NetoHigherOrderBernoulli.lean:97
Uses 1 · Used by 2
TheoremHigher-order Bernoulli numbers of order one#✓fully proved
For every natural number \(n\), the higher-order Bernoulli number \(B_n^{(1)}\) equals the Bernoulli number \(B_n\), with the convention \(B_1 = -1/2\).
public theorem higherOrderBernoulliNumber_one (n : ℕ) :
higherOrderBernoulliNumber 1 n = (bernoulli n : ℂ)
public noncomputable def higherOrderBernoulliNumber (α : ℂ) (n : ℕ) : ℂ :=
(n.factorial : ℂ) * PowerSeries.coeff n
(PowerSeries.subst (bernoulliPowerSeries ℂ - 1) (PowerSeries.binomialSeries ℂ α))
MathlibExt/NumberTheory/NetoHigherOrderBernoulli.lean:97
Uses 2
TheoremNeto's identity for higher-order Bernoulli numbers#✓fully proved
For every complex number \(x\) and all natural numbers \(n\) and \(\ell\),
where \(B_m^{(x)}\) is the higher-order Bernoulli number of order \(x\).
Proof: The substitution coefficient formula makes both sides polynomial in the order. For natural orders, the cited Bernoulli EGF gives the reflection identity, formal differentiation proves Neto's formula, and equality on all natural complex points finishes the polynomial argument.
public theorem neto_higher_order_bernoulli_identity (x : ℂ) (n ℓ : ℕ) :
(∑ k ∈ Finset.range (n + 1),
x ^ (n - k) * (Nat.choose n k : ℂ) * higherOrderBernoulliNumber x (ℓ + k)) =
(-1 : ℂ) ^ (n + ℓ) *
∑ k ∈ Finset.range (ℓ + 1),
x ^ (ℓ - k) * (Nat.choose ℓ k : ℂ) * higherOrderBernoulliNumber x (n + k)
public noncomputable def higherOrderBernoulliNumber (α : ℂ) (n : ℕ) : ℂ :=
(n.factorial : ℂ) * PowerSeries.coeff n
(PowerSeries.subst (bernoulliPowerSeries ℂ - 1) (PowerSeries.binomialSeries ℂ α))
MathlibExt/NumberTheory/NetoHigherOrderBernoulli.lean:97
Sources
- https://cs.uwaterloo.ca/journals/JIS/VOL23/Chellal/chellal7.tex, lines 662-670.
Uses 2
Loading dependency explorer…
The dependency explorer needs JavaScript and its graph data.