Cominciamo con qualcosa di facile. Prendiamo i numeri naturali:
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, …
Cancelliamo ora ogni secondo numero:
1, 3, 5, 7, 9, 11, 13, …
Facciamo infine le somme parziali:
1, 4, 9, 16, 25, 36, 49, …
Abbiamo così ottenuto la successione dei quadrati perfetti. Nulla di strano: sappiamo che i quadrati sono la somma dei numeri dispari, lo si vede disegnando gli gnomoni.

Ma ora andiamo più in là. Ripartiamo dai numeri naturali:
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, …
ma stavolta togliamo ogni terzo numero:
1, 2, 4, 5, 7, 8, 10, 11, 13, 14, 16, 17, 19, 20, …
Facciamo anche stavolta le somme parziali:
1, 3, 7, 12, 19, 27, 37, 48, 61, 75, 91, 108, 127, 147, …
Togliamo ogni secondo numero:
1, 7, 19, 37, 61, 91, 127, …
E rifacciamo infine le somme parziali:
1, 8, 27, 64, 125, 216, 343, …
Abbiamo ottenuto la successione dei cubi!

Come arrivare alle quarte potenze? Semplice: aggiungiamo come passo iniziale quello di eliminare ogni quarto numero prima di sommare.
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16 , …
1, 2, 3, 5, 6, 7, 9, 10, 11, 13, 14, 15, …
1, 3, 6, 11, 17, 24, 33, 43, 54, 67, 81, 96, …
Ora si fa il passaggio di cancellare ogni terzo numero e poi fare le somme parziali:
1, 3, 11, 17, 33, 43, 67, 81, …
1, 4, 15, 32, 65, 108, 175, 256, …
E finalmente si cancella ogni secondo numero e si fanno ancora una volta le somme parziali:
1, 15, 65, 175, …
1, 16, 81, 256, …
Et voilà!

Questa struttura aritmetica è stata trovata da Alfred Moessner, che nel 1951 scrisse un articolo intitolato “Eine Bemerkung über die Potenzen der natürlichen Zahlen” (un’osservazione delle potenze dei numeri naturali) e dimostrato lo stesso anno da Oskar Perron. Non so chi abbia deciso di chiamarlo “la magia di Moessner”, e non più banalmente il teorema di Moessner – nome comunque usato, come si vede su Wikipedia – ma secondo me un po’ di magia c’è davvero…

PS: se si eliminano i numeri triangolari, si ottiene invece la successione dei fattoriali!