Increasing subsequences of random walks

# Increasing subsequences of random walks

Omer Angel University of British Columbia, Vancouver BC, V6T 1Z2, Canada. E-mail: angel@math.ubc.ca. Supported in part by NSERC.    Richárd Balka University of British Columbia, and Pacific Institute for the Mathematical Sciences, Vancouver BC, V6T 1Z2, Canada. E-mail: balka@math.ubc.ca. Former affiliations: Department of Mathematics, University of Washington, and Alfréd Rényi Institute of Mathematics, Hungarian Academy of Sciences. Supported by the National Research, Development and Innovation Office-NKFIH, 104178.    Yuval Peres Microsoft Research, 1 Microsoft Way, Redmond, WA 98052, USA. E-mail: peres@microsoft.com.
###### Abstract

Given a sequence of real numbers , we consider the longest weakly increasing subsequence, namely with and maximal. When the elements are i.i.d. uniform random variables, Vershik and Kerov, and Logan and Shepp proved that .

We consider the case when is a random walk on with increments of mean zero and finite (positive) variance. In this case, it is well known (e.g., using record times) that the length of the longest increasing subsequence satisfies . Our main result is an upper bound , establishing the leading asymptotic behavior. If is a simple random walk on , we improve the lower bound by showing that .

We also show that if is a simple random walk in , then there is a subsequence of of expected length at least that is increasing in each coordinate. The above one-dimensional result yields an upper bound of . The problem of determining the correct exponent remains open.

2010 Mathematics Subject Classification: 60G17, 60G50.Keywords: random walk, restriction, monotone, increasing subsequence.

## 1 Introduction

For a function , its restriction to a subset of its domain is denoted . We say that is increasing if for all with . Define

 LIS(S|[0,n))=max{|A|:A⊂[0,n),S|A is % increasing}.

The main goal of this paper is to investigate when is a random walk. The simple random walk is the most natural case, but our results apply to walks with steps of mean zero and finite (positive) variance, that is, such that is an i.i.d. sequence with and . By normalising we may clearly assume that . We say that is the simple random walk if .

The famous Erdős-Szekeres Theorem  implies that must contain either an increasing or a decreasing subsequence of length at least . This is sharp for general sequences, and it is easy to see that there are even -step simple walks on for which the longest increasing subsequence has length of order . By symmetry, increasing and decreasing subsequences have the same length distribution, but this does not immediately imply that a similar bound holds in high probability.

In random settings, there have been extensive studies of the longest increasing subsequence in a uniformly random permutation initiated by Ulam . This is easily equivalent also to the case of a sequence of i.i.d. (non-atomic) random variables. A rich theory rose from the study of this question, which is closely related to last passage percolation and other models. It was proved by Vershik and Kerov  and by Logan and Shepp  that and in probability as . In this case, much more is known. Baik, Deift and Johansson  proved that the fluctuations of scaled by converge to the Tracy-Widom distribution, first arising in the study of the Gaussian Unitary Ensemble. We refer the reader to Romik’s book  for an excellent survey of this problem.

On the other hand, it appears that this problem has not been studied so far even for a simple random walk . The expected length of the longest strictly increasing subsequence of is at most the expected size of the range , hence is . Thus we consider (weakly) increasing subsequences. Taking the set of record times, or alternatively the zero set of both yield increasing subsequences of expected length . It is not immediate how to do any better. The largest level set of still has size . See Figure 1 for the longest increasing subsequence in one random walk instance. Note that the set of record times yields a similar lower bound for a general random walk with mean zero and finite variance.

On some reflection, one finds a number of arguments that yield the weaker bound for a simple random walk . For example, first one can show that with high probability, in every interval , each value is visited at most times. Assume is such that is increasing. For each define the interval , where is the first (and is the last) time such that . By monotonicity the intervals are disjoint. The length of the subsequence is then bounded by , where the number of intervals is at most . As with high probability, the Cauchy-Schwarz inequality gives the upper bound . However, going beyond the exponent requires more delicate arguments even in the case of a simple random walk.

