2. Generalized backward induction and the Shapiro Catalan triangle; plan for the proof of Theorem 1.1

We start by patiently wading through some backward induction steps, rewarding us with a recognized pattern. We will want to decide whether to stop or continue when u is below and near the Brownian boundary. Since we will be using the Brownian motion value function heavily for everything that follows, we switch back to using (u,b) for position, as Shepp did in his wonderful paper [Reference Shepp13] which inspired our theoretical work on bounds. We remark again that in our development from now on, u is real, not just an integer, even though the Chow–Robbins game itself has only integer values for positions. This is important because our analysis depends on approximating V(u,b) which is continuous in u. The graphs in Section 1.4 suggest the advantage of extending to the real case to do the analysis.

The famous backward induction principle of optimal stopping (see, for example, [Reference Chow, Robbins and Siegmund3]) applied to this simple random walk is

(2.1) \begin{align}&V(u,b) = \max \left\{ {\frac{u}{b},\frac{1}{2}V(u + 1,b + 1) + \frac{1}{2}V(u - 1,b + 1)} \right\}.\\&\text{Do not stop if }\frac{1}{2}V(u + 1,b + 1) + \frac{1}{2}V(u - 1,b + 1) > \frac{u}{b}; \quad\text{stop otherwise.} \nonumber\end{align}

This is the starting point for everything. Let us use it to get a preliminary result. Numerical evidence showed that if $\delta = \alpha \sqrt b \,\, - \,\,u$ is larger than 1/2 minus a hair (the hair being of the order of ${b^{ - 1/4}}$ ), you should not stop. Our bounds on V are not alone good enough to prove that without any backward steps, but let us see what we get that way. Using only our lower bound from Lemma 1.2 (assuming $b > 1600$ ), and a differential approximation ${V_W}(u,b) \ge \left( {1 + {\delta ^2}{b^{ - 1}}} \right)u/b$ (from (3.6) in the next section), we get $V(u,b) \ge {V_W}(u,b)\left( {1 - 0.43{b^{ - 1}}} \right) \ge \left( {1 + {\delta ^2}{b^{ - 1}}} \right)\left( {1 - 0.43{b^{ - 1}}} \right)u/b$ , and this is greater than $u/b$ if $\delta > 0.66$ . Well, that is something: continue if $\delta > 0.66$ . But just one step of backward induction with our lower bound will show how it begins to close in on 1/2. The distance of $u - 1$ from the Brownian boundary is $\alpha \sqrt {b + 1} \,\, - \,\,(u - 1) = \alpha \sqrt b + h - u + 1 = \delta + 1 + h$ , where $0 < h < \alpha {b^{ - 1/2}}/2$ . So

\begin{align*}\frac{1}{2}V(u + 1,b + 1) + \frac{1}{2}V(u - 1,b + 1) &\geqslant \frac{1}{2}\frac{{u + 1}}{{b + 1}} + \frac{1}{2}{V_W}(u - 1,b + 1)\left( {1 - \frac{{0.43}}{{b + 1}}} \right)\\ &\geqslant \frac{1}{2}\frac{{u + 1}}{{b + 1}} + \frac{1}{2}\frac{{u - 1}}{{b + 1}}\left( {1 + \frac{{{{(1 + \delta )}^2}}}{{b + 1}}} \right)\left( {1 - \frac{{0.43}}{{b + 1}}} \right)\\ &= \frac{u}{b} - \frac{u}{{b(b + 1)}} \\&+ \frac{1}{2}\frac{{u\left( {1 - \frac{1}{{\alpha \sqrt b - \delta }}} \right)}}{{{{\left( {b + 1} \right)}^2}}}\!\left(\! {{{(1 + \delta )}^2} - 0.43 - \frac{{0.43{{(1 + \delta )}^2}}}{{b + 1}}} \right).\end{align*}

For sufficiently large b, so that we can throw out small stuff, the condition for this to be greater than $u/b$ is clearly ${(1 + \delta )^2} - 0.43 > 2$ , or $\delta > 0.56$ . But to be concrete, assume $b > 1600$ and $u = \alpha \sqrt b \,\, - \delta \,\, > \alpha \sqrt b \,\, - 0.66$ (we already know to continue if $\delta $ greater than 0.66). Then some arithmetic shows that $\delta > 0.58$ is sufficient. That is an improvement. And we could go more steps back and get a better go bound.

Similarly, we can use our upper bound and close in on 1/2 from above, and it is convenient to do one step of that, to get a preliminary stop bound, to avoid annoyances later. Stop if $(V(u + 1,b + 1) + V(u - 1,b + 1))/2 < u/b$ . And $\alpha \sqrt {b + 1} \,\, - \,\,(u - 1) = \delta + 1 + h$ where h is as before. Looking ahead to (3.5) of the next section, ${V_W}(u,b) \le u/b + \alpha {\delta ^2}{b^{ - 3/2}}$ . We have

\begin{align*}\frac{1}{2}V(u + 1,b + 1) + \frac{1}{2}V(u - 1,b + 1) &\leqslant \frac{1}{2}\frac{{u + 1}}{{b + 1}} + \frac{1}{2}\left( {\frac{{u - 1}}{{b + 1}} + \alpha \frac{{{{(1 + \delta + h)}^2}}}{{{{\left( {b + 1} \right)}^{3/2}}}}} \right)\\ &= \frac{u}{b} - \frac{u}{{b(b + 1)}} + \frac{1}{2}\alpha \frac{{{{(1 + \delta + h)}^2}}}{{{{\left( {b + 1} \right)}^{3/2}}}}\\ &= \frac{u}{b} - \frac{{\alpha \sqrt b - \delta }}{{b(b + 1)}} + \frac{1}{2}\alpha \frac{{{{(1 + \delta + h)}^2}}}{{{{\left( {b + 1} \right)}^{3/2}}}}.\end{align*}

For sufficiently large b, the condition for this to be less than $u/b$ is clearly ${(1 + \delta )^2} < 2$ , or $\delta < 0.414$ . But again to be concrete, assuming $b > 1600$ , $\delta < 0.38$ can be shown to be sufficient.

With just one step of backward induction, we have already narrowed the range for the stopping boundary and proved the following lemma.

This will be useful later, so it is noted. But the goal is to get to $\alpha \sqrt b - 1/2 + c{b^{ - 1/4}} + o({b^{ - 1/4}})$ , by continuing way down the backward induction tree.

We proceed to systematize going backward. Continued backward induction leads to the tree in Figure 4, where further branching is stopped at $u + 1$ , creating leaves of the tree at those nodes, which are shown boxed. The V value at a node of the tree is greater than or equal to the average of the V values at its two parents. Figure 4 is a picture of eight rows of the backward induction tree. The meaning of the coefficients (weights) will be explained shortly, though it is perhaps already obvious from the way backward induction works for this simple symmetric random walk.

Figure 4.

Backward induction tree, with leaves boxed.

To explain the coefficients (the weights) displayed in Figure 4, consider the following succession of inequalities, using only the basic backward induction inequality (2.1):

\begin{align*}&V(u,b) \geqslant \frac{1}{2}V(u + 1,b + 1) + \frac{1}{2}V(u - 1,b + 1),\\&\frac{1}{2}V(u - 1,b + 1) \geqslant \frac{1}{4}V(u,b + 2) + \frac{1}{4}V(u - 2,b + 2),\\&\frac{1}{4}V(u,b + 2) + \frac{1}{4}V(u - 2,b + 2)\\&\quad\geqslant \frac{1}{8}V(u + 1,b + 3) + \frac{1}{8}V(u - 1,b + 3) + \frac{1}{8}V(u - 1,b + 3) + \frac{1}{8}V(u - 3,b + 3)\\&\quad= \frac{1}{8}V(u + 1,b + 3) + \frac{2}{8}V(u - 1,b + 3) + \frac{1}{8}V(u - 3,b + 3),\\&\frac{2}{8}V(u - 1,b + 3) + \frac{1}{8}V(u - 3,b + 3)\\&\quad\geqslant \frac{2}{{16}}V(u,b + 4) + \frac{2}{{16}}V(u - 2,b + 4) + \frac{1}{{16}}V(u - 2,b + 4) + \frac{1}{6}V(u - 4,b + 4)\\&\quad= \frac{2}{{16}}V(u,b + 4) + \frac{3}{{16}}V(u - 2,b + 4) + \frac{1}{6}V(u - 4,b + 4),\\\ldots.\\\end{align*}

That is how the weights are generated, recursively. Considering any row of the tree, V(u,b) is greater than or equal to the sum of the weights times values at the leaves at or above that row, plus the weights times values at the non-leaf nodes at that row. The picture shows the weights through row 7, with the initial node at row 0. We were patient enough to carry out the trivial calculations by hand through seven rows, at which point we recognized the leaf weight numerators as being the Catalan numbers.

Now formalize the notation and the recursion. First, look at the nodes that are not leaves. With m indicating row number starting from 0 and j column number starting from 0 in the tree, let T(m,j) be the coefficient of ${2^{ - m}}V(u - j,b + m),m \geqslant 0,0 \leqslant j \leqslant m$ , and $T(m,j) = 0$ outside this range. The initial condition is $T(0,0) = 1$ . The recursion $T(m,j) = T(m - 1,j - 1) + T(m - 1,j + 1)$ follows purely from backward induction. Using this recursion produces the table in Figure 5, through $m = 7$ .

Figure 5.

T(m,j), coefficient numerator in row m, col. j of tree in Figure 4.

In what follows we will only make use of the odd rows and columns of the tree in Figure 4, which correspond to the rows of the tree with leaves. Let $B(n,k) = T(2n - 1,2k - 1),n \ge 1,k \ge 1$ ; Figure 6 shows the first four rows of B. The recursion for T implies (in two steps) the following recursion for B: $B(n,k) = B(n - 1,k - 1) + 2B(n - 1,k) + B(n - 1,k + 1),n \ge 2,k \ge 2$ , with initial conditions $B(1,1) = 1$ and $B(1,k) = 0,k > 1$ . This is the recursion and initial conditions for the Shapiro Catalan triangle [Reference Shapiro12], one of the many arrays referred to in the literature as a Catalan triangle, or unfortunately sometimes as ‘the’ Catalan triangle. The columns of B (omitting the leading zeros) appear as every other diagonal in another popular ‘Catalan’s triangle’, the one that appears, for example, in the Wikipedia article of that name, which gives references to the history. But we need to reference our tree by row and column indices, so B is the right one for us. Fortunately, Miana and Romero [Reference Miana and Romero11] have the row moments for B (see (4.2)), which is exactly what we need for our purposes.

