This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Unit interval parking functions and the rr-Fubini numbers

S. Alex Bradt School of Mathematical and Statistical Sciences, Arizona State University, Tempe, AZ 85281 [email protected] Jennifer Elder Department of Computer Science, Math & Physics, Missouri Western State University, Saint Joseph, MO 64501 [email protected] Pamela E. Harris Department of Mathematical Sciences, University of Wisconsin-Milwaukee, Milwaukee, WI 53211 [email protected] Gordon Rojas Kirby Department of Mathematics and Statistics, San Diego State University, San Diego, CA 92182 [email protected] Eva Reutercrona Pacific Lutheran University, Tacoma, WA 98447 [email protected] Yuxuan (Susan) Wang Mount Holyoke College, South Hadley, MA 01075 [email protected]  and  Juliet Whidden Vassar College, Poughkeepsie, NY 12604 [email protected]
Abstract.

We recall that unit interval parking functions of length nn are a subset of parking functions in which every car parks in its preference or in the spot after its preference, and Fubini rankings of length nn are rankings of nn competitors allowing for ties. We present an independent proof of a result of Hadaway, which establishes that unit interval parking functions and Fubini rankings are in bijection. We also prove that the cardinality of these sets are given by Fubini numbers. In addition, we give a complete characterization of unit interval parking functions by determining when a rearrangement of a unit interval parking function is again a unit interval parking function. This yields an identity for the Fubini numbers as a sum of multinomials over compositions. Moreover, we introduce a generalization of Fubini rankings, which we call the rr-Fubini rankings of length n+rn+r. We show that this set is in bijection with unit interval parking functions of length n+rn+r where the first rr cars have distinct preferences. We conclude by establishing that these sets are enumerated by the rr-Fubini numbers.

Key words and phrases:
ennumerative combinatorics, parking function, fubini number, fubini ranking
1991 Mathematics Subject Classification:
Primary: 05A05; Secondary 05A19
P. E. Harris was supported through a Karen Uhlenbeck EDGE Fellowship.
This material is based upon work supported by the National Science Foundation under Grant No. DMS-1929284 while the authors were in residence at the Institute for Computational and Experimental Research in Mathematics in Providence, RI

1. Introduction

Throughout we let n={1,2,3,}n\in\mathbb{N}=\{1,2,3,\ldots\} and [n]={1,2,,n}[n]=\{1,2,\ldots,n\}. In this work, we are concerned with studying and enumerating families of nn-tuples in [n]n[n]^{n} satisfying certain properties and providing bijections between them. Our main objects of study are unit interval parking functions and Fubini rankings. We begin by defining each as subsets of [n]n[n]^{n} and recalling some known results from the literature.

Parking functions, which were introduced by Konheim and Weiss [9] in the context of hashing problems, can be defined as nn-tuples α=(a1,,an)[n]n\alpha=(a_{1},\dots,a_{n})\in[n]^{n} such that at least ii entries of α\alpha are at most ii, for all i[n]i\in[n]. We denote the set of parking functions of length nn by PFn\text{PF}_{n}, and it is well-known that |PFn|=(n+1)n1|\text{PF}_{n}|=(n+1)^{n-1}, [13].

We consider a special subset of parking functions called unit interval parking functions, as defined by Hadaway in [8]. In this case, one considers a queue of nn cars entering a one-way street with nn parking spots labeled increasingly from 11 to nn. The preferences of the nn cars are provided in a preference list α=(a1,a2,,an)[n]n\alpha=(a_{1},a_{2},\ldots,a_{n})\in[n]^{n}, where car ii prefers parking spot aia_{i}, for all i[n]i\in[n]. Car ii drives to its preference aia_{i} and if that parking spot is available it parks there. Otherwise, it attempts to park in spot ai+1a_{i}+1. If that spot is available, then parking is successful. Otherwise, car ii fails to park. If all cars are able to park using this unit interval parking rule, then we say that α\alpha is a unit interval parking function of length nn. For example, (1,1,2)(1,1,2) is a unit interval parking function while the rearrangement (2,1,1)(2,1,1) is not. We let UPFn\text{UPF}_{n} denote the set of unit interval parking functions of length nn.

Now consider the set of possible rankings for how nn competitors can rank in a competition allowing ties. In this setting, no ranks can be skipped, whenever two competitors tie for a rank they “cover” that rank and the next, and similar whenever more than two competitors tie they cover. For example, (2,3,5,3,1,5)(2,3,5,3,1,5) is one such ranking, while (1,1,4,4,5,6)(1,1,4,4,5,6) is not. These rankings were studied by Gross [7], who showed that the number of such rankings is given by a Fubini number111The numbers were named after Fubini by Comtet in [4]. We recall that the Fubini numbers, also known as the ordered Bell numbers[12, A000670], are defined by

(1) Fbn=k=0nj=0k(1)kj(kj)jn.\displaystyle\text{Fb}_{n}=\sum_{{k=0}}^{n}\sum_{{j=0}}^{k}(-1)^{{k-j}}{\binom{k}{j}}j^{n}.

Hadaway called these rankings “Fubini rankings” and established a bijection between them and unit interval parking functions [8, Theorem 5.12]. Her bijection depended on the content of a parking function, i.e. the multiset describing the values appearing in the tuple.

Motivated by the results of Hadaway [8], in this work we give a direct bijection between unit interval parking functions and Fubini rankings, which highlights their “block structure” (Definition 2.7). The block structure of these combinatorial objects plays a key role in our answer of Hadaway’s question of which rearrangements of a unit interval parking function is again a unit interval parking function. Moreover, this helped us uncover a new formula for the Fubini numbers. The block structure also allows us to fully generalize the result to a new family of combinatorial objects which we call rr-Fubini rankings.

This article is organized as follows. In Section 2, we study unit interval parking functions and we:

  1. (1)

    Give a complete characterization of unit interval parking functions by proving which rearrangements of unit interval parking functions are again unit interval parking functions (Theorem 2.9).

  2. (2)

    Establish a bijection between the set of unit interval parking functions of length nn and Fubini rankings of length nn (Theorem 2.5). This implies immediately unit interval parking functions of length nn are enumerated by the Fubini numbers, Equation (1). As a consequence, we give a new formula for the Fubini numbers (Corollary 2.13):

    Fbn=k=1n(c=(c1,c2,,ck)n(nc1,c2,,ck)),\text{Fb}_{n}=\displaystyle\sum_{k=1}^{n}\left(\sum_{\textbf{c}=(c_{1},c_{2},\ldots,c_{k})\models n}\binom{n}{c_{1},c_{2},\ldots,c_{k}}\right),

    where c=(c1,c2,,ck)n\textbf{c}=(c_{1},c_{2},\ldots,c_{k})\models n denotes that c is a composition of nn with kk parts.

