]> An IDEAL Group, CLC, Project

ON COMPARABILITY OF RANDOM PERMUTATIONS

DISSERTATIO N

Presented in Partial Fulfillment of the Requirements for

the Degree Doctor of Philosophy in the Graduate

School of the Ohio State University

By

Adam Hammett, B.S.

**** *

The Ohio State University

2007

Dissertation Committee:

133C1bClbl11111111lbbCC. Approved by

Dr. Boris Pittel, Advisor

Dr. Gerald Edgar

1 ċ.

Advisor

Dr. Akos Seress Graduate Program in

Graduate Program in Mathematics

ABSTRACT

Two permutations of [n] :={1, 2, ..., n} are comparable in the Bruhat order if one

can be obtained from the other by a sequence of transpositions decreasing the number

of inversions. We show that the total number of pairs of permutations (π, σ) with

π<σ is of order (n!)2/n2 at most. Equivalently, if π, a are chosen uniformly at

random and independently of each other, then P(π<σ) is of order n-2 at most. By

a direct probabilistic argument we prove P(π<σ) is of order (0.708)n at least, so

that there is currently a wide qualitative gap between the upper and lower bounds.

Next, emboldened by a connection with Ferrers diagrams and plane partitions

implicit in Bressoud's book [13], we return to the Bruhat order upper bound and show

that for n permutations π1, ..., πr selected independently and uniformly at random,

P(π1<...<πr)=O(n-r(r-1)),

thus providing an extension of our result for pairs of permutations to chains of length

r>2.

Turning to the related weak order "-"- when only adjacent transpositions are

admissible - we use a non-inversion set criterion to prove that Pn*:=P(πσ) is

submultiplicative, thus showing existence of ρ=limPn*n. We demonstrate that ρ

is 0.362 at most. Moreover, we prove the lower bound i=1n(H(i)/i) for Pn*, where

ii

H(i):=j=1i1/j. In light of numerical experiments, we conjecture that for each

order the upper bounds for permutation-pairs are qualitatively close to the actual

behavior. We believe that extensions to r-chains similar to that for the Bruhat order

upper bound can be made for our other bounds in each order, and are presently

working in this direction.

Finally, the weak order poset happens to be a lattice, and we study some properties

of its infimums and supremums. Namely, we prove that the number of r-tuples