Figure 6.

B(n,k), the odd rows and columns of T; it is the Shapiro Catalan triangle.

We now have a formula for the coefficients of the non-leaves in the odd rows of our backward induction tree in terms of well-known numbers:

(2.2) \begin{equation}{\textrm{coefficient of }}\,V(u - 2j + 1,b + 2m - 1) = {2^{ - 2m + 1}}B(m,j),\quad m \ge 1,j \ge 1.\end{equation}

There are simple known formulas for the entries in B:

\begin{align*}B(m,j) = \frac{j}{m}\left( {\begin{array}{*{20}{c}} {2m} \\ {m - j}\end{array}} \right) = \left( {\begin{array}{*{20}{c}} {2m - 1} \\ {m - j}\end{array}} \right) - \left( {\begin{array}{*{20}{c}} {2m - 1} \\ {m - j - 1}\end{array}} \right).\end{align*}

The first column is

\begin{align*}B(m,1) = \frac{1}{m}\left( {\begin{array}{*{20}{c}} {2m} \\ {m - 1}\end{array}} \right) = \frac{1}{{m + 1}}\left( {\begin{array}{*{20}{c}} {2m} \\ m\end{array}} \right) = \left( {\begin{array}{*{20}{c}} {2m - 1} \\ {m - 1}\end{array}} \right) - \left( {\begin{array}{*{20}{c}} {2m - 1} \\ {m - 2}\end{array}} \right) = {C_m}, \quad m \geqslant 1,\end{align*}

recognized as the mth Catalan number. ${C_0} = 1$ by definition.

For the leaves, let L(n) be the coefficient of ${2^{ - 2n + 1}}V(u + 1,b + 2n - 1),n \ge 1$ , where $n = 1,2,3,4,\ldots$ corresponds to rows $1,3,5,7,\ldots$ of the tree of Figure 4. The first four values of L(n) are 1,1,2,5, and from the way the tree is built, it is seen that $L(n) = T(2n - 2,0) = T(2n - 3,1),n \ge 2$ , with $L(1) = 1$ . But $T(2n - 3,1) = B(n - 1,1) = {C_{n - 1}}$ for $n \ge 2$ . So, we have the following formula for the leaf weights:

(2.3) \begin{equation}{\text{coefficient of }}\,V(u + 1,b + 2n - 1) = {2^{ - 2n + 1}}{C_{n - 1}},\quad n \geqslant 1.\end{equation}

Returning to the inequality we started with, for any row, V(u,b) is greater than or equal to the sum of the weights times values at the leaves at or above that row, plus the weights times values at the non-leaf nodes in that row. Looking at only odd rows of the tree, for $n \ge 1$ , let

\begin{align*}{S_L(n,u,b)} \,:\!= \sum\limits_{m = 0}^{n - 1} {{2^{ - 2m - 1}}{C_m}} V(u + 1,b + 2m + 1),\end{align*}

the sum of the leaves down through the row indexed by n; and

\begin{align*}{S_R}(n,u,b) \,:\!= {2^{ - 2n + 1}}\sum\limits_{j = 1}^n {B(n,j)V(u - 2j + 1,b + 2n - 1)},\end{align*}

the sum of the non-leaves across the row indexed by n (an odd row of the original tree). Define

(2.4) \begin{equation}TreeSum(n,u,b) \,:\!= S_L(n,u,b) + {S_R}(n,u,b).\end{equation}

Thus for any $n \ge 1$ , $V(u,b) \ge TreeSum(n,u,b)$ . This inequality is the summary of the chain of inequalities from repeated backward induction.

One property of TreeSum(n,u,b) is immediate from the basic backward induction principle (2.1): $TreeSum(n + 1,u,b) \le TreeSum(n,u,b),\,\,n \ge 1$ , which implies $TreeSum(n,u,b) \le TreeSum(1,u,b) = \left( {V(u + 1,b + 1) + V(u - 1,b + 1)} \right)/2$ . Thus for any n, $TreeSum(n,u,b) > u/b \Rightarrow TreeSum(1,u,b) > u/b \Rightarrow V(u,b) > u/b$ , so do not stop at (u,b). To prove the other direction, let $j \ge 0, k \ge 0$ . $V(u,b) > u/b \Rightarrow V(u-j,b) > (u-j)/b$ is obvious. Then $V(u-j,b) > (u-j)/b \Rightarrow u-j < {\beta}_b \Rightarrow u-b < {\beta}_{b+k}$ , from the monotonicity of Dvoretzky’s stop rules [Reference Dvoretzky6], so $V(u-j,b+k) > (u-j)/(b+k)$ . In other words, if you should not stop at (u,b), then you should not stop for a smaller u or a larger b. This implies that if $TreeSum(1,u,b) > u/b$ , then at every non-leaf node of the backward induction tree, the average of the values of its two children is greater than the ratio, so $TreeSum(n,u,b) = TreeSum(1,u,b)$ and $TreeSum(n,u,b) > u/b$ . This leads to the following lemma.

We can now explain the plan for the proof of Theorem 1.1, and give some indication of why it should work. The proof is accomplished in three stages.

Stage one, in Section 4, gets preliminary $O(b^{-1/4})$ bounds on the stop rule $\beta_b$ . For this stage, we consider ns of the form $c\sqrt b$ and u near enough to the boundary, that is, $\delta = \alpha \sqrt b \,\, - \,\,u$ small enough, so that $V(u + 1,b + 2m + 1) = (u + 1)/(b + 2m + 1)$ exactly for $m < n$ . With a binomial expansion of that ratio, the leaf sum can be approximated to any accuracy desired using simple formulas for $\sum_{m = 0}^{n - 1} {{2^{ - 2m - 1}}{m^k}{C_m}} $ , which we use for $k = 0,1,2$ . For the row sum, we use our simple upper and lower bounds on V in terms of ${V_W}$ , and approximate ${V_W}$ by a Taylor expansion about the boundary using four derivatives, and this leads to sums ${2^{ - 2n + 1}}\sum_{j = 1}^n {{j^k}B(n,j)} $ , which have simple known formulas. What is the good of letting n grow big? The point is that the larger n, the more of the weight of the tree sum is on the leaves, for which the value is exactly known as long as n does not get out of range, and the less is on the row sums, where we are limited in accuracy by our approximate bounds. This is our algebraic manifestation of Häggström and Wästlund’s idea to let the errors in the bounds wash out by moving the horizon back. In fact these sum formulas are all just simple expressions involving n and the central binomial ${2^{ - 2n}}\left( {\begin{array}{*{20}{c}} {2n} \\ n\end{array}} \right) \cong {(\pi n)^{ - 1/2}}$ , and n will be of order $\sqrt b $ , which hints at why ${b^{ - 1/4}}$ shows up in the answer. The first stage results in $O({b^{ - 1/4}})$ upper and lower bounds on the stop rule; this is Lemma 4.4, which in fact was our original goal. But it is not able to get the exact coefficient of ${b^{ - 1/4}}$ .

Stage two, in Section 5, feeds the result of stage one back into the Value approximations developed in stage one, to obtain a much sharper approximation for V near the boundary. It is perhaps the most conceptually tricky part of the proof, with a repeated feedback argument that will be better motivated when we get there. It shows that V is essentially piecewise linear near the boundary, below and tangent at integer points of $\delta $ , to the quadratic approximation to ${V_W}$ near the boundary. This is Theorem 5.1, perhaps of independent interest in showing the manner in which the Brownian Value overestimates V.

Finally, stage three, in Section 6, uses this improved estimate of V to go a bit further down the tree, by estimating the leaf values $V(u + 1,b + 2m + 1)$ for a range of m such that V is no longer just the ratio. By going just far enough down the tree, we are able to get the upper and lower bounds to come together, within an $o({b^{ - 1/4}})$ error, finding the exact coefficient of ${b^{ - 1/4}}$ and proving Theorem 1.1. The plan is straightforward except perhaps for Section 5, and uses standard approximations, but lots of them, so it looks worse than it is when all the approximations and sums are written out.

3. Approximating the Brownian motion value ${V_W}$

From formula (1.1), one may differentiate under the integral sign to see that all the derivatives with respect to u are positive for $u \le \alpha \sqrt b $ , so they all take their maximum value on the boundary; we use that several times below. To compute derivatives, and also for numerical work, it is best to write it in terms of standard functions. We have

\begin{align*}{V_W}(u,b) &= (1 - {\alpha ^2})\int_0^\infty {\exp \left( {\lambda u - \frac{{{\lambda ^2}b}}{2}} \right)d\lambda }\\&=(1 - {\alpha ^2}){b^{ - 1/2}}\exp \left( {\frac{{{u^2}}}{{2b}}} \right)\int_{ - \infty }^{u/\sqrt b } {\exp \left( { - \frac{{{w^2}}}{2}} \right)dw}\\ &= (1 - {\alpha ^2}){b^{ - 1/2}}G( {u/\sqrt b })/g( {u/\sqrt b })\\ &= (1 - {\alpha ^2}){b^{ - 1/2}}H( {u/\sqrt b }),\end{align*}

where G and g are the cumulative distribution function (CDF) and probability density function (PDF) of the standard normal distribution, respectively, and we have defined $H(x) \,:\!= G\left( x \right)/g\left( x \right)$ . This is related to the Mills ratio by $H(x)=1/g(x) - m(x),$ but we do not use that. On the boundary, ${u_0} = \alpha \sqrt b $ , ${V_W}({u_0},b) = {u_0}/b = \alpha {b^{ - 1/2}} = (1 - {\alpha ^2}){b^{ - 1/2}}H( {{u_0}/\sqrt b }) = (1 - {\alpha ^2}){b^{ - 1/2}}H\left( \alpha \right)$ , so we get the equation that alpha satisfies as $\alpha = (1 - {\alpha ^2})H\left( \alpha \right),$ which is useful for computing $\alpha $ using a library function for the Normal, or converted to an error function representation, to use that library function instead, if needed. We did that for the double-double precision numerical work for the extremely large number of positions we considered.

