A Magical Theorem Was Undiscovered for Thousands of Years
And is still kept hidden from students - why?
A couple of nights ago I was lying in bed watching some lectures, reading some articles, and so on, and at some point, I decided what I thought to be a wise decision at the time, namely that I would only watch one more math video before going to sleep.
I had absolutely no clue that the video I picked would keep me sleepless for hours. Not because it was long, but because I was shocked that I had never heard about the beautiful mathematics it contained.
Now, I have a Master’s degree in mathematics and I was so in love with the field (I still am) that I spent all my optional courses on subjects within mathematics as well. That amounts to about 40 different university-level mathematical disciplines under my belt and yet I had never heard of this theorem… NEVER!
This must mean that almost no one teaches it.
What is even more strange is the fact that the result itself went undiscovered for more than 2500 years. Not because it was too hard to prove (which is often the case), but because no one had discovered the pattern.
The discovery
The mathematics involved is so simple that school kids would be able to understand it yet there’s so much we still don’t know about it. Despite the extremely simple mathematics involved, the underlying truth stayed hidden from ancient times until 1951 when a mathematician by the name of Alfred Moessner made a miracle discovery.
To explain the magic let’s start with something very simple yet very pretty. Consider the partial sums of the odd natural numbers i.e.
1 = 1, 1+3 = 4, 1+3+5 = 9, 1+3+5+7 = 16, 1+3+5+7+9 = 25.
Do you spot the pattern? What we get are our square numbers i.e. 1=1², 4=2², 9=3²,… So partial sums of odd numbers are squares! The ancient Greeks knew of this and they even provided a genius visual proof of this fact.

In the above image, the colored balls represent the first 4 odd numbers and we clearly see that the sum of the first odd numbers gives a square number.
What Moessner did was simple, he took this fact and looked at it from a slightly different perspective. In particular, he framed this fact a bit more constructively. Almost algorithmically. His construction of the above is the following:
Imagine starting with all the natural numbers: 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,…
where the three dots mean “continue this pattern to infinity”. Then remove every second number. Then we get the odd numbers of course. So the new sequence of numbers is: 1, 3, 5, 7, 9, …
And take partial sums to get the squares 1, 4, 9,...
Well, yes we just saw that. But when we use this approach it is easier to generalize it (something we mathematicians like to do).
Now take all the natural numbers again but remove every third number to get the number sequence: 1, 2, 4, 5, 7, 8, 10, 11, 13,…
Take partial sums to get a new sequence: 1, 3, 7, 12, 19, 27, 37, 48, 61,…
Remove every second number to get: 1, 7, 19, 37, 61,…
Take partial sums again to get: 1, 8, 27, 64, 125,…
Do you recognize these numbers?… They are our cubes! 1=1³, 8=2³, 27=3³,…
This calculation can be structured in the following image:

This pattern continues. If we start by removing every n’th number from the natural numbers, take partial sums to form a new sequence, remove every n-1'th number from that sequence, and take partial sums, remove the n-2'th number, etc, we will eventually end up with the powers of n.
How come I had never heard of that? So simple. So beautiful.
But we are only just getting started because it turns out that there is much more to be found here. For example, if you remove other interesting sequences of numbers to begin with, what sequence will then be produced? Is there a general pattern?

In the above image, the triangular numbers are removed from the natural numbers, forming small triangles of numbers. By deleting numbers in the same way as before for each sequence, these triangles end up in the numbers 1, 2, 6, 24, 120, … and these are the factorials i.e. 2 = 2⋅1, 6 = 3⋅2⋅1, 24 = 4⋅3⋅2⋅1, etc.
Note that the triangular numbers are numbers of the form 1+2+3+⋅⋅⋅+n and they are converted into numbers of the form 1⋅2⋅3⋅⋅⋅n. Coincidence? No.
Also note that if we remove numbers of the form n, 2n, 3n, 4n, … where n > 1, then the output will be 1^n, 2^n, 3^n, 4^n, … so there is some kind of exponential relationship between the input and output here. It seems that in general if we multiply each number in the input sequence by a number m, the elements in the output sequence will all become raised to the power m.
Experimenting
To test this idea, I wrote a program in Python able to output these sequences based on an input sequence of numbers to be removed. I even considered other operations than addition.
If some of you would like to look for other interesting findings, then I have included the small code snippet and an example here:





