Identités de Newton

En mathématiques, et plus particulièrement en algèbre, les identités de Newton (connues également sous le nom de formules de Newton-Girard) sont des relations entre deux types de polynômes symétriques, les polynômes symétriques élémentaires, et les sommes de Newton, c'est-à-dire les sommes de puissances des indéterminées. Évaluées aux racines d'un polynôme P à une variable, ces identités permettent d'exprimer les sommes des k-ièmes puissances de toutes les racines de P (comptées avec leur multiplicité) en fonction des coefficients de P, sans qu'il soit nécessaire de déterminer ces racines. Ces formules furent redécouvertes par Isaac Newton vers 1666, apparemment sans avoir eu connaissance du travail analogue d'Albert Girard en 1629. Elles ont des applications dans de nombreux domaines mathématiques, tels que la théorie de Galois, la théorie des invariants, la théorie des groupes, la combinatoire, et même dans des domaines non mathématiques, comme en relativité générale.

Page de garde de Invention en l'Algèbre nouvelle, Albert Girard, Amsterdam, 1629

Énoncé mathématique

modifier

Formulation en fonction des polynômes symétriques

modifier

Soient x1,...,xn des variables ; pour tout entier naturel k non nul, on note pk(x1,...,xn) la somme des puissances k-ièmes :

  (appelée somme de Newton)

et pour k ≥ 0, on note ek(x1,...,xn) le polynôme symétrique élémentaire qui est la somme des produits distincts de k variables distinctes parmi les n ; ainsi en particulier

 

Les identités de Newton s'écrivent alors :

 

relations vraies pour tout  . On obtient ainsi pour les premières valeurs de k :

 

La forme de ces relations ne dépend pas du nombre n de variables (mais le membre de gauche devient nul après la n-ième identité), ce qui permet de les énoncer comme des identités dans l'anneau des polynômes symétriques. Dans cet anneau, on a :

 

et ainsi de suite ; dans ce cas, le membre de gauche ne s'annule jamais.

Ces équations permettent d'exprimer par récurrence les ei en fonction des pk ; inversement, les réécrivant sous la forme :

 

on peut exprimer les pi en fonction des ek.

Application aux racines d'un polynôme

modifier

Les xi étant pris comme paramètres et non comme variables, considérons le polynôme unitaire en t ayant pour racines x1,...,xn :

 