We will later use five derivatives of ${V_W}$ . Use ${D_u}$ for derivative operator with respect to the first argument. Then ${D_u}^{n}{V_W}(u,b) = (1 - {\alpha ^2}){b^{ - (n + 1)/2}}{H^{(n)}}( {u/\sqrt b } ),$ so we need to find the successive derivatives of H. It satisfies $H'(x) = xH(x) + 1$ , which makes this straightforward, and can easily be made systematic in terms of polynomials in x, similar to Hermite polynomials. One shows

(3.1) \begin{align}{H^{(n)}} = {P_n}H + {Q_n},\quad \text{with } {P_{n + 1}} = P'_n + x{P_n}, {Q_{n + 1}} = {Q_n}' + {P_n}, {P_0} = 1, {Q_0} = 0.\end{align}

So ${D_u}^{n}{V_W}(u,b) = (1 - {\alpha ^2}){b^{ - (n + 1)/2}}\left( {P_n}( {u/\sqrt b} )H( {u/\sqrt b } ) + {Q_n}( {u/\sqrt b } ) \right)\!.$

Remark. Just like for the Hermites, we have $P'_n = n{P_{n - 1}}$ (an Appell sequence), so just like for the Hermites, there is a computationally practical recurrence, not involving derivatives: ${P_{n + 1}}(x) = x{P_n}(x) + n{P_{n - 1}}(x)$ , the only difference from the Hermites being the positive sign. But we do not need to pursue this for our purposes, we just want a few derivatives.

On the boundary, $x = {u_0}/\sqrt b = \alpha \sqrt b /\sqrt b = \alpha $ and $H(\alpha ) = \alpha {(1 - {\alpha ^2})^{ - 1}},$ so

\begin{align*}{D_u}^{n}{V_W}({u_0},b) &= (1 - {\alpha ^2}){b^{ - n/2 + 1}}\left( {{P_n}\left( \alpha \right)H\left( \alpha \right) + {Q_n}\left( \alpha \right)} \right)\\&= {b^{ - (n + 1)/2}}\left( {\alpha {P_n}\left( \alpha \right) + (1 - {\alpha ^2}){Q_n}\left( \alpha \right)} \right).\end{align*}

To get the polynomials, turn the recursion crank in (3.1) five times, resulting in

(3.2) \begin{align} {P_0} &= 1,{Q_0} = 0;\quad {P_1} = x,{Q_1} = 1;\quad {P_2} = 1 + {x^2},{Q_2} = x;\quad {P_3} = 3x + {x^3},{Q_3} = 2 + {x^2};\nonumber \\ {P_4} &= 3 + 6{x^2} + {x^4},\,{Q_4} = 5x + {x^3};\quad {P_5} = 15x + 10{x^3} + {x^5},\,{Q_5} = 8 + 9{x^2} + {x^4}.\end{align}

To prove Lemma 1.2, we use three derivatives of H, which from (3.1) and (3.2) are

(3.3) \begin{align}{H^{(1)}}(x) &= xH(x) + 1,\nonumber\\{H^{(2)}}(x) &= (1 + {x^2})H(x) + x; \\{H^{(3)}}(x) &= (3x + {x^3})H(x) + 2 + {x^2}.\nonumber\end{align}

The first five derivatives of $V_W$ on the boundary are

(3.4) \begin{align}{D_u}{V_W}({u_0},b) &= {b^{ - 1}},\nonumber\\{D_u}^{2}{V_W}({u_0},b) &= {b^{ - 3/2}}2\alpha, \nonumber \\{D_u}^{3}{V_W}({u_0},b) &= {b^{ - 2}}\left( {2 + 2{\alpha ^2}} \right),\\{D_u}^{4}{V_W}({u_0},b) &= {b^{ - 5/2}}\left( {8\alpha + 2{\alpha ^3}} \right), \nonumber\\{D_u}^{5}{V_W}({u_0},b) &= {b^{ - 3}}\left( {8 + 16{\alpha ^2} + 2{\alpha ^4}} \right). \nonumber\end{align}

Let $\delta = {u_0} - u = \alpha \sqrt b - u$ . Using Taylor with just the first two derivatives,

\begin{align*}{V_W}(u,b) &= {V_W}({u_0},b) + (u - {u_0}){D_u}{V_W}({u_0},b) + \,{(u - {u_0})^2}{D_u}^{2}{V_W}({u_*},b)/2\\ &\le \alpha {b^{ - 1/2}} - \delta {b^{ - 1}} + {\delta ^2}{D_u}^{2}{V_W}({u_0},b)/2 = u{b^{ - 1}} + \alpha {\delta ^2}{b^{ - 3/2}}.\end{align*}

The inequality is because ${D_u}^{2}{V_W}(u,b)$ is an increasing function (recall all derivatives are positive), so is largest on the boundary. In summary,

(3.5) \begin{equation}{V_W}(u,b) \leqslant u{b^{ - 1}} + \alpha {\delta ^2}{b^{ - 3/2}}.\end{equation}

Using the third derivative with a similar monotonicity argument gives a lower bound, the one that we already used above in proving Lemma 2.1:

\begin{align*}{V_W}(u,b) &\ge u{b^{ - 1}} + \alpha {\delta ^2}{b^{ - 3/2}} - (1 + {\alpha ^2}){\delta ^3}{b^{ - 2}}/3 = u{b^{ - 1}} + \alpha \sqrt{b}{\delta ^2}{b^{ - 2}} - (1 + {\alpha ^2}){\delta ^3}{b^{ - 2}}/3\\ &= u{b^{ - 1}} + (u + \delta ){\delta ^2}{b^{ - 2}} - (1 + {\alpha ^2}){\delta ^3}{b^{ - 2}}/3 \quad\text{(recall }u+ \delta = \alpha \sqrt{b} \text{)} \\ &= u{b^{ - 1}}\left( {1 + {\delta ^2}{b^{ - 1}}} \right) + {\delta ^3}{b^{ - 2}}\left( {1 - (1 + {\alpha ^2})/3} \right) \ge u{b^{ - 1}}\left( {1 + {\delta ^2}{b^{ - 1}}} \right).\end{align*}

Thus

(3.6) \begin{equation}{V_W}(u,b) \geqslant u{b^{ - 1}}\left( {1 + {\delta ^2}{b^{ - 1}}} \right).\end{equation}

Using four derivatives, the monotonicity argument yields the upper bound

(3.7) \begin{equation}{V_W}(u,b) \leqslant u{b^{ - 1}} + \alpha {\delta ^2}{b^{ - 3/2}} - (1 + {\alpha ^2}){\delta ^3}{b^{ - 2}}/3 + (4\alpha + {\alpha ^3}){\delta ^4}{b^{ - 5/2}}/12.\end{equation}

With the fifth derivative term, one gets the lower bound

(3.8) \begin{eqnarray}{V_W}(u,b) \geqslant u{b^{ - 1}} + \alpha {\delta ^2}{b^{ - 3/2}} - (1 + {\alpha ^2}){\delta ^3}{b^{ - 2}}/3 + (4\alpha + {\alpha ^3}){\delta ^4}{b^{ - 5/2}}/12\\ - (4 + 8{\alpha ^2} + {\alpha ^4}){\delta ^5}{b^{ - 3}}/60.\nonumber\end{eqnarray}

These last two inequalities are used in proving Lemma 4.2.

4. Computing bounds on tree sums

This section is based on a lot of computation. We have resorted to big-O notation rather than getting specific constants, for which we apologized earlier. We will use big-O notation in inequalities, since we deal with approximate upper and lower bounds. The meaning will probably be clear from the context, but to be certain, we define our notation. Let f,g be real-valued functions of real variable b, with $g(b) > 0$ for all sufficiently large b. Define $f(b) \le O(g(b))$ to mean there exist $b_0$ and $K > 0$ such that $f(b) \le Kg(b), b \ge b_0$ . Usually only $f(b) = O(g(b))$ is defined in texts, requiring $\left| f(b) \right| \le Kg(b), b \ge b_0$ . That would not be convenient for our purposes. We will write expressions such as $f(b) \le U(b) + O\left( b^{ - 2} \right)$ , to mean there exist $b_0$ and $K > 0$ such that $f(b) \le U(b) + K{b^{ - 2}}$ for all $b \ge b_0$ . That is the same as $f(b) - U(b) \le O\left( {b^{ - 2}} \right)$ in our definition. U is an approximate upper bound for f, but f(b) could be arbitrarily far below U(b). This seems very natural for dealing with inequalities rather than equalities, but does not seem to be a standard notation. Similarly, still assuming $g(b) > 0$ for all sufficiently large b, define $f(b) \ge O\left( g(b) \right)$ to mean there exist $b_0$ and $K>0$ such that $f(b) \ge -Kg(b), b \ge b_0$ . For example, $f(b) \ge U(b) + O\left({b^{-2}} \right)$ means there exist $b_0$ and $K>0$ such that $f(b) \ge U(b) -K{b^{-2}}$ for all $b\ge b_0$ .

Here is a way to describe the situation using set language:

\begin{align*}O(g(b)) = \left\{ h\,:\text{ there exist }K > 0 \text{ and } b_0 \text{ such that }\left| h(b) \right| \le Kg(b)\text{ for all }b\ge b_0 \right\}.\end{align*}

Note $h\in O(g(b))$ implies $|h|$ , $-|h|$ , and $-h$ are in O(g(b)), and $O(g(b))=-O(g(b))$ . Now $f(b) \le O( g(b))$ means there exists $h \in O( g(b))$ such that $f(b) \le h(b) \text{ for all } b$ ; $f(b) = O( g(b))$ means there exists $h \in O(g(b))$ such that $f(b) = h(b) \text{ for all }b$ ; and $f(b) \ge O( g(b))$ means there exists $h\in O( g(b))$ such that $f(b) \ge h(b)\text{ for all } b$ . That unifies the settings: for some $h \in O(g(b))$ , the inequality or equality is true when h(b) replaces O(g(b)).

