Newton's Identities

lean4-prooflinear-algebravisualization
pk=i=1k1(1)i1eipki+(1)k1kekp_k = \sum_{i=1}^{k-1} (-1)^{i-1} e_i p_{k-i} + (-1)^{k-1} k e_k

Statement

Let x1,,xnx_1, \ldots, x_n be indeterminates over a commutative ring RR. Define the elementary symmetric polynomials ek=i1<<ikxi1xike_k = \sum_{i_1 < \cdots < i_k} x_{i_1} \cdots x_{i_k} and the power sums pk=i=1nxikp_k = \sum_{i=1}^n x_i^k.

Newton's identities state: for 1kn1 \le k \le n,

pke1pk1+e2pk2+(1)k1ek1p1+(1)kkek=0.p_k - e_1 p_{k-1} + e_2 p_{k-2} - \cdots + (-1)^{k-1} e_{k-1} p_1 + (-1)^k k e_k = 0.

Equivalently (for knk \le n):

pk=i=1k1(1)i1eipki+(1)k1kek.p_k = \sum_{i=1}^{k-1} (-1)^{i-1} e_i p_{k-i} + (-1)^{k-1} k e_k.

Visualization

Roots {1,2,3}\{1, 2, 3\} — compute eke_k and pkp_k directly, then verify the identity for k=3k = 3:

QuantityFormulaValue
e1e_11+2+31+2+366
e2e_212+13+231\cdot2 + 1\cdot3 + 2\cdot31111
e3e_31231\cdot2\cdot366
p1p_111+21+311^1+2^1+3^166
p2p_212+22+321^2+2^2+3^21414
p3p_313+23+331^3+2^3+3^33636

Check k=3k = 3: p3e1p2+e2p13e3=36614+11636=3684+6618=0p_3 - e_1 p_2 + e_2 p_1 - 3e_3 = 36 - 6\cdot14 + 11\cdot6 - 3\cdot6 = 36 - 84 + 66 - 18 = 0. Confirmed.

Proof Sketch

  1. Generating function. The identity follows from comparing coefficients in the logarithmic derivative of i(1xit)\prod_i (1 - x_i t), which simultaneously encodes both eke_k (via expansion) and pkp_k (via its derivative).
  2. Combinatorial proof (Zeilberger). Define a weight function on pairs (A,j)(A, j) where A[n]A \subseteq [n] and j[n]j \in [n], and build an involution on the set of such pairs of weight kk. The involution shows the alternating sum of weighted counts is zero — this zero-sum is exactly Newton's identity.
  3. Recursive recovery. The identity lets you recover eke_k from p1,,pkp_1, \ldots, p_k and vice versa, establishing that both families generate the same ring of symmetric polynomials.

Connections

Newton's identities connect to Vieta's FormulasVieta's Formulasx1+x2++xn=an1anx_1 + x_2 + \cdots + x_n = -\frac{a_{n-1}}{a_n}Relating the coefficients of a polynomial to the sums and products of its rootsRead more → (which express eke_k in terms of roots of a polynomial) and to Fundamental Theorem of AlgebraFundamental Theorem of Algebradegp1,  pC[z]    z0C:p(z0)=0\deg p \geq 1,\; p \in \mathbb{C}[z] \implies \exists z_0 \in \mathbb{C}: p(z_0) = 0Every 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