Fun with Newton polynomials

by t0rajir0u, Dec 10, 2007, 1:59 AM

The standard method of representing a polynomial $ P(x)$ is with coefficients attached to terms of the form $ x^k$, like so:

$ P(x) = a_0 + a_1 x + a_2 x^2 + ... + a_n x^n$.

Such a representation of a polynomial is very useful if one wants to do calculus or apply Vieta's formulas to find out something about the roots, but there are number-theoretic reasons we might prefer a different representation.

- A polynomial which returns integer values at integer arguments need not have integer coefficients. For example, $ P(x) = \frac {x(x - 1)}{2}$ is integer at integer values of $ x$.

- The forward difference operator $ \Delta P(x) = P(x + 1) - P(x)$ is useful for problems that only give values for $ P$ at integers, but it is messy to calculate using the above representation. (More generally, in an interpolation problem we do not necessarily want to invert a Vandermonde matrix.)

So what to do? In the language of vector spaces, we are using $ \{ 1, x, x^2, ... \}$ as a basis for the space of polynomials over some ring (the "standard basis"). This basis makes representing the derivative (as a linear operator) fairly straightforward. We would like a correspondingly straightforward representation of the forward difference operator for solving certain problems.

The answer is to use a different basis. For a given non-negative integer $ k$, denote $ {x \choose k} = \frac {x(x - 1)...(x - (k - 1))}{k!}$ (where $ {x \choose 0} = 1$, the empty product). These are binomial coefficients with a general argument (which are used in the General Binomial Theorem), or the Newton polynomials. We can use

$ \{ {x \choose 0}, {x \choose 1}, {x \choose 2}, ... \}$

as a polynomial basis, which addresses both of the concerns we had above.


Lemma 1: Given a polynomial $ P(x) = a_0 + a_1 {x \choose 1} + a_2 {x \choose 2} + ... + a_n {x \choose n}$ in "Newton form," $ P(x)$ is integer for integer values of $ x$ if and only if $ a_i$ is an integer for $ i = 0, 1, 2, ... n$.

Proof: The left direction is easy. If $ a_i$ is integer then so is $ a_i {x \choose i}$ (for integer $ x$). To prove the right direction, it suffices to consider $ x = 0, 1, 2, ... n$ and work inductively. We claim by strong induction that $ a_i$ is an integer for $ i = 0, 1, 2, ... n$. When $ i = 0$,

$ P(0) = a_0 \in \mathbb{Z}$

so our base case is true. Suppose we have shown that $ a_0, a_1, ... a_k$ are integers. Then

$ P(k + 1) = a_0 + a_1 {k + 1 \choose 1} + ... + a_k {k + 1 \choose k} + a_{k + 1} \in \mathbb{Z}$

and since the first $ k$ terms of the expression are integers, $ a_{k + 1}$ is an integer. This concludes the proof.

The next lemma makes dealing with polynomials that take integer values at integer arguments even easier.

Lemma 2: If $ P(x) = a_0 + a_1 {x \choose 1} + a_2 {x \choose 2} + ... + a_n {x \choose n}$, then

$ \Delta P(x) = a_1 + a_2 {x \choose 1} + a_3 {x \choose 2} + ... + a_n {x \choose n - 1}$

Proof: This is an immediate consequence of the combinatorial identity $ {x + 1 \choose k} - {x \choose k} = {x \choose k - 1}$, which can be proven for general $ x$ algebraically.


Now for some problems.

Problem 1: Find a closed formula for $ \sum_{i = 1}^{n} i^k$ for a particular small $ k$.

Solution: A demonstration for $ k = 4$. Other small cases follow similarly.

Our goal is to write $ i^4$ in the form $ a_0 + a_1 {i \choose 1} + a_2 {i \choose 2} + a_3 {i \choose 3} + a_4 {i \choose 4}$. To do so, it's enough to consider $ i = 0, 1, 2, 3, 4$ and apply a similar process as in the proof of Lemma 1.

$ 0^4 = a_0 \implies a_0 = 0$.
$ 1^4 = a_0 + a_1 \implies a_1 = 1$.
$ 2^4 = 16 = a_0 + 2a_1 + a_2 \implies a_2 = 14$.
$ 3^4 = 81 = a_0 + 3a_1 + 3a_2 + a_3 \implies a_3 = 36$.
$ 4^4 = 256 = a_0 + 4a_1 + 6a_2 + 4a_3 + a_4 \implies a_4 = 24$.

Hence

$ i^4 = {i \choose 1} + 14 {i \choose 2} + 36 {i \choose 3} + 24 {i \choose 4}$

and

$ \sum_{i = 1}^{n} i^4 = {n + 1 \choose 2} + 14 {n + 1 \choose 3} + 36 {n + 1 \choose 4} + 24 {n + 1 \choose 5}$

by the Hockeystick identity (which can be viewed as a corollary of Lemma 2).


Remark: In general, the coefficients of $ x^k$ in Newton form are given by

$ x^k = \sum_{j = 0}^{k} S(k, j) j! {x \choose j}$