Throughout this section, for integer $b>1600$ and real number u, let $\delta = \alpha \sqrt b - u$ . If $\delta < 0.38 $ or $\delta > 0.58$ , we already know whether to stop or keep going, by Lemma 2.1. But we will also want to estimate tree sums for a somewhat larger delta, to estimate leaf sums in proving Theorem 5.1. For this section, we assert that $\delta \le b^p$ , where $0 \le p \le 1/10$ . We also assert that $n = O(b^{1/2 + p})$ , where, as above, n indexes how far we go down the tree, looking at only the odd rows (the ones with leaves). We prove Lemma 4.3 with this generality. In our first application of Lemma 4.3, to get preliminary stop bounds, p will be 0. But in Section 5 we will use Lemma 4.3 with larger delta to get an improved estimate of Value, and then in Section 6 with larger n to finally prove Theorem 1.1.

Let ${G_n} \,:\!= {2^{ - 2n}}\left( {\begin{array}{*{20}{c}}{2n}\\n\end{array}} \right)$ , for which it is well known as a central binomial coefficient that ${G_n} \cong {\left( {\pi n} \right)^{ - 1/2}}$ ; to be more precise, ${\left( {\pi (n + 1/2)} \right)^{ - 1/2}} \le {G_n} \le {\left( {\pi n} \right)^{ - 1/2}}$ . The formulas that will be needed for the moments of Catalan numbers and Shapiro Catalan triangle rows can be simply expressed in terms of ${G_n}$ .

Lemma 4.1. (Catalan and Shapiro Catalan triangle moments.) We have

(4.1) \begin{align}&\sum\limits_{j = 0}^{n - 1} {{C_j}} {2^{ - 2j - 1}} = 1 - {G_n},\quad\sum\limits_{j = 0}^{n - 1} {{C_j}} {2^{ - 2j - 1}}j = n{G_n} + {G_n} - 1, \\&\sum\limits_{j = 0}^{n - 1} {{C_j}} {2^{ - 2j - 1}}{j^2} = \frac{1}{3}{n^2}{G_n} - \frac{4}{3}n{G_n} - {G_n} + 1; \nonumber\end{align}

and

(4.2) \begin{align}&{2^{ - 2n + 1}}\sum\limits_{j = 1}^n {B(n,j)} = {G_n},\quad{2^{ - 2n + 1}}\sum\limits_{j = 1}^n {jB(n,j)} = \frac{1}{2},\\&\,\,{2^{ - 2n + 1}}\sum\limits_{j = 1}^n {{j^2}B(n,j)} = n{G_n},\quad {2^{ - 2n + 1}}\sum\limits_{j = 1}^n {{j^3}B(n,j)} = \frac{{3n - 1}}{4}, \nonumber\\&{2^{ - 2n + 1}}\sum\limits_{j = 1}^n {{j^4}B(n,j)} = n(2n - 1){G_n},\quad{2^{ - 2n + 1}}\sum\limits_{j = 1}^n {{j^5}B(n,j)} = \frac{{15n(n - 1) + 2}}{8}. \nonumber\end{align}

When proving Lemma 4.4, we will be using the ratio for the value, on all the leaves, to compute the leaf sum. For getting an upper bound we will be restricting n to be small enough so that the ratio is in fact the correct value on all leaves, and for the lower bound we will be restricting n to be small enough so that the ratio (which is always a lower bound) will not be far off from the correct value. But in finally proving Theorem 1.1 in Section 6, n will be large enough so that the values at some of the leaves will be significantly larger than the ratio, and we have to take that into account. It will be convenient to split the leaf sum into two parts. For general u and b, define ${V_E}(u,b) = V(u,b) - u/b$ , the amount the value exceeds the ratio. Let

\begin{align*}{S_{LR}}(n,u,b) \,:\!= \sum\limits_{m = 0}^{n - 1} {{2^{ - 2m - 1}}{C_m}} \frac{{u + 1}}{{b + 2m + 1}},\end{align*}

the contribution from the ratio, and

\begin{align*}{S_{LE}}(n,u,b) \,:\!= \sum\limits_{m = 0}^{n - 1} {{2^{ - 2m - 1}}{C_m}} {V_E}(u + 1,b + 2m + 1),\end{align*}

from the excess over the ratio. So

\begin{align*}{S_L} = {S_{LR}} + {S_{LE}}.\end{align*}

We have

\begin{align*}{S_{LR}} &= \sum\limits_{m = 0}^{n - 1} {{2^{ - 2m - 1}}{C_m}} \frac{{u + 1}}{{b + 2m + 1}} = \frac{{u + 1}}{b}\sum\limits_{m = 0}^{n - 1} {{2^{ - 2m - 1}}{C_m}} {\left( {1 + \frac{{2m + 1}}{b}} \right)^{ - 1}}\\ &= \frac{{u + 1}}{b}\sum\limits_{m = 0}^{n - 1} {{2^{ - 2m - 1}}{C_m}} \left( {1 - \frac{{2m + 1}}{b} + \gamma (m)\frac{{{{\left( {2m + 1} \right)}^2}}}{{{b^2}}}} \right),\end{align*}

where $1 - \frac{{2m + 1}}{b} \le \gamma (m) \le 1$ . From (4.1),

\begin{align*}\sum\limits_{m = 0}^{n - 1} {{2^{ - 2m - 1}}{C_m}} {\left( {2m + 1} \right)^2} &= \frac{4}{3}{n^2}{G_n} - \frac{4}{3}n{G_n} - {G_n} + 1, \\\sum\limits_{m = 0}^{n - 1} {{2^{ - 2m - 1}}{C_m}} {\left( {2m + 1} \right)} &= 2n{G_n} - 1 + {G_n},\\\sum\limits_{m = 0}^{n - 1} {{2^{ - 2m - 1}}{C_m}} &= 1 - {G_n}.\end{align*}

Since $\gamma (m) \le 1,$

\begin{align*}{S_{LR}} &\le \dfrac{{u + 1}}{b}(1 - {G_n}) - \dfrac{{u + 1}}{{{b^2}}}\left( {2n{G_n} - 1 + {G_n} - \dfrac{4}{3}\dfrac{{{n^2}}}{b}{G_n} + \dfrac{4}{3}\dfrac{n}{b}{G_n} - \dfrac{1}{b} + \dfrac{{{G_n}}}{b}} \right)\\&= \dfrac{{u + 1}}{b}(1 - {G_n}) - \dfrac{{u + 1}}{{{b^2}}}\left( {2n{G_n} - 1 + {G_n} - \dfrac{4}{3}\dfrac{{{n^2}}}{b}{G_n} + O({b^{ - 1/2}})} \right),\end{align*}

where we used $ n = O({b^{1/2 + p}}), \frac{n}{b}{G_n} = O({b^{ - 3/4 + p/2}}) = O({b^{ - 1/2}})$ for our range of p.

For a lower bound, use $\gamma (m) \ge 1 -2{b^{ - 1/2 + p}}.$ Note ${b^{ - 1/2 + p}}({n^2}{G_n}/{b}) = O({b^{ - 3/4 + 5p/2}}) = O({b^{ - 1/2}})$ , from which it follows that the above upper bound for $S_{LR}$ is also a lower bound, to that order.

Anticipating combining $S_{LR}$ with the row sum $S_R$ and then getting the difference between the total and the ratio in terms of $\delta$ , write

\begin{align*}{S_{LR}} = \dfrac{u}{b} + \dfrac{1- (u + 1)G_n}{b} - \dfrac{\alpha}{b^{3/2}}\left(1+\dfrac{1 - \delta}{\alpha \sqrt{b}}\right)\left( {2n{G_n} - 1 + {G_n} - \dfrac{4}{3}\dfrac{{{n^2}}}{b}{G_n} + O({b^{ - 1/2}})} \right).\end{align*}

Using $\delta = O(b^p), n=O(p^{1/2+p})$ and $p \le 1/10$ , one can show $\dfrac{\delta{n^2}{G_n}}{b^{3/2}} = O(b^{-1/2}).$ It follows that

(4.3) \begin{align}{S_{LR}} = &\frac{u}{b} + \frac{1- (u + 1)G_n}{b}\nonumber \\& - \frac{\alpha}{b^{3/2}}\left( {2n{G_n} - 1 + {G_n} - \frac{4}{3}\frac{{{n^2}}}{b}{G_n}}+ \frac{2n{G_n}}{\alpha \sqrt{b}} - \delta \frac{2n{G_n}}{\alpha \sqrt{b}} + O({b^{ -1/2 + p}})\right).\end{align}

That takes care of the leaves; the row sum will be added to this. To get the row sum, we first replace V with $V_W$ and define

\begin{align*}{S_{RW}}(n,u,b) \,:\!= {2^{ - 2n + 1}}\sum\limits_{j = 1}^n {B(n,j)V_W(u - 2j + 1,b + 2n - 1)},\end{align*}

and use the Taylor approximation from Section 3 to compute it. This will be an upper bound for $S_R$ .

Lemma 4.2. (Row sum with $V_W$ ) We have

\begin{align*}{S_{RW}}(n,u,b) = \dfrac{(u + 1)G_n -1} {b} + \dfrac{\alpha }{{b^{3/2}}}X,\end{align*}

where

\begin{align*}X\,:\!=2n{G_n}-2+2{G_n}+ 2(\alpha + 1/\alpha) \dfrac{{n{G_n}}}{{\sqrt b }} - \dfrac{(4+\alpha ^2)}{3}\dfrac{{{n^2}{G_n}}}{b}\\ + 2\delta \left( {1 - {G_n} -( \alpha + 1/\alpha) \dfrac{{n{G_n}}}{{\sqrt b }}} \right) + {\delta ^2}{G_n} + O({b^{ - 1/2 + 2p}}).\end{align*}

The proof is along the lines of the above calculation of $S_{LR}$ except using (4.2) instead of (4.1). But it is very much longer, so we defer the proof to Appendix B.

To get a lower bound for $S_R$ , using Lemma 1.2 we can multiply ${S_{RW}}$ by

\begin{align*}1 - \dfrac{{5/12}}{{b + 2n - 1}}\left( {1 + \dfrac{1}{{\sqrt {b + 2n - 1} }}} \right) = 1 - \dfrac{5}{{12b}} + O\left( {{b^{ - 3/2 + p}}} \right).\end{align*}

