Angel Ivanov Raychev

A Generalization of Descent Polynomials

Angel Raychev

Research paperPublished in Discrete Mathematics · 2023arXiv:2109.04519Journal

PDF: September 2021 author preprint; not the publisher's typeset version.

Sole-author research in enumerative combinatorics. Written in 2021; published in Discrete Mathematics in 2023.

A permutation records an order. Its descent set records exactly where that order falls: the positions at which an entry is larger than the next one. Fix those positions and count all permutations that realize them. Although the permutations become longer as the alphabet grows, their number is eventually governed by a polynomial.

In this paper, I extended that counting problem to multiset permutations, where every label appears several times. The resulting family has two parameters: the number of distinct labels and their common multiplicity. I proved polynomiality in the first parameter, found the sharp point at which increasing the second stops changing the answer, and obtained an explicit formula for the stabilized polynomial. Its coefficients have a direct counting interpretation in a particular binomial basis, together with a sharp description of which shifts preserve nonnegativity.

The useful structural idea is that a prescribed descent set leaves only a bounded initial segment free. Everything after the last descent has to be sorted. That converts a count of increasingly long words into a finite collection of order patterns.

Fixing the positions of the descents

For a word w=(w1,…,wN)w=(w_1,\ldots,w_N), define

Des⁡(w)={i∈{1,…,N−1}:wi>wi+1}.\operatorname{Des}(w)=\{i\in\{1,\ldots,N-1\}:w_i>w_{i+1}\}.

For example,

w=(4,1,3,2)⟹Des⁡(w)={1,3}.w=(4,1,3,2) \quad\Longrightarrow\quad \operatorname{Des}(w)=\{1,3\}.

The positions are what matter, not the values of the entries at those positions. Equality is not a descent: in a word with repeated letters, 2,22,2 is a weakly increasing pair.

The ordinary descent polynomial counts permutations of {1,…,n}\{1,\ldots,n\} with one specified descent set II. This is different from an Eulerian polynomial, which, for a fixed nn, packages permutations according to their number of descents. Here the pattern II stays fixed while nn varies.

My generalization replaces the set by the multiset

{1m,2m,…,nm},\{1^m,2^m,\ldots,n^m\},

where the notation means mm copies of each label, not exponentiation. Copies of the same label are indistinguishable. If Sn(m)S_n^{(m)} denotes its distinct permutations, the counting function is

dm(I,n)=#{w∈Sn(m):Des⁡(w)=I}.\mathfrak d^m(I,n) =\#\{w\in S_n^{(m)}:\operatorname{Des}(w)=I\}.

Setting m=1m=1 recovers the ordinary descent polynomial. The empty descent set always gives dm(∅,n)=1\mathfrak d^m(\varnothing,n)=1: there is exactly one weakly increasing arrangement.

A complete example

Take n=3n=3, m=2m=2 and I={2}I=\{2\}. We want arrangements of 1,1,2,2,3,31,1,2,2,3,3 with a descent at position 22 and nowhere else. There are exactly five:

The five multiset permutations with descent set {2}: 121233, 131223, 221133, 231123 and 331122. The first two entries determine each word; the remaining four are sorted.
The five arrangements with exactly the required descent pattern. The shaded prefix contains the only choices; each remaining suffix is forced into weakly increasing order.

This example already exposes the mechanism behind the general results. Once the first two entries are chosen, the other copies have only one allowable order. We must still check the boundary between the prefix and the suffix: a sorted suffix alone does not guarantee the required final descent.

Two parameters and two kinds of structure

Let

I={α1<⋯<αt}≠∅,a=αt.I=\{\alpha_1<\cdots<\alpha_t\}\ne\varnothing, \qquad a=\alpha_t.

Let LL be the length of the longest consecutive run in II. For example, I={2,3,5,7,10,11,12}I=\{2,3,5,7,10,11,12\} has L=3L=3 because of the run 10,11,1210,11,12.

The main results describe different directions in this two-parameter family.

Polynomiality in the alphabet size. For every fixed positive multiplicity mm, the function dm(I,n)\mathfrak d^m(I,n) agrees with a polynomial in nn for n≥an\ge a.

Stabilization in the multiplicity. For n>an>a, the counts are weakly increasing as mm grows. They stop changing once

m≥M=a−t+1.m\ge M=a-t+1.

This threshold is sharp when n≥L+1n\ge L+1: MM is the smallest positive multiplicity from which the count stays constant. If M=1M=1, the family is already stable at its first allowed multiplicity. The condition on nn matters: a run of LL consecutive descents needs L+1L+1 distinct labels; with fewer labels, every count is zero.

The two statements should not be conflated. One says how the count changes when more labels become available. The other says that, beyond a precise point, more copies of the same labels add no new possibilities.

Seeing stabilization in a formula

For the single descent set I={2}I=\{2\}, the threshold is M=2M=2. The two regimes have simple forms:

d1({2},n)=(n2)−1,dm({2},n)=(n+12)−1(m≥2).\begin{aligned} \mathfrak d^1(\{2\},n)&=\binom n2-1,\\ \mathfrak d^m(\{2\},n)&=\binom{n+1}{2}-1 \qquad(m\ge2). \end{aligned}

Here n≥2n\ge2. The first formula chooses two distinct prefix labels; the second allows a repeated label as well. In each case, exactly one prefix fails to produce the boundary descent.

Distinct labels nnOne copy m=1m=1Two copies m=2m=2Three copies m=3m=3
2022
3255
4599
591414
6142020

The identical last two columns illustrate a theorem about all m≥2m\ge2, not a conjecture inferred from the displayed cases.

Why longer words eventually stop adding choices

After position a=max⁡Ia=\max I, no descent is permitted. Consequently, the suffix must be weakly increasing, and the first aa entries determine the entire word.

At a descent position, the entry cannot be 11. There are tt such positions, so at most a−ta-t copies of 11 can occur in the prefix. Once m≥a−t+1m\ge a-t+1, a copy of 11 is guaranteed to remain for the sorted suffix. Its first entry is therefore 11.

The last required descent now has the particularly simple condition

wa>1.w_a>1.

Also, any one label can appear at most a−t+1a-t+1 times in an admissible prefix. To see this, cut the prefix at every weak increase, partitioning it into maximal strictly decreasing runs. Its t−1t-1 internal descents leave exactly a−t+1a-t+1 such runs, and a fixed label appears at most once in each. Thus the multiplicity bound no longer excludes any admissible prefix once m≥Mm\ge M.

The stabilized problem is therefore a count of finite words v=(v1,…,va)v=(v_1,\ldots,v_a) satisfying

Des⁡(v)=I∖{a},vj∈{1,…,n},va>1.\operatorname{Des}(v)=I\setminus\{a\}, \qquad v_j\in\{1,\ldots,n\}, \qquad v_a>1.

The growing sorted suffix has disappeared from the problem. To prove sharpness, the paper constructs an admissible prefix containing exactly MM copies of one label. That prefix cannot exist at multiplicity M−1M-1; the assumption n≥L+1n\ge L+1 supplies enough labels for its strict drops.

Write d∞(I,n)d^\infty(I,n) for the resulting stable polynomial. In the polynomial range n≥an\ge a, it agrees with dm(I,n)\mathfrak d^m(I,n) for every m≥Mm\ge M. The superscript ∞\infty names this finite stabilization; it does not require infinite words or a limiting approximation.

Two proofs of polynomiality

The paper gives both an elementary counting proof and a symmetric-function proof. They explain the same phenomenon through different representations.

Count the order patterns in the prefix

Work in the polynomial range n≥an\ge a. Temporarily omit the requirement that position aa be a descent. Let Dm(I,n)D^m(I,n) count prefixes of length aa with descent set I∖{a}I\setminus\{a\}, using each label at most mm times.

A prefix uses some number r≤ar\le a of distinct labels. Record their multiplicities in increasing label order:

A=(a1,…,ar),a1+⋯+ar=a,1≤aj≤m.A=(a_1,\ldots,a_r),\qquad a_1+\cdots+a_r=a,\qquad 1\le a_j\le m.

Let C(A,I)C(A,I) be the number of such prefixes after their labels are normalized to 1,…,r1,\ldots,r. Choosing the actual labels contributes (nr)\binom nr, independently of their order pattern. Hence

