Friday, July 27, 2007

Newton's Identities: Proof of Newton's Formula

In a previous blog, I talked about Newton's Identities and the formula that Sir Isaac Newton used to generate these identities.

In today's blog, I will show a proof for Newton's formula. Interestingly, the proof revolves around a peculiar definition. The proof is taken from Jean-Pierre Tignol's Galois' Theory of Algebraic Equations and the definition in question is denoted as τ(a,b).

If you are reading Tignol's book, you may wonder why I am using same notation (σk and sk) in the opposite way. I do this in order to stay consistent with Harold Edwards who I referenced in my previous blogs.

τ(a,b) is defined as the sum of unique a-combinations of the n roots of a polynomial labeled x1, x2, ..., xn where one of the root is set to the power of b [For review of why a polynomial of degree n necessarily has n roots, see here for discussion on the fundamental theorem of algebra which guarantees this result]. This means that there are two different conditions that need to be considered. When b ≥ 2, regardless of which number is put to the power of b, there are n*C(n-1,a-1) terms in this equation. When b = 1, there are still C(n,a) possible terms in this equation.

If a ≥ 2, then n*C(n,a-1) is greater than C(n,a) since C(n,a) = (n/a)*C(n-1,a-1)

[Here's the reason why: (n/a)*C(n-1,a-1) = (n/a)*(n-1)!/[(n - 1 -[a-1])!(a-1)!] = n!/[(a)[n - 1 - a + 1]!(a-1)!] = n!/[a!][n-a] = C(n,a)]

With these two conditions in mind, we can now look at the definitions that are needed for the proof.

Definition 1: τ(a,b)

for b ≥ 2:

τ(a,b) = ∑ (i=1, n) ∑ (j=1, C(n-1,a-1)) xibxfj(1)*xfj(2)*...*xfj(a-1)

for b = 1:

τ(a,1) = ∑ (j=1,C(n,a)) xfj(1)*xfj(2)*...*xfj(k)

where:

C(n-1,a-1) = (n-1)!/{(a-1)!(n-1-[a-1])!}

fj(x) is a function where:
  • the range of fj(x) is {1, 2, ..., n}
  • fj is such that for all x,y if x ≠ y, then fj(x) ≠ fj(y)
  • each fj is a unique (a-1)-way selection of n-1 integers (not including i) so that we can assume that fj(1) is less than fj(2), etc. [since there is only one permutation of any selection that has this characteristic]
Example:

For n=3 with {x1, x2, x3 }

C(n-1,a-1) = C(2,1) = 2!/1!(2-1)! = 2

So there are n*C(n-1,a-1) =3*2 = 6 terms

τ(2,2) =
x12x2 + x12x3 + x22x1 +
x22x3 +
x32x1 +
x32x2

Definition 2: sk

sk = ∑ (i=1, n) xik

Example:

For n=3 with {x1, x2, x3},

s4 = x14 + x24 + x34

Lemma 1: τ(1,b) = sb

Proof:

(1) τ(1,b) = ∑ (i=1,n) xib using Definition 1 above.

(2) Then τ(1,b) = sb using Definition 2 above.

QED

Definition 3: σi

σk = ∑ (i=1,C(n,k)) xfi(1)*xfi(2)*...*xfi(k)

where:

C(n,k) = (n)!/{(k)!(n-k)!}

fi(x) is a function where:
  • the range of fi(x) is {1, 2, ..., n}
  • fi is such that for x,y if x ≠ y, then fi(x) ≠ fi(y)
  • each fi is a unique k-way selection of n integers so that we can assume that fi(1) is less than fi(2), etc. [since there is only one permutation of any selection that has this characteristic]
Example: σ2

For n=3 with {x1, x2, x3 }

C(n,k) = C(3,2) = 3!/2!(3-2)! = 3

So there are 3 terms

σ2 = x1*x2 + x1*x3 + x2*x3

Lemma 2: τ(a,1) = σa

Proof:

This is clear since for b=1, τ(a,1) and σa have the same definition.

QED

Lemma 3:

For a less than n and b greater than 2:

τ(a,b) = σasb-1 - τ(a+1,b-1)

Proof:

(1) τ(a,b) = ∑ (i=1, n) ∑ (j=1, C(n-1,a-1)) xibxfj(1)*xfj(2)*...*xfj(a-1) [See Definition 1 above]

= ∑ (i=1, n) ∑ (j=1, C(n-1,a-1)) xi(b-1)xixfj(1)*xfj(2)*...*xfj(a-1)

(2) σk = ∑ (i=1,C(n,k)) xfi(1)*xfi(2)*...*xfi(k) [See Definition 3 above]

(3) sb-1 = ∑ (i=1, n) xi(b-1)

(4) σasb-1 = [ ∑ (i=1,C(n,a)) xfi(1)*xfi(2)*...*xfi(a)][ ∑ (j=1, n) xj(b-1)]

(5) τ(a+1,b-1) = ∑ (i=1, n) ∑ (j=1, C(n-1,a)) xib-1xfj(1)*xfj(2)*...*xfj(a)

(6) Now each term in step #4 consists of the following form:

xjb-1*xfi(1)*...*xfi(a)

(7) Each term can then be categorized into two cases:

Case I: xj is the same as one of the values in xfi(1), xfi(2), ..., xfi(a)

Case II: xj is not the same as any of the values in xfi(1), xfi(2), ..., xfi(a)

(8) If C1 is the sum of all the terms that fall in Case I and C2 is the sum of all terms that fall into Case II, we can see that:

σasb-1 = C1 + C2

(9) Now, it is clear that for a less than n and b greater than 2:

C1 = τ(a,b) [See step #1 above]

(10) It is also clear that for a less than n and b greater than 2:

C2 = τ(a+1,b-1)

(11) Since C1 = σasb-1 - C2, it follows that:

τ(a,b) = σasb-1 - τ(a+1,b-1)

QED

Lemma 4:

if a is less than n:

τ(a,2) = σas1 - (a+1)σa+1

Proof:

(1) τ(a,2) = ∑ (i=1, n) ∑ (j=1, C(n-1,a-1)) xi2xfj(1)*xfj(2)*...*xfj(a-1) = = ∑ (i=1, n) ∑ (j=1, C(n-1,a-1)) [xixfj(1)*xfj(2)*...*xfj(a-1) ]*xi

(2) σa = ∑ (i=1,C(n,a)) xfi(1)*xfi(2)*...*xfi(a)

(3) s1 = ∑ (i=1, n) xi

(4) σas1 = [∑ (i=1,C(n,a)) xfi(1)*xfi(2)*...*xfi(a)][∑ (j=1, n) xj]

(5) Now each term in step #4 consists of the following form:

xj*xfi(1)*...*xfi(a)

(7) Each term can then be categorized into two cases:

Case I: xj is the same as one of the values in xfi(1), xfi(2), ..., xfi(a)

Case II: xj is not the same as any of the values in xfi(1), xfi(2), ..., xfi(a)

(8) If C1 is the sum of all the terms that fall in Case I and C2 is the sum of all terms that fall into Case II, we can see that:
σas1 = C1 + C2

(9) Now, it is clear that for a less than n:

C1 = τ(a,2) [See step #1 above]

(10) It is also clear that for a less than n:

C2 = (a+1)σa+1 since:

(a) σ(a+1) = ∑ (i=1,C(n,a+1)) xfi(1)*xfi(2)*...*xfi(a)*xfi(a+1)

(b) Now, removing all the terms that are case I, it is clear that there a+1 ways to reach each term in step (a).

For example, the combination x1*x2*...*xa*xa+1

This can be reached from:

x1 * (x2x3*...*xa+1) x2 * (x1x3*...*xa+1) ... xa+1*(x1x2*...*xa)

We can make the same argument for each possible term in C2

(11) Since C1 = σas1 - C2, it follows that:

τ(a,2) = σas1 - (a+1)σa+1

QED

Lemma 5:

if b ≥ 2, then:

τ(n,b) = σnsb-1

Proof:

(1) τ(n,b) = ∑ (i=1, n) ∑ (j=1, C(n-1,a-1)) xibxfj(1)*xfj(2)*...*xfj(a-1)

(2) Since b ≥ 2, we have:

τ(n,b) = x1bx2*...*xn + x2bx1*...*xn + ... + xnbx1*...*xn-1 =

= x1b-1x1x2*...*xn + x2b-1x1x2*...*xn + ... + xnb-1x1*...*xn =

= (x1b-1 + x2b-1 + ... + xnb-1)(x1x2*...xn)

(3) σn = x1*x2*...*xn

(4) sb-1 = x1b-1 + x2b-1 + ... + xnb-1

(5) It follows from step #2, that:

τ(n,b) = σnsb-1

QED

Theorem 6: Newton's Formula



Proof:

(1) τ(1,k) = sk [See Lemma 1 above]

(2) Using Lemma 3 above if k ≥ 3, we have:

τ(1,k) = σ1sk-1 - τ(2,k-1) = sk

(3) Assuming that k-1 ≥ 3, Lemma 3 gives us:

τ(2,k-1) = σ2sk-2 - τ(3,k-2)

(4) Following this pattern, if k ≤ n, we eventually get the following:

sk = σ1sk-1 - σ2sk-2 + ... + (-1)kτ(k-1,2)

(5) Using Lemma 4 above:

τ(k-1,2) = σk-1s1 - (k)σk

(6) Thus, we have:

sk = σ1sk-1 - σ2sk-2 + ... + (-1)kσk-1s1 + (-1)k+1k

(7) Subtracting sk from both sides and then multiplying -1 gives us:

0 = sk - σ1sk-1 + σ2sk-2 + ... + (-1)k-1σk-1s1 + (-1)kk

(8) If k is greater than n, then we eventually get to:

sk = σ1sk-1 - σ2sk-2 + ... + (-1)n+1τ(n,k+1-n)

(9) Using Lemma 5 above gives us:

sk = σ1sk-1 - σ2sk-2 + ... + (-1)n+1σnsk-n

(10) Subtracting sk from both sides and then multiplying by -1 gives us:

0 = sk - σ1sk-1 + σ2sk-2 + ... + (-1)nσnsk-n

(11) So, combining step #7 and step #10, if we assume that:

σ0 = 1 and if k is greater than n, σk=0, then:



QED

References

Monday, July 16, 2007

Newton's Identities: Newton's Formula

Newton's identities (see here for an introduction to Newton's identities) are relationships between the roots of a cubic polynomial and its coefficients. They were first presented by Albert Girard but were presented in even greater detail independently by Sir Isaac Newton.

In today's blog, I will present Newton's formula and show how it leads to Newton's identities. In a future blog, I will provide the proof for the formula.

Definition 1: Newton's Formula for Newton's Identities

For a polynomial of degree n, with roots { r1, ..., rn }:

For any integer i, let si be the sum based on the roots:

si = r1i + r2i + ... + rni.

For any integer j, Let σj be the elementary symmetric polynomial for j in n [see here for details if needed on elementary symmetric polynomials]

Then, Newton's formula is:



where σ0 = 1 and if k is greater than n, σk=0.

We can now use it to build some equations for si and σj.

If k=1, then we have:

s1 = σ1

If k ≥ 2, then we have:

sk = ∑ (i=1, k-1) (-1)i+1sk-iσi + (-1)k+1k

Using the above formula gives us:

s1 = σ1

s2 = s1σ1 - 2σ2

s3 = s2σ1 - s1σ2 + 3σ3

s4 = s3σ1 - s2σ2 + s1σ3 - 4σ4

s5 = s4σ1 - s3σ2 + s2σ3 - s1σ4 + 5σ5

...

Now, building each formula based on the previous formula gives us the following formulas for sk in term so of σi:

s1 = σ1

s2 = s1σ1 - 2σ2 = (σ11 - 2σ2 = σ12 - 2σ2

s3 = s2σ1 - s1σ2 + 3σ3 = (σ12 - 2σ21 - (σ12 + 3σ3 = σ13 -3σ1σ2 + 3σ3

s4 = s3σ1 - s2σ2 + s1σ3 - 4σ4 = σ1( σ13 -3σ1σ2 + 3σ3) - σ212 - 2σ2) + σ31) - 4σ4 =

= σ14 - 3σ12σ2 + 3σ1σ3 - σ12σ2 + 2σ22 + σ1σ3 - 4σ4 =

= σ14 -4σ12σ2 + 4σ1σ3 + 2σ22 - 4σ4


...

Further, we can use these same formulas for each σi so that:

σ1 = s1

σ2 = (1/2)[s1σ1 - s2]

σ3 = (1/3)[s3 - s2σ1 + s1σ2]

σ4 = (1/4)[-s4 + s3σ1 - s2σ2 + s1σ1

σ5 = (1/5)[s5 - s4σ1 + s3σ2 - s2σ3 + s1σ4]


Now, I will show these equations can be used to derive Newton's identities for cubic polynomials (that is, where n = 3). [See here for review of Newton's identities] where I am assuming an equation of the following form:

x3 + bx2 + cx + d = 0

Here are the justifications for each formula presented previously.

Identity 1:
r + s + t = -b

Proof:

σ1 = r + s + t
[See Definition 1, here]

σ1 = (-1)1(b) = -b
[See Lemma 1, here]

QED

Identity 2: r2 + s2 + t2 = b2 - 2c

Proof:

s2 = r2 + s2 + t2 [See Definition 1 above]

s2 = σ12 - 2σ2 [See formula above]

σ12 - 2σ2 = ((-1)1b)2 - 2(-1)2(c) = b2 - 2c.

QED

Identity 3: r3 + s3 + t3 = -b3 + 3bc - 3d

Proof:

s3 = r3 + s3 + t3 [See Definition 1 above]

s3 = σ13 -3σ1σ2 + 3σ3 [See formula above]

σ13 -3σ1σ2 + 3σ3 = [(-1)1b]3 - 3(-1)1b(-1)2c + 3(-1)3d =

= -b3 +3bc -3d.

QED

Identity 4: rs + rt + st = c

Proof:

σ2 = rs + rt + sr + st +tr + ts [See Definition 1, here]

σ2 = (-1)2c = c

QED

Identity 5: r2s + r2t + s2r + s2t + t2r + t2s = -bc + 3d

Proof:

(rs + rt + sr + st +tr + ts)(r + s + t) - 3rst = r2s + r2t + s2r + s2t + t2r + t2s

(rs + rt + sr + st +tr + ts)(r + s + t) - 3rst = σ2σ1 - 3σ3 = (-1)1b(-1)2c - 3*(-1)3d = -bc + 3d.

QED

Identity 6: r3s + r3t + s3r + s3t + t3r + t3s = b2c - 2c2 - bd

Proof:

(r3 + s3 + t3)(r + s + t) - (r4 + s4 + t4) = r3s + r3t + s3r + s3t + t3r + t3s

(r3 + s3 + t3)(r + s + t) - (r4 + s4 + t4) =( s3)(σ1) - s4

( s3)(σ1) - s4 = (σ13 -3σ1σ2 + 3σ3)(σ1) - (σ14 - 4σ12σ2 + 4σ1σ3 + 2σ22 - 4σ4) =

= σ12σ2 - σ1σ3 - 2σ22 + 4σ4 = [(-1)b]2(-1)2c - 2[(-1)2c]2 - [(-1)b(-1)3d] =

= b2c - 2c2 - bd + 0 = b2c - 2c2 - bd

QED

Identity 7: r2s2 + r2t2 + s2t2 = c2 - 2bd

Proof:

r2s2 + r2t2 + s2t2 = (1/2)(r2 + s2 + t2)(r2 + s2 + t2) - (1/2)(r4 + s4 + t4) - =

= (1/2)s2*s2 - (1/2)s4 =

= (1/2)(σ12 - 2σ2)(σ12 - 2σ2) - (1/2)( σ14 - 4σ12σ2 + 4σ1σ3 + 2σ22 - 4σ4) =

= (1/2)σ14 - 2σ12σ2 + 2σ22 - (1/2)σ14 + 2σ12σ2 - 2σ1σ3 - σ22 + 2σ4 =

= σ22 - 2σ1σ3 + 2σ4 =

= (c)2 - 2*(-1)(b)(-1)(d) + 2*(0) = c2 -2bd.

QED

Identity 8: r3s2 + r3t2 + s3r2 + s3t2 + t3r2 + t3s2 = -bc2 + 2b2d + cd

Proof:

(1) (r2s2 + r2t2 + s2t2)(r + s + t) - (rst)(rs + rt + st) =
r3s2 + r3t2 + s3r2 + s3t2 + t3r2 + t3s2

(2) (r2s2 + r2t2 + s2t2)(r + s + t) - (rst)(rs + rt + st) =

=(σ22 - 2σ1σ3 + 2σ4 )1) - (σ3)(σ2) =

= σ1σ22 - 2σ12σ3 + 2σ1σ4 - σ2σ3 =

= (-1)b(c)2 - 2(b)2(-d) + 2(-b)(0) - (c)(-d) =

= -bc2 + 2b2d + cd

QED


Identity 9: r3s3 + r3t3 + s3t3 = c3 - 3bcd + 3d2

Proof:

(1) r3s3 + r3t3 + s3t3 = (r2s2 + r2t2 + s2t2)(rs + rt + st) - (r3s2t + r3st2 + s3r2t + s3t2r + t3r2s + t3s2r )

(2) (r2s2 + r2t2 + s2t2)(rs + rt + st) - (r3s2t + r3st2 + s3r2t + s3t2r + t3r2s + t3s2r ) = (c2 - 2bd)(σ2) - (bcd - 3d2) =

=
(c2 - 2bd)(c) - (bcd - 3d2) = c3 - 2bcd -bcd + 3d2 =

= c3 - 3bcd +3d2

QED

Identity 10: rst = -d

Proof:

σ3 = rst

σ3 = (-1)3d = -d

QED

Identity 11: r2st + s2rt + t2rs = bd

Proof:

r2st + s2rt + t2rs = rst(r + s + t) = σ31 = (-1)(d)(-1)(b) = bd

QED


Identity 12: r3st + s3rt + t3rs = -b2d + 2cd

Proof:

(rst)(r2 + s2 + t2) = r3st + s3rt + t3rs

(rst)(r2 + s2 + t2) = σ3s2 =

= σ3(
σ12 - 2σ2) = (-d)(b2 -2c) = -b2d + 2cd.

QED

Identity 13: r2s2t + r2st2 + rs2t2 = -cd

Proof:

(rst)(rt + rs + st) = r2s2t + r2st2 + rs2t2

(rst)(rt + rs + st) = σ3σ2 = (-d)(c) = -cd.

QED
Identity 14: r3s2t + r3st2 + s3r2t + s3t2r + t3r2s + t3s2r = bcd - 3d2

Proof:

(rst)[(rs + rt + st)(r + s + t) - 3rst] = r3s2t + r3st2 + s3r2t + s3t2r + t3r2s + t3s2r
(rst)[(rs + rt + st)(r + s + t) - 3rst] = σ3[(σ2)(σ1) - 3σ3] =
σ3[(σ2)(σ1) - 3σ3] = (-1)d[c(-b) - 3(-1)d] = (-d)[-bc + 3d] = bcd - 3d2

QED

Identity 15: r3s3t + r3t3s + s3t3r = -c2d + 2bd2

Proof:

(rst)[(1/2)(r2 + s2 + t2)(r2 + s2 + t2) - (1/2)(r4 + s4 + t4) ] = r3s3t + r3t3s + s3t3r

(rst)[(1/2)(r2 + s2 + t2)(r2 + s2 + t2) - (1/2)(r4 + s4 + t4) ] =

= σ3[(1/2)(s2)(s2) - (1/2)s4 ] =

= σ3[(1/2)(σ12 - 2σ2)(σ12 - 2σ2) - (1/2)( σ14 - 4σ12σ2 + 4σ1σ3 + 2σ22 - 4σ4)] =

= σ3[(1/2)σ14 - 2σ12σ2 + 2σ22 - (1/2)σ14 + 2σ12σ2 - 2σ1σ3 - 2σ22 + 2*0)] =

= σ322 -2σ1σ3] = (-1)d[(c)2 - 2*(-1)b(-1)d] =

= (-d)[c2 - 2bd] = -c2d + bd2.

QED

Identity 16: r2s2t2 = d2

Proof:

r2s2t2 = (rst)2 = (σ3)2 = [(-1)3d]2 = d2

QED

Identity 17: r3s2t2 + s3r2t2 + t3r2s2 = -bd2

Proof:

(rst)(rst)(r + s + t) = r3s2t2 + s3r2t2 + t3r2s2

(rst)(rst)(r + s + t) = σ321 = (-b)d2= -bd2

QED

Identity 18: r3s3t2 + r3t3s2 + s3t3r2 = cd2

Proof:

(rs + rt + st)(rst)(rst) =
r3s3t2 + r3t3s2 + s3t3r2

(rs + rt + st)(rst)(rst) = σ232 = (c)(d)2 = cd2

QED

Identity 19: r3s3t3 = -d3

Proof:

r3s3t3 = (rst)3 = (σ3)3 = (-d)3 = -d3

QED

References