Now

\begin{align*}\dfrac{5}{{12b}}S_{RW} =\dfrac{5}{{12b}}\dfrac{uG_{n}+1}{b} + O(b^{-2})=\dfrac{\alpha }{{b^{3/2}}}\dfrac{5G_{n}}{12} + O(b^{-2+p}),\end{align*}

so

\begin{align*}S_{RW} - \dfrac{\alpha }{{b^{3/2}}}\dfrac{5G_{n}}{12} + O(b^{-2+p}) \le S_R.\end{align*}

We get some simplification when we add

\begin{align*}S_{LR} + S_{RW} = \dfrac{u}{b} + \dfrac{\alpha }{{b^{3/2}}}Y,\end{align*}

where

\begin{align*}Y\,:\!=-1+{G_n}+ 2\alpha \dfrac{{n{G_n}}}{{\sqrt b }} - \dfrac{{{\alpha ^2}}}{3}\dfrac{{{n^2}{G_n}}}{b} + 2\delta \left( {1 - {G_n} - \alpha \dfrac{{n{G_n}}}{{\sqrt b }}} \right) + {\delta ^2}{G_n} + O({b^{ - 1/2 + 2p}}).\end{align*}

Thus

\begin{align*}\dfrac{\alpha }{{b^{3/2}}}Y' \le S_{LR} + S_{R}- \dfrac{u}{b} \le \dfrac{\alpha }{{b^{3/2}}}Y,\end{align*}

where $Y'=Y-\dfrac{5}{12}G_n.$ We have proved the following result.

Lemma 4.3. (Bounds on tree sum.) Let $0 \le p \le 1/10$ , $n = O({b^{1/2 + p}})$ and $\delta \le {b^p}$ , where $\delta = \alpha \sqrt b - u$ . Then

\begin{align*}\frac{\alpha }{b^{3/2}}\left( C' + \delta B + {\delta ^2}A \right) \le TreeSum(n,u,b) - S_{LE} - \frac{u}{b} \le \frac{\alpha }{b^{3/2}}\left( C + \delta B + {\delta ^2}A \right),\end{align*}

where

\begin{align*}&C \,:\!= - (1 - {C_1}{G_n})+ O(b^{ - 1/2 + 2p}),\quad C' \,:\!= - (1 - C_1'{G_n})+ O(b^{ - 1/2 + 2p}),\quad A\,:\!=G_n,\\&B\,:\!=2(1- {B_1}{G_n}),\quad\text{with }C_1\,:\!=1+2\frac{{\alpha}n}{\sqrt b}-\frac{1}{3}\frac{\alpha^2 n^2}{b},C_1'\,:\!=C_1- 5/12, B_1\,:\!=1+\frac{\alpha n}{\sqrt b}.\end{align*}

We can use this immediately to get preliminary $O({b^{ - 1/4}})$ bounds on the stop rule. If we do not go too far down the tree, ${S_{LE}}$ will be zero. By Lemma 2.1, the value $V(u + 1,b + 2n - 1)$ at the leaf will be just the ratio if $\alpha \sqrt {b + 2n - 1} - (u + 1) < 0.38$ . But $\alpha \sqrt {b + 2n - 1} - (u + 1) = \alpha \sqrt b {\left( {1 + \frac{{2n - 1}}{b}} \right)^{1/2}} - (\alpha \sqrt b - \delta + 1) < \alpha n/\sqrt{b} + \delta - 1 \le 0.38$ if $\alpha n / {\sqrt b} \le 0.8$ and $\delta \le 0.58$ , and then ${S_{LE}} = 0$ . We make note of this for use later:

(4.4) \begin{equation}\alpha n/\sqrt b \le 0.8 \text{ and } \delta \le 0.58 \text{ implies } S_{LE} = 0.\end{equation}

If ${S_{LE}} = 0$ and $C + \delta B + {{\delta ^2}}A \le 0$ , the tree sum will not exceed the ratio, and if $C' + \delta B + {{\delta ^2}}A \ge 0$ , the ratio will not exceed the tree sum. Let ${\delta _0}$ and $\delta _0'$ be the positive numbers satisfying $C + {\delta _0}B + {\delta _0^2}A = 0$ and $C' + \delta _0'B + {\delta _0'^2}A = 0$ , respectively. From Lemma 2.2, it follows that $\alpha \sqrt b \,\, - \delta _0' \le {\beta _b} \le \alpha \sqrt b \,\, - {\delta _0}\,$ . In (a) of the next lemma, we will estimate ${\delta _0}$ and $\delta _0'$ using some convenient ns to get preliminary stop bounds from this. In addition, in (b) below we will estimate ${\delta _0}$ and $\delta _0'$ for a value of n that will be used throughout Section 5 for convenience.

Lemma 4.4. (Preliminary bounds on stop rule)

\begin{align*}\text{(a) } \alpha \sqrt b -1/2+\frac{0.231}{\sqrt \pi}b^{-1/4}+O(b^{-1/2}) \le \beta _b \le \alpha \sqrt b -1/2+\frac{0.429}{\sqrt \pi}b^{-1/4}+O(b^{-1/2}).&&\end{align*}

(b) Let $n=\lfloor{ \sqrt b /2}\rfloor$ , and let $\delta _0$ and $\delta _0'$ be the positive numbers satisfying $C + {\delta _0}B + {\delta _0^2}A = 0$ and $C' + \delta _0'B + {\delta _0'^2}A = 0$ , respectively, where A, B, C, and C’ are as in Lemma 5.2. Then there exists $b_0$ such that

\begin{align*}1/2 - \frac{0.44}{\sqrt \pi}b^{-1/4} \le \delta_0 \le {\delta_0}' \le 1/2 - \frac{0.13}{\sqrt \pi}b^{-1/4},\quad b \ge b_0.\end{align*}

Proof. We will have $n=\Theta(b^{-1/2})$ in all cases here, so $A=G_n=O(n^{-1/2})=O(b^{-1/4})$ . For $\delta _0$ , approximate the solution to the quadratic equation

\begin{align*}\delta _0 &= \dfrac{-B+\sqrt {B^2 - 4AC}}{2A}\\&=\dfrac{-B+B\left(1-2ACB^{-2}-2{A^2}{C^2}{B^{-4}}+O(A^3)\right)}{2A}\\&=-\dfrac{C}{B}-\dfrac{AC^2}{{B^3}}+O(b^{-1/2})\end{align*}

since $A=O(b^{-1/4}).$ Now

\begin{align*}-\dfrac{C}{B}-\dfrac{AC^2}{{B^3}}&=\dfrac{1-C_1G_n+O(b^{-1/2+2p})}{2(1-B_1G_n)}-G_n\dfrac{\left(1-C_1G_n+O(b^{-1/2+2p})\right)^2}{8(1-B_1G_n)^3}\\&=\dfrac{1}{2}(1-C_1G_n)(1+B_1G_n)-\dfrac{G_n}{8}+O(b^{-1/2+2p})\\&=\dfrac{1}{2}-\dfrac{1}{2}\left(C_1-B_1+\dfrac{1}{4}\right)G_n+O(b^{-1/2+2p}).\end{align*}

This gives

(4.5) \begin{equation}\delta_0=\frac{1}{2}-\frac{1}{2}\left(\frac{\alpha n}{\sqrt b} - \frac{1}{3} \frac{{\alpha^2}{n^2}}{b}+\frac{1}{4}\right)G_n +O(b^{-1/2+2p}).\end{equation}

Changing C to C’ in the above calculation gives

(4.6) \begin{equation}\delta_0'=\frac{1}{2}-\frac{1}{2}\left(\frac{\alpha n}{\sqrt b} - \frac{1}{3} \frac{{\alpha^2}{n^2}}{b}-\frac{1}{6}\right)G_n +O(b^{-1/2+2p}).\end{equation}

To get an upper bound on the stop rule, we may assume $\delta \le 0.58$ by Lemma 2.1. Take $n = \lfloor {0.8 \sqrt b / \alpha } \rfloor $ , so ${\alpha n / \sqrt b} \le 0.8$ and ${S_{LE}} = 0$ , and ${\alpha n / \sqrt b } = 0.8 + O({b^{ - 1/2}})$ , and $p = 0$ . Then

\begin{align*}A = {G_n} = \dfrac{1}{{\sqrt {\pi n} }}\left( {1 + O({n^{ - 1}})} \right) = \sqrt {\dfrac{\alpha }{{0.8}}} \dfrac{{{b^{ - 1/4}}}}{{\sqrt \pi }}\left( {1 + O({b^{ - 1/2}})} \right).\end{align*}

By (4.5),

\begin{align*}\delta _0 &= \dfrac{1}{2} - \dfrac{1}{2}\left( {0.8 - \dfrac{1}{3}0.64 + \dfrac{1}{4}} \right){G_n} + O({b^{ - 1/2}})\\ &= \dfrac{1}{2} - .4184\sqrt {\dfrac{\alpha }{{0.8}}} \dfrac{{{b^{ - 1/4}}}}{{\sqrt \pi }} + O({b^{ - 1/2}})\\ &\ge \dfrac{1}{2} - \dfrac{{0.429}}{{\sqrt \pi }}{b^{ - 1/4}} + O({b^{ - 1/2}}),\end{align*}

which proves the upper bound.

For the lower bound, we decide to take $n = \lfloor 1.1{\sqrt b / \alpha } \rfloor $ so $\alpha n / \sqrt b = 1.1 + O({b^{ - 1/2}})$ , to get a slightly better result. For this n, it is not quite true that ${S_{LE}} = 0$ , but the expression is still a lower bound. With that choice, (4.6) gives

\begin{align*}\delta _0' &= \dfrac{1}{2} - \dfrac{1}{2}\left( {1.1 - \dfrac{1}{3}1.21 - \dfrac{1}{6}} \right)\sqrt {\dfrac{\alpha }{{1.1}}} \dfrac{{{b^{ - 1/4}}}}{{\sqrt \pi }} + O({b^{ - 1/2}})\\ &\le \dfrac{1}{2} - \dfrac{{0.231}}{{\sqrt \pi }}{b^{ - 1/4}}+ O({b^{ - 1/2}}),\end{align*}