In Section 3, we introduce the rr-Fubini rankings of length n+rn+r (Definition 3.2). Our main results establish that they are enumerated by the rr-Fubini numbers. We recall that the rr-Stirling numbers of second kind222The 22-Stirling numbers of the second kind are [12, A143494]. are defined by

{n+rk+r}r=1(kr)!i=0kr(1)k+i+r(kri)(i+r)nr,\left\{\begin{matrix}n+r\\ k+r\end{matrix}\right\}_{r}=\frac{1}{(k-r)!}\sum_{i=0}^{k-r}(-1)^{k+i+r}\binom{k-r}{i}(i+r)^{n-r},

and the rr-Fubini numbers333Note that the 11-Fubini numbers are the Fubini numbers defined in (1). are a generalization of the Fubini numbers defined in [1, 11] by

Fbnr=k=0n(k+r)!{n+rk+r}r.\text{Fb}^{r}_{n}=\sum_{k=0}^{n}(k+r)!\left\{\begin{matrix}n+r\\ k+r\end{matrix}\right\}_{r}.
  1. (3)

    We establish that rr-Fubini rankings are enumerated by the rr-Fubini numbers (Theorem 3.6) by giving a bijection between rr-Fubini rankings and unit interval parking functions starting with rr distinct values (Theorem 3.5).

We end in Section 4 detailing some directions for future research.

2. Bijection between unit interval parking functions and Fubini rankings

In this section, we establish the bijection between UPFn\text{UPF}_{n} and the set of Fubini rankings of length nn, denoted by FRn\text{FR}_{n} (Theorem 2.5). We remark that this result was originally established by Hadaway in [8]. We begin by giving a formal definition of Fubini rankings.

Definition 2.1.

The nn-tuple β=(b1,b2,,bn)[n]n\beta=(b_{1},b_{2},\ldots,b_{n})\in[n]^{n} is a Fubini ranking of length nn if the following holds: For all x[n]x\in[n], if k>0k>0 entries of β\beta are equal to xx, then the next largest value in β\beta is x+kx+k. Denote the set of Fubini rankings of length nn by FRn\text{FR}_{n}.

Note that (1,2,3)(1,2,3) is a Fubini ranking, while (1,3,3)(1,3,3) is not. The following result is due to Cayley [2].

Proposition 2.2.

For n1n\geq 1, |FRn|=Fbn|\text{FR}_{n}|=\text{Fb}_{n}.

We remark that any rearrangement of a Fubini ranking is again a Fubini ranking, as a Fubini ranking only depends on the ranks and not who places in which rank. Recall that this is also true for parking functions as any rearrangements of a parking function is also a parking function. However, this is not true for unit interval parking functions. For example, (1,1,2)(1,1,2) is a unit interval parking function while (1,2,1)(1,2,1) is not. We soon address the question of when a rearrangement of a unit interval parking function is again a unit interval parking function (Theorem 2.9). Before doing so, we consider Fubini rankings as parking preferences and establish that Fubini rankings are always parking functions.

Lemma 2.3.

For all n1n\geq 1, FRnPFn\text{FR}_{n}\subseteq\text{PF}_{n}.

Proof.

Suppose that β=(b1,,bn)\beta=(b_{1},\dots,b_{n}) is a Fubini ranking. Then by our previous remark we know that any rearrangement is also a Fubini ranking. Thus consider the weakly increasing rearrangement of β\beta which we call β=(d1,d2,,dn)\beta^{\uparrow}=(d_{1},d_{2},\ldots,d_{n}) which satisfies that didi+1d_{i}\leq d_{i+1} for all i[n1]i\in[n-1]. To show that βPFn\beta\in\text{PF}_{n} it suffices to check that diid_{i}\leq i for all i[n]i\in[n]. By Definition 2.1, there exists j[n]j\in[n] values d1=d2==dj=1d_{1}=d_{2}=\ldots=d_{j}=1, thus, values diid_{i}\leq i for all i[j]i\in[j]. Then, again by definition, the next entry is dj+1=j+1d_{j+1}=j+1, which satisfies dj+1j+1d_{j+1}\leq j+1. Iterating this process, for k>j+1k>j+1 if there are ties then dk<kd_{k}<k or there is not tie and the value at dkd_{k} is a new rank, hence dk=kkd_{k}=k\leq k. Thus, diid_{i}\leq i, for all i[n]i\in[n], and βPFn\beta\in\text{PF}_{n}. ∎

Observation 2.4.

Consider β=(b1,,bn)FRn[n]n\beta=(b_{1},\dots,b_{n})\in\text{FR}_{n}\subseteq[n]^{n} as a parking preference. Observe that cars at indices corresponding to a repeated value jj in β\beta do not interact with any car cic_{i} such that bijb_{i}\neq j. Thus, if cic_{i} has parking preference jj, i.e. bi=jb_{i}=j, then cic_{i} parks in spot j+kj+k, where kk denotes the number of occurrences of jj in β\beta up to car cic_{i}.

Next we provide an independent proof of the following result [8, Theorem 5.12].

Theorem 2.5.

If n1n\geq 1, then the sets UPFn\text{UPF}_{n} and FRn\text{FR}_{n} are in bijection, and |UPFn|=|FRn|=Fbn|\text{UPF}_{n}|=|\text{FR}_{n}|=\text{Fb}_{n}.

To make our bijection precise we begin by defining the displacement of a parking function: Given a parking function α=(a1,a2,,an)\alpha=(a_{1},a_{2},\ldots,a_{n}), if car ii parks in spot sis_{i}, then siais_{i}-a_{i} is the displacement of car ii, and the sum d(α)=i=1n(siai)d(\alpha)=\sum_{i=1}^{n}(s_{i}-a_{i}) is called the total displacement of α\alpha. We note that there are interesting results related to the study of parking functions and their displacement. This includes a bijection that takes the total displacement of parking functions to the number of inversions of labeled trees [10]. In our work, displacement is utilized as a way to test whether a parking function is a unit interval parking function. Namely, we know that a parking functions is a unit interval parking function if and only if each car has displacement at most 1. In what follows, we say that whenever the displacement of a car is equal to zero, then the car is lucky. Lucky cars and tree inversions were studied in [6, 15].

To prove Theorem 2.5 we define a bijection from Fubini rankings (which we know are parking functions by Lemma 2.3) to unit interval parking functions that modifies the entries of the Fubini ranking. This map will do this in such a way so as to ensure that each car still parks in the same space, but is displaced at most one, which makes it a unit interval parking function.

Lemma 2.6.

