Quickest Online Selection of an Increasing Subsequence of Specified Size*


 Curtis Townsend
 6 months ago
 Views:
Transcription
1 Quickest Online Selection of an Increasing Subsequence of Specified Size* Alessandro Arlotto, Elchanan Mossel, 2 J. Michael Steele 3 The Fuqua School of Business, Duke University, 00 Fuqua Drive, Durham, North Carolina 27708; 2 Wharton School, Department of Statistics, Huntsman Hall 459, University of Pennsylvania, Philadelphia, Pennsylvania 904 and Department of Statistics, University of California, Berkeley, 40 Evans Hall, Berkeley, California 94720; 3 Wharton School, Department of Statistics, Huntsman Hall 447, University of Pennsylvania, Philadelphia, Pennsylvania 904; Received 26 December 204; accepted 5 May 205 Published online 5 February 206 in Wiley Online Library (wileyonlinelibrary.com). DOI 0.002/rsa ABSTRACT: Given a sequence of independent random variables with a common continuous distribution, we consider the online decision problem where one seeks to minimize the expected value of the time that is needed to complete the selection of a monotone increasing subsequence of a prespecified length n. This problem is dual to some online decision problems that have been considered earlier, and this dual problem has some notable advantages. In particular, the recursions and equations of optimality lead with relative ease to asymptotic formulas for mean and variance of the minimal selection time. 206 Wiley Periodicals, Inc. Random Struct. Alg., 49, , 206 Keywords: increasing subsequence problem; online selection; sequential selection; timefocused decision problem; dynamic programming; Markov decision problem. INCREASING SUBSEQUENCES AND TIMEFOCUSED SELECTION If X, X 2,... is a sequence of independent random variables with a common continuous distribution F, then L n = max{k : X i < X i2 < < X ik, where i < i 2 < < i k n} Correspondence to: J. Michael Steele *E.M. acknowledges the support of NSF grants DMS and CCF 32005, ONR grant number N and grant from the Simons Foundation. 206 Wiley Periodicals, Inc. 235
2 236 ARLOTTO, MOSSEL, AND STEELE represents the length of the longest monotone increasing subsequence in the sample {X, X 2,..., X n }. This random variable has been subject to a remarkable sequence of investigations beginning with Ulam [7] and culminating with Bail et al. [5] where it was found that n /6 (L n 2 n) converges in distribution to the TracyWidom law, a new universal law that was introduced just a few years earlier in Tracy and Widom [6]. The review of Aldous and Diaconis [] and the monograph of Romik [4] draw rich connections between this increasing subsequence problem and topics as diverse as card sorting, triangulation of Riemann surfaces, and especially the theory of integer partitions. Here we consider a new kind of decision problem where one seeks to select as quickly as possible an increasing subsequence of a prespecified length n. More precisely, at time i, when the decision maker is first presented with X i, a decision must be made either to accept X i a member of the selected subsequence or else to reject X i forever. The decision at time i is assumed to be a deterministic function of the observations {X, X 2,..., X i },so the times τ <τ 2 < <τ n of affirmative selections give us a strictly increasing sequence stopping times that are adapted to the sequence of σ fields F i = σ {X, X 2,..., X i }, i <. The quantity of most interest is β(n) := min E[τ n ] () π where the minimum is over all sequences π = (τ, τ 2,..., τ n ) of stopping times such that τ <τ 2 < <τ n and X τ < X τ2 < < X τ n. Such a sequence π will be called a selection policy, and the set of all such selection policies with E[τ n ] < will be denoted by (n). It is useful to note that the value of β(n) is not changed if we replace each X i with F (X i ), so we may as well assume from the beginning that the X i s are all uniformly distributed on [0, ]. Our main results concern the behavior of n β(n) and the structure of the policy that attains the minimum (). Theorem. The function n β(n) = min π (n) E[τ n] is convex, β() =, and for all n 2 one has the bounds 2 n2 β(n) 2 n2 + n log n. (2) We can also add some precision to this result at no additional cost. We just need to establish that there is always an optimal policy within the subclass of threshold policies π = (τ, τ 2,..., τ n ) (n) that are determined by real valued sequences {t i [0, ] : i n} and the corresponding recursions τ k+ = min{i >τ k : X i [X τk, X τk + t n k ( X τk )]}, 0 k < n. (3) Here, the recursion begins with τ 0 = 0 and X 0 = 0, and one can think of t n k as the threshold parameter that specifies the maximum fraction that one would be willing to spend from a residual budget ( X τk ) in order to accept a value that arrives after the time τ k when the k th selection was made. Random Structures and Algorithms DOI 0.002/rsa
3 QUICKEST ONLINE SELECTION OF AN INCREASING SUBSEQUENCE 237 Theorem 2. one has There is a unique threshold policy π = (τ, τ 2,..., τ n ) (n) for which β(n) = min E[τ n]=e[τ n ], (4) π (n) and for this optimal policy π one has for all α>2that Var[τ n ]= 3 n3 + O(n 2 log α n) as n. (5) In the next section, we prove the existence and uniqueness of an optimal threshold policy, and in Section 3 we complete the proof of Theorem after deriving some recursions that permit the exact computation of the optimal threshold values. Section 4 deals with the asymptotics of the variance and completes the proof of Theorem 2. In Sections 5 and 6, we develop the relationship between Theorems and 2 and the more traditional sizefocused online selection problem which was first studied in Samuels and Steele [5] and then studied much more extensively by Bruss and Robertson [9], Gnedin [], Bruss and Delbaen [7,8], and Arlotto et al. [3]. On an intuitive level, the timefocused selection problem and the sizefocused selection problems are dual to each other, and it is curious to consider the extent to which rigorous relationships that can be developed between the two. Finally, in Section 7 we underscore some open issues, and, in particular, we note that there are several other selection problems where it may be beneficial to explore the possibility of time versus size duality. 2. THRESHOLD POLICIES: EXISTENCE AND OPTIMALITY A beneficial feature of the timefocused monotone selection problem is that there is a natural similarity between the problems of size n and size n. This similarity entails a scaled regeneration and leads to a useful recursion for β(n). Lemma (Variational Beta Recursion). For all n =, 2,...we have [ β(n) = inf E τ + ] β(n ), (6) τ X τ where the minimum is over all stopping times τ and where we initialize the recursion by setting β(0) = 0. Proof. Suppose that the first value X τ has been selected. When we condition on the value of X τ, we then face a problem that amounts to the selection of n values from a thinned sequence of uniformly distributed random variables on [X τ,]. The thinned sequence may be obtained from the original sequence of uniform random variables on [0, ] by eliminating those elements that fall in [0, X τ ]. The remaining thinned observations are then separated by geometric time gaps that have mean /( X τ ), so conditional on X τ, the optimal expected time needed for the remaining selections is β(n )/( X τ ). Optimization over τ then yields (6). This quick proof of Lemma is essentially complete, but there is a variation on this argument that adds further information and perspective. If we first note that β() =, then Random Structures and Algorithms DOI 0.002/rsa
4 238 ARLOTTO, MOSSEL, AND STEELE one can confirm (6) for n = simply by taking τ =. One can then argue by induction. In particular, we take n 2 and consider a arbitrary selection policy π = (τ, τ 2,..., τ n ). If we set π = (τ 2 τ, τ 3 τ,..., τ n τ ), then one can view π as a selection policy for the sequence (X, X 2,...) = (X +τ, X 2+τ,...) where one can only make selections from those values that fall in the interval [X τ,]. As before, remaining thinned observations are then separated by geometric time gaps that have mean /( X τ ). Thus, conditional on X τ, the optimal expected time needed for the remaining selections is β(n )/( X τ ). By the putative suboptimality of the selection times τ, τ 2,..., τ n one has β(n ) ( X τ ) E[τ n τ τ, X τ ]= τ + E[τ n τ, X τ ]. Moreover, we note that the sum β(n ) τ + ( X τ ) is the conditional expectation given τ of a strategy that first selects X τ and then proceeds optimally with the selection of n further values. The suboptimality of this strategy and the further suboptimality of the policy π = (τ, τ 2,..., τ n ) give us β(n ) β(n) E[τ + ( X τ ) ] E[τ n]. (7) The bounds of (7) may now look obvious, but since they are valid for any strategy one gains at least a bit of new information. At a minimum, one gets (6) immediately just by taking the infimum in (7) over all π in (n). Thus, in a sense, the bound (7) generalizes the recursion (6) which has several uses. In particular, (6) helps one to show that there is a unique threshold policy that achieves the minimal expectation β(n). Lemma 2 (Existence and Uniqueness of an Optimal Threshold Policy). For i n, there are constants 0 t i such that the threshold policy π (n) defined by (3) is the unique optimal policy. That is, for π = (τ, τ 2,..., τ n ) one has β(n) = min E[τ n]=e[τ n ], π (n) and π is the only policy in (n) that achieves this minimum. Proof. The proof again proceeds by induction. The case n = is trivial since the only optimal policy is to take any element which presented. This corresponds to the threshold policy with t = and β() =. For the moment, we consider an arbitrary policy π = (τ, τ 2,..., τ n ) (n). Wehave E[τ ] <, and we introduce a parameter t by setting t = (E[τ ]). Next, we define a new, threshold stopping time τ by setting τ = min{i : X i < t}, and we note that this construction gives us E[τ ]=E[τ ]=/t. Fors [0, t], we also have the trivial inequality (X τ < s) τ i= Random Structures and Algorithms DOI 0.002/rsa (X i < s),
5 QUICKEST ONLINE SELECTION OF AN INCREASING SUBSEQUENCE 239 so by Wald s equation we also have P(X τ < s) E[ τ i= (X i < s)] =se[τ ]=s/t. (8) The definition of τ implies that X τ is uniformly distributed on [0, t], so we further have P(X τ < s) = min{, s/t}. Comparison with (8) then gives us the domination relation P(X τ < s) P(X τ < s) for all 0 s. (9) From (9) and the monotonicity of x ( x), we have by integration that β(n ) β(n ) E[ )] E[ ], (0) X τ X τ and one has a strict inequality in (9) and (0) unless τ = τ with probability one. If we now add E[τ ]=E[τ ] to the corresponding sides of (0) and take the infimum over all τ, then the beta recursion (6) gives us [ β(n ) ] E[τ ]+E ) inf X τ τ { β(n ) E[τ ]+E[ ] X τ } = β(n). In other words, the first selection of an optimal policy is given by uniquely by a threshold rule. To see that all subsequent selections must be made by threshold rules, we just need to note that given the time τ and value X τ = x of the first selection, one is left with a selection problem of size n from the smaller set {X i : i >τ and X i > x}. The induction hypothesis applies to this problem of size n, so we conclude that there is a unique threshold policy (τ2, τ 3,..., τ n ) that is optimal for these selections. Taken as a whole, we have a unique threshold policy (τ, τ 2,..., τ n ) (n) for the problem of selecting an increasing subsequence of size n in minimal time. Lemma 2 completes the proof of the first assertion (4) of Theorem 2. After we develop a little more information on the behavior of the mean, we will return to the proof of the second assertion (5) of Theorem LOWER AND UPPER BOUNDS FOR THE MEAN The recursion (6) for β(n) is informative, but to determine its asymptotic behavior we need more concrete and more structured recursions. The key relations are summarized in the next lemma. Lemma 3 (Recursions for β(n) and the Optimal Thresholds). t (0, ) we let g(x, t) = t + x t log ( t Random Structures and Algorithms DOI 0.002/rsa For each x and ), G(x) = min g(x, t), and () 0<t< H(x) = arg min g(x, t). (2) 0<t<
6 240 ARLOTTO, MOSSEL, AND STEELE We then have β() =, and we have the recursion β(n + ) = G(β(n)) for all n. (3) Moreover, if the deterministic sequence t, t 2,...is defined by the recursion t = and t n+ = H(β(n)) for all n, (4) then the minimum in the defining equation () for β(n) is uniquely achieved by the sequence of stopping times given by the threshold recursion (3). Proof. An optimal first selection time has the form τ = min{i : X i < t}, so we can rewrite the recursion (6) as { } { ) β(n ) t } β(n) = min + E[β(n ] = min + 0<t< t X τ 0<t< t t 0 s ds = min g(t, β(n )) G(β(n )). (5) 0<t< The selection rule for the first element is given by τ = min{i : X i < t n } so by (5) and the definitions of g and H we have t n = H(β(n )). Lemma 3 already gives us enough to prove the first assertion of Theorem which states that the map n β(n) is convex. Lemma 4. The map n (n) := β(n + ) β(n) is an increasing function. Proof. One can give a variational characterization of that makes this evident. First, by the defining relations () and the recursion (3) we have β(n + ) β(n) = G(β(n)) β(n) = min 0<t< { t + β(n) [ t log ( ) ] }, t so if we set ĝ(x, t) = [ ( ) ] t + x t log, t then we have (n) = β(n + ) β(n) = min ĝ(β(n), t). (6) 0<t< Now, for 0 x y and t (0, ) we then have [ ( ) ] ĝ(x, t) ĝ(y, t) = (x y) t log = (x y) t so from the monotonicity β(n) β(n + ),weget ĝ(β(n), t) ĝ(β(n + ), t) for all 0 < t <. k=2 k tk 0; When we minimize over t (0, ), we see that (6) gives us (n) (n + ). Random Structures and Algorithms DOI 0.002/rsa
7 QUICKEST ONLINE SELECTION OF AN INCREASING SUBSEQUENCE 24 We next show that the two definitions in () can be used to give an a priori lower bound on G. An induction argument using the recursion (3) can then be used to obtain the lower half of (2). Lemma 5 (Lower Bounding G Recursion). For the function x G(x) defined by (), we have ( ) 2 (x + )2 G 2 x2 for all x. (7) Proof. To prove (7), we first note that by () it suffices to show that one has ( ) δ(x, t) = (x + ) 2 t 2 x 2 log 0 (8) t for all x and t (0, ).Forx the map t δ(x, t) is twicecontinuous differentiable and concave in t. Hence there is a unique value t (0, ) such that t = argmax 0<t< δ(x, t), and such that t = t (x) satisfies the first order condition Solving this equation gives us (x + ) 2 ( t ) x 2 = 0. t = 2x + ( (x + ), and δ(x, 2 t ) = + 2x 2x 2 log + ), x so the familiar bound x ( 2x log + ) 2 x for x, gives us ( δ(x, t) δ(x, t ) + 2x 2x 2 x ) = 0, 2x 2 and this is just what we needed to complete the proof of (8). Now, to argue by induction, we consider the hypothesis that one has 2 n2 β(n). (9) This holds for n = since β() =, and, if it holds for some n, then by the monotonicity of G we have G(n 2 /2) G(β(n)). Now, by (7) and (3) we have ( ) 2 (n + )2 G 2 n2 G(β(n)) = β(n + ), and this completes our induction step from (9). Random Structures and Algorithms DOI 0.002/rsa
8 242 ARLOTTO, MOSSEL, AND STEELE 3.. An Alternative Optimization Argument From the definition of the threshold policy (3) we also know that for 0 k < n the difference τ k+ τ k of the selection times is a geometrically distributed random variable with parameter λ k = t n k ( X τk ), so one has the representation Moreover, we also have n [ ] β(n) = E. (20) λ k k=0 E[X τk+ X τk X τk ]= 2 λ k so, when we take expectations and sum, we see that telescoping gives us n E[λ k ]=E[X τ 2 n], (2) k=0 where the last inequality is a trivial consequence of the fact that X τ n [0, ]. The relations (20) and (2) now suggest a natural optimization problem. Specifically, for n 2 we let {l k :0 k n } denote a sequence of random variables, and we consider the following convex program minimize subject to n k=0 E[l k ] n E[l k ] 2 k=0 E[l k ] 0, k = 0, 2,..., n. (22) This is a convex program in the n real variables E[l k ], k = 0,,..., n, and, by routine methods, one can confirm that its unique optimal solution is given by taking E[l k ]=2/n for all 0 k < n. The corresponding optimal value is obviously equal to n 2 /2. By the definition of the random variables {λ k : 0 k < n}, we have the bounds 0 λ k, and by (2) we see that their expectations satisfy the first constraint of the optimization problem (22), so {E[λ k ] :0 k < n} is a feasible solution of (22). Feasibility of the values {E[λ k ] :0 k < n} then implies that n 2 n2 k=0 ], (23) E [λ k and, by convexity of x /x on [0, ], we obtain from Jensen s inequality that n k=0 n ] E [λ k k=0 [ ] E = β(n). (24) λ k This subsection is based on the kind suggestions of an anonymous referee. Random Structures and Algorithms DOI 0.002/rsa
9 QUICKEST ONLINE SELECTION OF AN INCREASING SUBSEQUENCE 243 Taken together (23) and (24) tell us that n 2 /2 β(n), so we have a second proof of (9), or the lower bound (2) of Theorem. Here we should note that Gnedin [] used a similar convex programming argument to prove a tight upper bound in the corresponding sizefocused selection problem. Thus, we have here our first instance of the kind of duality that is discussed more fully in Section Completion of the Proof of Theorem We now complete the proof of Theorem by proving the upper half of (2). The argument again depends on an a priori bound on G. The proof is brief but delicate. Lemma 6 (Upper Bounding G Recursion). For the function x G(x) defined by () one has G( 2 x2 + x log(x)) 2 (x + )2 + (x + ) log(x + ) for all x. Proof. If we set f (x) := x 2 /2 + x log(x), then we need to show that G(f (x)) f (x + ). If we take t = 2/(x + 2) then the defining relation () for G tells us that G(f (x)) g(f (x), t ) = x x Next, for any y 0 integration over (0, y) of the inequality ( log + 2 ) f (x). (25) x u + u2 + 2u + 2 y(y + 2) implies the bound log( + y) 2(u + ) 2 2(y + ). If we now set y = 2/x and substitute this last bound in (25), we obtain G(f (x)) x ( ) f (x) x = f (x + ) + + (x + ){log(x) log(x + )} f (x + ), 2 just as needed to complete the proof of the lemma. One can now use Lemma 6 and induction to prove that for all n 2 one has β(n) 2 n2 + n log(n). (26) Since β(2) = G() = min 0<t< g(, t) <3.5 and 2( + log(2)) 3.39, one has (26) for n = 2. Now, for n 2, the monotonicity of G gives us that β(n) ( ) 2 n2 + n log(n) implies G(β(n)) G 2 n2 + n log(n). Random Structures and Algorithms DOI 0.002/rsa
10 244 ARLOTTO, MOSSEL, AND STEELE Finally, by the recursion (3) and Lemma 6 we have ( ) β(n + ) = G(β(n)) G 2 n2 + n log(n) 2 (n + )2 + (n + ) log(n + ). This completes the induction step and establishes (26) for all n 2. This also completes last part of the proof of Theorem. 4. ASYMPTOTICS FOR THE VARIANCE To complete the proof of Theorem 2, we only need to prove that one has the asymptotic formula (5) for Var[τn ]. This will first require an understanding of the size of the threshold t n, and we can get this from our bounds on β(n) once we have an asymptotic formula for H. The next lemma gives us what we need. Lemma 7. that For x G(x) and x H(x) defined by () and (2), we have for x G(x) = (x /2 + 2 /2 ) 2 { + O(/x)} and (27) 2 H(x) = x { + O(x /2 )}. (28) Proof. For any fixed x wehaveg(t, x) when t 0ort, so the minimum of g(t, x) is obtained at an interior point 0 < t <. Computing the tderivative g t (t, x) gives us g t (t, x) = t x 2 t log( 2 t ) + x t( t), so the firstorder condition g t (t, x) = 0 implies that at the minimum we have the equality t = x 2 t log( 2 t ) + x t( t). Writing this more informatively as = log( t) + t x t = t2 2 + i=3 i t i, (29) i we see the righthand side is monotone in t, so there is a unique value t = t (x) that solves (29) for t. The last sum on the righthand side of (29) tells us that 2 t2 or, equivalently, x and when we use these bounds in (29) we have t 2 x, t 2 2 x t2 2 + ( ) 2 i/2 t2 x 2 + O(x 3/2 ). i=3 Random Structures and Algorithms DOI 0.002/rsa
11 QUICKEST ONLINE SELECTION OF AN INCREASING SUBSEQUENCE 245 Solving these inequalities for t, we then have by the definition (2) of H that 2 H(x) = t = x { + O(x /2 )}. (30) Finally, to confirm the approximation (27), we substitute H(x) = t into the definition () of G and use the asymptotic formula (30) for H(x) to compute G(x) = g(x, t ) = + x ( ) ( ) log = t t i + x t t t i i= ( xt + + xt2 2 = { + xt + xt2 t 2 + O(xt 3 )}= t x = x O() = ( x /2 + 2 /2) 2 { + O(/x)}, ) + O() and this completes the proof of the lemma. The recursion (4) tells us that t n = H(β(n)) and the upper and lower bounds (2) of Theorem tell us that β(n) = n 2 /2 + O(n log n), so by the asymptotic formula (28) for H we have t n = 2 n + O(n 2 log n). (3) To make good use of this formula we only need two more tools. First, we need to note that the stopping time τn satisfies a natural distributional identity. This will lead in turn to a recursion from which we can extract the required asymptotic formula for v(n) = Var(τn ). If t n is the threshold value defined by the recursion (4), we let γ(t n ) denote a geometric random variable of parameter p = t n and we let U(t n ) denote a random variable with the uniform distribution on the interval [0, t n ]. Now, if we take the random variables γ(t n ), U(t n ), and τ n to be independent, then we have the distributional identity, τ d τn n = γ(t n ) + ( U(t n )), (32) and this leads to a useful recursion for the variance of τn. To set this up, we first put R(t) = ( U(t)) where U(t) is uniformly distributed on [0, t], and we note E[R(t)] = t log( t) = + t/2 + t 2 /3 + O(t 3 ); moreover, since E[R 2 (t)] =( t) = + t + t 2 + O(t 3 ), we also have Var[R(t)] =( t) t 2 log 2 ( t) = t2 2 + O(t3 ). (33) Lemma 8 (Approximate Variance Recursion). For the variance v(n) := Var(τn ) one has the approximate recursion { v(n) = + 2 } n + O(n 2 log n) v(n ) + 3 n2 + O(n log n). (34) Random Structures and Algorithms DOI 0.002/rsa
12 246 ARLOTTO, MOSSEL, AND STEELE Proof. By independence of the random variables on the right side of (32), we have v(n) = Var(γ (t n )) + Var[R(t n )τ n ]. (35) From (3) we have t n = 2/n + O(n 2 log n), so for the first summand we have Var(γ (t n )) = t 2 n t n = 4 n2 + O(n log n). (36) To estimate the second summand, we first use the complete variance formula and independence to get Var(R(t n )τ n ) = E[R2 (t n )]E[(τ n )2 ] E[R(t n )] 2 (E[τ n ])2 Now from (3) and (33) we have = E[R 2 (t n )]v(n ) + Var[R(t n )](E[τ n ])2. (37) Var[R(t n )]= 3 n 2 + O(n 3 log n), and from (2) we have { } 2 (E[τ n ])2 = 2 (n )2 + O(n log n) = 4 n4 + O(n 3 log n), so from (35), (36) and (37) we get v(n) ={ + 2 n + O(n 2 log n)}v(n ) + 3 n2 + O(n log n), and this completes the proof of (34). To conclude the proof of Theorem 2, it only remains to show that the approximate recursion (34) implies the asymptotic formula v(n) = 3 n2 (n + ) + O(n 2 log α n) for α>2. (38) If we define r(n) by setting v(n) = 3 n 2 (n + ) + r(n), then substitution of v(n) into (34) gives us a recursion for r(n), { r(n) = + 2 } n + O(n 2 log n) r(n ) + O(n log n). We then consider the normalized values ˆr(n) = r(n)/(n 2 log α (n)), and we note they satisfy the recursion ˆr(n) ={ + O(n 2 )}ˆr(n ) + O(n log α n). This is a recursion of the form ˆr(n + ) = ρ nˆr(n) + ɛ n, and one finds by induction its solution has the representation n ˆr(n) =ˆr(0)ρ 0 ρ ρ n + ɛ k ρ k+ ρ n. Here, the product of the evolution factors ρ n is convergent and the sum of the impulse terms ɛ n is finite, so the sequence ˆr(n) is bounded, and, consequently, (34) gives us (38). This completes the proof of the last part of Theorem 2. Random Structures and Algorithms DOI 0.002/rsa k=0
13 QUICKEST ONLINE SELECTION OF AN INCREASING SUBSEQUENCE SUBOPTIMAL POLICIES AND A BLOCKING INEQUALITY Several inequalities for β(n) can be obtained through the construction of suboptimal policies. The next lemma illustrates this method with an inequality that leads to an alternative proof of (9), the uniform lower bound for β(n). Lemma 9 (Blocking Inequality). For nonnegative integers n and m one has the inequality β(mn) min{m 2 β(n), n 2 β(m)}. (39) Proof. First, we fix n and we consider a policy π that achieves the minimal expectation β(n). The idea is to use π to build a suboptimal π policy for the selection of an increasing subsequence of length mn. We take X i, i =, 2,...to be a sequence of independent random variables with the uniform distribution on [0, ], and we partition [0, ] into the subintervals I =[0, /m), I 2 =[/m,2/m),..., I m =[(m )/m,]. We then define π by three rules: i. Beginning with i = wesayx i is feasible value if X i I.IfX i is feasible, we accept X i if X i would be accepted by the policy π applied to the sequence of feasible values after we rescale those values to be uniform in [0, ]. We continue this way until the time τ when we have selected n values. ii. Next, beginning with i = τ, we follow the previous rule except that now we say X i is feasible value if X i I 2. We continue in this way until time τ + τ 2 when n additional increasing values have been selected. iii. We repeat this process m 2 more times for the successive intervals I 3, I 4,..., I m. At time τ + τ 2 + +τ m, the policy π will have selected nm increasing values. For each j m we have E[τ j ]=mβ(n), so by suboptimality of π we have β(mn) E[τ + τ 2 + +τ m ]=m2 β(n). We can now interchange the roles of m and n, so the proof of the lemma is complete. The blocking inequality (39) implies that even the crude asymptotic relation β(n) = n 2 /2 + o(n 2 ) is strong enough to imply the uniform lower bound n 2 /2 β(n). Specifically, one simply notes from (39) and β(n) = n 2 /2 + o(n 2 ) that β(mn) (mn) 2 β(n) β(mn) and lim n 2 m (mn) = 2 2. This third derivation of the uniform bound n 2 /2 β(n) seems to have almost nothing in common with the proof by induction that was used in the proof of Lemma 6. Still, it does require the bootstrap bound β(n) = n 2 /2 + o(n 2 ), and this does require at least some of the machinery of Lemma DUALITY AND THE SIZEFOCUSED SELECTION PROBLEM In the online sizefocused selection problem one considers a set of policies s (n) that depend on the size n of a sample {X, X 2,..., X n }, and the goal is to make sequential Random Structures and Algorithms DOI 0.002/rsa
14 248 ARLOTTO, MOSSEL, AND STEELE selections in order to maximize the expected size of the selected increasing subsequence. More precisely, a policy π n s (n) is determined by stopping times τ i, i =, 2,...such that τ <τ 2 < <τ k n and X τ X τ2 X τk. The random variable of interest is L o n (π n) = max{k : X τ < X τ2 < < X τk where τ <τ 2 < <τ k n}, and most previous analyses have focused on the asymptotic behavior of l(n) := max πn s(n) E[Lo n (π n)]. For example, Samuels and Steele [5] found that l(n) 2n, but now a number of refinements of this are known. Our goal here is to show how some of these refinements follow from the preceding theory. 6.. Uniform Upper Bound for l(n) via Duality Perhaps the most elegant refinement of l(n) 2n is the following uniform upper bound that follows independently from related analyses of Bruss and Robertson [9] and Gnedin []. Bruss and Robertson [9] obtained the upper bound (40) by exploiting the equivalence of the sizefocused monotone subsequence problem with a special knapsack problem that we will briefly revisit in Section 7. On the other hand, Gnedin [] showed that (40) can be proved by an optimization argument like the one used to prove (22). Here, we give a third argument with the curious feature that one gets a sharp bound that holds for all n from an apparently looser relation that only holds approximately for large n. Proposition (Uniform Upper Bound). For all n, one has l(n) 2n. (40) As we noted above, this proposition is now well understood, but it still seems instructive to see how it can be derived from β(n) = (/2)n 2 + O(n log n). The basic idea is that one exploits duality with a suboptimality argument like the one used in Section 5, although, in this case a bit more work is required. We fix n, and, for a much larger integer k, weset N k = (k 2k 2/3 )l(n) and r k = k k 2/3. (4) The idea of the proof is to give an algorithm that is guaranteed to select from {X, X 2,...} an increasing subsequence of length N k.ift k is the number of the elements that the algorithm inspects before returning the increasing subsequence, then by the definition of β( ) we have β(n k ) E[T k ]. One then argues that (40) follows from this relation. We now consider [0, ] and for i r k, we consider the disjoint intervals I i = [(i )/k, i/k) and a final reserve interval I =[r k /k,] that is added to complete the partition of [0, ]. Next, we let ν() be the first integer such that S :={X, X 2,..., X ν() } I Random Structures and Algorithms DOI 0.002/rsa
15 QUICKEST ONLINE SELECTION OF AN INCREASING SUBSEQUENCE 249 has cardinality n, and for each i > we define ν(i) to be least integer greater ν(i ) for which the set S i :={X ν(i )+, X ν(i )+2,..., X ν(i) } I i has cardinality n. By Wald s lemma and (4) we have E[ν(r k )]=nkr k where r k = k k 2/3. (42) Now, for each i n, we run the optimal sizefocused sequential selection algorithm on S i, and we let L(n, i) be the length of the subsequence that one obtains. The random variables L(n, i), i r k, are independent, identically distributed, and each has mean equal to l(n). We then set L(n, r k ) = L(n,) + L(n,2) + +L(n, r k ), and we note that if L(n, r k ) N k, for N k as defined in (4), then we have extracted an increasing subsequence of length at least N k ; in this case, we halt the procedure. On the other hand if L(n, r k )<N k, we need to send in the reserves. Specifically, we recall that I =[r k /k,] and we consider the postν(r k ) reserve subsequence S :={X i : i >ν(r k ) and X i I }. We now rescale the elements of S to the unit interval, and we run the optimal timefocused algorithm on S until we get an increasing sequence of length N k.ifweletr(n, k) denote the number of observations from S that are examined in this case, then we have E[R(n, k)] =β(n k ) by the definition of β. Finally, since I has length at least k /3, the expected number of elements of {X i : i >ν(r k )} that need to be inspected before we have selected our increasing subsequence of length N k is bounded above by k /3 β(n k ). The second phase of our procedure may seem wasteful, but one rarely needs to use the reserve subsequence. In any event, our procedure does guarantee that we find an increasing subsequence of length N k in a finite amount of time T k. By (42) and the upper bound k /3 β(n k ) on the incremental cost when one needs to use the reserve subsequence, we have β(n k ) E[T k ] knr k +{knr k + k /3 β(n k )}P(L(n, r k )<N k ), (43) where, as noted earlier, the first inequality comes from the definition of β. The summands of L(n, r k ) are uniformly bounded by n and E[L(n, r k )]=r k l(n), soby the definition (4) of N k and r k we see from Hoeffding s inequality that P(L(n, r k )<N k ) P ( L(n, r k ) E[L(n, r k )] < (k 2/3 )l(n) ) (44) exp{ A n k /3 }, for constants A n, K n, and all k K n. The exponential bound (44) tells us that for each n there is a constant C n such that the last summand in (43) is bounded by C n for all k.bythe bounds (2) of Theorem we have β(n k ) = (/2)N 2 k + O(N k log N k ), and by (4) we have N k = (k 2k 2/3 )l(n) + O(), r k = k k /3 + O(), so, in the end, our estimate (43) tell us 2 l2 (n){k 2 2k 5/3 + 4k 4/3 } k 2 n + o n (k 2 ). When we divide by k 2 and let k,wefindl(n) 2n, just as we hoped. Random Structures and Algorithms DOI 0.002/rsa
16 250 ARLOTTO, MOSSEL, AND STEELE 6.2. Lower Bounds for l(n) and the Duality Gap One can use the timefocused tools to get a lower bound for l(n), but in this case the slippage, or duality gap, is substantial. To sketch the argument, we first let T r denote the time required by the optimal timefocused selection policy to select r values. We then follow the rtarget timefocused policy. Naturally, we stop if we have selected r values, but if we have not selected r values by time n, then we quit, no matter how many values we have selected. This suboptimal strategy gives us the bound rp(t r n) l(n), and from this bound and Chebyshev s inequality, we then have r{ Var(T r )/(n E[T r ]) 2 } l(n). (45) If we then use the estimates (2) and (5) for E[T r ] and Var[T r ] and optimize over r, then (45) gives us the lower bound (2n) /2 O(n /3 ). However in this case the timefocused bounds and the duality argument leave a big gap. Earlier, by different methods and for different reasons Rhee and Talagrand [3] and Gnedin [] obtained the lower bound (2n) /2 O(n /4 ) l(n). Subsequently, Bruss and Delbaen [7] studied a continuous time interpretation of the online increasing subsequence problem where the observations are presented to the decision maker at the arrival times of a unitrate Poisson process on the time interval [0, t), and, in this new formulation, they found the stunning lower bound 2t O(log t). Much later, Arlotto et al. [3] showed by a depoissonization argument that the lower bound of Bruss and Delbaen [7] can be used to obtain 2n O(log n) l(n) for all n under the traditional discretetime model for sequential selection. Duality estimates such as (45) seem unlikely to recapture this bound. 7. OBSERVATIONS, CONNECTIONS, AND PROBLEMS A notable challenge that remains is that of understanding the asymptotic distribution of τn, the time at which one completes the selection of n increasing values by following the unique optimal policy π that minimizes the expected time E[τn ]=β(n). By Theorems and 2 we know the behavior of the mean E[τn ] and of the variance Var[τ n ] for large n, but the asymptotic behavior of higher moments or characteristic functions does not seem amenable to the methods used here. Some second order asymptotic distribution theory should be available for τn, but it is not likely to be easy. Still, some encouragement can be drawn from the central limit theorem for the sizefocused selection problem that was obtained by Bruss and Delbaen [8] under the Poisson model introduced in Bruss and Delbaen [7]. Also, by quite different means, Arlotto et al. [3] obtained the central limit theorem for the sizefocused selection problem under the traditional model discussed in Section 6. Each of these developments depends on martingale arguments that do not readily suggest any natural analog in the case of timefocused selection. There are also several other sizefocused selection problems where it may be informative to investigate the corresponding timefocused problem. The most natural of these is probably the unimodal subsequence selection problems where one considers the random variable Random Structures and Algorithms DOI 0.002/rsa
17 QUICKEST ONLINE SELECTION OF AN INCREASING SUBSEQUENCE 25 U o n (π n) = max{k : X τ < X τ2 < < X τ t > X τt+ > > X τk, where τ <τ 2 < <τ k n}, and where, as usual, each τ k is a stopping time. Arlotto and Steele [4] found that max E[U o n (π n)] 2 n, πn asn, and an analogous timefocused selection problem is easy to pose. Unfortunately, it does not seem easy to frame a useful analog of Lemma 3, so the timefocused analysis of the unimodal selection problem r s open. For another noteworthy possibility one can consider the timefocused problem in a multidimensional setting. Here one would take the random variables X i, i =, 2,...,to have the uniform distribution on [0, ] d, and the decision maker s task is then to select as quickly as possible a subsequence that is monotone increasing in each coordinate. The dual, sizefocused, version of this problem was studied by Baryshnikov and Gnedin [6] who characterized the asymptotic behavior of the optimal mean. Finally, a further possibility is the analysis of a timefocused knapsack problem. Here one considers a sequence X, X 2,... of independent item sizes with common continuous distribution F, and the decision maker s task is to select as quickly as possible a set of n items for which the total size does not exceed a given knapsack capacity c. The analytical goal is then to estimate the size of the objective function [ φ F,c (n) = min E τ n : τ,...,τn n k= ] X τk c, where, just as before, the decision variables τ < τ 2 <..., < τ n are stopping times. The optimal mean φ F,c (n) now depends on the model distribution F and on the capacity parameter c. Never the less, if F represents the uniform distribution on [0, ] and if one takes c =, then the timefocused knapsack problem is equivalent to the timefocused increasing subsequence problem; specifically, one has φ Unif, (n) = β(n), for all n. The dual, sizefocused, version of this problem was first studied by Coffman et al. [0], Rhee and Talagrand [3], and Bruss and Robertson [9] who proved asymptotic estimates for the mean. While the structure of the optimal policy is now well understood (see, e.g. Papastavrou et al. [2]), only little is known about the limiting distribution of the optimal number of sizefocused knapsack selections. From Arlotto et al. [2] one has a suggestive upper bound on the variance, but at this point one does not know for sure that this bound gives the principle term of the variance for large n. ACKNOWLEDGMENTS The authors thank the referees for their useful suggestions, including a clarification of the proof of Lemma and the alternative optimization proof of (9) given in Section 3. Random Structures and Algorithms DOI 0.002/rsa
18 252 ARLOTTO, MOSSEL, AND STEELE REFERENCES [] D. Aldous and P. Diaconis, Longest increasing subsequences: from patience sorting to the BaikDeiftJohansson theorem, Bull Amer Math Soc (NS) 36 (999), [2] A. Arlotto, N. Gans, and J. M. Steele, Markov decision problems where means bound variances, Oper Res 62 (204), [3] A. Arlotto, V. V. Nguyen, and J. M. Steele, Optimal online selection of a monotone subsequence: A central limit theorem, Stoch Process Appl 25 (205), [4] A. Arlotto and J. M. Steele, Optimal sequential selection of a unimodal subsequence of a random sequence, Combin Probab Comput 20 (20), [5] J. Baik, P. Deift, and K. Johansson, On the distribution of the length of the longest increasing subsequence of random permutations, J Amer Math Soc 2 (999), [6] Y. M. Baryshnikov and A. V. Gnedin, Sequential selection of an increasing sequence from a multidimensional random sample, Ann Appl Probab 0 (2000), [7] F. T. Bruss and F. Delbaen, Optimal rules for the sequential selection of monotone subsequences of maximum expected length, Stoch Process Appl 96 (200), [8] F. T. Bruss and F. Delbaen, A central limit theorem for the optimal selection process for monotone subsequences of maximum expected length, Stoch Process Appl 4 (2004), [9] F. T. Bruss and J. B. Robertson, Wald s lemma for sums of order statistics of i.i.d. random variables, Adv Appl Probab 23 (99), [0] E. G. Coffman Jr., L. Flatto, and R. R. Weber, Optimal selection of stochastic intervals under a sum constraint, Adv Appl Probab 9 (987), [] A. V. Gnedin, Sequential selection of an increasing subsequence from a sample of random size, J Appl Probab 36 (999), [2] J. D. Papastavrou, S. Rajagopalan, and A. J. Kleywegt, The dynamic and stochastic knapsack problem with deadlines, Manage Sci 42 (996), [3] W. Rhee and M. Talagrand, A note on the selection of random variables under a sum constraint, J Appl Probab 28 (99), [4] D. Romik, The surprising mathematics of longest increasing subsequences, Cambridge University Press, Cambridge, 205. [5] S. M. Samuels and J. M. Steele, Optimal sequential selection of a monotone sequence from a random sample, Ann Probab 9 (98), [6] C. A. Tracy and H. Widom, Levelspacing distributions and the Airy kernel, Comm Math Phys 59 (994), [7] S. M. Ulam, Monte Carlo calculations in problems of mathematical physics, Modern mathematics for the engineer: Second series, McGrawHill, New York, 96, pp Random Structures and Algorithms DOI 0.002/rsa