Alexander Bass

Plus Nine Divide Five

Here’s a fun little math problem: take any number and add nine, and divide by five. We’ll use 1515 as a starting example

(15+9)5=4.8 \frac{(15 + 9)}{5} = 4.8

Then add nine, and divide by five

(4.8+9)5=2.76 \frac{(4.8+ 9)}{5} = 2.76

repeat

(2.76+9)5=2.352 \frac{(2.76+ 9)}{5} = 2.352 (2.352+9)5=2.2704 \frac{(2.352+ 9)}{5} = 2.2704

[5 more skipped]

2.2500013 2.2500013\dots

The value seems to converge to 2.252.25. Let’s try another starting number: 58-58 and repeat the process. And let’s be a bit more formal with our notation. We can describe each result of adding nine and dividing by five to be a step in a sequence. So we can describe the next in a sequence given the last in the sequence

Sn+1=(Sn+9)5 S_{n+1} = \frac{(S_n + 9)}{5}

And we’ll define S0=58S_0 = -58

S0=58S1=9.8S2=0.16 S3=1.768S4=2.1536S5=2.23072S6=2.246144S7=2.249S8=2.250 \begin{aligned} S_0 &=-58& S_1 &= -9.8 & S_2 &= -0.16 \\\ S_3 &= 1.768& S_4&=2.1536&S_5&=2.23072 \\ S_6&= 2.246144&S_7&=2.249\dots&S_8&=2.250\dots \end{aligned}

Funnily enough, even with a different starting value, the process converges to the same value of 2.252.25. Try a few starting numbers yourself; you’ll see that they all converge to 2.252.25.

What’s the big deal with 2.252.25? Well, let’s make it a fraction

94 \frac{9}{4}

Nine fourths is pretty similar to the rule used for the problem add 9 divide by 5. Could it be that the value converges to a fraction of the value you add ( 99) over one less than the value you divide by ( 515-1)?

Now, we can redefine the problem: add 7 divide by 2. If the hypothesis is correct, we should see all starting values converge to 71\frac{7}{1}. We redefine the sequence

Sn+1=(Sn+7)2 S_{n+1} = \frac{(S_n+ 7)}{2}

Then try it out

S0=78S1=42.5S2=24.75S3=15.875S4=11.4375S5=9.21875S6=8.109375S7=7.55469S8=7.27734S9=7.13867S10=7.06934S11=7.03467 \begin{aligned} S_0 &=78& S_1 &= 42.5 & S_2 &= 24.75 \\ S_3 &= 15.875& S_4&=11.4375&S_5&=9.21875 \\ S_6&= 8.109375&S_7&=7.55469\dots&S_8&=7.27734\dots \\ S_9 &= 7.13867\dots&S_{10}&=7.06934\dots &S_{11}&= 7.03467\dots \end{aligned}

If you were to continue on you would see that the value converges to 77.

Let’s notate the problem to not depend on any 77’s or 22’s

Sn+1=Sn+ab S_{n+1} = \frac{S_n+a}{b}

And with what we learned from the two examples above, we can state our hypothesis

Sn+1ab1 S_{n+1} \approx \frac{a}{b-1}

The next iteration in the sequence is approximately equal to aa over bb minus one. This approximation will become more accurate with more steps. We could even be fancy and write it with limits

limnSn=ab1 \lim_{n \to \infty} S_n = \frac{a}{b-1}

But how do we prove it?

Into the weeds

Instead of evaluating each step in the sequence with numbers, what if we were to do it with algebra?

Setting everything up we have the definition of the sequence

Sn+1=Sn+abS_{n+1}=\frac{S_{n}+a}{b}

Given that S0S_0 is some unimportant constant

Lets evaluate the first few steps with algebra. Of course S0=S0S_0 = S_0 but the following ones are more interesting

S1=S0+ab=S0b+abS2=S1+ab= S0+ab+ab=S0b2+ab2+abS3=S2+ab= S0+ab+ab+ab=S0b3+ab3+ab2+ab \begin{aligned} S_1 &= \frac{S_0 + a}{b} &&&= \frac{S_0}{b} + \frac{a}{b} \\ S_2 &= \frac{S_1 + a}{b} &= &~\frac{\frac{S_0 + a}{b} + a}{b} &= \frac{S_0}{b^2} + \frac{a}{b^2} + \frac{a}{b} \\ S_3 &= \frac{S_2 + a}{b} &= &~\frac{\frac{\frac{S_0+a}{b}+a}{b}+a}{b} &= \frac{S_0}{b^3} + \frac{a}{b^3} + \frac{a}{b^2} + \frac{a}{b} \end{aligned}

After dividing through all the repeated fractions, we get a relatively clean pattern where, the starting value is divided by a certain amount of bb’s which is added to the sum of aa divided by increasing powers of bb. Lets factor out the form of the equation into two parts

Sn=S0bn+a(1b1+1b2++1bn1+1bn) S_n = \frac{S_0}{b^n} +a\left(\frac{1}{b^{1}} + \frac{1}{b^{2}} + \dots + \frac{1}{b^{n-1}} + \frac{1}{b^n} \right)

We can simplify (or complicate depending on how you look at it) the sum of reciprocal powers of bb by using a the summation notation

Sn=S0bn+ak=1n1bk S_{n}=\frac{S_0}{b^n}+ a\sum^n_{k=1}{\frac{1}{b^k}}

