Binomial convolutions of rational power series

DefinitionBinomial (Hurwitz) product of power series#✓fully proved

Let \(K\) be a commutative ring. The binomial product of formal power series \(A = \sum_n a_n X^n\) and \(B = \sum_n b_n X^n\) in \(K[[X]]\) is the power series \(\sum_n c_n X^n\) with \(c_n = \sum_{k=0}^{n} \binom{n}{k} a_k b_{n-k}\).

public noncomputable def binomialProduct {K : Type*} [CommRing K] (A B : PowerSeries K) :
    PowerSeries K :=
  PowerSeries.mk fun n => ∑ k ∈ Finset.range (n + 1),
    (n.choose k : K) * PowerSeries.coeff k A * PowerSeries.coeff (n - k) B
Sources
  • Ira M. Gessel and Ishan Kar, Binomial Convolutions for Rational Power Series, Journal of Integer Sequences 27 (2024), Article 24.1.3, equation (1).
Used by 9

TheoremCoefficients of the binomial product#✓fully proved

Let \(K\) be a commutative ring and let \(A = \sum_n a_n X^n\) and \(B = \sum_n b_n X^n\) be formal power series over \(K\). For every natural number \(n\), the coefficient of \(X^n\) in the binomial product of \(A\) and \(B\) is \(\sum_{k=0}^{n} \binom{n}{k} a_k b_{n-k}\).

public theorem coeff_binomialProduct {K : Type*} [CommRing K] (A B : PowerSeries K) (n : ℕ) :
    PowerSeries.coeff n (binomialProduct A B) = ∑ k ∈ Finset.range (n + 1),
      (n.choose k : K) * PowerSeries.coeff k A * PowerSeries.coeff (n - k) B
Sources
  • Ira M. Gessel and Ishan Kar, Binomial Convolutions for Rational Power Series, Journal of Integer Sequences 27 (2024), Article 24.1.3, equation (1).
Uses 1

TheoremConstant coefficient of the binomial product#✓fully proved

Let \(K\) be a commutative ring. For formal power series \(A\) and \(B\) over \(K\), the constant coefficient of the binomial product of \(A\) and \(B\) is the product of the constant coefficients of \(A\) and \(B\).

public theorem constantCoeff_binomialProduct {K : Type*} [CommRing K] (A B : PowerSeries K) :
    PowerSeries.coeff 0 (binomialProduct A B) = PowerSeries.coeff 0 A * PowerSeries.coeff 0 B
Sources
  • Ira M. Gessel and Ishan Kar, Binomial Convolutions for Rational Power Series, Journal of Integer Sequences 27 (2024), Article 24.1.3, equation (1).
Uses 1

TheoremRationality of power series descends from field extensions#●ready to prove

Let \(F\) be a field, let \(G\) be a field extension of \(F\) (a field with an \(F\)-algebra structure), and let \(A\) be a power series over \(F\). If there are polynomials \(p\), \(q\) over \(G\) with \(q \ne 0\) and \(A q = p\), where \(A\) is viewed as a power series over \(G\), then there are polynomials \(P\), \(Q\) over \(F\) with \(Q \ne 0\) and \(A Q = P\).

public theorem_wanted exists_polynomial_quotient_of_extension {F G : Type*} [Field F] [Field G]
    [Algebra F G] (A : PowerSeries F) (p q : Polynomial G) (hq : q ≠ 0)
    (h : PowerSeries.map (algebraMap F G) A * (q : PowerSeries G) = p) :
    ∃ P Q : Polynomial F, Q ≠ 0 ∧ A * (Q : PowerSeries F) = P
Sources
  • Ira M. Gessel and Ishan Kar, Binomial Convolutions for Rational Power Series, Journal of Integer Sequences 27 (2024), Article 24.1.3, Theorem 1 (jis-2024-24.1.3-thm-1).

TheoremOrthogonal vectors descend from field extensions#●ready to prove

