### A better estimate of Kempner’s series

#### by Willie Wong

The Kempner series recently regained some notoriety due to a Saturday Morning Breakfast Cereal comic (the last panel). The observation first appeared in a 1914 American Mathematical Monthly article, in which it was shown that the series consisting of the usual harmonic series

but with all the terms, whose decimal expansion includes the digit ‘9’, removed, in fact converges to some number below 80. The original proof is given in the Wikipedia article linked above, so I will not repeat it. But to make it easier to see the idea: let us first think about the case where the number is expressed in base 2. In base 2, all the positive integers has the leading binary bit being 1 (since it cannot be zero). Therefore there are *no* binary positive numbers without the bit ‘1’ in its expansion. So the corresponding series converges trivially to zero. How about the case of the bit ‘0’? The only binary numbers without any ‘0’ bits are

.

So the corresponding series actually becomes

So somewhere from the heavily divergent harmonic series, we pick up a rapidly converging geometric series. So what’s at work here? Among all the n-bit binary numbers, exactly 1 has all bits not being 0. So the *density* of these kinds of numbers decays rather quickly: in base 2, there are numbers that are exactly n-bit long. So if a number has a binary representation that is exactly n bits long (which means that ), the chances that it is one of the special type of numbers is . This probability we can treat then as a *density*: replacing the discrete sum by the integral (calculus students may recognize this as the germ of the “integral test”) and replacing the by the density , we get the estimate

.

Doing the same thing with the original Kempner series gives that the chances a n-digit number does not contain the digit nine to be

The length of the decimal expansion of a natural number is basically . So the density we are interested in becomes

From this we can do an integral estimate

The integral can be computed using that

to get

Notice that this estimate is much closer to the currently known value of roughly 22.92 than to the original upper bound of 80 computed by Kempner.

Kempner’s estimate is a heavy overestimate because he performed a summation replacing every n-digit long number that does not contain the digit 9 by ; this number can be many times (up to 9) times smaller than the original number. Our estimate is low because among the n-digit long numbers, the numbers that do not contain the digit 9 are not evenly distributed: they tend to crowd in the front rather than in the back (in fact, we do not allow them to crowd in the back because none of the numbers that start with the digit 9 is admissible). So if in the original question we had asked for numbers that do not contain the digit 1, then our computation will give an overestimate instead since these numbers tend to crowd to the back.

Thanks for the post! The binary case does make the convergence seem much more intuitive, and the evaluation of that integral is pretty slick.