Dm(I,n)=∑A∈Comp⁡(a,m)C(A,I)(nℓ(A)),D^m(I,n)= \sum_{A\in\operatorname{Comp}(a,m)} C(A,I)\binom n{\ell(A)},

where Comp⁡(a,m)\operatorname{Comp}(a,m) is the finite set of compositions of aa with parts at most mm, and ℓ(A)\ell(A) is the number of parts. Every summand is a polynomial in nn, and all remaining data depend only on II and mm.

Completing the sorted suffix produces either the desired last descent or no descent at that boundary. Thus

Dm(I,n)=dm(I,n)+dm(I∖{a},n).D^m(I,n)= \mathfrak d^m(I,n)+ \mathfrak d^m(I\setminus\{a\},n).

Induction on the size of the descent set proves polynomiality. The proof turns a global arrangement problem into the choice of labels for finitely many normalized patterns.

Encode the words as ribbon tableaux

There is also a geometric representation. Break a word at its descents. Each weakly increasing run becomes a row of a ribbon-shaped tableau; strict drops become the strict inequalities in the overlapping columns. A ribbon is a connected skew Young diagram containing no 2×22\times2 square.

The word 1,2,1,2,3,3 is represented by a ribbon with a bottom row containing 1,2 and an upper row containing 1,2,3,3. The overlap places 2 below 1, recording the strict descent at position 2.
The same example as a ribbon tableau. Read the bottom row left to right, then the row above. Weak increases along rows and the strict drop at the shared column encode exactly the word's descent pattern.

For the appropriate shape λ/μ\lambda/\mu, the count becomes a coefficient of a skew Schur function:

dm(I,n)=[x1m⋯xnm]sλ/μ.\mathfrak d^m(I,n) =[x_1^m\cdots x_n^m]s_{\lambda/\mu}.

The Jacobi–Trudi identity expresses this Schur function as a determinant of complete homogeneous symmetric functions. The determinant’s dimension is fixed by II. As nn grows, only the ribbon’s long final row grows; each determinant term has one large factor hnm−uh_{nm-u} and finitely many factors of fixed total degree uu.

Extracting the equal-multiplicity coefficient again reduces to choosing labels for a fixed pattern. Each term is eventually polynomial in nn, so their finite signed sum is too. This proof connects the counting family to symmetric functions, rather than treating the polynomiality as an isolated numerical property.

An explicit formula for the stable polynomial

In the stable regime, no multiplicity constraint remains in the prefix. Its weakly increasing runs can be counted by stars and bars, and inclusion–exclusion enforces a strict drop at each prescribed internal boundary.

Define the gaps between the descent positions by

β1=α1,βj=αj−αj−1(2≤j≤t).\beta_1=\alpha_1, \qquad \beta_j=\alpha_j-\alpha_{j-1} \quad(2\le j\le t).

For a composition A=(a1,…,as)A=(a_1,\ldots,a_s) of tt, group consecutive gaps into ss blocks; let σj\sigma_j be the sum of the gaps in the jjth block. For example, A=(2,1)A=(2,1) groups (β1,β2,β3)(\beta_1,\beta_2,\beta_3) into (β1+β2,β3)(\beta_1+\beta_2,\beta_3).

Put

Hr(n)=(n+r−1r).H_r(n)=\binom{n+r-1}{r}.

This is the number of weakly increasing words of length rr on nn labels. The paper’s formula can then be written compactly as

d∞(I,n)=∑A=(a1,…,as)∈Comp⁡(t)(−1)t−s(∏j=1s−1Hσj(n))(Hσs(n)−1).\boxed{ d^\infty(I,n)= \sum_{A=(a_1,\ldots,a_s)\in\operatorname{Comp}(t)} (-1)^{t-s} \left(\prod_{j=1}^{s-1}H_{\sigma_j}(n)\right) \bigl(H_{\sigma_s}(n)-1\bigr). }

Merging adjacent runs removes a required descent, producing the inclusion–exclusion sign. The final subtraction removes the all-11 final block, the only weakly increasing block that would end at 11 rather than above it.

For I={1,3}I=\{1,3\}, the gap sequence is (1,2)(1,2). There are only two compositions of 22, so