Let \(F\) be a field, let \(G\) be a field extension of \(F\) (a field with an \(F\)-algebra structure), let \(n \in \mathbb{N}\) and let \(S \subseteq F^n\). If some nonzero \(t \in G^n\) satisfies \(\sum_{i=1}^n s_i t_i = 0\) in \(G\) for every \(s \in S\), then some nonzero \(u \in F^n\) satisfies \(\sum_{i=1}^n s_i u_i = 0\) for every \(s \in S\).

public theorem_wanted exists_orthogonal_of_orthogonal_over_extension {F G : Type*} [Field F]
    [Field G] [Algebra F G] {n : ℕ} (S : Set (Fin n → F)) (t : Fin n → G) (ht : t ≠ 0)
    (horth : ∀ s ∈ S, ∑ i, algebraMap F G (s i) * t i = 0) :
    ∃ u : Fin n → F, u ≠ 0 ∧ ∀ s ∈ S, ∑ i, s i * u i = 0
Sources
  • Ira M. Gessel and Ishan Kar, Binomial Convolutions for Rational Power Series, Journal of Integer Sequences 27 (2024), Article 24.1.3, Lemma 3 (jis-2024-24.1.3-lem-3).

TheoremBinomial product of shifted powers of geometric series#●ready to prove

For every field \(K\) of characteristic \(0\), all \(j, k \in \mathbb{N}\) and all \(\alpha, \beta \in K\),

\[\frac{x^j}{(1 - \alpha x)^{j+1}} \odot \frac{x^k}{(1 - \beta x)^{k+1}} = \binom{j + k}{j} \frac{x^{j+k}}{(1 - (\alpha + \beta) x)^{j+k+1}},\]

where \(\odot\) is the binomial product.

public theorem_wanted binomialProduct_X_pow_mul_inv_pow {K : Type*} [Field K] [CharZero K]
    (j k : ℕ) (α β : K) :
    binomialProduct (X ^ j * ((1 - C α * X)⁻¹) ^ (j + 1))
        (X ^ k * ((1 - C β * X)⁻¹) ^ (k + 1)) =
      ((j + k).choose j : PowerSeries K) * X ^ (j + k) * ((1 - C (α + β) * X)⁻¹) ^ (j + k + 1)
Sources
  • Ira M. Gessel and Ishan Kar, Binomial Convolutions for Rational Power Series, Journal of Integer Sequences 27 (2024), Article 24.1.3, Lemma 4 (jis-2024-24.1.3-lem-4).
Uses 1

TheoremBinomial product with a geometric series#●ready to prove

For every field \(K\) of characteristic \(0\), every power series \(A(x)\) over \(K\) and every \(\beta \in K\),

\[A(x) \odot \frac{1}{1 - \beta x} = \frac{1}{1 - \beta x} A\left(\frac{x}{1 - \beta x}\right),\]

where \(\odot\) is the binomial product and \(A(x/(1 - \beta x))\) is the substitution of the series \(x/(1 - \beta x)\), which has zero constant term, into \(A\).

public theorem_wanted binomialProduct_inv_one_sub {K : Type*} [Field K] [CharZero K]
    (A : PowerSeries K) (β : K) :
    binomialProduct A (1 - C β * X)⁻¹ =
      (1 - C β * X)⁻¹ * PowerSeries.subst (X * (1 - C β * X)⁻¹) A
Sources
  • Ira M. Gessel and Ishan Kar, Binomial Convolutions for Rational Power Series, Journal of Integer Sequences 27 (2024), Article 24.1.3, Corollary 5 (jis-2024-24.1.3-cor-5).
Uses 1

TheoremBinomial products of proper rational series#●ready to prove

Let \(K\) be a field of characteristic \(0\), let \(\alpha_1, \dots, \alpha_m, \beta_1, \dots, \beta_n \in K\), not necessarily distinct or nonzero, and let \(R\), \(S\) be polynomials over \(K\) with \(\deg R < m\) and \(\deg S < n\). If \(A\), \(B\) are power series over \(K\) with \(A(x) \prod_{i=1}^m (1 - \alpha_i x) = R(x)\) and \(B(x) \prod_{j=1}^n (1 - \beta_j x) = S(x)\), then there is a polynomial \(T\) over \(K\) with \(\deg T < mn\) such that