Although Ulam’s problem and the problem of this paper are superficially similar, it turns out that their structure is different, and finding monotone subsequences in random walks is more closely related in nature to restriction theorems for continuous functions. In the continuous setting the upper Minkowski dimension plays the role of counting. Balka and Peres  showed that a Brownian motion is not monotone on any set of Hausdorff dimension greater than . However, working with the Minkowski dimension requires understanding the structure of a set at a specific scale, which is not needed for the Hausdorff dimension. The proof of  is based on Kaufman’s uniform dimension doubling theorem for two-dimensional Brownian motion. A key fact used there is that a closed set in of Hausdorff dimension strictly greater than 1/2 intersects the zero set of a Brownian motion with positive probability. Hausdorff dimension cannot be replaced here by upper Minkowski dimension. Therefore, the methods developed in  are not powerful enough to prove the leading term of the upper bound even in the case of a simple random walk. However, Balka, Máthé and Peres analyzed the behaviour of some self-affine functions motivated by  and questions of Kahane and Katznelson , which had a large impact on this paper. Later on, Máthé and the authors of this paper  proved restriction theorems for fractional Brownian motion using similar methods to this paper’s ones, and the case of self-affine functions was also handled in . In particular, they proved that a Brownian motion is not monotone on any set of upper Minkowski dimension greater than . Hence the interaction between this paper and  played an important role in both directions. For other restriction theorems in the continuous, deterministic case, and for generic functions in the sense of Baire category see Elekes , Kahane and Katznelson , and Máthé .

The main goal of this paper is to prove the following theorem.

###### Theorem 1.

Let be a random walk with i.i.d. steps satisfying and . For all large enough, for all we have

 P(LIS(S|[0,n))≥ℓ√n25√[b]log4n(log4log4n))≤(2log4n)2ℓ.

Moreover, if for some then we have

 P(LIS(S|[0,n))≥ℓ√n25√[b]log4n(log4log4n))≤n−ℓ/2.

In the following corollaries let be a random walk as in Theorem 1.

###### Corollary 1.1.

For all with probability we have

 LIS(S|[0,n))≤n1/2+ε.
###### Corollary 1.2.

For all and large enough

 ELIS(S|[0,n))≤n1/2+ε.

In the other direction, we show that in the case of a simple random walk, with high probability there are increasing subsequences somewhat longer than the trivially found ones.

###### Theorem 2.

Let be a simple random walk. For every for all large enough

 P(LIS(S|[0,n))<ε√nlog2n)≤250ε.

Consequently, for all large enough we have

 E(LIS(S|[0,n)))≥(1/1000)√nlog2n.

In Section 4 we consider higher dimensional random walks. Let and let . We say that is increasing on a set if all the coordinate functions of are non-decreasing, i.e.  is increasing with respect to the coordinate-wise partial order on . Generalizing , we define

 LIS(S|[0,n))=max{|A|:A⊂[0,n),S|A is % increasing}.

Since the restriction of a random walk to a single coordinate is again a random walk, if is a -dimensional random walk with mean and bounded second moment then Corollary 1.1 implies that with probability . For a large class of two-dimensional random walks we are able to prove a lower bound as well. However, the problem of determining the correct exponent remains open.

###### Theorem 3.

Let be a two-dimensional random walk with steps for which

• the mean is the zero vector,

• the covariance matrix is the identity matrix,

• the coordinates of have finite moments for some .

Then there is a constant such that for eevery and

 P(LIS(S|[0,n))<εn1/3)≤cε.

Consequently, for all we have

 ELIS(S|[0,n))≥14cn1/3.

Finally, in Section 5 we state some open questions.

### Acknowledgments