Given that our goal is to find the value SnS_n converges to with an ever-growing nn, let’s take the limit

limn(S0bn+ak=1n1bk) \lim_{n\to\infty} \left( \frac{S_0}{b^n}+ a\sum^n_{k=1}{\frac{1}{b^k}}\right)

With an ever increasing nn we can see that the term S0bn\frac{S_0}{b^n} will become ever smaller and smaller and with that, we’ve proven that the starting value doesn’t matter. The term which has S0S_0 will go to zero as nn goes to infinity.

We still haven’t proven our hypothesis of

limnSn=ab1 \lim_{n \to \infty} S_n = \frac{a}{b-1}

But we can simplify the equation a bit by removing the term which approaches zero, and putting the limit into the upper bound of the sum

ak=11bk a\sum^{\infty}_{k=1}{\frac{1}{b^k}}

Unless we can untangle that infinite sum of a increasing powers, we can’t yet prove it. Through this method, and please forgive me, but this route is somewhat of a dead end road (I’ll explain why I say somewhat later).

The simple way

So, here’s the problem in case you forgot it

Sn+1=Sn+ab S_{n+1} = \frac{S_n +a}{b}

If this sequence does converge to one specific value, then there should be an SnS_n such that Sn+1=SnS_{n+1} = S_n . Or in plainer words, there should be a value which we can add aa to, then divide by bb and get that same value back. We can setup the equation to say such

Sn=Sn+abS_n = \frac{S_n+a}{b}

Multiply both sides by bb

Sn×b=Sn+a S_n \times b = S_n + a

subtract SnS_n from both sides

Sn×(b1)=a S_n\times(b-1) = a

Solve for SnS_n

Sn=ab1 S_n = \frac{a}{b-1}

So now (if we are willing to assume that there is only one value the process converges on) we have proven that

limnSn=ab1 \lim_{n \to \infty} S_n = \frac{a}{b-1}

Or, the process will converge to be the fraction aa over bb minus one.

But we have proven much more than that! Allow me to restate some equations which we have seen to be valid

limnSn=ak=1(1b)klimnSn=ab1 \begin{aligned} \lim_{n \to \infty} S_n &= a\sum^{\infty}_{k=1}{\left(\frac{1}{b}\right)^k}\\ \lim_{n \to \infty} S_n &= \frac{a}{b-1}\\ \end{aligned}

Then we can set them to be equal to each other

ab1=ak=1(1b)k \frac{a}{b-1} = a\sum^{\infty}_{k=1}{\left(\frac{1}{b}\right)^k}

We can then divide out the aa

1b1=k=1(1b)k \frac{1}{b-1} = \sum^{\infty}_{k=1}{\left(\frac{1}{b}\right)^k}

And now − all from playing a silly number game − we have a simple identity for what a complicated infinite sum evaluates to. I hope you now understand why I took the time to go down the ‘dead end road.’

Closing

The infinite sum identity discovered in this article is nothing new. It’s well known and called the Geometric Series. I’ve been curious about add nine, and divide by five for more than a year at this point − bouncing around in the back of my mind. Only in the past few months did I solve the problem. I first made a computer program that brute forced plenty of different values just to get an idea of the problem, then I went on to solve it in a very similar way to how I’ve written it out today. I’m still a bit fascinated that a problem as simple as add nine, and divide by five could reveal a seemingly obscure identity for an infinite sum.

Technical notes

(bonus stuff for nerds and pedants)

The values of aa and bb within the following statement are assumed to be positive and non-zero

Sn+1=Sn+ab S_{n+1} = \frac{S_n+a}{b}

And in general, any time you see a fraction in this article, it is assumed that division by zero is not considered.


The series I called the geometric series is not quite a real geometric series

n=1(1b)n \sum^{\infty}_{n=1}{\left(\frac{1}{b}\right)^n}

The real one has a lower bound of n=0n=0 and is typically not expressed as a fraction

n=0pn \sum^{\infty}_{n=0}p^{n}

We define p=1bp = \tfrac{1}{b}

These two sums are equivalent except for the fact that the real geometric series has one more term at n=0n=0. We can find that term by evaluating p0p^0 which equals 11. Therefor, we can state that the ‘fake’ geometric series is one less than the real one, and then set them to be equal to each other.

n=0pn=n=1(1b)n+1 \sum^{\infty}_{n=0}p^{n} = \sum^{\infty}_{n=1}\left(\frac{1}{b}\right)^n +1

And as established earlier, we know the solution to the ‘fake’ geometric sum, so we can then solve for the real one.

n=0pn=1b1+1 \sum^{\infty}_{n=0}p^{n} = \frac{1}{b-1}+1

Let’s substitute out bb in favor of pp to make the equation clearer and then crunch some algebra

n=0pn=11p1+1=p1p+1=p+1p1p=11p \sum^{\infty}_{n=0} p^{n} = \frac{1}{\frac{1}{p} -1}+1 = \frac{p}{1-p}+1 = \frac{p+1-p}{1-p} = \frac{1}{1-p}

Which yields us the standard form you may be familiar with

n=0pn=11p \sum^{\infty}_{n=0} p^{n} = \frac{1}{1-p}

If you catch any mistakes in my reasoning, math, or notation: please let me know.