\[(A \odot B)(x) \prod_{i=1}^m \prod_{j=1}^n \bigl(1 - (\alpha_i + \beta_j) x\bigr) = T(x),\]

where \(\odot\) is the binomial product.

Here \(\deg 0 = -\infty\), so \(T = 0\) when \(mn = 0\).

public theorem_wanted binomialProduct_rational {K : Type*} [Field K] [CharZero K] (m n : ℕ)
    (α : Fin m → K) (β : Fin n → K) (R S : Polynomial K) (hR : R.degree < m) (hS : S.degree < n)
    (A B : PowerSeries K) (hA : A * ∏ i, (1 - C (α i) * X) = (R : PowerSeries K))
    (hB : B * ∏ j, (1 - C (β j) * X) = (S : PowerSeries K)) :
    ∃ T : Polynomial K, T.degree < (m * n : ℕ) ∧
      binomialProduct A B * ∏ i, ∏ j, (1 - C (α i + β j) * X) = (T : PowerSeries K)
Sources
  • Ira M. Gessel and Ishan Kar, Binomial Convolutions for Rational Power Series, Journal of Integer Sequences 27 (2024), Article 24.1.3, Theorem 6 (jis-2024-24.1.3-thm-6).
Uses 1

TheoremBinomial products of arbitrary rational series#●ready to prove

Let \(K\) be a field of characteristic \(0\), let \(\alpha_1, \dots, \alpha_m\) and \(\beta_1, \dots, \beta_n\) be nonzero elements of \(K\), and let \(R\), \(S\) be polynomials over \(K\). Let \(A\), \(B\) be power series over \(K\) with \(A(x) \prod_{i=1}^m (1 - \alpha_i x) = R(x)\) and \(B(x) \prod_{j=1}^n (1 - \beta_j x) = S(x)\). Put \(u = \max(d_R + 1 - m, 0)\) and \(v = \max(d_S + 1 - n, 0)\), where \(d_R\) and \(d_S\) are the degrees of \(R\) and \(S\), the zero polynomial counting as degree \(0\). Then there is a polynomial \(T\) over \(K\) of degree less than \((u + m)(v + n)\) such that

\[(A \odot B)(x) \Big(\prod_{i=1}^m (1 - \alpha_i x)\Big)^{v} \Big(\prod_{j=1}^n (1 - \beta_j x)\Big)^{u} \prod_{i=1}^m \prod_{j=1}^n (1 - (\alpha_i + \beta_j) x) = T(x),\]

where \(\odot\) is the binomial product.

public theorem_wanted binomialProduct_rational_general {K : Type*} [Field K] [CharZero K]
    (m n : ℕ) (α : Fin m → K) (β : Fin n → K) (hα : ∀ i, α i ≠ 0) (hβ : ∀ j, β j ≠ 0)
    (R S : Polynomial K) (A B : PowerSeries K) (hA : A * ∏ i, (1 - C (α i) * X) = (R : PowerSeries K))
    (hB : B * ∏ j, (1 - C (β j) * X) = (S : PowerSeries K)) :
    ∃ T : Polynomial K, T.degree < (((R.natDegree + 1 - m) + m) * ((S.natDegree + 1 - n) + n) : ℕ) ∧
      binomialProduct A B * ((∏ i, (1 - C (α i) * X)) ^ (S.natDegree + 1 - n) *
        (∏ j, (1 - C (β j) * X)) ^ (R.natDegree + 1 - m) *
        ∏ i, ∏ j, (1 - C (α i + β j) * X)) = (T : PowerSeries K)
Sources
  • Ira M. Gessel and Ishan Kar, Binomial Convolutions for Rational Power Series, Journal of Integer Sequences 27 (2024), Article 24.1.3, Corollary 7 (jis-2024-24.1.3-cor-7).
Uses 1

TheoremBinomial product with a monomial#●ready to prove