completing the proof of (a).

For (b), we take $n=\lfloor \sqrt b / 2 \rfloor$ , to be used in the next section with larger $\delta$ , where p is not zero. Then (4.5) and (4.6) with this n yield

\begin{align*}1/2 - \dfrac{{0.433}}{{\sqrt \pi }}{b^{ - 1/4}} + O({b^{ - 1/2 + 2p}}) \le {\delta _0} \le \delta _0' \le 1/2 - \dfrac{{0.137}}{{\sqrt \pi }}{b^{ - 1/4}} + O({b^{ - 1/2 + 2p}}),\end{align*}

which for some $b_0$ implies the result in (b).

5. Improving estimate of V near the boundary, using feedback

Our next goal is to use the above to get an improved estimate of V that is very accurate when not too far from the boundary. This will allow us to go further down the backward induction tree by estimating the leaf values $V(u + 1,b + 2m - 1)$ for larger m, for which $u + 1$ is more than 1/2 below the boundary and the value is no longer just the ratio. This is used in Section 6 to obtain the correct coefficient for the ${b^{ - 1/4}}$ term. We shall show V is approximately piecewise linear near the boundary: that is Theorem 5.1. As a first step, Lemma 4.4 showed that $V(u,b) = u/b$ if $\delta \le 1/2\,\, - \,0.43{b^{ - 1/4}}/\sqrt \pi $ for b sufficiently large, already an improvement over our previous upper bound $V(u,b) \le u/b + \alpha {\delta ^2}{b^{ - 3/2}}$ coming from the Brownian motion case. The Brownian upper bound is too big by about $\alpha {b^{-3/2}}/4$ when delta is near 1/2.

We introduce a bit more notation. In this section and the next, we will be doing some estimates, in a chained fashion, an unbounded number of times, which could cause a problem if just using big-O notation. The book [Reference Graham, Knuth and Pataschnik7] gives notation (attributed to de Bruijn) which will be convenient. We let $f(b) = L(g(b))$ mean that $\left| {f(b)} \right| \le \left| {g(b)} \right|$ for all b. In set language, $L(g) = \left\{ {f:\left| {f(b)} \right| \le \left| {g(b)} \right|\text{ for all } b} \right\}$ . Also, let ${L^ + }(g) = \left\{ {f \in L(g):f(b) \ge 0\text{ for all }b} \right\}$ . Similarly, we let ${O^ + }(g)$ be the non-negative members of O(g).

We assert that throughout Section 5, $n = \lfloor {\sqrt b /2} \rfloor $ . Lemma 4.4 (b) showed that using this n, $1/2 - 0.44{b^{ - 1/4}}/\sqrt\pi \le {\delta _0} \le \delta _0' \le 1/2 - 0.13{b^{ - 1/4}}/\sqrt\pi, b \ge {b_0}$ , where ${\delta _0}$ satisfies $C + {\delta _0}B + \delta _0^2A = 0$ , and $\delta _0'$ satisfies $C' + \delta _0'B + {\delta _0'^2}A = 0$ . A,B,C,C’ are defined in Lemma 4.3. We assert ${b_0} > 1600$ . Now $n \ge \sqrt b /2 - 1 = \sqrt b /2\left( {1 - 2/\sqrt b } \right)$ , so ${G_n} \le 1/\sqrt {\pi n} \le \sqrt {2.14/\pi } {b^{ - 1/4}} \le 0.83{b^{ - 1/4}}$ for $b \ge {b_0}$ , a bound that will be used several times in the following. In going further away from boundary, we will use the following lemma. Recall that ${V_E}(u,b) = V(u,b) - u/b$ .

Lemma 5.1. For $1/2 \le \delta \le {b^p}$ , with $p \le 1/10$ , $n = \lfloor {\sqrt b /2} \rfloor $ , and $b \ge {b_0}$ ,

\begin{align*}{V_E}(u,b) = \frac{\alpha }{{{b^{3/2}}}}\left( {2(\delta - 1/2) + L\left( {\max \{ {\delta ^2},1\} {b^{ - 1/4}}} \right)} \right) + {S_{LE}}(n,u,b).\end{align*}

Proof. Let $1/2 \le \delta \le {b^p}$ . The Value is the same as the tree sum because $\delta \ge \delta _0'$ , so u is below the stop boundary. $\delta _0' = 1/2 - \gamma '{b^{ - 1/4}}$ for some $0.13/\sqrt \pi \le \gamma ' \le 0.44/\sqrt \pi $ . Let $x = \delta - \delta _0' = \delta - 1/2 + \gamma '{b^{ - 1/4}}$ . From Lemma 4.3 in the lower bound case,

\begin{align*}{V_E}(u,b) - {S_{LE}} &\ge \alpha {b^{-3/2}}\left( {C' + B\delta + A{\delta ^2}} \right)\\ &= \alpha {b^{-3/2}}\left( {C' + B\delta _0' + A{\delta _0'}^2 + Bx + A{x^2} + 2Ax\delta _0'} \right)\\ &= \alpha{b^{-3/2}}\left( {Bx + A{x^2} + 2Ax\delta _0'} \right).\end{align*}