Given β=(b1,,bn)FRn\beta=(b_{1},\ldots,b_{n})\in\text{FR}_{n}, define ϕ:FRnUPFn\phi:\text{FR}_{n}\to\text{UPF}_{n} by ϕ(β)=(a1,,an)UPFn\phi(\beta)=(a_{1},\dots,a_{n})\in\text{UPF}_{n}, where

ai={biif |{1ji:bj=bi}|=1bi+k2if |{1ji:bj=bi}|=k>1.a_{i}=\begin{cases}b_{i}&\text{if $|\{1\leq j\leq i:b_{j}=b_{i}\}|=1$}\\ b_{i}+k-2&\text{if $|\{1\leq j\leq i:b_{j}=b_{i}\}|=k>1$.}\end{cases}

The map ϕ\phi is well defined.

Proof.

We show that ϕ(β)\phi(\beta) is a unit interval parking function whenever β\beta is a Fubini ranking. To do so we induct on mm, the number of distinct values in β\beta.

If m=1m=1, then β=(1,1,1)\beta=(1,1\ldots,1) and ϕ(β)=(1,1,2,3,4,,n1)\phi(\beta)=(1,1,2,3,4,\ldots,n-1). As a parking function, ϕ(β)\phi(\beta) parks the cars in order 1,2,,n1,2,\ldots,n, and they all park within one spot of their preference. Thus ϕ(β)\phi(\beta) a unit interval parking function.

Suppose that if β\beta has mm distinct values, then ϕ(β)\phi(\beta) is a unit interval parking function. Now consider β\beta having m+1m+1 distinct values.

Let I[n]I\subseteq[n] be the set of indices where the smallest mm distinct values of β\beta occur, and let J=[n]IJ=[n]\setminus I be the set of indices where the max value of β\beta occurs. Note that by assumption, the entries in ϕ(β)\phi(\beta) indexed by II park in the first |I||I| parking spots, and they do so by parking at most one away from their preference.

Now consider the entries in β\beta indexed by JJ, which contain the values max(β)\max(\beta). Since β\beta is a Fubini ranking, we know that the next available rank is |I|+1|I|+1. Moreover note that |I|+1=max(β)|I|+1=\max(\beta), as otherwise the rank |I|+1|I|+1 would be unearned and β\beta would not have been a Fubini ranking. To finish building ϕ(β)\phi(\beta), we need only consider the values in index set JJ. By definition, those values in β\beta were all max(β)\max(\beta) and in ϕ(β)\phi(\beta) they have now been replaced with the values max(β),max(β),max(β)+1,,max(β)+|J|2\max(\beta),\max(\beta),\max(\beta)+1,\ldots,\max(\beta)+|J|-2, appearing in this relative order at the indices in JJ. We conclude by now noting that under ϕ(β)\phi(\beta) the cars indexed by JJ park in spots |I|+1,|I|+2,,n=max(β),max(β)+1,,n|I|+1,|I|+2,\ldots,n=\max(\beta),\max(\beta)+1,\ldots,n. Thus, the cars indexed by JJ parking under ϕ(β)\phi(\beta) are displaced at most one from their preference, as desired. ∎

In [8, Lemma 5.6] Hadaway established that if αUPFn\alpha\in\text{UPF}_{n}, then a value aa may appear at most twice. Below, we characterize when a rearrangement of a unit interval parking function is again a unit interval parking function.

To make our approach precise, we begin with the following general definition on parking functions, which in Theorem 2.9 below we specialize to unit interval parking functions.

Definition 2.7.

Given α=(a1,,an)PFn\alpha=(a_{1},\ldots,a_{n})\in\text{PF}_{n}, let α=(a1,,an)\alpha^{\uparrow}=(a^{\prime}_{1},\ldots,a^{\prime}_{n}) be the weakly increasing rearrangement of α\alpha. The partition of α\alpha^{\uparrow} as a concatenation denoted α=π1|π2||πm\alpha^{\prime}=\pi_{1}|\pi_{2}|\cdots|\pi_{m}, where πj\pi_{j} begins at (and includes) the jjth entry α\alpha^{\uparrow} satisfying ai=ia_{i}^{\prime}=i, is the block structure of α\alpha. Each πi\pi_{i} is called a block of α\alpha.

For example, if α=(8,1,5,5,1,2,4,7)PF8\alpha=(8,1,5,5,1,2,4,7)\in\text{PF}_{8}, then α=(1,1,2,4,5,5,7,8)\alpha^{\uparrow}=(1,1,2,4,5,5,7,8) and α=112|4|55|7|8\alpha^{\prime}=112|4|55|7|8.

Observation 2.8.

Note that if α\alpha is a unit interval parking function with block structure π1|π2||πm\pi_{1}|\pi_{2}|\cdots|\pi_{m}, then as cars in each block park, they park without interacting with cars in other blocks. To make this precise we let (πi)\ell(\pi_{i}) denote the length of πi\pi_{i} for all i[m]i\in[m]. Then cars in π1\pi_{1} park and occupy spots 11 through (π1)\ell(\pi_{1}). By definition of π2\pi_{2} car c(π1)+1c_{\ell(\pi_{1})+1} has preference (π1)+1\ell(\pi_{1})+1, and hence parks precisely in spot (π1)+1\ell(\pi_{1})+1. As all other elements in π2\pi_{2} are in weakly increasing order the remaining cars of π2\pi_{2} park in spots (π1)+2,(π1)+3,,(π1)+(π2)\ell(\pi_{1})+2,\ell(\pi_{1})+3,\ldots,\ell(\pi_{1})+\ell(\pi_{2}). This shows that the cars in π1\pi_{1} and π2\pi_{2} park without affecting each other and this fact extends to all πj\pi_{j} with 2<jm2<j\leq m.

Lastly, we note that each block of α\alpha corresponds to a block of repeated values in a Fubini ranking, as described in Observation 2.4.

Theorem 2.9.

Given α=(a1,,an)UPFn\alpha=(a_{1},\ldots,a_{n})\in\text{UPF}_{n}, let α\alpha^{\uparrow} be its weakly increasing rearrangement and α=π1|π2||πm\alpha^{\prime}=\pi_{1}|\pi_{2}|\dots|\pi_{m} be the block structure of α\alpha (as in Definition 2.7).

  1. (1)

    There are

    (2) (n|π1|,,|πm|)\displaystyle\binom{n}{|\pi_{1}|,\ldots,|\pi_{m}|}

    possible rearrangements σ\sigma of α\alpha such that σ\sigma is still a unit interval parking function.

  2. (2)

    A rearrangement σ\sigma of α\alpha is in UPFn\text{UPF}_{n} if and only if the entries in σ\sigma respect the relative order of the entries in each of the blocks π1,π2,,πm\pi_{1},\pi_{2},\ldots,\pi_{m}.

Proof.

To establish Part (1) we prove the following claims by inducting on the number of blocks in α\alpha^{\prime}.

  • The weakly increasing rearrangement of an element α\alpha in UPFn\text{UPF}_{n} is also in UPFn\text{UPF}_{n}.

  • If πj\pi_{j} is a weakly increasing block of α\alpha^{\prime}, then those elements appear in weakly increasing order in α\alpha.

  • Any two blocks πi\pi_{i} and πj\pi_{j} are independent of each other, i.e. as sets πiπj=\pi_{i}\cap\pi_{j}=\emptyset for all iji\neq j.

  • The formula in (2) holds.

Suppose m=1m=1. Then α=(a1,,an)\alpha^{\uparrow}=(a_{1}^{\prime},\ldots,a_{n}^{\prime}) is such that 1aiai+11\leq a_{i}^{\prime}\leq a_{i+1}^{\prime}, aiia_{i}^{\prime}\leq i since it is a parking function, and aiia_{i}^{\prime}\neq i for i>1i>1 by assumption. Thus, a1=1a_{1}^{\prime}=1, which implies aa2<2a\leq a_{2}^{\prime}<2. That is a2=1a_{2}^{\prime}=1 and car 2 parks in spot 2 and is displaced 1 spot from its preference. For car 3 to have displacement at most 1 we must have 2<a32<a_{3}^{\prime} since spots 1 and 2 are occupied so that 2a3<32\leq a_{3}^{\prime}<3, i.e. a3=2a_{3}^{\prime}=2. Continuing in this fashion establishes α=(1,1,2,3,n1)\alpha^{\uparrow}=(1,1,2,3\ldots,n-1). Next we establish that α=α\alpha=\alpha^{\uparrow} by supposing that α\alpha is out of weakly increasing order to obtain a contradiction.

Let ii be the smallest index such that ai>i1a_{i}>i-1 and there exists some k>ik>i satisfying ai>aja_{i}>a_{j}. Let jj be the largest such kk. Parking under α\alpha, cars 1,2,3,,i11,2,3,\ldots,i-1 park in spots 1,2,,i11,2,\ldots,i-1, respectively. Car ii with preference ai>i1a_{i}>i-1 parks in spot aia_{i} leaving spots i,i+1,,ai1i,i+1,\ldots,a_{i}-1 unoccupied after it parks. Note that of the cars numbered i+1,i+2,,aiaj+ii+1,i+2,\ldots,a_{i}-a_{j}+i, those with preferences in the set {i,i+1,,ai1}\{i,i+1,\ldots,a_{i}-1\} park without displacement and those with preference smaller that aja_{j} get displaced one unit. Also, one of the cars with preference smaller than aja_{j} has parked in spot aja_{j}. Further, any car with preference larger than aja_{j} parks with a displacement at most one. This is implied because α\alpha consists of the numbers 1,1,2,3,,n11,1,2,3,\ldots,n-1 so that all numbers not equal to one appear exactly once in α\alpha. This then ensures that spots i,i+1,ai1i,i+1,\ldots a_{i}-1 are occupied by cars arriving and parking before car jj. Thus, the displacement of car jj is at best aiaj+1>1+12a_{i}-a_{j}+1>1+1\geq 2, which contradicts that αUPFn\alpha\in UPF_{n}. Therefore, the unique unit interval parking function consisting of a single block is (1,1,2,,n1)(1,1,2,\ldots,n-1), and it is in weakly increasing order. Furthermore, there are (n|π1|)=(nn)=1\binom{n}{|\pi_{1}|}=\binom{n}{n}=1 possible rearrangements that produce a unit parking function.

Suppose that if αUPFn\alpha\in\text{UPF}_{n} is such that α\alpha^{\prime} has mm distinct blocks, then those blocks appear in α\alpha as weakly increasing strings, and there are (n|π1|,,|πm|)\binom{n}{|\pi_{1}|,\ldots,|\pi_{m}|}, ways to reorder α\alpha so that the result is again in UPFn\text{UPF}_{n}.

Now consider αUPFn\alpha\in\text{UPF}_{n} having m+1m+1 distinct blocks in α\alpha^{\prime}. We may consider the block structure of α\alpha as the concatenation α=βγ\alpha^{\prime}=\beta^{\prime}\gamma^{\prime}, where the block structure of β\beta is given by β=π1|π2||πm\beta^{\prime}=\pi_{1}|\pi_{2}|\dots|\pi_{m}, and the block structure of γ\gamma is given by γ=πm+1\gamma^{\prime}=\pi_{m+1}. As in Observation 2.8 we may now consider both β\beta and γ\gamma as unit parking functions on their own, also noting that max(β)+1=min(γ)\max(\beta)+1=\min(\gamma), otherwise we would not have a parking function.

By induction, βUPFn|πm+1|\beta\in\text{UPF}_{n-|\pi_{m+1}|}, and there are (n|πm+1||π1||πm|)\binom{n-|\pi_{m+1}|}{|\pi_{1}|\ldots|\pi_{m}|} ways to rewrite it as a unit parking function, by induction hypothesis. Separate from this, γ\gamma has the form of an element in UPF|πm+1|\text{UPF}_{|\pi_{m+1}|}, and there is one way to arrange these elements. Because γ\gamma fills parking spots n|πm+1|+1,n|πm+1|+2,,nn-|\pi_{m+1}|+1,n-|\pi_{m+1}|+2,\ldots,n, it does not interfere with β\beta being a unit parking function, and thus β\beta and γ\gamma can be parked “without interactions.” That is, cars can be given preferences from a valid rearrangement of β\beta at the same time as cars are given preferences from γ\gamma. Then if β\beta is a unit parking function produced by rearranging β\beta, there are

(nn|πm+1|)(|πm+1||πm+1|)=(nn|πm+1|,|πm+1|)\binom{n}{n-|\pi_{m+1}|}\binom{|\pi_{m+1}|}{|\pi_{m+1}|}=\binom{n}{n-|\pi_{m+1}|,|\pi_{m+1}|}

ways that we could intermix the preferences from β\beta and γ\gamma. By induction, there are (n|πm+1||π1|,,|πm|)\binom{n-|\pi_{m+1}|}{|\pi_{1}|,\ldots,|\pi_{m}|} ways to rearrange β\beta as a unit parking function, thus there are

(n|πm+1||π1|,,|πm|)(|πm+1||πm+1|)=(n|π1|,,|πm|,|πm+1|)\binom{n-|\pi_{m+1}|}{|\pi_{1}|,\ldots,|\pi_{m}|}\binom{|\pi_{m+1}|}{|\pi_{m+1}|}=\binom{n}{|\pi_{1}|,\ldots,|\pi_{m}|,|\pi_{m+1}|}

ways to rewrite α\alpha as a unit parking function, as desired.

To prove Part (2) first note that the above proof of Part (1) implies that if σ\sigma respects the relative orders of the entries in the blocks, then σUPFn\sigma\in\text{UPF}_{n}. It suffices to prove the reverse direction, that given α\alpha a unit interval parking function, then it respects the relative order of the entries in each block. To show this, note that the set of cars with preferences within a block of α\alpha park independently of each other. We know that every block of α\alpha must be a unit interval parking function in its own right, if α\alpha is a unit interval parking function. From the base case of Part (1), we established that a single block is a unit interval parking function if and only if its entries are in weakly increasing order. Thus, a unit interval parking function must have all of the values within its blocks appearing in weakly increasing order, as claimed. ∎

Remark 2.10.

Chaves Meyles, Harris, Jordaan, Rojas Kirby, Sehayek, and Spingarn give a very interesting and surprising connection between unit interval parking functions and the permutohedron [3]. They also give an independent proof of Part (2) of Theorem 2.9 using a “prime decomposition” of parking functions.

We are ready to introduce the inverse to the map in Lemma 2.6.

Lemma 2.11.

Given α=(a1,,an)UPFn\alpha=(a_{1},\ldots,a_{n})\in\text{UPF}_{n}, let α\alpha^{\uparrow} be its weakly increasing rearrangement and α=π1|π2||πm\alpha^{\prime}=\pi_{1}|\pi_{2}|\dots|\pi_{m} be the block structure of α\alpha (as in Definition 2.7). Define ψ:UPFnFRn\psi:\text{UPF}_{n}\mapsto\text{FR}_{n} by ψ(α)=(b1,bn)\psi(\alpha)=(b_{1},\dots b_{n}), where, for all i[n]i\in[n],

bi=min{aπj:aiπj}.b_{i}=\min\{a\in\pi_{j}:a_{i}\in\pi_{j}\}.

The map ψ\psi is well defined.

Before establishing Lemma 2.11 we provide an example to illustrate the map ψ\psi.

Example 2.12.

Let n=10n=10 and consider α=(2,4,7,4,1,5,7,8,2,9)\alpha=(2,4,7,4,1,5,7,8,2,9). We leave it to the reader to verify that αUPFn\alpha\in\text{UPF}_{n}. The weakly increasing rearrangement alphaalpha^{\uparrow} and the block structure of α\alpha, denoted α\alpha^{\prime}, are as follows:

α\displaystyle\alpha^{\uparrow} =(1,2,2,4,4,5,7,7,8,9)\displaystyle=(1,2,2,4,4,5,7,7,8,9)
α\displaystyle\alpha^{\prime} =1π1|22π2|445π3|7789π4.\displaystyle=\underbrace{1}_{\pi_{1}}|\underbrace{22}_{\pi_{2}}|\underbrace{445}_{\pi_{3}}|\underbrace{7789}_{\pi_{4}}.

For each i[10]i\in[10], bib_{i} is the smallest element in the block containing the value aia_{i}; that is, β=(2,4,7,4,1,4,7,7,2,7)\beta=(2,4,7,4,1,4,7,7,2,7). One can verify that β\beta is a Fubini ranking.

Proof of Lemma 2.11.

Let α=(a1,a2,,an)UPFn\alpha=(a_{1},a_{2},\ldots,a_{n})\in\text{UPF}_{n} and let α=(a1,a2,,an)\alpha^{\uparrow}=(a_{1}^{\prime},a_{2}^{\prime},\ldots,a_{n}^{\prime}) be the weakly increasing rearrangement of α\alpha.

Let I[n]I\subseteq[n] be the set of indices satisfying ai=ia_{i}^{\prime}=i and let m=|I|m=|I|. Then the block structure of α\alpha can be partitioned as the concatenation α=π1|π2||πm\alpha^{\prime}=\pi_{1}|\pi_{2}|\dots|\pi_{m}, where a new block πj\pi_{j} begins at (and includes) each aia^{\prime}_{i} satisfying ai=ia^{\prime}_{i}=i.

Because α\alpha is a unit interval parking function, and thus α\alpha\uparrow is a parking function, we note that ahha^{\prime}_{h}\leq h for all h[n]h\in[n]. This means that if ah=ha^{\prime}_{h}=h, and is the minimum element of a new block, this is the first occurrence of the number hh in α\alpha^{\uparrow}. Therefore, there are h1h-1 elements in α\alpha^{\uparrow} less than hh, all appearing in the blocks before whichever block contains hh. Let ah=hπja^{\prime}_{h}=h\in\pi_{j}. This integer, hh, may then be partitioned as follows:

h=1+(h1)=1+l=1j1|πl|.h=1+(h-1)=1+\sum_{l=1}^{j-1}|\pi_{l}|.

Therefore, in general, we have that

min(πj)=1+l=1j1|πl|.\min(\pi_{j})=1+\sum_{l=1}^{j-1}|\pi_{l}|.

We also note that the elements min(πj)\min(\pi_{j}) and min(πj+1)\min(\pi_{j+1}) then have the relationship min(πj)+|πj|=min(πj+1)\min(\pi_{j})+|\pi_{j}|=\min(\pi_{j+1}). Moreover, by definition of bib_{i}, there are exactly |πj||\pi_{j}| copies of the value min(πj)\min(\pi_{j}) in ψ(α)\psi(\alpha). Also, note that by definition α\alpha^{\uparrow} begins with 1, and so min(π1)=1\min(\pi_{1})=1, and hence there exists a bj=1b_{j}=1 in ψ(α)\psi(\alpha). These facts together imply that ψ(α)\psi(\alpha) is a Fubini ranking. ∎

The following result gives a formula for the Fubini numbers, which we did not find in the literature.

Corollary 2.13.

If n1n\geq 1 and c=(c1,c2,,ck)n\textbf{c}=(c_{1},c_{2},\ldots,c_{k})\vDash n denotes a composition of nn with kk parts, then

Fbn=k=1n(c=(c1,c2,,ck)n(nc1,c2,,ck)).\text{Fb}_{n}=\displaystyle\sum_{k=1}^{n}\left(\sum_{\textbf{c}=(c_{1},c_{2},\ldots,c_{k})\models n}\binom{n}{c_{1},c_{2},\ldots,c_{k}}\right).
Proof.

We count the number of Fubini rankings with kk blocks (1kn1\leq k\leq n), where blocks are defined as in Definition 2.7. From Theorem 2.9, we know that, for each composition (c1,c2,,ck)(c_{1},c_{2},\ldots,c_{k}) of nn, the number of Fubini rankings with cic_{i} elements in the iith block is given by the multinomial (nc1,c2,,ck)\binom{n}{c_{1},c_{2},\ldots,c_{k}}. ∎

We are now ready to prove our main result.

Proof of Theorem 2.5.

It suffices to show that ϕ\phi and ψ\psi are inverses of each other.

First, we prove that for every βFRn\beta\in\text{FR}_{n}, ψ(ϕ(β))=β\psi(\phi(\beta))=\beta. Suppose β\beta has mm distinct values, x1<x2<<xmx_{1}<x_{2}<\cdots<x_{m}.

For each i[m]i\in[m], let I(xi):={j[n]:bj=xi}I(x_{i}):=\{j\in[n]:b_{j}=x_{i}\} be the set of indices containing the value xix_{i} in β\beta. Since β\beta is a Fubini ranking, we have that 1+j=1i1|I(xj)|=xi1+\sum_{j=1}^{i-1}|I(x_{j})|=x_{i} for all i[m]i\in[m]. (In particular, x1=1x_{1}=1.) For each i[m]i\in[m], ϕ\phi takes the values xix_{i} at index set I(xi)I(x_{i}) in β\beta and replaces them by the values xi,xi,xi+1,xi+2,,xi+|I(xi)|2x_{i},x_{i},x_{i}+1,x_{i}+2,\ldots,x_{i}+|I(x_{i})|-2 at index set I(xi)I(x_{i}) (in that relative order) in ϕ(β)\phi(\beta).

To apply ψ\psi to ϕ(β)\phi(\beta), we first rearrange ϕ(β)\phi(\beta) into weakly increasing order, which is given by

ϕ(β)=(x1,x1,x1+1,,x1+|I(x1)|2|I(x1)|,,xm,xm,xm+1,,xm+|I(xm)|2|I(xm)|).\phi(\beta)^{\uparrow}=(\underbrace{x_{1},x_{1},x_{1}+1,\ldots,x_{1}+|I(x_{1})|-2}_{|I(x_{1})|},\ldots,\underbrace{x_{m},x_{m},x_{m}+1,\ldots,x_{m}+|I(x_{m})|-2}_{|I(x_{m})|}).

To partition the block structure of ϕ(β)\phi(\beta) we begin by noting that, for all j[m]j\in[m], since xj=1+k=1j1|I(xk)|x_{j}=1+\sum_{k=1}^{j-1}|I(x_{k})| and the first xjx_{j} appears in index 1+k=1j1|I(xk)|1+\sum_{k=1}^{j-1}|I(x_{k})|, the first instance of xjx_{j} begins a new block. Thus

ϕ(β)=π1|π2||πm\phi(\beta)^{\prime}=\pi_{1}|\pi_{2}|\cdots|\pi_{m}

where, for all i[m]i\in[m], we have that

(3) πj=xj,xj,xj+1,,xj+|I(xj)|2|I(xj)|.\displaystyle\pi_{j}=\underbrace{x_{j},x_{j},x_{j}+1,\ldots,x_{j}+|I(x_{j})|-2}_{|I(x_{j})|}.

Now observe that the entries in ψ(ϕ(β))\psi(\phi(\beta)) are defined by

ci=min{aπj:aiπj}.c_{i}=\min\{a\in\pi_{j}:a_{i}\in\pi_{j}\}.

Claim: ci=bic_{i}=b_{i} for all i[n]i\in[n].

To prove the claim we consider the values in the index set I(xi)I(x_{i}). In β\beta, these are xix_{i}’s. In ϕ(β)\phi(\beta), these are the values in πi\pi_{i} which are xi,xi,xi+1,,xi+|I(xi)|2x_{i},x_{i},x_{i}+1,\ldots,x_{i}+|I(x_{i})|-2 (still occurring at the indices I(xi)I(x_{i})). In ψ(ϕ(β))\psi(\phi(\beta)) for any iI(xi)i\in I(x_{i}), the value is min{aπi:aiπi}=xi\min\{a\in\pi_{i}:a_{i}\in\pi_{i}\}=x_{i}. So for every iI(xi)i\in I(x_{i}), ψ(ϕ(β))=xi\psi(\phi(\beta))=x_{i} which is precisely the value of β\beta at every iI(xi)i\in I(x_{i}). Since this holds for all i[m]i\in[m] and since i=1mI(xi)=[n]\cup_{i=1}^{m}I(x_{i})=[n], we have established that ψ(ϕ(β))i=bi\psi(\phi(\beta))_{i}=b_{i} (showing the lists agree at every entry). Hence ψ(ϕ(β))=β\psi(\phi(\beta))=\beta.

For the reverse composition, we begin with a unit interval parking function α=(a1,a2,,an)UPFn\alpha=(a_{1},a_{2},\ldots,a_{n})\in\text{UPF}_{n} and let α=(a1,a2,,an)\alpha^{\uparrow}=(a_{1}^{\prime},a_{2}^{\prime},\ldots,a_{n}^{\prime}) be the weakly increasing rearrangement of α\alpha and we partition the block structure of α\alpha as

α=π1|π2||πm\alpha^{\prime}=\pi_{1}|\pi_{2}|\cdots|\pi_{m}

where a new block πj\pi_{j} begins and includes each aia_{i}^{\prime} satisfying ai=ia_{i}^{\prime}=i. For each i[m]i\in[m], let I(πj)={i[n]:aiπj}I(\pi_{j})=\{i\in[n]:a_{i}\in\pi_{j}\}. Then use Lemma 2.11 to find ψ(α)=(b1,b2,,bn)\psi(\alpha)=(b_{1},b_{2},\ldots,b_{n}) where bi=min{aπj:aiπj}b_{i}=\min\{a\in\pi_{j}:a_{i}\in\pi_{j}\}. Note ψ(α)\psi(\alpha) is a Fubini ranking. Now consider ϕ(ψ(α))=(c1,c2,,cn)\phi(\psi(\alpha))=(c_{1},c_{2},\ldots,c_{n}).

Claim: ci=aic_{i}=a_{i} for all i[n]i\in[n].

To prove the claim we consider the values cic_{i} and aia_{i} where iI(πj)i\in I(\pi_{j}) for an arbitrary j[m]j\in[m]. Fix j[m]j\in[m] and let I(πj)={k1<k2<<k}I(\pi_{j})=\{k_{1}<k_{2}<\ldots<k_{\ell}\} where =|I(πj)|\ell=|I(\pi_{j})|. Let ψ(α)|I(πj)=πj,k1πj,k2πj,k\psi(\alpha)|_{I(\pi_{j})}=\pi_{j,k_{1}}\pi_{j,k_{2}}\cdots\pi_{j,k_{\ell}} denote the substring of ψ(α)\psi(\alpha) made up of the values at the indices in I(πj)I(\pi_{j}). Note that by definition of ψ\psi, πj,kb=min(πj)\pi_{j,k_{b}}=\min(\pi_{j}) for all 1b1\leq b\leq\ell. Let ϕ(ψ(α))|I(πj)=ck1ck2ck\phi(\psi(\alpha))|_{I(\pi_{j})}=c_{k_{1}}c_{k_{2}}\cdots c_{k_{\ell}} denote the substring of ϕ(ψ(α))\phi(\psi(\alpha)) made up of the values at the indices in I(πj)I(\pi_{j}). By definition of ϕ\phi,

ckb={min(πj) if b=1min(πj)+b2 if b>1.c_{k_{b}}=\begin{cases}\min(\pi_{j})&\mbox{ if $b=1$}\\ \min(\pi_{j})+b-2&\mbox{ if $b>1$.}\end{cases}

By construction/definition ak1ak2ak=πja_{k_{1}}^{\prime}a_{k_{2}}^{\prime}\cdots a_{k_{\ell}}^{\prime}=\pi_{j} and aki=ckia_{k_{i}}^{\prime}=c_{k_{i}} for all i[]i\in[\ell]. By Theorem 2.9, these values appear in α\alpha in the exact same relative order. Thus, aki=ckia_{k_{i}}=c_{k_{i}} for all i[]i\in[\ell]. ∎

3. Bijection between unit interval parking functions starting with rr distinct values and rr-Fubini rankings

We begin this section by defining rr-Fubini numbers and show that they enumerate rr-Fubini rankings. We also establish that the rr-Fubini rankings are in bijection with the set of unit interval parking functions, which begin with rr distinct values.

Definition 3.1.

The rr-Fubini numbers444The 22-Fubini numbers are defined in [12, A232472]. are given by

Fbnr=k=0n(k+r)!{n+rk+r}r\text{Fb}^{r}_{n}=\sum_{k=0}^{n}(k+r)!\left\{\begin{matrix}n+r\\ k+r\end{matrix}\right\}_{r}

where {n+rk+r}r\left\{\begin{matrix}n+r\\ k+r\end{matrix}\right\}_{r} denotes the rr-Stirling numbers of the second kind555The 22-Stirling numbers are defined in [12, A008277]., which enumerate set partitions of length n+rn+r into k+rk+r blocks, in which the first rr values appear in distinct blocks.

For ease of reference, we provide Table 1, where we give the values of these numbers for 1r81\leq r\leq 8, 1m81\leq m\leq 8, and m=n+rm=n+r.

m=n+rrm=n+r\setminus r 1 2 3 4 5 6 7 8
1 1 0 0 0 0 0 0 0
2 3 2 0 0 0 0 0 0
3 13 10 6 0 0 0 0 0
4 75 62 42 24 0 0 0 0
5 541 466 342 216 120 0 0 0
6 4683 4142 3210 2184 1320 720 0 0
7 47293 42610 34326 24696 15960 9360 5040 0
8 545835 498542 413322 310344 211560 131760 75600 40320
Table 1. The rr-Fubini numbers for 1r81\leq r\leq 8, and 1m81\leq m\leq 8, where m=n+rm=n+r.
Definition 3.2.

An 𝒓\boldsymbol{r}-Fubini ranking of length n+rn+r is a Fubini ranking of length n+rn+r whose first rr values are distinct. More precisely, α=(a1,a2,,ar,ar+1,,an+r)Fbn+r\alpha=(a_{1},a_{2},\ldots,a_{r},a_{r+1},\ldots,a_{n+r})\in\text{Fb}_{n+r} is an rr-Fubini ranking if |{a1,a2,,ar}|=r|\{a_{1},a_{2},\ldots,a_{r}\}|=r. We let FRn+rr\text{FR}^{r}_{n+r} denote the set of rr-Fubini rankings of length n+rn+r.

Example 3.3.

There are ten 22-Fubini rankings of length three with the first two values being distinct. These are: (1,3,1),(3,1,1),(1,2,2),(2,1,2),(1,2,3),(1,3,2),(2,1,3),(2,3,1),(3,1,2),(3,2,1)(1,3,1),(3,1,1),(1,2,2),(2,1,2),(1,2,3),(1,3,2),(2,1,3),(2,3,1),(3,1,2),(3,2,1). Note that this is precisely the value in Table 1 with m=3m=3 and r=2r=2.

Note that when r=1r=1, the 11-Fubini rankings allow the repetition of any of the values, and thus we recover the definition of Fubini rankings given by Hadaway [8] and which we studied in Section 2. Hence, the 11-Fubini rankings are enumerated by the Fubini numbers. Moreover, the rr-Fubini rankings are nested, meaning that

𝔖m=FRmmFRmm1FRm2FRm1=FRm,\mathfrak{S}_{m}=\text{FR}^{m}_{m}\subseteq\text{FR}^{m-1}_{m}\subseteq\cdots\subseteq\text{FR}^{2}_{m}\subseteq\text{FR}^{1}_{m}=\text{FR}_{m},

where 𝔖m\mathfrak{S}_{m} denotes the set of permutations on the set [m][m].

In what follows we let UPFmrUPFm\text{UPF}^{r}_{m}\subseteq\text{UPF}_{m} denote the set of unit interval parking functions of length mn+rm\coloneqq n+r in which the first rr values are distinct. As with rr-Fubini rankings, the sets UPFmr\text{UPF}^{r}_{m} are also nested:

𝔖m=UPFmmUPFmm1UPFm2UPFm1=UPFm.\mathfrak{S}_{m}=\text{UPF}^{m}_{m}\subseteq\text{UPF}^{m-1}_{m}\subseteq\cdots\subseteq\text{UPF}^{2}_{m}\subseteq\text{UPF}^{1}_{m}=\text{UPF}_{m}.

Even though the sets share similar properties, such as this nesting, we now provide an example which illustrates that the sets UPFn+rr\text{UPF}^{r}_{n+r} and FRn+rr\text{FR}^{r}_{n+r} are not the same.

Example 3.4.

Note (3,2,1,4,4,4)(3,2,1,4,4,4) is in FR64\text{FR}^{4}_{6} but not in UPF64\text{UPF}^{4}_{6}, since car 66 would not park within a unit interval. Note (3,2,1,4,4,5)(3,2,1,4,4,5) is in UPF64\text{UPF}^{4}_{6} and not in FR64\text{FR}^{4}_{6}, because because the tie at rank 4 would disallow rank 5 from appearing.

We are now ready to state and prove our main result of this section.

Theorem 3.5.

For r1r\geq 1 and n0n\geq 0, if m=n+rm=n+r, then the sets UPFmr\text{UPF}^{r}_{m} and FRnr\text{FR}^{r}_{n} are in bijection.

Proof.

Given a UPFm\text{UPF}_{m} in which the first rr values are distinct, we use the bijection in Theorem 2.5 to find a unique FRm\text{FR}_{m}. We know that each of the first rr cars will park in its preferred spot, so each of those values will be starting a new π\pi block and will thus be mapped to itself in the Fubini ranking. So the Fubini ranking will start with the same rr distinct values as the unit interval parking function. Thus, |UPFmr|=|FRnr||\text{UPF}^{r}_{m}|=|\text{FR}^{r}_{n}|. ∎

Next we enumerate the elements of UPFn+rr\text{UPF}^{r}_{n+r}.

Theorem 3.6.

If r1r\geq 1 and n0n\geq 0, then |UPFn+rr|=Fbnr|\text{UPF}^{r}_{n+r}|=\text{Fb}^{r}_{n}.

Proof.

Let α=(a1,a2,,ar,ar+1,,an+r)UPFn+rr\alpha=(a_{1},a_{2},\ldots,a_{r},a_{r+1},\ldots,a_{n+r})\in\text{UPF}^{r}_{n+r} be weakly increasing. As UPFn+rrUPFn+r\text{UPF}^{r}_{n+r}\subset\text{UPF}_{n+r}, by Theorem 2.9 we know the structure of α\alpha, meaning that there exists r+kr+k, with k[n]k\in[n], indices of α\alpha such that ai=ia_{i}=i. These values always begin a block πi\pi_{i}. By Theorem 2.9, we know that we can construct all elements of UPFn+rr\text{UPF}^{r}_{n+r} by permuting the blocks while keeping the relative order of the elements of each block. The number of ways to do this, for a fixed k[n]k\in[n], is given by

(k+r)!{n+rk+r}r(k+r)!\left\{\begin{matrix}n+r\\ k+r\end{matrix}\right\}_{r}

where (as in Definition 3.1) the definition of the rr-Stirling number ensures that the first k+rk+r values appear in distinct blocks. The result follows by taking the sum over all possible k[n]k\in[n]. Therefore,

|UPFn+rr|=k=0n(k+r)!{n+rk+r}r=Fbnr.|\text{UPF}^{r}_{n+r}|=\sum_{k=0}^{n}(k+r)!\left\{\begin{matrix}n+r\\ k+r\end{matrix}\right\}_{r}=\text{Fb}_{n}^{r}.\qed

Together Theorems 3.5 and 3.6 imply the following.

Corollary 3.7.

If r1r\geq 1 and n0n\geq 0, then |FRn+rr|=Fbnr|\text{FR}^{r}_{n+r}|=\text{Fb}^{r}_{n}.

4. Future Work

To begin, we wonder if it is possible to construct a simple proof of Part 2 of Theorem 2.5, that does not rely on Part 1, as that would eliminate some of the complexity in our argument. Moreover, there is much to be discovered about rr-Fubini rankings and we now provide some directions for further study:

  1. (1)

    We remark that the intersection between Fubini rankings and unit interval parking functions is non-trivial. In fact, computationally the sequence for |FRnUPFn||\text{FR}_{n}\cap\text{UPF}_{n}| is given by [12, A080599] whose terms for 1n71\leq n\leq 7 are:

    1,3,12,66,450,3690,35280.1,3,12,66,450,3690,35280.

    Stanley notes that this sequence gives the number of intervals in the weak (Bruhat) order of 𝔖n\mathfrak{S}_{n} that are Boolean algebras. A new bijective proof has now appeared [5].

  2. (2)

    One could extend the above to the case of rr-Fubini rankings with r>1r>1. Namely, enumerate the cardinalities of the set FRnrUPFn+rr\text{FR}^{r}_{n}\cap\text{UPF}^{r}_{n+r}, where we fix one of the parameters while varying the other.

  3. (3)

    The study of statistics of permutations is a well-researched mathematical area. Schumacher has results related to ascents, descents, and ties for parking functions [14] and it would be of interest to study statistics for unit interval parking functions and Fubini rankings, as well as their generalizations UPFn+rr\text{UPF}^{r}_{n+r} and FRn+rr\text{FR}^{r}_{n+r}. We hope to follow up this paper with results in this direction.

References

  • [1] Leonard Carlitz. Weighted stirling numbers of the first and second kind–ii. Fibonacci Quart, 18(3):242–257, 1980.
  • [2] Arthur Cayley. On the analytical forms called Trees, with application to the theory of chemical combinations, volume 9 of Cambridge Library Collection - Mathematics, page 427–460. Cambridge University Press, 2009.
  • [3] Lucas Chaves Meyles, Pamela E. Harris, Richter Jordaan, Gordon Rojas Kirby, Sam Sehayek, and Ethan Spingarn. Unit-interval parking functions and the permutohedron, 2023. Arxiv:2305.15554.
  • [4] Louis Comtet. Advanced combinatorics. D. Reidel Publishing Co., Dordrecht, enlarged edition, 1974. The art of finite and infinite expansions.
  • [5] Jennifer Elder, Pamela E. Harris, Jan Kretschmann, and J. Carlos Martínez Mori. Boolean intervals in the weak order of 𝔖n\mathfrak{S}_{n}, 2023. Arxiv:2201.12887, accepted at Journal of Combinatorics.
  • [6] Ira M. Gessel and Seunghyun Seo. A refinement of Cayley’s formula for trees. Electron. J. Combin., 11(2):Research Paper 27, 23, 2004/06.
  • [7] Oliver A. Gross. Preferential arrangements. Amer. Math. Monthly, 69:4–8, 1962.
  • [8] Kimberly P. Hadaway. On combinatorical problems of generalized parking functions. Williams College, Honors Thesis, 2022.
  • [9] Alan G. Konheim and Benjamin Weiss. An occupancy discipline and applications. SIAM Journal on Applied Mathematics, 14(6):1266–1274, 1966.
  • [10] Germain Kreweras. Une famille de polynômes ayant plusieurs propriétés énumeratives. Period. Math. Hungar., 11(4):309–320, 1980.
  • [11] István Mező. Periodicity of the last digits of some combinatorial sequences. Journal of Integer Sequences [electronic only], 17, 08 2013.
  • [12] OEIS Foundation Inc. The On-Line Encyclopedia of Integer Sequences, 2023. Published electronically at http://oeis.org.
  • [13] John Riordan. Ballots and trees. J. Combinatorial Theory, 6:408–411, 1969.
  • [14] Paul R. F. Schumacher. Descents in parking functions. J. Integer Seq., 21(2):Art. 18.2.3, 8, 2018.
  • [15] Seunghyun Seo and Heesung Shin. A generalized enumeration of labeled trees and reverse Prüfer algorithm. J. Combin. Theory Ser. A, 114(7):1357–1361, 2007.