For every field \(K\) of characteristic \(0\), every \(m \in \mathbb{N}\) and every power series \(A(x)\) over \(K\),

\[x^m \odot A(x) = \frac{x^m}{m!} \frac{d^m}{dx^m}\bigl(x^m A(x)\bigr),\]

where \(\odot\) is the binomial product and \(d/dx\) is the formal derivative.

public theorem_wanted binomialProduct_X_pow {K : Type*} [Field K] [CharZero K] (m : ℕ)
    (A : PowerSeries K) :
    binomialProduct (X ^ m) A =
      C ((m.factorial : K)⁻¹) * X ^ m * (⇑(PowerSeries.derivative (R := K)))^[m] (X ^ m * A)
Sources
  • Ira M. Gessel and Ishan Kar, Binomial Convolutions for Rational Power Series, Journal of Integer Sequences 27 (2024), Article 24.1.3, Lemma 8 (jis-2024-24.1.3-lem-8).
Uses 1

TheoremResultant identity for pairwise sums#●ready to prove

Let \(K\) be a field of characteristic \(0\), let \(m, n \ge 0\) be integers, and let \(\alpha_1, \dots, \alpha_m\) and \(\beta_1, \dots, \beta_n\) be elements of \(K\). Regard \(f(y) = \prod_{i=1}^{m} (1 - y - \alpha_i x)\) and \(g(y) = \prod_{j=1}^{n} (y - \beta_j x)\) as polynomials in \(y\) with coefficients in \(K[x]\), and let \(\operatorname{Res}_{m,n}(f, g)\) be the determinant of their Sylvester matrix formed with the formal degrees \(m\) and \(n\). Then in \(K[x]\)

\[ \prod_{i=1}^{m} \prod_{j=1}^{n} \bigl(1 - (\alpha_i + \beta_j)x\bigr) = (-1)^{mn} \operatorname{Res}_{m,n}(f, g). \]

With \(U(x) = \prod_{i=1}^{m} (1 - \alpha_i x)\) and \(V(x) = \prod_{j=1}^{n} (1 - \beta_j x)\), these are \(f(y) = (1 - y)^m U(x / (1 - y))\) and \(g(y) = y^n V(x / y)\), the form used in the source.

public theorem_wanted prod_one_sub_add_eq_resultant {K : Type*} [Field K] [CharZero K] (m n : ℕ)
    (α : Fin m → K) (β : Fin n → K) :
    ∏ i, ∏ j, (1 - Polynomial.C (α i + β j) * Polynomial.X) =
      (-1) ^ (m * n) * Polynomial.resultant
        (∏ i, (1 - Polynomial.X - Polynomial.C (Polynomial.C (α i) * Polynomial.X)))
        (∏ j, (Polynomial.X - Polynomial.C (Polynomial.C (β j) * Polynomial.X))) m n
Sources
  • Ira M. Gessel and Ishan Kar, Binomial Convolutions for Rational Power Series, Journal of Integer Sequences 27 (2024), Article 24.1.3, Theorem 10 (jis-2024-24.1.3-thm-10).

TheoremBinomial square of the Lucas generating function#●ready to prove

Let \(K\) be a field of characteristic \(0\), let \(a, b \in K\), and let \(g(x) = \sum_{n \ge 0} G_n x^n\) be a power series over \(K\) with \(g(x)(1 - ax - bx^2) = 2 - ax\). Then

\[g \odot g = \frac{2}{1 - ax} + g(2x),\]

where \(\odot\) is the binomial product, and \(\sum_{k=0}^{n} \binom{n}{k} G_k G_{n-k} = 2a^n + 2^n G_n\) for every \(n \ge 0\).

public theorem_wanted binomialProduct_lucasV_generatingFunction {K : Type*} [Field K]
    [CharZero K] (a b : K) (g : PowerSeries K)
    (hg : g * (1 - C a * X - C b * X ^ 2) = 2 - C a * X) :
    binomialProduct g g = 2 * (1 - C a * X)⁻¹ + PowerSeries.rescale 2 g ∧
      ∀ n : ℕ, ∑ k ∈ Finset.range (n + 1),
          (n.choose k : K) * PowerSeries.coeff k g * PowerSeries.coeff (n - k) g =
        2 * a ^ n + 2 ^ n * PowerSeries.coeff n g