Now $B = 2\left( {1 - {B_1}{G_n}} \right) = 2 - 2\left( {1 + \alpha n/\sqrt b} \right){G_n} \ge 2 - \left( {2 + \alpha } \right){G_n}$ and $A = {G_n}$ , so $Bx + A{x^2} + 2Ax\delta _0'\ge 2(\delta - 1/2) + 2\gamma '{b^{ - 1/4}} +\left( { - \left( {2 + \alpha } \right)x + {x^2} + 2x\delta _0'} \right){G_n} = 2(\delta - 1/2) + 2\gamma '{b^{ - 1/4}} + x\left( {x - \alpha - 1 - 2\gamma '} \right){G_n}.$ But $x\left( {x - \alpha - 1 - 2\gamma '} \right)$ has a minimum of $ - {\left( {\alpha + 1 + 2\gamma '} \right)^2}/4 > - 1$ , so $x\left( {x - \alpha - 1 - 2\gamma '} \right){G_n} > - {b^{ - 1/4}}$ , so ${V_E}(u,b) - {S_{LE}} \ge \alpha {b^{-3/2}}\left( {2(\delta - 1/2) - {b^{ - 1/4}}} \right)$ .

Next, let $x = \delta - {\delta _0} = \delta - 1/2 + \gamma {b^{ - 1/4}}$ , where $0.13/\sqrt \pi \le \gamma \le 0.44/\sqrt \pi $ . From Lemma 4.3 in the upper bound case, similar to lower bound case, we get

\begin{align*}{V_E}(u,b) - {S_{LE}} \le \alpha {b^{-3/2}}( {Bx + A{x^2} + 2Ax{\delta _0}}).\end{align*}

Now $B \le 2 - 2{G_n}$ , so

\begin{align*}Bx + A{x^2} + 2Ax{\delta _0} &\le 2(\delta - 1/2) + 2\gamma {b^{ - 1/4}} - 2x{G_n} + ( {{x^2} + x}){G_n}\\ &= 2(\delta - 1/2) + 2\gamma {b^{ - 1/4}} + ({x^2} - x){G_n}. \end{align*}

If $\delta \le 1$ , then $x \le 1$ , so $2\gamma {b^{ - 1/4}} + ({x^2} - x){G_n} \le 2\gamma {b^{ - 1/4}} < {b^{ - 1/4}}$ . If $\delta > 1$ , then ${x^2} - x=(\delta - \delta_0)^2 - (\delta - \delta_0)=\delta^2 - \delta -\delta_0 (2\delta - \delta_0 -1) \le \delta^2 - 1.$ So $2\gamma {b^{ - 1/4}} + ({x^2} - x){G_n} \le \left( {.88/\sqrt \pi } \right){b^{ - 1/4}} + ({\delta ^2} - 1)\left( {\sqrt {2.14} /\sqrt \pi } \right){b^{ - 1/4}} \le {\delta ^2}{b^{ - 1/4}}$ .

Now we look at $S_{LE},$ first considering a few small ranges of delta, to establish a pattern. If $0 \le \delta \le 1/2$ , then from (4.4), ${S_{LE}} = 0$ . If $\delta \le \,{\delta _0}$ , $V(u,b) = u/b$ . Now assume ${\delta _0} \le \delta \le 1/2$ . Let $x = \delta - {\delta _0}$ , so $x \le \gamma {b^{ - 1/4}}$ . From Lemma 4.3, proceeding as in the above proof, ${V_E}(u,b) \le \alpha {b^{ - 3/2}}\left( {C + B\delta + A{\delta ^2}} \right) = \alpha {b^{ - 3/2}}\left( {Bx + A{x^2} + 2Ax{\delta _0}} \right).$ With this x, it is easy to see that $Bx + A{x^2} + 2Ax{\delta _0} \le {b^{ - 1/4}};$ we do not care to be any more precise than that. We get

(5.1) \begin{align}0 \le \delta \le 1/2 \text{ implies }{V_E}(u,b) = \frac{\alpha }{{{b^{3/2}}}}\left( {{L^ + }({b^{ - 1/4}})} \right),\quad b \ge b_0.\end{align}

When $1/2 \le \delta $ , the value V(u,b) is the tree sum, bigger than the ratio. The leaf sum part involves $V(u + 1,b + 2m + 1)$ , where $m \le n - 1$ . The distance of $u + 1$ below the boundary is $d = \max \left\{ {\alpha \sqrt {b + 2m + 1} \,\, - (u + 1),0} \right\}$ . But $\alpha \sqrt {b + 2m + 1} \,\, - (u + 1) = \alpha \sqrt b {\left( {1 + (2m + 1)/b} \right)^{1/2}} - \alpha \sqrt b + \delta - 1$ . Standard estimates give $1 + m{b^{ - 1}} \le {\left( {1 + (2m + 1)/b} \right)^{1/2}} \le 1 + m{b^{ - 1}} + {b^{ - 1}}/2$ , so $\alpha m{b^{ - 1/2}} \le \alpha \sqrt b {\left( {1 + (2m + 1)/b} \right)^{1/2}} - \alpha \sqrt b \le \alpha m{b^{ - 1/2}} + \alpha {b^{ - 1/2}}/2$ . Thus for $m \le n-1$ ,

(5.2) \begin{align}d = \left( {\alpha \sqrt {b + 2m + 1} \,\, - (u + 1)} \right)^+ = \left( {\frac{{\alpha m}}{{\sqrt b }} + \delta - 1 + {L^ + }\left( {\alpha {b^{ - 1/2}}/2} \right)} \right)^+.\end{align}

We will use this in several places below.

Consider $1/2 \le \delta \le 1$ . Since $\delta \le 1$ , $d \le \alpha m/\sqrt b + {L^ + }\left( {\alpha {b^{ - 1/2}}/2} \right)$ . But $m \le \sqrt b /2$ , so $d \le \alpha /2 + {L^ + }\left( {\alpha {b^{ - 1/2}}/2} \right)$ , and computation shows this is less than $1/2 - 0.44{b^{ - 1/4}}/\sqrt \pi $ if $b \ge {b_0}$ , assumed. So ${V_E}(u + 1,b + 2m + 1) = 0$ (that is, the value is just the ratio), and ${S_{LE}} = 0$ . Then Lemma 5.1 gives

(5.3) \begin{align}1/2 \le \delta \le 1 \text{ implies }{V_E}(u,b) = \frac{\alpha }{{{b^{3/2}}}}\left( {2(\delta - 1/2) + L({b^{ - 1/4}})} \right), b \ge {b_0}.\end{align}

Our improved estimate so far for $V_E(u,b)$ , $0 \le \delta \le 1$ , is piecewise linear plus an $O({b^{ - 1/4}})$ correction; zero for delta from 0 to 1/2, then slope $2 \alpha b^{-3/2}$ from 1/2 to 1. The piecewise-linear part is accurate to order ${b^{ - 1/4}}$ , relative to $b^{-3/2}$ . Compare this to our previous upper bound estimate from the Brownian value, $V_E \le {V_W}(u,b) - u{b^{ - 1}} \le \alpha {b^{ - 3/2}}{\delta ^2}$ , and about equal to that for delta small compared to $\sqrt b $ ; that is, quadratic for delta not too big. Our piecewise-linear function matches the quadratic one at $\delta = 0$ and $\delta = 1$ , and is tangent to the parabola at those points. At $\delta = 1/2$ , the piecewise-linear function is below the parabola by 1/4. That is where the upper bound from the Brownian is worst.

Now we show that same linear piece continues up to 3/2. Let $1 \le \delta \le 3/2$ . The distance of $u + 1$ from the boundary $\alpha \sqrt {b + 2m + 1}$ is, from (5.2), $d = \alpha m{b^{ - 1/2}} + \delta - 1 + {L^ + }\left( {\alpha {b^{ - 1/2}}/2} \right) \le \alpha /2 + 1/2 + {L^ + }\left( {\alpha {b^{ - 1/2}}/2} \right) \le 1$ , for $b \ge {b_0}$ . If also $d \ge 1/2$ , we can feed this into (5.3) with d in the role of $\delta $ and $b + 2m + 1$ in the role of b. This is the key feedback idea that will be used in proving Theorem 5.1. If $d < 1/2$ , (5.3) still holds if we replace $d - 1/2$ by $(d - 1/2)^+$ . So ${V_E}(u + 1,b + 2m + 1) = \alpha {\left( {b + 2m + 1} \right)^{ - 3/2}}\left( 2(d - 1/2)^+ + L({b^{ - 1/4}}) \right)$ . Note $d - 1/2 \le \alpha m{b^{ - 1/2}} + \alpha {b^{ - 1/2}}/2$ . Using (4.1),

\begin{align*}{S_{LE}} &\le \sum\limits_{m = 0}^{n - 1} {{2^{ - 2m - 1}}{C_m}\alpha {{\left( {b + 2m + 1} \right)}^{ - 3/2}}\left( {2\alpha m{b^{ - 1/2}} + \alpha {b^{ - 1/2}} + {b^{ - 1/4}}} \right)}\\ &\le \alpha {b^{ - 3/2}}\left( {n{G_n}2\alpha {b^{ - 1/2}} + \alpha {b^{ - 1/2}} + {b^{ - 1/4}}} \right) \le \alpha {b^{ - 3/2}}\left( {2{b^{ - 1/4}}} \right),\quad \mbox{for }b \ge {b_0}.\end{align*}

Using Lemma 5.1 and adding in this estimate of $S_{LE}$ gives the same answer as (5.3), except for increasing the error bound, which then covers (5.3) as well:

(5.4) \begin{align}1/2 \le \delta \le 3/2 \text{ implies }{V_E}(u,b) = \frac{\alpha }{{{b^{3/2}}}}\left( {2(\delta - 1/2) + L(5{b^{ - 1/4}})} \right),\quad b \ge {b_0}.\end{align}

Now that we see that the piecewise-linear function below and tangent to the quadratic at integers gets us this far, we guess that this pattern continues (up to some error), and Theorem 5.1 will prove this. By stepping along one 1/2 unit at a time, and feeding the result back into the previous step the way we did to get to 3/2, we will get our piecewise-linear estimate that will be good enough for our purpose.

Theorem 5.1. (Value near the boundary.) Assume $b \ge {b_0}$ . Then

\begin{align*}{V_E}&(u,b) =\frac{\alpha }{{{b^{3/2}}}}\left\{ {2j(\delta - 1/2) - j(j - 1) + L\left( {{M_j}{b^{ - 1/4}}} \right)} \right\},\quad j - 1/2 \le \delta \le j + 1/2,\\&0 \le j \le {b^{1/10}},\end{align*}

where j is integral, $M_0 = 1$ , and $M_j \le 5{j^3}, j \ge 1$ .

Figure 7 is a graph of the piecewise-linear function $2j(\delta - 1/2) - j(j - 1)$ for $0 \le \delta \le 2.5$ , compared to ${\delta ^2}$ . which is shown dashed. The j changes at the half-integer points. It is tangent to the parabola at integers, and 1/4 below the parabola at half-integers, with a straight-line graph between half-integer points.

Figure 7.

Graph of piecewise-linear approximation to scaled excess value, ${\alpha ^{-1}}{b^{3/2}}V_E$ , as a function of distance $\delta$ from boundary. The dashed curve is $\delta ^2$ , the quadratic approximation to the scaled excess value that comes from the Brownian upper bound.

The induction proof of this theorem uses the same feedback logic above that established the case $j=1$ at (5.4), but it is long and tedious, so we put it in Appendix C.

6. Proof of Theorem 1.1

In this section, let ${n_0} = \lfloor {\sqrt b /\alpha } \rfloor $ , so $\alpha {n_0}/\sqrt{b} = 1 - {L^ + }(\alpha {b^{ - 1/2}})$ and $\alpha /\sqrt{b} = 1/{n_0} - {L^ + }(\alpha {b^{ - 1}})$ ; we will see that replacing $\sqrt b /\alpha $ by its integer part is not going to matter, to the order of interest. Let $0 < p \le 1/10$ , to be decided later (in the end, it will be 1/12). Let $J + 1 = \left\lfloor {{b^p}} \right\rfloor $ , so $J \sim {b^p}$ . In this section, we are going to let $n = (J + 1){n_0}$ , so n will go to infinity faster than $\sqrt b $ , but only slightly. This larger n will cause some of the leaf values to be more than just the ratio. By dividing n into J stretches of size ${n_0}$ , we will get a linear approximation to the extra part of the leaf values, ${V_E}(u + 1,b + 2m + 1),$ on each stretch, via Theorem 5.1; that is Lemma 6.1.

Assume throughout this section that $b \ge {b_0}$ and $\delta = \alpha \sqrt b - u.$ In this section, where we seek the stop rule, we only consider $\delta$ in the range $1/2 - \,0.44{\pi ^{ - 1/2}}{b^{ - 1/4}} \le \delta \le 1/2 - \,0.13{\pi ^{ - 1/2}}{b^{ - 1/4}}$ , because we already know that the exact stop value occurs for some delta in this range. So $\delta = 1/2 - \gamma {b^{ - 1/4}}$ , for some $0.13{\pi ^{ - 1/2}} \le \gamma \le 0.44{\pi ^{ - 1/2}}$ .

We will estimate $V_E(u+1,b+2m+1)$ by using results from Section 5, when $u+1$ is further away from $\alpha \sqrt{b+2m+1}$ , but in this section we will use d instead of $\delta$ to denote that distance, since $\delta$ is fixed throughout this section to be $\alpha \sqrt b - u$ . To avoid confusion, note that we will be using results from Section 5 with d in the role of $\delta$ there, and $b+2m+1$ in the role of b.

By letting J go to infinity with b, at just the right rate, we can make the upper and lower bounds from Lemma 4.3 come together, to order $o({b^{ - 1/4}})$ , as we will see later.

Lemma 6.1. There exists K such that for $\delta = 1/2 - {L^ + }({b^{ - 1/4}})$ and for $j = 1,\ldots,J$ ,

\begin{align*}{V_E}(u + 1,b + 2m + 1) = \frac{\alpha }{{{b^{3/2}}}}\left\{ {2j\frac{m}{{{n_0}}} - j(j + 1) + L\left( {K{j^3}{b^{ - 1/4}}} \right)} \right\},\quad j{n_0} \le m \le (j + 1){n_0}.\end{align*}

To prove the main theorem, we have to put in the extra contribution ${S_{LE}}$ to the leaf sum that comes from the leaf values that exceed the ratio, when going down to $n=(J + 1){n_0}.$

Lemma 6.2.

\begin{align*}{S_{LE}} = \frac{\alpha }{{{b^{3/2}}}}\left\{ {{C_{LE}}{G_n} + O\left( {{J^{5/2}}{b^{ - 1/2}}} \right)} \right\},\end{align*}

where

\begin{align*}{C_{LE}} \,:\!= \frac{1}{3}{J^2} - \frac{1}{3}J - 4\zeta (\! -1/2)\sqrt {J + 1} + O(1).\end{align*}

Proof. Let ${S_{LE}}(j)$ denote the contribution to ${S_{LE}}$ from summing ${{2^{ - 2m - 1}}{C_m}} {V_E}(u + 1,b + 2m + 1)$ over m in the range $j{n_0} \le m < (j + 1){n_0}$ . Let $W(j) = {S_{LE}}(j){b^{3/2}}/\alpha$ . $W(0) = 0$ . Now let $j \ge 1$ . By Lemmas 6.1 and 4.1,

\begin{align*}W(j) &= \frac{{2j}}{{{n_0}}}\sum\limits_{m = j{n_0}}^{(j + 1){n_0} - 1} {m{2^{ - 2m - 1}}{C_m}} - \left( {j(j + 1) + L\left( {K{j^3}{b^{ - 1/4}}} \right)} \right)\sum\limits_{m = j{n_0}}^{(j + 1){n_0} - 1} {{2^{ - 2m - 1}}{C_m}} \\&= \frac{{2j}}{{{n_0}}}\left( {(j + 1){n_0}{G_{(j + 1){n_0}}} - j{n_0}{G_{j{n_0}}}} \right) - \left( {\frac{{2j}}{{{n_0}}} + j(j + 1) + L\left( {K{j^3}{b^{ - 1/4}}} \right)} \right)\left( {{G_{j{n_0}}} - {G_{(j + 1){n_0}}}} \right)\\&= 3j(j + 1){G_{(j + 1){n_0}}} - j(3j + 1){G_{j{n_0}}} + L\left( {(K + 1){j^3}{b^{ - 1/4}}} \right)\left( {{G_{j{n_0}}} - {G_{(j + 1){n_0}}}} \right).\end{align*}

Now

\begin{align*}{G_{j{n_0}}} - {G_{(j + 1){n_0}}} \le \dfrac{1}{{\sqrt \pi }}\left( {\dfrac{1}{{\sqrt {{n_0}j} }} - \dfrac{1}{{\sqrt {{n_0}(j + 1) + 1/2} }}} \right) \le \dfrac{{{j^{ - 3/2}}}}{{2\sqrt {\pi {n_0}} }} \le \dfrac{{{j^{ - 3/2}}{b^{ - 1/4}}}}{{2\sqrt \pi }},\end{align*}

and

\begin{align*}\dfrac{1}{{\sqrt {k\pi } }}\left( {1 - \dfrac{1}{{4k}}} \right) \le {G_k} \le \dfrac{1}{{\sqrt {k\pi } }} \implies \dfrac{{{G_{j{n_0}}}}}{{{G_{(J + 1){n_0}}}}} = \dfrac{{\sqrt {J + 1} }}{{\sqrt j }}\left( {1 + L\left( {\dfrac{1}{{4j{n_0}}}} \right)} \right).\end{align*}

Using the latter to express $G_{jn_0}$ and $G_{(j+1)n_0}$ in terms of $G_{(J+1)n_0}$ ,

\begin{align*}W(j) &= \left( \begin{array}{l}3j\sqrt {j + 1} \left( {1 + L\left( {{{(4(j + 1){n_0})}^{ - 1}}} \right)} \right)\\ - \sqrt j (3j + 1)\left( {1 + L\left( {{{(4j{n_0})}^{ - 1}}} \right)} \right)\end{array} \right)\sqrt {J + 1} {G_{(J + 1){n_0}}} + L\left( {(K + 1){j^{3/2}}{b^{ - 1/2}}} \right)\\&= \left( {3j\sqrt {(j + 1)} - 3(j - 1)\sqrt j - 4\sqrt j } \right)\sqrt {J + 1} {G_n} + L\left( {(K + 2){j^{3/2}}{b^{ - 1/2}}} \right).\end{align*}

We wrote it that way so that telescoping occurs for the first two terms when summing

\begin{align*}\sum\limits_{j=1}^{J}W(j)&={\left( {\sum\limits_{j = 1}^J {\left(3j\sqrt {(j + 1)} - 3(j - 1)\sqrt {j}\right)} - 4\sum\limits_{j = 1}^J {\sqrt j } } \right)\sqrt {J + 1} {G_n} + O\left( {J^{5/2}}{ b^{ - 1/2}} \right)}\\&= {\left( {3J\sqrt {(J + 1)} - 4\sum\limits_{j = 1}^J {\sqrt j } } \right)\sqrt {J + 1} {G_n} + O\left( {{J^{5/2}}{b^{ - 1/2}}} \right)}.\end{align*}

At this point we are through summing over an unbounded range, so we replaced the big-L with big-O notation, with no harm. There is a well-known asymptotic formula for the sum of square roots as a generalized harmonic number: ${H_J}^{( -1/2)} = \sum\limits_{j = 1}^J {\sqrt j } = 2{J^{3/2}}/3 + {J^{1/2}}/2 + \zeta ( -1/2) + O\left( {{J^{ - 1/2}}} \right)$ [Reference Graham, Knuth and Pataschnik7, p. 594]. And $\sqrt {J + 1} = {J^{1/2}} + {J^{ - 1/2}}/2 + O\left( {{J^{ - 3/2}}} \right)$ , so $4{H_J}^{( -1/2)}\sqrt {J + 1} = 8{J^2}/3 + 10J/3 + 4\zeta ( -1/2)\sqrt {J + 1} + O(1),$ and $3J(J + 1) - 4{H_J}^{( - 1/2)}\sqrt {J + 1} = {J^2}/3 - J/3 - 4\zeta ( -1/2)\sqrt {J + 1} + O(1).$ So ${S_{LE}} = \alpha {b^{ - 3/2}}\left\{ {\left( {{J^2}/3 - J/3 - 4\zeta ( -1/2)\sqrt {J + 1} } +O(1)\right){G_n} + O\left( {{J^{5/2}}{b^{ - 1/2}}} \right)} \right\}.$

We can now conclude the proof of Theorem 1.1.

Proof of Theorem 1.1. Go back to Lemma 4.3, to get the bounds on the tree sum using this n, adding in the ${S_{LE}}$ term from Lemma 6.2. Our assumption about J implies $n = O({b^{1/2 + p}})$ with $p \le 1/10$ , so Lemma 4.3 applies. For the upper bound, $TreeSum(n,u,b) - u/b \le \alpha {b^{ - 3/2}}\left\{ {C + \delta B + {\delta ^2}A} \right\} + {S_{LE}} = \alpha {b^{ - 3/2}}\left\{ {C^* + \delta B + {\delta ^2}A} \right\}$ , where ${C^*} = - \left( {1 - ({C_1} + {C_{LE}}){G_n}} \right) + O({J^{^{5/2}}}{b^{ - 1/2}}).$ Proceeding as done in the proof of Lemma 4.4,

\begin{align*}{\delta _0} = \dfrac{{ - B + B\left( {1 - 2A{C^*}B^{-2} - 2{A^2}{C^{*2}}B^{-4} - O({A^3})} \right)}}{{2A}} = - \dfrac{{{C^*}}}{B} - \dfrac{{{G_n}{C^{*2}}}}{{{B^3}}} + O({b^{ - 1/2 - p}}),\end{align*}

since $A={G_n} = O({b^{ - 1/4 - p/2}}).$ As in the earlier proof,

\begin{align*}- \dfrac{{{C^*}}}{B} - \dfrac{{G_n{C^{*2}}}}{{{B^3}}}= \dfrac{1}{2} - \dfrac{1}{2}\left( {{C_1} + {C_{LE}} - {B_1} + \dfrac{1}{4}} \right){G_n} + O({J^{5/2}}{b^{ - 1/2}}).\end{align*}

But $\dfrac{{\alpha n}}{{\sqrt b }}=J + 1 + O\left( {(J + 1){b^{ - 1/2}}} \right),$ so ${{C_1} + {C_{LE}} - {B_1} + \dfrac{1}{4}}={C_{LE}} + \dfrac{{\alpha n}}{{\sqrt b }} - \dfrac{{{\alpha ^2}{n^2}}}{{3b}} + \dfrac{1}{4}=\dfrac{1}{3}{J^2} - \dfrac{1}{3}J - 4\zeta ( -1/2)\sqrt {J + 1} - \dfrac{1}{3}{J^2} + \dfrac{1}{3}J + O(1) = - 4\zeta ( -1/2)\sqrt {J + 1} + O(1).$ The higher-power terms miraculously canceled, exposing the dominant zeta term. So

\begin{align*}{\delta _0}&=\frac{1}{2}- \frac{1}{2}\left( { - 4\zeta (\! -1/2)\sqrt {J + 1} + O(1)} \right)\frac{1}{{\sqrt {\pi (J + 1){n_0}} }}+ O\left( {{J^{5/2}}{b^{ - 1/2}}} \right)\\&= \frac{1}{2} - \frac{1}{2}\left( { - 4\zeta (\! -1/2) + O({J^{ - 1/2}})} \right)\frac{{\sqrt \alpha }}{{\sqrt \pi }}{b^{ - 1/4}} + O({J^{5/2}}{b^{ - 1/2}})\\&= \frac{1}{2} - \frac{{\left( { - 2\zeta (\! -1/2)} \right)\sqrt \alpha }}{{\sqrt \pi }}{b^{ - 1/4}} + O({J^{ - 1/2}}{b^{ - 1/4}}) + O({J^{5/2}}{b^{ - 1/2}}).\end{align*}

Since $J \sim {b^p}$ , the best we can do with this is to let $p = 1/12$ . We get

\begin{align*}{\delta _0} = \frac{1}{2} - \frac{{\left( { - 2\zeta (\! -1/2)} \right)\sqrt \alpha }}{{\sqrt \pi }}{b^{ - 1/4}} + O({b^{ - 7/24}}).\end{align*}

For the lower bound, to get $\delta _0'$ , the only change is that ${C_1}$ is replaced by $C_1'={C_1}-5/12,$ which does not affect the calculation to order $O({b^{ - 7/24}})$ . This completes the proof of Theorem 1.1.

As a final remark, we could say that we expected those higher terms to cancel, based on the idea that going further down the tree leads to less weight on the row and more weight on the leaves, where, at least for a while, the errors are quite small thanks to Theorem 5.1. But honestly, when that happened with just the right choice of p, we thanked Tyche rather than crediting our insight, since we did not really know if it would happen before doing it.