(π1 , . . . ' πr) of n-permutations with minimal infimum, 12 ... n, asymptotically equals

-(n!)rhr(z*)(z*)n+1, r>2, n . (1)

Here, z*=z*(r)(1, 2) is the unique (positive) root of the equation

hr(z):=j0(-1)j(j!)rzj=0

within the disk |z|<2. Moreover, (1) is also the asymptotic number of r-tuples with

maximal supremum, n(n -1) ...1.

iii

To my wife, Rachael, for her unending suppport.

iv

ACKNOWLEDGMENTS

First and foremost, I wish to thank my advisor, Boris Pittel. Our collaboration was

quite intensive on this research project, and without his countless suggestions and

endless encouragement, none of these results would have materialized. He is truly an

outstanding example of what a mentor should be.

This work was inspired in large part by a thought-provoking talk Mark Skandera

gave at the MIT Combinatorics Conference honoring Richard Stanley's 60th birthday.

We are grateful to Mark for an enlightening follow-up discussion of comparability

criteria for the Bruhat order. We thank Sergey Fomin for encouragement and for

introducing us to an instrumental notion of the permutation-induced poset, and also

for many insightful suggestions regarding the history of Bruhat order. Without Ed

Overman's invaluable guidance we would not have been able to obtain our numerical

results. Craig Lennon gave us an idea for proving an exponential lower bound in

the case of Bruhat order. We thank Mikl6s B6na for his interest in this work and

for drawing our attention to the lower bound for the number of linear extensions in

Richard Stanley's book.

v

VITA

April 28, 1979. . . . . . . . . . . . . . . . . . . . . . . . . Born - Santa Maria, CA

1997-2001 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Undergraduate,

Westmont College - Santa Barbara, CA

2001 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . B.S. in Mathematics,

Westmont College

2001-Present . . . . . . . . . . . . . . . . . . . . . . . . . . Graduate Teaching Associate,

The Ohio State University

FIELDS OF STUDY

Major Field: Mathematics

Specialization: Combinatorial Probability

vi

TABLE OF CONTENTS

Abstract ii

Dedication . 1v

Acknowledgments . v

Vita vi

List of Figures x

List of Tables xi

CHAPTER PAGE

1 Introduction 1

1.1 Bruhat Order, Preliminaries 1

1.2 Main Results Related to the Bruhat Order 5

1.3 Weak Order, Preliminaries 7

1.4 Main Results Related to the Weak Order 8

2The Proof of the Bruhat Order Upper Bound 14

2.1 A Necessary Condition for Bruhat Comparability 14

2.2 A Probabilistic Simplification 21

2.3 Asymptotics 23

3The Proof of the Bruhat Order Lower Bound . 26

3.1 A Sufficient Condition for Bruhat Comparability 26

3.2 A Reduction to Uniforms 29

3.3 An Algorithm to Maximize the Bound 34

vii

4

An Extension to Chains in Bruhat Order

39

4.1 The General Setup . 39

4.2 A Tractable Necessary Condition 41

4.3 The Core Counting Problem 42

4.4 A Probabilistic Simplification 56

4.5 Asymptotics 57

5Some Properties of the Weak Ordering 62

5.1 A Criterion for Weak Comparability 62

5.2 Submultiplicativity of Pn* 65

6 The Proof of the Weak Order Upper Bound 67

6.1 A Necessary Condition for Weak Comparability 67

6.2 A Reduction to Uniforms 69

6.3 An Algorithm to Minimize the Bound 71

7The Proof of the Weak Order Lower Bound 75

7.1 A Formula for P(Ei(π)Ei(σ)) . . 75

7.2 Positive Correlation of the Events {Ei(π)Ei(σ)} 77

7.3 Linear Extensions of Arbitrary Finite Posets 81

8Numerics . 83

8.1 Bruhat Order Numerics 83

8.2 Weak Order Numerics 85

9On Infs and Sups in the Weak Order Lattice 89

9.1 A Connection with Complete, Directed, Acyclic Graphs 89

9.2 Computing Infs and Sups in the Weak Order Lattice 93

9.3 Some Equivalent Criteria for inf{π1 , . . . ' πr}=12 ... n 96

9.4 Submultiplicativity Again 99

9.5 Sharp Asymptotics of Pn(r) . 1Ot

9.5.1 Step 1: An Exact Formula for Pn(r) lai

9.5.2 Step 2: A Generating Function for Pn(r).. 105

9.5.3 Step 3: Asymptotics 108

Vlll

10

Open Problems

114

10.1 The Problems 114

Bibliography 115

IX

LIST OF FIGURES

FIGURE PAGE

1.1 The Bruhat order on S3 and S4. 2

1.2 Superimposing M(π) and M(σ) to form M(π, σ). 5

1.3 The permutation-induced poset 7(2143). 10

2.1 Finding a necessary condition for π<σ. 16

2.2 Selection of first m=m1 ( n/2-m resp.) rows (columns resp.) in

corners to support M(π). 21

3.1 G5 and an emboldened subgrid C. 27

3.2 Deletion of 2 largest elements of π, σ, and its affect on C. . 28

4.1 Finding a necessary condition for π1<... <πr. 43

4.2 A non-intersecting nest of 7 lattice paths. 46

4.3 An intersecting nest of 7 lattice paths. 47

4.4 "Swapping" the tails in Figure 4.3. 49

8.1 Experimental determination of the exponent -a in the asymptotic

equation Pncn-a. 85

8.2 Experimental determination of the ratio ρ in the asymptotic equation

Pn*cρn. 86

x

LIST OF TABLES

TABLE PAGE

3.1 Exact computation of Nk for smallish k. 37

6.1 Exact computation of Nk* for smallish k. 74

8.1 Computer simulation data for Pn. 84

8.2 Exact computation of Pn for smallish n. 84

8.3 Computer simulation data for Pn*. 87

8.4 Exact computation of Pn* for smallish n. 87

8.5 Our theoretical lower bound for Pn* applied for smallish n. 88

Xi

(JIAPTER 1

INTRODUCTION

In this chapter, we present the fundamental ideas necessary to understand the ma-

terial in subsequent chapters. We also state the main results to be proved in this

dissertation. It is assumed that the reader has some familiarity with probability the-

ory and combinatorics. A good graduate-level introduction to these subjects can be

found, for instance, in Billingsley [4] and Stanley's volumes on enumerative combina-

torics, [46] and [47].

1.1 Bruhat Order, Preliminaries

Let n>1 be an integer. Two permutations of [n] :={1, ..., n} are comparable in the

Bruhat order if one can be obtained from the other by a sequence of transpositions of

pairs of elements forming an inversion. Here is a precise definition of the Bruhat order

on the set of permutations Sn (see [46, p. 172, ex. 75. a.], Humphreys [30, p. 119]).

If ω =ω(1) ... ω(n) Sn, then a reduction of ω is a permutation obtained from ω

by interchanging some ω (i) with some ω (j) provided i<j and ω (i)>ω (j). We say

that π<σ in the Bruhat order if there is a chain σ=ω1ω2...ωs=π, where

each ωt is a reduction of ωt-1. The number of inversions in ωt strictly decreases with t.

1

4321

321

312

132

123

1234

Figure 1.1: The Bruhat order on S3 and S4.

Indeed, one can show that if ω2 is a reduction of ω1 via the interchange ω1(i)ω1(j),

i<j, then

inv(ω1)=inv(ω2)+2N(ω1)+1,

N(ω1):=|{k : i<k<j, ω1(i)>ω1(k)>ω1(j)}|;

here inv(ω1), say, is the number of inversions in ω1 (see Björner and Brenti [6]). Figure

1.1 illustrates this poset on S3 and S4. The Bruhat order notion can be extended

to other Coxeter groups (see Björner [5], Deodhar [20], and [6, p. 63] for historical

background), but we will be dealing with the symmetric group Sn only.

The definition of the Bruhat order is very transparent, and yet deciding for given π, a

2

whether π<σ from the definition is computationally difficult, even for smallish n.

Fortunately, there exist efficient algorithms for checking Bruhat comparability, which

can all be traced back to an algorithmic comparability criterion due to Ehresmann

(1934) [22] (see also Knuth [34], Björner and Brenti [6]). The Ehresmann "tableau

criterion" states that π<σ if and only if πi,j<σi,j for all 1<i<j<n -1, where

πi,j and σi,j are the i-th entry in the increasing rearrangement of π(1) , . . . , π(j)

and of a(1) , . . . ' σ(j). These arrangements form two staircase tableaux, hence the

term tableau criterion" For example, 41523 > 21534 is verified by element-wise

comparisons of the two tableaux

12451235

14 5 1 2 5

1412

4 2

Also, it is well-known that Ehresmann's criterion is equivalent to the (0, 1)-matrix

criterion. It involves comparing the number of Fs contained in certain submatrices

of the (0, 1)-permutation matrices representing π and a (see B6na[10], [6]). Later,

Björner and Brenti [7] were able to improve on the result of [22], giving a tableau

criterion that requires fewer operations. Very recently, Drake, Gerrish and Skandera

[21] have found two new comparability criteria, involving totally nonnegative poly-

nomials and the Schur functions respectively. We are aware of other criteria (see [5],

Fulton [25, pp. 173-177], Lascoux and Schiitzenberger [36], [20] ), but we found the

(0, 1)-matrix and Ehresmann criteria most amenable to probabilistic study.

3

The (0, 1)-matrix criterion for Bruhat order on Sn says that for π, σSn, π<σ if

and only if for all i, j<n, the number of π(1) , . . . , π(i) that are at most j exceeds

(or equals) the number of a(1), . . . ' a(i) that are at most j (see [11] for this version).

It is referred to as the (0, 1)-matrix criterion because of the following recasting of this

condition: let M (π), M(σ) be the permutation matrices corresponding to π, σ, so

that for instance the (i, j)-entry of M(π) is 1 if π(j)=i and 0 otherwise. Here, we

are labeling columns 1, 2, . . . ' n when reading from left to right, and rows are labeled

1, 2, . . . ' n when reading from bottom to top so that this interpretation is like placing

ones at points (i, π(i)) of the n ×n integer lattice and zeroes elsewhere. Denoting

submatrices of M (ċ) corresponding to rows I and columns J by M (ċ)I,J, this criterion

says that π<σ if and only if for all i, j<n, the number of ones in M (π)[i],[j] is at

least the number of ones in M (σ)[i],[j] (see [21] for this version).

An effective way of visualizing this criterion is to imagine the matrices M (π) and

M (σ) as being superimposed on one another into a single matrix, M (π, σ), with the

ones for M (π) represented by ×'s("crosses"), the ones for M (σ) by 's ("balls")

and the zeroes for both by empty entries. Note that some entries of a (π, σ) may

be occupied by both a cross and a ball. Then the (0, 1)-matrix criterion says that

π<σ if and only if every southwest submatrix of M (π, σ) contains at least as many

crosses as balls. Here, in the notation above, a southwest submatrix is a submatrix

M (π, σ)[i],[j] of M (π, σ) for some i, j<n. It is clear that we could also check

π<σ by checking that crosses are at least as numerous as balls in every northeast

submatrix of M (π, σ). Likewise, π<σ if and only if balls are at least as numerous

as crosses in every northwest submatrix of M (π, σ), or similarly balls are at least

4

M(0)

Figure 1.2: Superimposing M(π) and M(σ) to form M(π, σ).

as numerous as crosses in every southeast submatrix of M (π, σ). Parts of all four

of these equivalent conditions will be used in our proofs. As a quick example, with

π= 21534 and σ=41523, π<σ is checked by examining southwest submatrices of

M (π, σ) in Figure 1.2. Also, the superimposing of M (π) with M (σ) to form M (π, σ)

is illustrated in this figure.

1.2 Main Results Related to the Bruhat Order

In this dissertation, we use the (0, 1)-matrix and the Ehresmann criteria to obtain

upper and lower bounds for the number of pairs (π, σ) with π<σ.

5

Theorem 1.2.1. Let n>1 be an integer, and let π, σSn be selected independently

and uniformly at random. Then there exist universal constants c1, c2>0 such that

c1(0.708)n<P(π<σ)<c2/n2

Equivalently, the number of pairs (π, σ) with π<σ is sandwiched between the

counts c1(0.708)n(n!)2 and c2n-2(n!)2. The lower bound follows from a sufficient

condition derived from the (0, 1)-matrix criterion, and a computer-aided tabulation

of an attendant function of a smallish integer argument. Empirical estimates based on

generating pairs of random permutations suggest that P(π<σ) is of order n-(2+δ),

for δ close to 0.5 from above. So apparently it is the upper bound which comes close

to the true proportion P(π<σ). It is certain that the constant 0.708 can be further

improved, but we do not know if our method could be extended to deliver a lower

bound (1 -o(1))n. A lower bound n-a, a qualitative match of the upper bound, seems

out of sight presently.

A deeper insight reveals a more general result, related to chains of length r in Bruhat

order, once we realize some connections with MacMahon's formula [13] for counting

plane partitions contained in an r×s×t box. Without going into much unnecessary

detail here, one can visualize a plane partition as stacks of unit cubes pushed into a

corner. The k-th Ehresmann condition contains a clear connection between Bruhat

order on permutations and counting combinatorial objects related to plane partitions,

namely non-intersecting lattice paths, a notion we will make precise later on. A closer

look at our methods for permutation-pairs in the spirit of Gessel and Viennot's work

6

[26] implies an extension of Theorem 1.2.1, upper bound, from pairs of permutations

to r-tuples:

Theorem 1.2.2. Let π1, ..., πrSn be selected independently and uniformly at

random. Then there exists a uniform constant c >0 such that

P(π1<...<πr)<c/nr(r-1).

Note that this result implies that there are at most cn-r(r-1) (n!)r length r chains in

the Bruhat order poset.

1.3 Weak Order, Preliminaries

Then we turn to the modified order on Sn, the rveak order "-" Here πσ if there

is a chain σ=ω1ω2...ωs=π, where each ωt is a simple reduction of ωt-1,

i.e. obtained from ωt-1 by transposing two adjacent elements ωt-1(i), ωt-1(i+1) with

ωt-1(i)>ωt-1(i+1). Since at each step the number of inversions decreases by 1,

all chains connecting a and π have the same length. Alternatively, there is a simple

non-inversion (resp. inversion) set criterion, contained in Berge [3], we can use to

check πσ. Indeed, given ωSn introduce the set of non-inversions of ω:

E(ω)={(i, j) : i<j, ω-1(i)<ω-1(j)}.

Similarly, for ωSn we introduce the set of inversions of ω:

7

E*(ω)={(i, j) : i>j, ω-1(i)<ω-1(j)}.

Then, for given π, σSn, we have πσ if and only if E(π)E(σ) (equivalently

E*(π)E*(σ)). Note that ωSn is uniquely determined by its E(ω) (resp. its

E*(ω)).

It turns out that the poset (Sn,) is a lattice (see [3]). Indeed, given π1, ..., πrSn,

there is an efficient way to compute E(inf{π1, ..., πr}) (resp. E*(sup{π1 , . . . , πr}))

from the set i=1rE(πi) (resp. i=1rE " ( πi )). We will see precisely how to do this later.

1.4 Main Results Related to the Weak Order

We prove the following probabilistic result for weak order comparability:

Theorem 1.4.1. Let π, σSn be selected independently and uniformly at random,

and let Pn*:=P(πσ). Then Pn* is submultiplicative, i.e. Pn1+n2*<Pn1*Pn2*. Con-

sequently there exists ρ=limPn*n. Furthermore, there exists an absolute constant

c>0 such that

i=1n(H(i)/i)<Pn*<c(0.362)n

where H(i):=j=1i1/j. Consequently, ρ< 0.362.

The proof of the upper bound is parallel to that of Theorem 1.2.1, lower bound,

while the lower bound follows from the non-inversion (resp. inversion) set criterion

8

described last section. Empirical estimates indicate that ρ is close to 0.3. So here too,

as in Theorem 1.2.1, the upper bound seems to be qualitatively close to the actual

probability Pn*. And our lower bound, though superior to the trivial bound 1/n!, is

decreasing superexponentially fast with n, which makes us believe that there ought

to be a way to vastly improve it.

Paradoxically, it is the lower bound that required a deeper combinatorial insight.

Clearly the number of π's below (or equal to) a equals e(P), the total number of

linear extensions of 7=P(σ), the poset induced by σ. (The important notion of

7(σ) was brought to our attention by Sergey Fomin [24].) We prove that for any

poset 7 of cardinality n,

e(P)>n!/iPd(i), (1.1)

where d(i):=| { jP : j<i in 7} |. (This bound is an exact value of e(P) if the

Hasse diagram is a (directed) rooted tree, Knuth [34, sect. 5.1.4, ex. 20], or a forest

of such trees, Björner and Wachs [8].) The bound (1.1) for e(P(σ)) together with the

independence of sequential ranks in the uniform permutation were the key ingredients

in the proof of Theorem 1.4.1, lower bound.

Mikl6s B6na[12] has informed us that this general lower bound for e(P) had been

stated by Richard Stanley as a level 2 exercise in [46, p. 312, ex. 1] without a

solution. We have decided to keep the proof in the dissertation, since we could

not find a published proof anywhere either. The classic hook formula provides an

example of a poset 7 for which (1.1) is markedly below e(P). It remains to be seen

whether (1.1) can be strengthened in general, or at least for 7(σ). As an illustration,

9

1 2

Figure 1.3: The permutation-induced poset 7(2143).

7=P(2143) has the Hasse diagram appearing in Figure 1.3. Then e(P) =4, but

(1.1) delivers only

e(P)>24/9 e(P)>3.

Regarding the lattice properties of (Sn,), note that the identity permutation 12 ... n

is the unique minimum, and n(n-1) ...1 is the unique maximum. Let π1, ..., πrSn

be selected independently and uniformly at random. It is natural to ask: "How

likely is it that the infimum (resp. supremum) of {π1 , . . . ' πr} is the unique mini-

mal (resp. maximal) element in the weak order lattice?" Equivalently, what is the

asymptotic number of r-tuples (π1 , . . . ' πr) such that inf{π1, . . . , πr}=12 ... n (resp.

sup{π1, ..., πr}=n(n -1) ...1), n ? It turns out that the answer is the same

whether we consider infs or sups, which allows us to focus only on infimums. We

prove the following:

Theorem 1.4.2. Let Pn(r)=P (inf{π1, ..., πr}=12 ...n). Then

1. As a function of n, Pn(r) is submultiplicative. Hence, there exists

to

p(r)=limPn(r)n= inf Pn(r)n, r>1.

2. For each r>1, put hr(z)=j0(-1)j(j!)rzj and Hr(z)=(hr(z))-1 Then, letting

P0(r)=1, we have

Hr(z)=n0Pn(r)zn,

from which we obtain (Darboux theorem [2])

Pn(r)-1z*hr(z*)1(z*)n, n .

Here, z*=z*(r)(1, 2) is the unique simple root of hr(z)=0 in the disc

|z|<2. Consequently, p(r)=1/z*

Unlike our results about comparability in the Bruhat and weak orderings, here we

have a sharp asymptotic formula. The key to the proof of this theorem is establishing

the exact formula

Pn(r)=k=0n-1(-1)k b1+...+bn-k=nb1,...,bn-k1,(b1!)r...(bn-k!)r1,

which follows from the principle of inclusion-exclusion. This formula for Pn(r) is, in

some sense, an r-analog of that for the Eulerian numbers ( B6na[11], Knuth [34]).

Indeed, it turns out that Pn(r) is the probability that the uniform, independent permu-

tations π1-1, ..., πr-1 have no common descents. Introduce the random variable Sn(r),

11

the number of these common descents, so that Pn(r)=P(Sn(r)=0). Another natural

question here is:

Problem. What is the limiting distribution of Sn(r)l.?

We believe that the answer here is "Gaussian", as it is in the case of the number

of descents in a single uniformly random permutation (Sachkov [44]). Our feeling

is that the proof will involve use of the bivariate generating function Fr(x, y)=

n1xnE[(1+y)Sn(r)], which we prove has the simple form

Fr(x, y)=xfr(xy)1-xfr(xy), fr(z).=j0zj(j+1)!r .

Interestingly, this generating function is a special case of a more general result proved

by Richard Stanley [45], although he was probably unaware of the connections his

work had with the weak ordering.

In conclusion we mention several papers that are in the same spirit of this dissertation.

First, [39] and [40] (both by Boris Pittel), where the "probability-of-comparability"

problems were solved for the poset of integer partitions of n under dominance order,

and for the poset of set partitions of [n] ordered by refinement. Also, [41] (again

by B. Pittel), where the "inf