Sources
  • Ira M. Gessel and Ishan Kar, Binomial Convolutions for Rational Power Series, Journal of Integer Sequences 27 (2024), Article 24.1.3, Theorem 11 (jis-2024-24.1.3-thm-11).
Uses 1

TheoremBinomial products of series with a cubic denominator#●ready to prove

Let \(K\) be a field of characteristic \(0\), let \(A, B, C \in K\) with \(C \ne 0\), and let \(D(x) = 1 + Ax + Bx^2 + Cx^3\). Suppose that \(D(x) = \prod_{i=1}^{3} (1 - \alpha_i x)\) with \(\alpha_1, \alpha_2, \alpha_3 \in K\) and that \(2\alpha_i \ne \alpha_j + \alpha_k\) whenever \(i, j, k \in \{1, 2, 3\}\) are not all equal. Let \(r\) and \(s\) be power series over \(K\) such that \(r D\) and \(s D\) are polynomials of degree less than \(3\). Then there are polynomials \(u\), \(v\) over \(K\) of degree at most \(2\) such that

\[r \odot s = \frac{u(x)}{D(2x)} + \frac{1}{1 + Ax} \odot \frac{v(x)}{D(-x)},\]

where \(\odot\) is the binomial product.

The second term uses the factor \(1/(1 + Ax)\); with \(1/(1 - Ax)\) in its place the statement already fails for the tribonacci series \(x/(1 - x - x^2 - x^3)\), where \(A = -1\).