OA thanks the organizers of the probability, combinatorics and geometry meeting at the Bellairs Institute, as well as several participants, in particular Simon Griffiths who proposed this problem, and Louigi Addario-Berry, Guillaume Chapuy, Luc Devroye, Gábor Lugosi and Neil Olver for useful discussions. The present collaboration took place mainly during visits of OA and RB to Microsoft Research. We are indebted to András Máthé and Boris Solomyak for useful suggestions.

## 2 Upper bound

To simplify notations in the proof, it is convenient to assume the length of the random walk is a power of . By monotonicity in of , we can interpolate for all other . Our main result Theorem 1 follows from the following by this monotonicity and the substitution .

###### Theorem 2.1.

Let be a random walk with and . For all large enough for all we have

 P(LIS(S|[0,4n))≥ℓ2n+4√nlog2n)≤(2n2)ℓ.

Moreover, if for some then we have

 P(LIS(S|[0,4n))≥ℓ2n+4√nlog2n)≤2−ℓn.

The main goal of this section is to prove Theorem 2.1. The key is a multi-scale argument, the time up to is split into intervals. We consider the number of these intervals that intersect our set , as well as the sizes of intersections. Repeating this allows us to get (inductively) better and better bounds. The dependence on the randomness of the walk is done through some estimates on the local time, which we derive in the following subsection.

Throughout this section, fix a random walk with and . Various constants below depend only on the law of . We will use the following theorems in this section.

###### Theorem 2.2 (Petrov, ).

There is a constant such that for all and we have

 supx∈RP(x≤S(n)≤x+λ)≤cλ+1√n.

For the following theorems see [8, Thm. A.2.5] and its corollaries.

###### Theorem 2.3.

For all and we have

 P(max0≤i≤n|S(i)|≥λ√n)≤1λ2.
###### Theorem 2.4.

Assume that for some . Then there is a constant such that for all and we have

 P(max0≤i≤n|S(i)|≥λ√n)≤e−cλ2.

### 2.1 Scaled local time estimates

###### Definition 2.5.

Let and . A time interval of order is of the form

 Im,p=[p4m,(p+1)4m)⊂N.

A value interval of order is of the form

 Jm,q=[q2m,(q+1)2m)⊂R.

Note that a time interval is a subset of , while a value interval is a real interval. For all let be the set of time intervals of order contained in . Clearly .

###### Definition 2.6.

The scaled local time is the number of order intervals in in which takes at least one value in :

 Sm,k,p,q=∣∣{I∈Im,k,p:∃x∈I, S(x)∈Jm−k,q}∣∣.

Our intermediate goal is to prove the following uniform estimate on scaled local times.

###### Proposition 2.7.

There is a such that for all large enough

 P(Sm,k,p,q≤γn2k for all k≤m≤n,p<4n−m,q∈Z)≥1−2−(n+1).

We begin with an estimate on the expectation of a single scaled local time. Let and denote the expectation and probability for a random walk started at .

###### Lemma 2.8.

For some absolute constant and for all we have

 ExSm,k,p,q≤c2k.
###### Proof.

The proof uses the idea that conditioned on the event that some time interval contributes to , with probability bounded from , the random walk is still nearby at the end of the time interval.

The strong Markov property of the walk and translation invariance imply that it is enough to consider the case . Since the central limit theorem yields , there exist and such that for all we have

 P0(0≤S(j)≤N√j)≥c1.

For each let be the event that and let be the event

that counts towards . By the above inequality, for all we have

 Px(Ai|Bi)≥c1.

By Theorem 2.2 there exists a constant such that for all and we have

 Px(Ai)≤c2(N+1)2m−k√i4m−k=c2(N+1)√i.

The above two inequalities imply that

 ExSm,k,0,0=4k∑i=1Px(Bi)≤4k∑i=1c2(N+1)c1√i≤c2k,

where . The proof is complete. ∎

Next we estimate the tail of a single scaled local time.

###### Lemma 2.9.

There is an absolute constant such that for all and we have

 Px(Sm,k,p,q≥Cℓ2k)≤2−ℓ.
###### Proof.