d∞({1,3},n)=H1(n)(H2(n)−1)−(H3(n)−1)=n((n+12)−1)−((n+23)−1)=n3−4n+33.\begin{aligned} d^\infty(\{1,3\},n) &=H_1(n)\bigl(H_2(n)-1\bigr) -\bigl(H_3(n)-1\bigr)\\ &=n\left(\binom{n+1}{2}-1\right) -\left(\binom{n+2}{3}-1\right)\\ &=\frac{n^3-4n+3}{3}. \end{aligned}

This is a degree-33 polynomial arising from words whose lengths keep growing. The last descent position, rather than the eventual word length, governs that degree.

Coefficients that count normalized patterns

A polynomial can have positive or negative coefficients depending on its basis. For the stabilized family, the basis

(n−10),(n−11),(n−12),…\binom{n-1}{0},\binom{n-1}{1},\binom{n-1}{2},\ldots

has a particularly natural meaning. Write

d∞(I,n)=∑i=0abi(I)(n−1i).d^\infty(I,n)=\sum_{i=0}^{a}b_i(I)\binom{n-1}{i}.

The coefficient bi(I)b_i(I) counts normalized prefixes vv of length aa such that:

  • Des⁡(v)=I∖{a}\operatorname{Des}(v)=I\setminus\{a\};
  • all entries lie in {1,…,i+1}\{1,\ldots,i+1\};
  • each of 2,…,i+12,\ldots,i+1 occurs at least once;
  • the final entry is greater than 11.

The label 11 is optional. The other ii labels are mandatory, and choosing their actual values from {2,…,n}\{2,\ldots,n\} contributes the factor (n−1i)\binom{n-1}{i}. The count therefore separates into a choice of labels and a normalized pattern.

The coefficient theorem is more precise than nonnegativity:

bi(I)>0⟺L≤i≤a,b_i(I)>0\quad\Longleftrightarrow\quad L\le i\le a,

and every other coefficient is zero. In particular, d∞(I,n)d^\infty(I,n) has degree exactly aa.

For I={2}I=\{2\}, the two patterns counted by b1b_1 are (1,2)(1,2) and (2,2)(2,2); the one pattern counted by b2b_2 is (2,3)(2,3). Thus

d∞({2},n)=2(n−11)+(n−12).d^\infty(\{2\},n) =2\binom{n-1}{1}+\binom{n-1}{2}.

In the unshifted basis, the same polynomial is

d∞({2},n)=−(n0)+(n1)+(n2).d^\infty(\{2\},n) =-\binom n0+\binom n1+\binom n2.

The negative constant is not a defect in the count. It shows that the choice of basis can reveal or obscure its direct combinatorial structure.

The shift boundary is sharp

The paper classifies every integer shift. In the expansion in ((n+ki))i≥0\bigl(\binom{n+k}{i}\bigr)_{i\ge0}, all coefficients are integers, and:

  • if k≤−1k\le-1, every coefficient is nonnegative;
  • if k≥0k\ge0, at least one coefficient is negative.

For the unshifted basis, write its coefficients as ci(I)c_i(I). There is an exact initial segment:

ci(I)=(−1)i+t(0≤i≤L).c_i(I)=(-1)^{i+t} \qquad(0\le i\le L).

Pascal’s identity transfers nonnegativity toward smaller shifts. Together with the alternating initial coefficients at shift 00, this identifies (n−1i)\binom{n-1}{i} as the boundary within this family of shifted binomial bases: it retains a direct nonnegative counting interpretation, while every shift to its right introduces a negative coefficient.

Paper and publication

I carried out most of this work at the Research Science Institute in summer 2021, mentored by Pakawut Jiradilok, who proposed the topic and guided the research. The paper is sole-authored.

Angel Raychev. “A generalization of descent polynomials.” Discrete Mathematics 346(1), Article 113105, 2023. DOI: 10.1016/j.disc.2022.113105.

The downloadable paper is the original September 2021 arXiv preprint, not the publisher’s typeset version. It contains the full proofs, seven figures, and references to the descent-polynomial and symmetric-function literature.

The central outcome is a structured family of counting polynomials: an elementary prefix decomposition explains polynomiality, ribbon tableaux place it in the language of symmetric functions, and stabilization exposes coefficients that themselves count finite order patterns.

Resources