public theorem_wanted binomialProduct_cubic_denominator {K : Type*} [Field K] [CharZero K]
    (A B C' : K) (hC : C' ≠ 0) (α : Fin 3 → K)
    (hD : 1 + C A * X + C B * X ^ 2 + C C' * X ^ 3 = ∏ i, (1 - C (α i) * X))
    (hα : ∀ i j k, ¬(i = j ∧ j = k) → 2 * α i ≠ α j + α k)
    (r s : PowerSeries K) (pr ps : Polynomial K) (hpr : pr.degree < 3) (hps : ps.degree < 3)
    (hr : r * (1 + C A * X + C B * X ^ 2 + C C' * X ^ 3) = (pr : PowerSeries K))
    (hs : s * (1 + C A * X + C B * X ^ 2 + C C' * X ^ 3) = (ps : PowerSeries K)) :
    ∃ u v : Polynomial K, u.natDegree ≤ 2 ∧ v.natDegree ≤ 2 ∧
      binomialProduct r s =
        (u : PowerSeries K) * (1 + C (2 * A) * X + C (4 * B) * X ^ 2 + C (8 * C') * X ^ 3)⁻¹ +
          binomialProduct (1 + C A * X)⁻¹
            ((v : PowerSeries K) * (1 - C A * X + C B * X ^ 2 - C C' * X ^ 3)⁻¹)
Sources
  • Ira M. Gessel and Ishan Kar, Binomial Convolutions for Rational Power Series, Journal of Integer Sequences 27 (2024), Article 24.1.3, Theorem 12 (jis-2024-24.1.3-thm-12).
Uses 1

TheoremHadamard products of rational series#●ready to prove

Let \(K\) be a field of characteristic \(0\), let \(\alpha_1, \dots, \alpha_m, \beta_1, \dots, \beta_n \in K\), and let \(A = \sum_k a_k x^k\) and \(B = \sum_k b_k x^k\) be power series over \(K\) such that \(A(x) \prod_{i=1}^m (1 - \alpha_i x)\) and \(B(x) \prod_{j=1}^n (1 - \beta_j x)\) are polynomials. Then their Hadamard product \(A * B = \sum_k a_k b_k x^k\) has the property that \((A * B)(x) \prod_{i=1}^m \prod_{j=1}^n (1 - \alpha_i \beta_j x)\) is a polynomial. Moreover, call a power series \(f\) over \(K\) proper if \(f Q = P\) for some polynomials \(P\), \(Q\) over \(K\) with \(Q(0) \ne 0\) and \(\deg P < \deg Q\). If \(A\) and \(B\) are proper, then so is \(A * B\).

public theorem_wanted hadamardProduct_rational {K : Type*} [Field K] [CharZero K] (m n : ℕ)
    (α : Fin m → K) (β : Fin n → K) (R S : Polynomial K) (A B : PowerSeries K)
    (hA : A * ∏ i, (1 - C (α i) * X) = (R : PowerSeries K)) (hB : B * ∏ j, (1 - C (β j) * X) = (S : PowerSeries K)) :
    (∃ T : Polynomial K, PowerSeries.mk (fun k => PowerSeries.coeff k A * PowerSeries.coeff k B) *
        ∏ i, ∏ j, (1 - C (α i * β j) * X) = (T : PowerSeries K)) ∧
      ((∃ P Q : Polynomial K, Q.coeff 0 ≠ 0 ∧ P.degree < Q.degree ∧ A * (Q : PowerSeries K) = (P : PowerSeries K)) →
        (∃ P Q : Polynomial K, Q.coeff 0 ≠ 0 ∧ P.degree < Q.degree ∧
          B * (Q : PowerSeries K) = (P : PowerSeries K)) →
        ∃ P Q : Polynomial K, Q.coeff 0 ≠ 0 ∧ P.degree < Q.degree ∧
          PowerSeries.mk (fun k => PowerSeries.coeff k A * PowerSeries.coeff k B) *
            (Q : PowerSeries K) = (P : PowerSeries K))
Sources
  • Ira M. Gessel and Ishan Kar, Binomial Convolutions for Rational Power Series, Journal of Integer Sequences 27 (2024), Article 24.1.3, Theorem 13 (jis-2024-24.1.3-thm-13).

TheoremResultant identity for pairwise products#●ready to prove

Let \(K\) be a field of characteristic \(0\), let \(m, n \ge 0\) be integers, and let \(\alpha_1, \dots, \alpha_m\) and \(\beta_1, \dots, \beta_n\) be elements of \(K\). Regard \(f(y) = \prod_{i=1}^{m} (1 - \alpha_i y)\) and \(g(y) = \prod_{j=1}^{n} (y - \beta_j x)\) as polynomials in \(y\) with coefficients in \(K[x]\), and let \(\operatorname{Res}_{m,n}(f, g)\) be the determinant of their Sylvester matrix formed with the formal degrees \(m\) and \(n\). Then in \(K[x]\)

\[ \prod_{i=1}^{m} \prod_{j=1}^{n} (1 - \alpha_i \beta_j x) = (-1)^{mn} \operatorname{Res}_{m,n}(f, g). \]

With \(U(x) = \prod_{i=1}^{m} (1 - \alpha_i x)\) and \(V(x) = \prod_{j=1}^{n} (1 - \beta_j x)\), these are \(f(y) = U(y)\) and \(g(y) = y^n V(x / y)\), the form used in the source.

public theorem_wanted prod_one_sub_mul_eq_resultant {K : Type*} [Field K] [CharZero K] (m n : ℕ)
    (α : Fin m → K) (β : Fin n → K) :
    ∏ i, ∏ j, (1 - Polynomial.C (α i * β j) * Polynomial.X) =
      (-1) ^ (m * n) * Polynomial.resultant
        (∏ i, (1 - Polynomial.C (Polynomial.C (α i)) * Polynomial.X))
        (∏ j, (Polynomial.X - Polynomial.C (Polynomial.C (β j) * Polynomial.X))) m n
Sources
  • Ira M. Gessel and Ishan Kar, Binomial Convolutions for Rational Power Series, Journal of Integer Sequences 27 (2024), Article 24.1.3, Theorem 14 (jis-2024-24.1.3-thm-14).
3 fully proved12 ready to prove

The dependency explorer needs JavaScript and its graph data.

Download the complete graph data.