Let , where is the constant of Lemma 2.8 and denotes rounding up. By Markov’s inequality we have , establishing the claim for . We proceed inductively: Assume that the claim holds for some . Observe the walk starting at time either until we reach time or until sub-intervals of order contribute to . The latter happens with probability at most . By the strong Markov property the conditional probability that there are additional sub-intervals contributing to is at most , proving the claim for . ∎

###### Proof of Proposition 2.7.

Let , where is the constant of Lemma 2.9. We apply Lemma 2.9 with to each of the relevant . Since , there are choices for each of and . As , there are at most options for . If then with have scaled local time . This is likely, as Theorem 2.3 yields that . These imply that for all large enough we have

 P(∃m,k,p,q, s.t. Sm,k,p,q>7Cn2k) ≤(n+1)24n22n+22−7n+2−2n ≤2⋅2−2n≤2−(n+1).

Clearly we may also require . ∎

### 2.2 No long increasing subsequence

Next, we use Proposition 2.7 to rule out the existence of very long increasing subsequences in the random walk. We need the following definition.

###### Definition 2.10.

Let be a function and let be a finite set such that . The variation of restricted to is defined as

 V1(S|A)=k−1∑i=1|S(ai+1)−S(ai)|.

Note that if is increasing then equals the diameter of . The upper bound of Theorem 2.1 follows from the following proposition.

###### Proposition 2.11.

Fix with . Assume that a walk is such that

1. the event of Proposition 2.7 occurs,

2. .

Then we have

 LIS(S|[0,4n))≤(γn2b+1)a+1.
###### Proof.

Let be a set such that is increasing. For let

 Dℓ={I∈Iab,ℓb,0:I∩A≠∅}% anddℓ=|Dℓ|(γn2b+1)ℓ

be the set of intervals of order that intersect , and its size with a convenient normalization. Clearly and . In order to prove the claim we prove inductively bounds on .

Let and index the elements of , and suppose that interval contains intervals in , so that . By assumption (1) for all we have that is visited in at most sub-intervals of . It follows that if then must visit at least value intervals of order . The diameter of the union of these visited intervals is at least . This leads to a variation bound

 V1(S|A∩Ii)≥(piγn2b−2)2(a−ℓ)b. (2.1)

Assumption (2) yields that . Thus

 |Dℓ−1|∑i=1V1(S|A∩Ii)≤γn2n. (2.2)

Inequalities (2.1) and (2.2) and imply that

 γn2n≥|Dℓ−1|∑i=1(piγn2b−2)2(a−ℓ)b=2n−ℓb⎛⎝|Dℓ−1|∑i=1piγn2b−2|Dℓ−1|⎞⎠,

and therefore

 |Dℓ−1|∑i=1piγn2b−2|Dℓ−1|≤γn2ℓb. (2.3)

Using and dividing (2.3) by yields

 dℓ−dℓ−1≤γn2b2(2γn)ℓ−1≤γn2b−ℓ,

where we have used that . As , the above inequality implies

 dℓ≤1+ℓ∑i=1γn2b−i≤γn2b+1

for every . In particular we get for

Finally, we use Proposition 2.11 to derive an estimate on the likelihood of long increasing subsequences in a random walk.

###### Proof of Theorem 2.1.

First we prove the theorem for . Let

 a=⌈√n/log2n⌉andb=⌈√nlog2n⌉,

where denotes rounding up. Note that . We consider up to time . For large enough with probability the event of Proposition 2.7 occurs for . Moreover, Theorem 2.3 implies that

 P(maxi≤4ab{|S(i)|}≥ab2ab)≤1(ab)2≤1n2.

Thus with probability at least the conditions and conclusion of Proposition 2.11 hold for .

Suppose additionally that for some . Then Theorem 2.4 yields that for some constant and for all large enough

 P(maxi≤4ab{|S(i)|}≥ab2ab)≤e−c(ab)2≤2−(n+1).

