Newton's Identities
Statement
Let be indeterminates over a commutative ring . Define the elementary symmetric polynomials and the power sums .
Newton's identities state: for ,
Equivalently (for ):
Visualization
Roots — compute and directly, then verify the identity for :
| Quantity | Formula | Value |
|---|---|---|
Check : . Confirmed.
Proof Sketch
- Generating function. The identity follows from comparing coefficients in the logarithmic derivative of , which simultaneously encodes both (via expansion) and (via its derivative).
- Combinatorial proof (Zeilberger). Define a weight function on pairs where and , and build an involution on the set of such pairs of weight . The involution shows the alternating sum of weighted counts is zero — this zero-sum is exactly Newton's identity.
- Recursive recovery. The identity lets you recover from and vice versa, establishing that both families generate the same ring of symmetric polynomials.
Connections
Newton's identities connect to Vieta's FormulasVieta's FormulasRelating the coefficients of a polynomial to the sums and products of its rootsRead more → (which express in terms of roots of a polynomial) and to Fundamental Theorem of AlgebraFundamental Theorem of AlgebraEvery non-constant polynomial with complex coefficients has at least one root in ℂRead more → (where symmetric polynomials arise as invariants under permutation of roots). The combinatorial involution technique also appears in Inclusion–Exclusion Principle.
Lean4 Proof
-- Mathlib: Mathlib.RingTheory.MvPolynomial.Symmetric.NewtonIdentities
-- The main theorem is `MvPolynomial.mul_esymm_eq_sum` (Newton's identity in sum form)
-- and `MvPolynomial.psum_eq_mul_esymm_sub_sum` (solved for p_k).
-- We alias the Mathlib statement directly.
open MvPolynomial in
theorem newton_identities_mul_esymm (σ R : Type*) [CommRing R] [Fintype σ]
[DecidableEq σ] (k : ℕ) :
k * esymm σ R k =
∑ a ∈ Finset.antidiagonal k with a.1 ∈ Set.Ioo 0 k,
(-1 : MvPolynomial σ R) ^ (a.fst + 1) * esymm σ R a.fst * psum σ R a.snd +
(-1) ^ k * psum σ R k :=
mul_esymm_eq_sum σ R k