09.12.2012 Views

Concrete mathematics : a foundation for computer science

Concrete mathematics : a foundation for computer science

Concrete mathematics : a foundation for computer science

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

448 ASYMPTOTICS<br />

We might as well do Ez now, since it is so easy:<br />

=2 = x(&&T)<br />

k>l<br />

This is the telescoping series (1 -;)+(;-$)+($-&)+... =l.<br />

Finally, X1 gives us the leading term of S,, the coefficient of Inn in (9.53):<br />

=1 = x k($2 - &).<br />

k>l<br />

Thisis (l-i)+(i-$)+(G-&)+... = $+$+$+- =HE’ =7r2/6. (If<br />

we hadn’t applied summation by parts earlier, we would have seen directly<br />

that S, N xk3,(lnn)/k2, because H,t-, -H,tmlP1 N Inn; so summation by<br />

parts didn’t help us to evaluate the leading term, although it did make some<br />

of our other work easier.)<br />

Now we have evaluated each of the E’s in (g.53), so we can put everything<br />

together and get the answer to Golomb’s problem:<br />

S, = glnn+,-&+0(h),<br />

Notice that this grows more slowly than our original hand-wavy estimate of<br />

C(logn)‘. Sometimes a discrete sum fails to obey a continuous intuition.<br />

Problem 6: Big Phi.<br />

Near the end of Chapter 4, we observed that the number of fractions in<br />

the Farey series 3,, is 1 + (#J (n) , where<br />

O(n) = q(l) +(p(2) +...+cP(n);<br />

and we showed in (4.62) that<br />

@(n) = i 1 p(k) ln/k1 11 + n/k1 .<br />

k21<br />

(9.55)<br />

Let us now try to estimate cD(n) when n is large. (It was sums like this that<br />

led Bachmann to invent O-notation in the first place.)<br />

Thinking BIG tells us that Q(n) will probably be proportional to n2.<br />

For if the final factor were just Ln/k] instead of 11 + n/k], we would have<br />

(0(n)( < i xka, [n/k]’ 6 i xk>,(n/k)2 = $n2, because the Mobius function<br />

p(k) is either -1, 0, or +l: The additional ‘1 + ’ in that final factor<br />

adds xka, p(k) Ln/k] ; but this is zero <strong>for</strong> k > n, so it cannot be more than<br />

nH, = O(nlog n) in absolute value.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!