Perfect powers of successive naturals: Can you always reach a constant difference?

I was thinking about what happens if you take a sequence of consecutive squares, for example 1,4,9, 16. Taking the differences gives you another sequence, 7,5,3. And taking the differences between those numbers, you get 2,2--a constant. Through elementary first-semester algebra, you can easily verify that this works no matter where you start your series of squares. If you use 9801, 10000, 10201, you will get the same result--after the second round of subtraction, you'll end up with a final constant difference of 2. The same thing works for cubes with the final constant difference being 6, although an additional round of subtraction is needed, in turn requiring one more number in the original sequence. Similarly as with the squares, it's not difficult to show that it will work for any sequence of natural number cubes.

My question is this: Is there a more general principle at work here? Suppose I take a series of sequential naturals n, n+1, n+2, n+3,... , and raise them to any given positive integral power p. Can it then be shown that if I take the differences repeatedly, proceeding as above, that I will eventually reach a constant difference? And if so, does this principle have anything to do with "difference engine" computing?

Solution 1:

Yes.

The first difference of a polynomial of degree $d$ is a polynomial of degree $d-1$.

By induction, the $m$-th difference of a polynomial of degree $d$, when $m \le d$, is a polynomial of degree $d-m$.

Setting $m = d$, the $d-th$ difference of a polynomial of degree $d$ is a constant, which is exactly what you discovered.

You next step is a proof.

Solution 2:

If you start with the sequence of $n$-th powers, the $n$-th differences will all be $n!$.

For a function $f$ defined on the integers define $(\Delta f)(x)=f(x+1)-f(x)$. $\Delta$ is a linear operator on such functions: you can easily check that $\Delta\big(af(x)+bg(x)\big)=a\Delta f(x)+b\Delta g(x)$ for any such functions $f$ and $g$ and constants $a$ and $b$.

Now suppose that $\Delta^k(x^k)$ is the constant function $k!$ for each $k<n$. For any constant function $f$ the difference function $\Delta f$ is identically $0$, so $\Delta\big((\Delta^k(x^k)\big)=\Delta(k!)=0$ for each $k<n$, and it follows that $\Delta^m(x^k)=0$ whenever $k<n$ and $m>k$. In particular, $\Delta^{n-1}(x^k)=0$ for $k<n-1$. Thus,

$$\begin{align*} \Delta^n(x^n)&=\Delta^{n-1}\big((x+1)^n-x^n\big)\\ &=\Delta^{n-1}\left(\sum_{k=0}^n\binom{n}kx^k-x^n\right)\\ &=\Delta^{n-1}\left(\sum_{k=0}^{n-1}\binom{n}kx^k\right)\\ &=\sum_{k=0}^{n-1}\binom{n}k\Delta^{n-1}(x^k)\\ &=\binom{n}{n-1}\Delta^{n-1}(x^{n-1})\\ &=n(n-1)!\\ &=n!\;, \end{align*}$$

and by induction $\Delta^n(x^n)=n!$ for all $n$. (You can easily check the basis step of the induction.)

These so-called forward differences are a special case of divided differences, after which Babbage’s Difference Engine was named.