où les coefficients ak sont donnés par les polynômes symétriques élémentaires des racines : ak =ek(x1,…,xn). Les sommes des puissances des racines (qu'on appelle encore des sommes de Newton),   peuvent alors être exprimées en fonction des coefficients du polynôme en appliquant les identités de Newton de proche en proche :

 

Application au polynôme caractéristique d'une matrice

modifier

Lorsque le polynôme du paragraphe précédent est le polynôme caractéristique d'une matrice A, les racines xi sont les valeurs propres de A (comptées avec leur multiplicité algébrique). Pour tout entier k, la matrice Ak a pour valeurs propres les xk
i
(avec les multiplicités correspondantes). Les coefficients du polynôme caractéristique de Ak sont alors donnés par les polynômes symétriques élémentaires en ces puissances xk
i
. En particulier, la somme des xk
i
est donnée par la trace de Ak : sk = tr(Ak).

Les identités de Newton relient donc les traces des Ak aux coefficients du polynôme caractéristique de A. Réciproquement, les utilisant pour calculer les polynômes symétriques élémentaires à partir des sommes de puissances, elles permettent donc de déterminer le polynôme caractéristique de A en calculant uniquement les puissances de A et leurs traces, puis en résolvant un système triangulaire d'équations. Le théorème de Cayley-Hamilton permet ensuite d'en déduire simplement la matrice inverse de A.

L'ensemble de ces calculs (réécrits sous une forme efficace) constitue l'algorithme de Faddeev-Leverrier, datant de 1840 ; une implémentation parallèle rapide de cet algorithme est due à Lazslo Csanky[1], en 1976. Il demande cependant des divisions (par des entiers), et n'est donc utilisable en général que dans des corps de caractéristique 0. L'implémentation de Csanky montre que ces différents calculs sont (dans ce cas) dans la classe de complexité NC.

Relation à la théorie de Galois

modifier

Pour un n donné, les polynômes symétriques élémentaires ek(x1, …, xn) pour k = 1, …, n forment une « base algébrique » de l'algèbre commutative des polynômes symétriques en x1,...,xn, c'est-à-dire que tout polynôme symétrique s'exprime en fonction polynomiale de ces polynômes symétriques élémentaires, et que cette expression est unique. Ce résultat général, connu sous le nom de théorème fondamental des polynômes symétriques, peut s'expliciter (à l'aide des identités de Newton) dans le cas des sommes de Newton. Ainsi, appliqué au polynôme unitaire   dont les coefficients ak sont considérés comme des paramètres, cela signifie que toute fonction polynomiale S(x1, …, xn) de ses racines peut s'écrire comme fonction polynomiale P(a1, …, an) de ses coefficients seuls, c'est-à-dire sans qu'il soit nécessaire de calculer ces racines. Cela peut également se déduire de la théorie de Galois, en voyant les ak comme éléments d'un corps de base ; les racines sont alors dans une extension de ce corps, et le groupe de Galois de cette extension les permute suivant le groupe symétrique entier ; le corps invariant par tous les éléments de ce groupe de Galois est le corps de base.

Réciproquement, les identités de Newton permettent d'exprimer les polynômes symétriques élémentaires en fonction des sommes de Newton, et de démontrer que ces sommes forment également une « base algébrique » de l'algèbre commutative des polynômes symétriques.

Démonstrations

modifier

Chaque identité peut être vérifiée directement par un simple calcul algébrique, mais le cas général demande une démonstration. Voici quelques pistes possibles :

À partir du cas particulier n = k

modifier

La k-ième identité de Newton en k variables s'obtient par substitution dans la formule définissant les ek :

 

Substituant xj à t, on a :

 

Sommant sur l'ensemble des j, on obtient :

 

(Les termes en i = 0 sont séparés de la somme parce que p0 n'est en général pas défini.)

Cette équation donne immédiatement l'identité cherchée. Les identités correspondant à n < k variables s'en déduisent en annulant les k − n variables restantes ; le cas où n > k se traite en remarquant que chaque monôme ne contient pas plus de k variables, et que ce monôme ne change pas si on annule les n − k autres variables ; il suffit alors d'utiliser l'identité de Newton correspondant à ces k variables.

Par identification de séries formelles

modifier

On peut également obtenir les identités de Newton à l'aide de manipulations formelles basées sur l'identité

 

reliant les racines et les coefficients d'un polynôme unitaire. Pour faciliter les calculs, on commence par substituer 1/t à t, et on multiplie les deux membres par tn, obtenant :

 

Remplaçant les ai par les polynômes symétriques, on obtient l'identité

 

Après différentiation par rapport à t, et multiplication par t, il vient :

 

où le polynôme de droite a d'abord été réécrit comme une fonction rationnelle, puis celle-ci développée en série formelle en t, et les coefficients de chaque t j regroupés pour obtenir une somme de puissances (la convergence de la série est en fait assurée pour t assez proche de zéro, mais comme on ne s'intéresse qu'aux coefficients des t j, cette question de convergence n'a pas d'importance réelle). Comparant les coefficients de tk des deux membres, on obtient

 ,

ce qui est la k-ième identité de Newton.

Comme somme télescopique d'identités

modifier

La dérivation suivante[2] est formulée dans l'anneau des polynômes symétriques, car alors les identités ne dépendent pas du nombre de variables. Pour un k > 0 fixé, on définit (pour 2 ≤ i ≤ k) la fonction symétrique r(i) comme la somme des monômes distincts de degré k obtenus en multipliant une variable élevée à la puissance i par k − i autres variables distinctes. En particulier, r(k) = pk ; le cas r(1) est exclu car les monômes n'auraient plus de variable jouant un rôle spécial. Tous les produits pieki peuvent être exprimés en fonction des r(j) (les cas extrêmes i=1 et i=k mis à part). On obtient  , puisque chaque produit des termes de gauche mettant en jeu des variables distinctes contribue à r(i), alors que ceux où les variables de pi figurent déjà parmi les variables du terme correspondant de eki contribuent à r(i + 1), et que tous les termes de droite sont ainsi obtenus une fois et une seule. Pour i = k, en multipliant par e0 = 1, on obtient trivialement  . Enfin le produit p1ek−1 (correspondant à i = 1) apporte des contributions à r(i + 1) = r(2) comme pour les autres valeurs de i < k, mais les contributions restantes sont égales à k fois chaque monôme de ek, puisque chacune des variables peut provenir du facteur p1 ; ainsi  .

La k-ième identité de Newton est alors obtenue en prenant la somme alternée de ces équations, laquelle est une somme télescopique : tous les termes de la forme r(i) disparaissent.

Identités analogues

modifier

De nombreuses familles d'identités analogues et étroitement reliées aux identités de Newton existent.

Utilisation des polynômes symétriques complètement homogènes

modifier

Notant hk le polynôme symétrique complètement homogène (en), c'est-à-dire la somme de tous les monômes de degré k, les sommes de puissances satisfont des identités semblables à celles de Newton, mais dont les termes sont tous positifs. Dans l'anneau des polynômes symétriques, elles s'écrivent  , pour tout k ≥ 1. Contrairement aux identités de Newton, le membre de gauche ne s'annule pas pour k assez grand, et le nombre de termes non nuls des membres de droite grandit indéfiniment. Pour les premières valeurs de k, on a

 

Ces relations peuvent se démontrer par un argument analogue à celui donné plus haut utilisant des séries formelles, mais en utilisant l'identité entre fonctions génératrices :

 .

En revanche, les autres démonstrations données précédemment ne peuvent s'adapter aisément à ces nouvelles identités.

Expression des polynômes symétriques élémentaires en fonction de sommes de Newton

modifier

Comme on l'a dit, les identités de Newton permettent d'exprimer par récurrence les polynômes symétriques élémentaires en fonction de sommes de Newton. Cela demande des divisions par des entiers, et ne peut donc être fait que dans l'anneau ΛQ des polynômes symétriques à coefficients rationnels :

 

et ainsi de suite. Pour un polynôme unitaire, ces formules expriment les coefficients en fonction des sommes de puissances des racines, en remplaçant chaque ei par ai et chaque pk par sk.

Expression des polynômes symétriques complètement homogènes en fonction de sommes de Newton

modifier

Les relations analogues concernant les polynômes symétriques complètement homogènes se développent de même, amenant aux équations :

 

où tous les termes sont positifs. Ces expressions correspondent exactement aux indicateurs de cycles des polynômes des groupes symétriques, en interprétant les sommes de Newton pi comme des indéterminées : le coefficient d'un monôme p1m1p2m2plml dans l'expression de hk est égal à la proportion de toutes les permutations de k ayant m1 points fixes, m2 cycles de longueur 2, …, et ml cycles de longueur l. Plus précisément, ce coefficient peut s'écrire 1/N, avec   ; N est le nombre de permutations avec une permutation fixée π ayant le type de cycle correspondant. Les expressions correspondant aux fonctions symétriques élémentaires ont des coefficients ayant les mêmes valeurs absolues, mais un signe égal à la signature de π, c'est-à-dire à (−1)m2+m4+….

Expression des sommes de Newton

modifier

Inversement, on peut exprimer les sommes de Newton en fonction des polynômes symétriques élémentaires, et ces expressions ont des coefficients entiers :

 

mais ces expressions ne semblent pas suivre de règle explicite. On voit cependant que le coefficient d'un monôme   dans l'expression de pk a le même signe que le coefficient du produit correspondant   dans l'expression de ek donnée plus haut, c'est-à-dire (−1)m2+m4+…. De plus, la valeur absolue du coefficient de M est la somme, prise sur l'ensemble des suites de polynômes symétrique élémentaires dont le produit est M, de l'indice du dernier polynôme de chaque suite : ainsi, le coefficient de e15e3e43 dans l'expression donnant p20 est  , puisque parmi les ordres distincts de cinq facteurs e1, un facteur e3 et trois facteurs e4, il y en a 280 se terminant par e1, 56 se terminant par e3, et 168 se terminant par e4.

Enfin, les identités concernant les polynômes complètement homogènes peuvent de même être inversées, amenant à :

 

Ces identités ont exactement la même forme que les précédentes, au signe près : le signe du monôme   est à présent −(−1)m1+m2+m3+….

Expressions sous forme de déterminants

modifier

Les expressions précédentes correspondant à la résolution de systèmes d'équations linéaires, elles peuvent se formuler explicitement à l'aide de déterminants, en utilisant la règle de Cramer. Par exemple, écrivant les identités de Newton sous la forme :

 

et en prenant  ,  ,  , …,  et   comme inconnues, on obtient pour cette dernière :

 

Les calculs sont analogues (mais un peu plus compliqués) pour les  , ou pour les expressions en fonction des polynômes symétriques complètement homogènes ; on obtient finalement[3] :

 

On peut remarquer[4] que la formule pour hn s'obtient en prenant le permanent de la matrice pour en au lieu de son déterminant, et plus généralement qu'une expression pour un polynôme de Schur quelconque peut être obtenue en prenant l'immanant correspondant de cette matrice.

Notes et références

modifier
(en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Newton's identities » (voir la liste des auteurs).
  1. (en) L. Csanky, « Fast parallel matrix inversion algorithms », SIAM J. Comput., vol. 5, no 4,‎ (lire en ligne [PDF]).
  2. (en) D. G. Mead, « Newton's Identities », The American Mathematical Monthly, Mathematical Association of America, vol. 99-8,‎ , p. 749–751 (DOI 10.2307/2324242, JSTOR 2324242).
  3. (en) Ian G. Macdonald, Symmetric functions and Hall polynomials, Oxford, The Clarendon Press, Oxford University Press, coll. « Oxford Mathematical Monographs », , viii+180 (ISBN 0-19-853530-9, MR 84g:05003), p. 20.
  4. (en) Dudley E. Littlewood, The theory of group characters and matrix representations of groups, Oxford, Oxford University Press, , viii+310 (ISBN 0-8218-4067-3), p. 84.

Voir aussi

modifier

Articles connexes

modifier

Bibliographie

modifier

Liens externes

modifier