Thus with probability at least the conditions and conclusion of Proposition 2.11 hold for .

Let be such that Proposition 2.11 holds for . Since is increasing in , we obtain that for large enough

 LIS(S|[0,4n))≤LIS(S|[0,4ab)) ≤(γab2b+1)a+1=2ab2b(2γab)a+1 ≤2n+3√nlog2n+O(√n/log2n) <2n+4√nlog2n−1.

This proves Theorem 2.1 if . For the general case fix , it is enough to prove that for all we have

 (2.4)

then setting concludes the proof, where denotes rounding down. Let . If is already defined then let be the minimal integer so that . Since increases by at most when incrementing , we actually have . By the strong Markov property at , we see that are i.i.d. copies of . However, requires for all , with probability at most . This implies (2.4), and the proof is complete. ∎

## 3 Lower bound for a simple random walk

The goal of this section is to prove Theorem 2. For simplicity, we present our argument only for the simple random walk on . However, it seems that the argument should extend with minor changes to any random walk with bounded integer steps of 0 mean, and finite variance. The construction relies on values appearing multiple times in the walk, and fails more fundamentally if the walk is not supported on multiples of some .

###### Definition 3.1.

Let denote the hitting time of by the simple random walk. Let be the 2-order of , that is, the number of times it is divisible by .

###### Lemma 3.2.

Consider a simple random walk from conditioned to hit before returning to , and stopped when it reaches . Let be the times of the first and last visits to . Then:

1. The number of visits to is geometric with mean .

2. The walk on is a walk conditioned to hit before returning to , and stopped when it reaches .

3. The walk on is a walk from conditioned to hit without returning to , and stopped when it reaches .

4. The two sub-walks and the geometric variable are independent.

###### Proof.

In order to prove the first statement we first consider a simple random walk from up to the time when it reaches either or . This walk has probability of returning to without hitting , at which time another excursion from begins. Therefore the number of visits to on is geometric with mean . Moreover, the number of visits to is independent of whether the walk hits or , so when conditioning on hitting the distribution is still geometric with mean , which proves the first statement.

Now we return to our original walk from . Excursions from either return to , or hit , or hit . The partition into excursions around does not give any information on the trajectory within each excursion, except for its type, and the other claims follow. ∎

###### Lemma 3.3.

Let be a simple random walk. For all we have

 ELIS(S|[0,τ2n))≥n2n−1,

and for any ,

 P(LIS(S|[0,τ2n))<(1−ε)n2n−1)≤2ε2n2.
###### Proof.

We construct an increasing subsequence of as follows. Informally, we take some times to be in our index set, greedily in decreasing order of the 2-order of .

For each integer we construct an interval . The intervals are such that if then . Given such intervals, we have that is increasing along , where

 A=2n−1⋃x=1{i∈Ix:S(i)=x}.

We start by setting and , where is the last visit to before . Let and be such that and assume by induction that are already defined for all for which . Now we define . Let and , then clearly . Thus and are already defined by the inductive hypothesis. Let , where is the first hitting time of after and is the time of the last visit to before . See Figure 2 for an example. Figure 2: The increasing subsequence constructed for Lemma 3.3 in a simple random walk stopped at 16. All visits to 8 are used, then all compatible visits to 4,12, followed by 2,6,10,14 and a single visit to each odd value. There exist longer subsequences of length 42 in this case.

Assuming , we show that the law of restricted to is that of a simple random walk started at conditioned to hit before returning to , and stopped when hitting . This is seen inductively using Lemma 3.2, and since the walk after the last visit to before cannot return to .

From the above, we deduce that for all with the number of visits to in is geometric with mean , and these geometric variables are all independent. Since for each there are integers with , we get that

 E|A|=n−1∑k=02n−k−12k=n2n−1.

As any geometric satisfies and our geometric random variables are all independent, we obtain that

 Var|A|=n−1∑k=02n−k−12k(2k−1)≤n−1∑k=02n+k−1≤22n−1.

