Generalized Eulerian polynomials and gamma-positivity
DefinitionNumber of double excedances#✓fully proved
For a permutation \(w\) of \(\{0, 1, \dots, n-1\}\), a double excedance of \(w\) is a point \(i\) with \(w^{-1}(i) < i < w(i)\), and \(\mathrm{dexc}(w)\) is the number of double excedances of \(w\).
def doubleExcedanceCount {n : ℕ} (w : Equiv.Perm (Fin n)) : ℕ :=
(Finset.univ.filter fun i ↦ w⁻¹ i < i ∧ i < w i).card
Sources
- Pasquale Petrullo and Domenico Senato, Generalized Eulerian Polynomials with a Nonnegative Gamma Vector, Journal of Integer Sequences 27 (2024), Article 24.1.5, Section 4.
TheoremThe identity has no double excedances#✓fully proved
For every \(n\), the identity permutation of \(\{0, 1, \dots, n-1\}\) has no double excedances.
theorem doubleExcedanceCount_one (n : ℕ) :
doubleExcedanceCount (1 : Equiv.Perm (Fin n)) = 0
def doubleExcedanceCount {n : ℕ} (w : Equiv.Perm (Fin n)) : ℕ :=
(Finset.univ.filter fun i ↦ w⁻¹ i < i ∧ i < w i).card
WantedExt/Combinatorics/Enumerative/GeneralizedEulerianGammaWanted.lean:46
Sources
- Pasquale Petrullo and Domenico Senato, Generalized Eulerian Polynomials with a Nonnegative Gamma Vector, Journal of Integer Sequences 27 (2024), Article 24.1.5, Section 4.
Uses 1
DefinitionGenerating function of the generalized Eulerian polynomials#✓fully proved
For rational numbers \(q\), \(t\), \(x\), \(y\) with \(x \ne y\), this is the formal power series in \(z\) over \(\mathbb{Q}\)
Here the base has constant term \(1\), and its \(q\)-th power is the binomial series \((1 + u)^q = \sum_{n \ge 0} \binom{q}{n} u^n\) with \(u\) the base minus \(1\). For \(x = y\) the denominator is zero and the formula above does not apply.
noncomputable def generalizedEulerianSeries (q t x y : ℚ) : PowerSeries ℚ :=
(binomialSeries ℚ q).subst
(C (x - y) * (C x * rescale y (exp ℚ) - C y * rescale x (exp ℚ))⁻¹ - 1) *
rescale ((x + y) * q * t) (exp ℚ)
Sources
- Pasquale Petrullo and Domenico Senato, Generalized Eulerian Polynomials with a Nonnegative Gamma Vector, Journal of Integer Sequences 27 (2024), Article 24.1.5, Equation (2).
DefinitionGeneralized Eulerian polynomials#✓fully proved
For rational numbers \(q, t, x, y\) and \(n \ge 0\), the generalized Eulerian polynomial \(P_n(q, t; x, y)\) is \(n!\) times the coefficient of \(z^n\) in the formal power series over \(\mathbb{Q}\)
where the \(q\)-th power is the binomial series \((1 + u)^q = \sum_{k \ge 0} \binom{q}{k} u^k\) applied to \(u = \frac{x - y}{x e^{yz} - y e^{xz}} - 1\). When \(x \ne y\) the quotient has constant term \(1\), so \(u\) has constant term \(0\); for \(x = y\) the formula is degenerate and the value of \(P_n\) is a formal convention.
noncomputable def generalizedEulerianP (q t x y : ℚ) (n : ℕ) : ℚ :=
(n.factorial : ℚ) * coeff n (generalizedEulerianSeries q t x y)
noncomputable def generalizedEulerianSeries (q t x y : ℚ) : PowerSeries ℚ :=
(binomialSeries ℚ q).subst
(C (x - y) * (C x * rescale y (exp ℚ) - C y * rescale x (exp ℚ))⁻¹ - 1) *
rescale ((x + y) * q * t) (exp ℚ)
WantedExt/Combinatorics/Enumerative/GeneralizedEulerianGammaWanted.lean:64
Sources
- Pasquale Petrullo and Domenico Senato, Generalized Eulerian Polynomials with a Nonnegative Gamma Vector, Journal of Integer Sequences 27 (2024), Article 24.1.5, Equation (3).
Uses 1 · Used by 6
TheoremGeneralized Eulerian polynomials with zero cycle weight#✓fully proved
For all rational numbers \(t, x, y\) and every \(n \ge 0\), the generalized Eulerian polynomial satisfies \(P_0(0, t; x, y) = 1\) and \(P_n(0, t; x, y) = 0\) for \(n \ge 1\).
Indeed, for \(q = 0\) the generating function of the polynomials \(P_n(0, t; x, y)\) is \(1\).
theorem generalizedEulerianP_zero_left (t x y : ℚ) (n : ℕ) :
generalizedEulerianP 0 t x y n = if n = 0 then 1 else 0
noncomputable def generalizedEulerianP (q t x y : ℚ) (n : ℕ) : ℚ :=
(n.factorial : ℚ) * coeff n (generalizedEulerianSeries q t x y)
WantedExt/Combinatorics/Enumerative/GeneralizedEulerianGammaWanted.lean:75
Sources
- Pasquale Petrullo and Domenico Senato, Generalized Eulerian Polynomials with a Nonnegative Gamma Vector, Journal of Integer Sequences 27 (2024), Article 24.1.5, Equation (2).
Uses 2
TheoremExponential generating function of excedances over cyclic permutations#●ready to prove
Let \(x \ne 1\) be a rational number. For \(n \ge 1\) let \(C_n(x) = \sum_{m=0}^{n} C_{n,m} x^m\), where the circular Eulerian number \(C_{n,m}\) is the number of permutations \(w\) of \(\{0, 1, \dots, n-1\}\) whose points all lie in a single cycle and which have exactly \(m\) excedances, that is, points \(i\) with \(w(i) > i\). Then, as formal power series in \(z\) over \(\mathbb{Q}\),
where the quotient has constant term \(1\) and \(\log\) is the formal logarithm of a power series with constant term \(1\).
public theorem_wanted cyclicExcedanceEGF_eq (x : ℚ) (hx : x ≠ 1) :
PowerSeries.mk (fun n ↦ if n = 0 then 0 else
(∑ m ∈ Finset.range (n + 1), (circularEulerian n m : ℚ) * x ^ m) / (n.factorial : ℚ)) =
X + logOf (C (x - 1) * (C x * exp ℚ - rescale x (exp ℚ))⁻¹)
def circularEulerian (n m : ℕ) : ℕ :=
(Finset.univ.filter (fun σ : Equiv.Perm (Fin n) =>
(0 < n ∧ ∀ i j : Fin n, σ.SameCycle i j) ∧
excedanceCount n σ = m)).card
MathlibExt/Combinatorics/AssociatedStirlingEulerian.lean:92
Sources
- Pasquale Petrullo and Domenico Senato, Generalized Eulerian Polynomials with a Nonnegative Gamma Vector, Journal of Integer Sequences 27 (2024), Article 24.1.5, Proposition 1 (
jis-2024-24.1.5-prop-1).
Uses 1
TheoremPermutation expansion of the generalized Eulerian polynomials#●ready to prove
Let \(P_n(q, t; x, y)\) be the generalized Eulerian polynomial and let \(S_n\) be the set of permutations of \(\{0, 1, \dots, n-1\}\). For \(w \in S_n\) write \(\mathrm{cyc}(w)\) for the number of cycles of \(w\) (fixed points counted as cycles), \(\mathrm{fix}(w)\) for the number of fixed points of \(w\), and \(\mathrm{exc}(w)\) for the number of excedances of \(w\), the points \(i\) with \(w(i) > i\). Then for all rational numbers \(q, t, x, y\) with \(x \ne y\) and every \(n \ge 0\),
public theorem_wanted generalizedEulerianP_eq_sum_perm (q t x y : ℚ) (hxy : x ≠ y) (n : ℕ) :
generalizedEulerianP q t x y n =
∑ w : Equiv.Perm (Fin n),
q ^ cycleQuotientCount w * t ^ Fintype.card (Function.fixedPoints w) *
x ^ excedanceCount n w * y ^ excedanceCount n w⁻¹ *
(x + y) ^ Fintype.card (Function.fixedPoints w)
noncomputable def generalizedEulerianP (q t x y : ℚ) (n : ℕ) : ℚ :=
(n.factorial : ℚ) * coeff n (generalizedEulerianSeries q t x y)
WantedExt/Combinatorics/Enumerative/GeneralizedEulerianGammaWanted.lean:75
public noncomputable def cycleQuotientCount {α : Type*} [Fintype α]
(σ : Equiv.Perm α) : ℕ :=
Nat.card (Quotient (Equiv.Perm.SameCycle.setoid σ))
MathlibExt/Combinatorics/Enumerative/Stirling.lean:27
def excedanceCount (n : ℕ) (σ : Equiv.Perm (Fin n)) : ℕ :=
(Finset.univ.filter (fun i : Fin n => i < σ i)).card
MathlibExt/Combinatorics/AssociatedStirlingEulerian.lean:26
Sources
- Pasquale Petrullo and Domenico Senato, Generalized Eulerian Polynomials with a Nonnegative Gamma Vector, Journal of Integer Sequences 27 (2024), Article 24.1.5, Theorem 2 (
jis-2024-24.1.5-thm-2).
Uses 3
TheoremThe gamma coefficients of the generalized Eulerian polynomials have nonnegative integer coefficients#●ready to prove
The gamma coefficients of the generalized Eulerian polynomials have nonnegative integer coefficients. Fix \(n \ge 0\) and let \(P_n(q, t; x, y)\) be the generalized Eulerian polynomial. There are natural numbers \(c_{k,i,j}\), for \(k, i, j \ge 0\), such that, writing \(\gamma_{n,k}(q, t) = \sum_{i=0}^{n} \sum_{j=0}^{n} c_{k,i,j} q^i t^j\), for all rational numbers \(q, t, x, y\) with \(x \ne y\),
That is, the gamma coefficients \(\gamma_{n,k}\) are polynomials in \(q\) and \(t\) with nonnegative integer coefficients.
public theorem_wanted generalizedEulerianP_gamma_nonneg (n : ℕ) :
∃ c : ℕ → ℕ → ℕ → ℕ, ∀ q t x y : ℚ, x ≠ y →
generalizedEulerianP q t x y n =
∑ k ∈ Finset.range (n / 2 + 1),
(∑ i ∈ Finset.range (n + 1), ∑ j ∈ Finset.range (n + 1),
(c k i j : ℚ) * q ^ i * t ^ j) * x ^ k * y ^ k * (x + y) ^ (n - 2 * k)
noncomputable def generalizedEulerianP (q t x y : ℚ) (n : ℕ) : ℚ :=
(n.factorial : ℚ) * coeff n (generalizedEulerianSeries q t x y)
WantedExt/Combinatorics/Enumerative/GeneralizedEulerianGammaWanted.lean:75
Sources
- Pasquale Petrullo and Domenico Senato, Generalized Eulerian Polynomials with a Nonnegative Gamma Vector, Journal of Integer Sequences 27 (2024), Article 24.1.5, Theorem 4 (
jis-2024-24.1.5-thm-4).
Uses 1
TheoremUnimodality of the generalized Eulerian polynomials#●ready to prove
Let \(P_n(q, t; x, y)\) be the generalized Eulerian polynomial. For all natural numbers \(q\), \(t\) and \(n\) there are a polynomial \(p\) with rational coefficients \(p_0, p_1, p_2, \dots\) and an index \(m\) such that \(p(x) = P_n(q, t; x, 1)\) for every rational \(x \ne 1\), \(p_i \le p_{i+1}\) for every \(i < m\), and \(p_{i+1} \le p_i\) for every \(i \ge m\). That is, \(P_n(q, t; x, 1)\) is a unimodal polynomial in \(x\).
Only unimodality is asserted, not symmetry of the coefficients.
public theorem_wanted generalizedEulerianP_unimodal (q t n : ℕ) :
∃ p : Polynomial ℚ, (∀ x : ℚ, x ≠ 1 → p.eval x = generalizedEulerianP q t x 1 n) ∧
∃ m : ℕ, (∀ i < m, p.coeff i ≤ p.coeff (i + 1)) ∧
∀ i, m ≤ i → p.coeff (i + 1) ≤ p.coeff i
noncomputable def generalizedEulerianP (q t x y : ℚ) (n : ℕ) : ℚ :=
(n.factorial : ℚ) * coeff n (generalizedEulerianSeries q t x y)
WantedExt/Combinatorics/Enumerative/GeneralizedEulerianGammaWanted.lean:75
Sources
- Pasquale Petrullo and Domenico Senato, Generalized Eulerian Polynomials with a Nonnegative Gamma Vector, Journal of Integer Sequences 27 (2024), Article 24.1.5, Corollary 5 (
jis-2024-24.1.5-cor-5).
Uses 1
TheoremThe gamma coefficients of the generalized Eulerian polynomials count permutations#●ready to prove
Fix \(n \ge 0\) and let \(P_n(q, t; x, y)\) be the generalized Eulerian polynomial. For a permutation \(w\) of \(\{0, 1, \dots, n-1\}\) write \(\mathrm{cyc}(w)\) for its number of cycles (fixed points counted as cycles), \(\mathrm{fix}(w)\) for its number of fixed points, \(\mathrm{exc}(w)\) for its number of excedances (points \(i\) with \(w(i) > i\)) and \(\mathrm{dexc}(w)\) for its number of double excedances. For \(k \ge 0\) and rational \(q, t\) let \(\gamma_{n,k}(q, t)\) be the sum of \(q^{\mathrm{cyc}(w)} t^{\mathrm{fix}(w)}\) over the permutations \(w\) of \(\{0, 1, \dots, n-1\}\) with \(\mathrm{exc}(w) = k\) and \(\mathrm{dexc}(w) = 0\). Then for all rational \(q, t, x, y\) with \(x \ne y\),
and for every \(k \ge 0\) and all rational \(q, t\), \(\gamma_{n,k}(q, t) = \sum_{i=0}^{n} \sum_{j=0}^{n} \gamma_{n,k,i,j} q^i t^j\), where \(\gamma_{n,k,i,j}\) is the number of permutations \(w\) of \(\{0, 1, \dots, n-1\}\) with \(\mathrm{exc}(w) = k\), \(\mathrm{dexc}(w) = 0\), \(\mathrm{cyc}(w) = i\) and \(\mathrm{fix}(w) = j\).
MetaMathlibExt.generalizedEulerianP_gamma_eq · 17 lines
public theorem_wanted generalizedEulerianP_gamma_eq (n : ℕ) :
(∀ q t x y : ℚ, x ≠ y →
generalizedEulerianP q t x y n =
∑ k ∈ Finset.range (n / 2 + 1),
(∑ w ∈ Finset.univ.filter (fun w : Equiv.Perm (Fin n) ↦
excedanceCount n w = k ∧ doubleExcedanceCount w = 0),
q ^ cycleQuotientCount w * t ^ Fintype.card (Function.fixedPoints w)) *
x ^ k * y ^ k * (x + y) ^ (n - 2 * k)) ∧
∀ (k : ℕ) (q t : ℚ),
∑ w ∈ Finset.univ.filter (fun w : Equiv.Perm (Fin n) ↦
excedanceCount n w = k ∧ doubleExcedanceCount w = 0),
q ^ cycleQuotientCount w * t ^ Fintype.card (Function.fixedPoints w) =
∑ i ∈ Finset.range (n + 1), ∑ j ∈ Finset.range (n + 1),
((Finset.univ.filter (fun w : Equiv.Perm (Fin n) ↦ excedanceCount n w = k ∧
doubleExcedanceCount w = 0 ∧
cycleQuotientCount w = i ∧ Fintype.card (Function.fixedPoints w) = j)).card : ℚ) *
q ^ i * t ^ j
noncomputable def generalizedEulerianP (q t x y : ℚ) (n : ℕ) : ℚ :=
(n.factorial : ℚ) * coeff n (generalizedEulerianSeries q t x y)
WantedExt/Combinatorics/Enumerative/GeneralizedEulerianGammaWanted.lean:75
def excedanceCount (n : ℕ) (σ : Equiv.Perm (Fin n)) : ℕ :=
(Finset.univ.filter (fun i : Fin n => i < σ i)).card
MathlibExt/Combinatorics/AssociatedStirlingEulerian.lean:26
def doubleExcedanceCount {n : ℕ} (w : Equiv.Perm (Fin n)) : ℕ :=
(Finset.univ.filter fun i ↦ w⁻¹ i < i ∧ i < w i).card
WantedExt/Combinatorics/Enumerative/GeneralizedEulerianGammaWanted.lean:46
public noncomputable def cycleQuotientCount {α : Type*} [Fintype α]
(σ : Equiv.Perm α) : ℕ :=
Nat.card (Quotient (Equiv.Perm.SameCycle.setoid σ))
MathlibExt/Combinatorics/Enumerative/Stirling.lean:27
Sources
- Pasquale Petrullo and Domenico Senato, Generalized Eulerian Polynomials with a Nonnegative Gamma Vector, Journal of Integer Sequences 27 (2024), Article 24.1.5, Theorem 7 (
jis-2024-24.1.5-thm-7).
Uses 4
TheoremGeneralized Eulerian polynomials at \(x = -1\), \(y = 1\)#●ready to prove
Let \(P_m(q, t; x, y)\) be the generalized Eulerian polynomial. For a permutation \(w\) write \(\mathrm{cyc}(w)\) for its number of cycles (fixed points counted as cycles), \(\mathrm{exc}(w)\) for its number of excedances (points \(i\) with \(w(i) > i\)) and \(\mathrm{dexc}(w)\) for its number of double excedances, and let \(D_{2n}\) be the set of derangements (permutations without fixed points) of \(\{0, 1, \dots, 2n-1\}\). Then for all rational numbers \(q, t\) and every \(n \ge 0\), \(P_{2n+1}(q, t; -1, 1) = 0\) and
where the exponent \(n - \mathrm{exc}(w)\) is an integer, possibly negative.
public theorem_wanted generalizedEulerianP_neg_one (q t : ℚ) (n : ℕ) :
generalizedEulerianP q t (-1) 1 (2 * n + 1) = 0 ∧
(-1) ^ n * generalizedEulerianP q t (-1) 1 (2 * n) =
∑ w ∈ Finset.univ.filter (fun w : Equiv.Perm (Fin (2 * n)) ↦
Fintype.card (Function.fixedPoints w) = 0),
q ^ cycleQuotientCount w * (-1 : ℚ) ^ ((n : ℤ) - excedanceCount (2 * n) w) ∧
∑ w ∈ Finset.univ.filter (fun w : Equiv.Perm (Fin (2 * n)) ↦
Fintype.card (Function.fixedPoints w) = 0),
q ^ cycleQuotientCount w * (-1 : ℚ) ^ ((n : ℤ) - excedanceCount (2 * n) w) =
∑ w ∈ Finset.univ.filter (fun w : Equiv.Perm (Fin (2 * n)) ↦
Fintype.card (Function.fixedPoints w) = 0 ∧ excedanceCount (2 * n) w = n ∧
doubleExcedanceCount w = 0), q ^ cycleQuotientCount w
noncomputable def generalizedEulerianP (q t x y : ℚ) (n : ℕ) : ℚ :=
(n.factorial : ℚ) * coeff n (generalizedEulerianSeries q t x y)
WantedExt/Combinatorics/Enumerative/GeneralizedEulerianGammaWanted.lean:75
public noncomputable def cycleQuotientCount {α : Type*} [Fintype α]
(σ : Equiv.Perm α) : ℕ :=
Nat.card (Quotient (Equiv.Perm.SameCycle.setoid σ))
MathlibExt/Combinatorics/Enumerative/Stirling.lean:27
def excedanceCount (n : ℕ) (σ : Equiv.Perm (Fin n)) : ℕ :=
(Finset.univ.filter (fun i : Fin n => i < σ i)).card
MathlibExt/Combinatorics/AssociatedStirlingEulerian.lean:26
def doubleExcedanceCount {n : ℕ} (w : Equiv.Perm (Fin n)) : ℕ :=
(Finset.univ.filter fun i ↦ w⁻¹ i < i ∧ i < w i).card
WantedExt/Combinatorics/Enumerative/GeneralizedEulerianGammaWanted.lean:46
Sources
- Pasquale Petrullo and Domenico Senato, Generalized Eulerian Polynomials with a Nonnegative Gamma Vector, Journal of Integer Sequences 27 (2024), Article 24.1.5, Corollary 8 (
jis-2024-24.1.5-cor-8).
Uses 4
Loading dependency explorer…
The dependency explorer needs JavaScript and its graph data.