where $ S(k, j)$ are the Stirling numbers of the second kind. This gives the identity

$ \sum_{i = 1}^{n} i^k = \sum_{j = 0}^{k} S(k, j) j! {n + 1 \choose j + 1}$

which is given in Engel; compare with Faulhaber's formula.


Problem 2: A polynomial $ P(x)$ of degree $ 2007$ satisfies $ P(k) = 2^k$ for $ k = 0, 1, 2, ... 2007$. Find $ P(2008)$.

Solution: If we try to calculate the coefficients of $ P(x)$ in Newton form, we might recall the identity $ {n \choose 0} + {n \choose 1} + ... + {n \choose n - 1} + {n \choose n} = 2^n$. Indeed, verify that

$ P(x) = {x \choose 0} + {x \choose 1} + {x \choose 2} + ... + {x \choose 2007}$

is the unique polynomial that satisfies the problem conditions. (When $ x < 2007$, the later terms cancel to zero.) Thus

$ P(2008) = {2008 \choose 0} + {2008 \choose 1} + ... + {2008 \choose 2007} = \boxed{ 2^{2008} - 1 }$.

Alternately, if we attempt to generally find a polynomial $ P_n(x)$ such that $ P_n(k) = 2^k, k = 0, 1, 2, ... n$, noting that $ \Delta P_n(x) = P_{n - 1}(x)$ makes a proof by induction quite simple.


Remark: All we really did in this problem was interpolate. It was the choice of Newton polynomials rather than the standard basis that made the interpolation easy to calculate. Note that Lagrange interpolation would've taken much longer to simplify.


Practice Problem (USAMO 1984): A polynomial $ P(x)$ of degree $ 3n$ satisfies

$ P(0) = P(3) = ... = P(3n) = 2$
$ P(1) = P(4) = ... = P(3n - 2) = 1$
$ P(2) = P(5) = ... = P(3n - 1) = 0$

and $ P(3n + 1) = 730$. Find all possible values of $ n$.

Comment

2 Comments

The post below has been deleted. Click to close.
This post has been deleted. Click here to see post.
I like the thinking going on here...I have looked at this similarly [but not as explicitly]. $ p(x) = a_n\prod (x - r_i)$ is another nice way to view polynomials.

The problem can also be solved with finite differences...

http://mathworld.wolfram.com/FiniteDifference.html

I think you hinted at this, but it is very explicitly stated that you find the coefficients in the 'leftmost' slanted column of the finite differences table.

The table is very easy to do with $ \{2^k\}_k$....


BTW congrats on getting into MIT. If anybody that I know deserves to go there, it is you.

by Altheman, Dec 15, 2007, 8:43 PM

The post below has been deleted. Click to close.
This post has been deleted. Click here to see post.
"In the language of vector spaces, we are using $ \{1, x, x^2, \dots \}$ as a basis for the space of polynomials over some ring"... It should be Field, not Ring. Because, you would then speak in the language of modules. :P

by Anonymous, Jun 20, 2008, 5:32 PM

Various mathematical thoughts, ranging from problem-solving techniques to attempts to tie together related concepts.

avatar

t0rajir0u
Archives
+ March 2009
+ October 2007
+ May 2007
Shouts
Submit
  • orz $~~~~$

    by clarkculus, Jan 10, 2025, 4:13 PM

  • Insanely auraful

    by centslordm, Jan 1, 2025, 11:17 PM

  • Fly High :(

    by Siddharthmaybe, Oct 22, 2024, 8:34 PM

  • Dang it he is gone :(( /revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive/revive.

    by 799786, Aug 4, 2022, 1:56 PM

  • annoying precision

    by centslordm, May 16, 2021, 7:34 PM

  • rip t0rajir0u

    by OlympusHero, Dec 5, 2020, 9:29 PM

  • Shoutbox bump xD

    by DuoDuoling0, Oct 4, 2020, 2:25 AM

  • dang hes gone :(

    by OlympusHero, Jul 28, 2020, 3:52 AM

  • First shout in July

    by smartguy888, Jul 20, 2020, 3:08 PM

  • https://artofproblemsolving.com/community/c2448

    has more.

    -πφ

    by piphi, Jun 12, 2020, 8:20 PM

  • wait hold up 310,000
    people visited this man?!?!??

    by srisainandan6, May 29, 2020, 5:16 PM

  • first shout in 2020

    by OlympusHero, Apr 4, 2020, 1:15 AM

  • in his latest post he says he moved to wordpress

    by MelonGirl, Nov 16, 2019, 2:43 AM

  • Please revive!

    by AopsUser101, Oct 30, 2019, 7:10 PM

  • first shout in october fj9odiais

    by bulbasaur., Oct 14, 2019, 1:14 AM

128 shouts
Tags
About Owner
  • Posts: 12167
  • Joined: Nov 20, 2005
Blog Stats
  • Blog created: Dec 5, 2006
  • Total entries: 48
  • Total visits: 321831
  • Total comments: 202
Search Blog