The second claim now follows by Chebyshev’s inequality. ∎

###### Proof of Theorem 2.

Fix . For large enough let be an integer such that

 19m2m≤ε√nlog2n<14m2m. (3.1)

Then we have that

 P(LIS(S|[0,n))<ε√nlog2n)≤P(LIS(S|[0,τ2m))<12m2m−1)+P(τ2m>n). (3.2)

Applying Lemma 3.3 for this we get

 P(LIS(S|[0,τ2m))<12m2m−1)≤8m2.

Moreover, [9, Thm. 2.17] and (3.1) imply that

 P(τ2m>n)≤12⋅2m√n≤108εlog2nm.

Since and as , plugging the previous bounds in (3.2) gives for large enough

 P(LIS(S|[0,n))<ε√nlog2n)≤8m2+108εlog2nm≤250ε.

Finally, applying the above inequality for implies that

 E(LIS(S|[0,n)))≥√nlog2n1000.\qed

## 4 Random walks in higher dimensions

The main goal of this section is to prove Theorem 3. As noted, the upper bound in the one-dimensional case holds trivially in every dimension. For sequences we use the notation if as . The lower bound is based on the following estimate by Denisov and Wachtel, see [4, Example 2] and see there the history of similar estimates for Brownian motion and random walks.

###### Theorem 4.1.

Let be a two-dimensional random walk satisfying the conditions of Theorem 3. Let be the hitting time of the positive quadrant: . Then there is some so that

 P(τ>n)∼cn−1/3.

More generally, for a higher dimensional random walk , define the hitting time

 τ=inf{n>0:S(n)∈Rd+}.

Denisov and Wachtel [4, Theorem 1] proved that for some and , where is the exponent corresponding to Brownian motion staying outside a quadrant up to time (assuming again that the walk is normalized so that and , and that for some ). Consequently, the following lemma completes the proof of Theorem 3 (with ), and gives a similar lower bound for random walks in higher dimensions.

###### Lemma 4.2.

Let be a random walk in , and let be such that

 P(τ>n)=O(n−α).

Then there is a constant such that for all and

 P(LIS(S|[0,n))<εnα)≤cε. (4.1)

Consequently, for all we have

 ELIS(S|[0,n))≥14cnα.
###### Proof.

Fix . Define the greedy increasing subsequence with time indices given by the recursion

 a0=0,ai+1=min{a>ai:S(a)−S(ai)∈Rd+}.

Setting , we see that if then . This gives a set with i.i.d. increments with the law of .

Choose such that for all

 P(τ>n)≤c1n−α,

and define the truncated variables . Then

 P(∃i≤kn:Xi≠Yi)≤knP(τ>n)≤c1ε. (4.2)

The also form an i.i.d. sequence and satisfy

 EYi≤n∑m=0P(τ>m)≤c2n1−α,

where depends only on and . By Markov’s inequality,

 P(kn∑i=1Yi≥n)≤knEY1n≤c2ε. (4.3)

Combining (4.2) and (4.3) we obtain

 P(kn∑i=1Xi≥n) ≤P(∃i≤kn:Xi≠Yi)+P(kn∑i=1Yi≥n) ≤(c1+c2)ε.

As noted, this is a bound on . Hence (4.1) holds with . Applying (4.1) with yields the second claim. ∎

## 5 Open Questions

There are many potential extensions of our results. Two central open problems are to reduce the gap between the lower and upper bounds in dimension one, and to determine the right order of magnitude in higher dimensions. Moreover, our lower bound in Theorem 2 is specific to the simple random walks, and our proof does not work for general random walks.

###### Question 5.1.

Let be a random walk with zero mean and finite (positive) variance. Is there a constant such that, with probability ,

 LIS(S|[0,n))≤√nlogan?

Does this upper bound hold at least when is a simple random walk?

###### Question 5.2.

Let and let be a -dimensional simple random walk. What is the order of magnitude of