Polynôme de Tutte

Un article de Wikipédia, l'encyclopédie libre.

Le polynôme de Tutte, aussi appelé polynôme dichromatique ou polynôme de Tutte–Whitney, est un polynôme invariant de graphes dont les valeurs expriment des propriétés d'un graphe. C'est un polynôme en deux variables qui joue un rôle important en théorie des graphes et en combinatoire. Il est défini pour tout graphe non orienté et contient des informations liées à ses propriétés de connexité.

L'importance de ce polynôme provient des informations qu'il contient sur le graphe . Étudié au départ dans le cadre de la théorie algébrique des graphes comme une généralisation des problèmes d'énumération liés à la coloration de graphes, il contient diverses spécialisations à d'autres disciplines, comme le polynôme de Jones en théorie des nœuds, les fonctions de partitions liées au modèle de Potts (en) en physique statistique, le polynôme énumérateur des poids en théorie des codes, le polynôme de fiabilité en théorie des réseaux. Tous peuvent être exprimés comme des spécialisations du polynôme de Tutte. Il est aussi à la source de divers problèmes algorithmiques en informatique théorique. L'interprétation combinatoire des polynômes de Tutte est en étroite relation avec l’énumération d'objets combinatoires par des méthodes de langages formels et séries formelles non commutatives[1]

Les polynômes de Tutte ont plusieurs nom et définitions équivalents. Un polynôme de Tutte est équivalent au rang polynomial de Whitney, au polynôme dichromatique de Tutte et au random cluster model de Fortuin–Kasteleyn par des transformations simples. C'est essentiellement une série génératrice comptant les ensembles d'arêtes d'une taille de composantes connexes donnés, avec une généralisation naturelle aux matroïdes. C'est également l'invariant de graphes le plus général définissable par une récurrence de type suppression-contraction. Plusieurs livres de théorie des graphes ou de matroïdes consacrent des chapitres entiers aux polynômes de Tutte[2],[3],[4].

Le polynôme est le polynôme de Tutte du graphe taureau. La ligne rouge indique l'intersection avec le plan , équivalent au polynôme chromatique.

Définitions et exemples[modifier | modifier le code]

Définitions équivalentes[modifier | modifier le code]

Soit un graphe non orienté, avec ensemble de sommets et ensemble d'arêtes . Pour , on note le graphe partiel qui a les mêmes sommets que , et ses arêtes dans .

Definition.- Le polynôme de Tutte d'un graphe est le polynôme défini par :

,

est le nombre de composantes connexes du graphe . L'exposant du deuxième facteur est le nombre cyclomatique de .

On peut donner une formulation légèrement différente de la même définition en posant

,

est le rang de Whitney du graphe . La fonction génératrice des rangs de Whitney est définie par :

.

Les deux fonctions sont équivalentes par un simple changement de variables. On a en effet :

.

Le polynôme dichromatique de Tutte résulte d'une autre transformation simple. On a :

.

Définition originale.- La définition originale de de Tutte est équivalente, mais moins facile à formuler. Pour un graphe connexe , on a

est le nombre d'arbres couvrants d'« activité interne » et d' « activité interne » [5].

Contraction-suppression.- Une troisième définition est par récurrence sur la suppression d'arêtes. La contraction d'arête d'un graphe produit le graphe obtenu en fusionnant les sommets et et en supprimant l'arête . On note le graphe où l'on a seulement supprimé l'arête , sans fusionner ses deux extrémités. Le polynôme de Tutte est défini par récurrence par

selon que est boucle un isthme, ou ni l'un ni l'autre, avec la condition initiale

si ne contient pas d'arête.

Le « random cluster model » de la mécanique statistique dû à Fortuin et Kasteleyn 1972 fournit encore une autre définition équivalente[6] : Le polynôme

est lié à par la transformation [7] :

Propriétés[modifier | modifier le code]

Le polynôme de Tutte se factorise pour les composantes connexes : Si est l'union disjointe de graphes et , alors

.

Si est un graphe planaire et si dénote son graphe dual, alors

.

En particulier, le polynôme chromatique d'un graphe planaire est le polynôme de flot de son dual.

Exemples[modifier | modifier le code]

Deux graphes isomorphes ont même polynôme de Tutte, mais la réciproque est fausse. Par exemple, le polynôme de Tutte de tout arbre ayant arêtes est .

Un polynôme de Tutte peut être donné sous forme de table, en présentant dans un tableau le coefficient de en ligne et colonne . Par exemple, le polynôme de Tutte du graphe de Petersen, qui est :

est donné par la table suivante :

0 36 84 75 35 9 1
36 168 171 65 10
120 240 105 15
180 170 30
170 70
114 12
56
21
6
1

Un autre exemple, est le polynôme de Tutte du graphe octaédrique et :

Note historique[modifier | modifier le code]

L'intérêt de W. T. Tutte pour la formule de contraction-suppression remonte à ses études undergraduate au Trinity College de Cambridge, motivé par les rectangles parfaits (en) et les arbres couvrants. Il a utilisé souvent la formule dans ses recherches et se demandait « s'il existait d'autres invariants de graphes intéressants, aussi invariants par isomorphisme, avec des formules semblables »[8]. Ronald M. Foster avait déjà observé que le polynôme chromatique était de cette nature, et Tutte commençait à en découvrir d'autres. Il appelait au début W-function, et V-function dans le cas de fonctions multiplicatives sur les composantes connexes, les invariants de graphes vérifiant la formule de contraction-suppression. Tutte écrit : « En jouant avec mes W-functions j'ai obtenu un polynôme en deux variables dont on peut déduire soit le polynôme chromatique, soit le polynôme de flots, en annulant l'une ou l'autre des variables et en ajustant les signes »[8] Tutte appelle cette fonction la dichromate, parce qu'il la concevait comme une généralisation du polynôme chromatique à deux variables mais elle est maintenant appelée polynôme de Tutte. Comme dit Tutte : « Ce n'est pas correct envers Hassler Whitney qui connaissait et utilisait ces coefficients analogues sans leur accoler deux variables ». Il y a une « confusion notable »[9] entre les termes « dichromate » et « dichromatic polynomial », introduit par Tutte dans un autre article, et qui n'en diffère que légèrement. La généralisation des polynômes de Tutte aux matroïdes a été publié en premier par Crapo, même si elle figure déjà dans la thèse de Tutte[10].

Indépendamment des recherches en théorie algébrique des graphes, Potts a commencé l'étude des fonctions de partition de certains modèles de mécanique statistique en 1952. Le travail de Fortuin et Kasteleyn[11] sur les random cluster model (en), est une généralisation du modèle de Potts, et fournit une expression unifiée qui montre leur relation avec les polynômes de Tutte[10].

Spécialisations[modifier | modifier le code]

En faisant varier les valeurs de , les polynômes de Tutte prennent des formes qui ont déjà été étudiées pour elles-mêmes dans divers domaines des mathématiques ou de physique. L'un des attraits des polynômes de Tutte réside dans leur propriété à fournir un cadre unifié pour l'analyse de ces quantités.

Polynôme chromatique[modifier | modifier le code]

Le polynôme chromatique dans le plan de Tutte.

Pour , le polynôme de Tutte se spécialises en le polynôme chromatique :

est le nombre de composantes connexes de .

Pour des valeurs entières de , la valeur du polynôme chromatique est égale au nombre de colorations du graphe en utilisant couleurs. Clairement ne dépend pas de l’ensemble de couleurs. Ce qui est moins clair, c'est que c'est la valeur en d'un polynôme à coefficients entiers. Pour voir cela, on observe que :

  1. si G a n sommets et pas d'arêtes, alors .
  2. si G contient une boucle, alors .
  3. si e est une arête qui n'est pas une boucle, alors

Ces trois conditions permettent de calculer en appliquant une suite de contractions-suppressions, mais elles ne garantissent pas que des séquences différentes de suppressions et contractions conduisent au même résultat ; ceci provient du fait que compte des propriétés combinatoires, indépendamment de la récurrence. En particulier,

donne le nombre d'orientations acycliques du graphe.

Polynôme de Jones[modifier | modifier le code]

Polynômes de Jones dans le plan de Tutte.

Sur l'hyperbole , le polynôme de Tutte d'un graphe planaire se spécialise en le polynôme de Jones d'un nœud alternant associé.

Valeur en des points particuliers[modifier | modifier le code]

(2,1)[modifier | modifier le code]

compte le nombre de forêts, c'est-à-dire le nombre des sous-ensembles d'arêtes sans cycle.

(1,1)[modifier | modifier le code]

compte le nombre de forêts couvrantes (c'est le nombre d'ensembles d'arêtes sans cycle et qui ont même nombre de composants connexes que G). Si le graphe est connexe, compte le nombre d'arbres couvrants.

(1,2)[modifier | modifier le code]

compte le nombre de sous-graphes couvrants (ensemble d'arêtes qui ont même nombre de composantes connexes que G).

(2,0 )[modifier | modifier le code]

compte le nombre d'orientations acycliques (en) de G[12],[13].

(0,2)[modifier | modifier le code]

compte le nombre d'orientations fortement connexes de G[14].

(0,−2)[modifier | modifier le code]

Si G est un graphe 4-régulier, alors

compte le nombre d'orientations eulériennes de G. Ici, est le nombre de composantes connexes de G[12].

(3,3)[modifier | modifier le code]

Si G est un graphe grille de paramètres , alors compte le nombre de manière de paver un rectangle de largeur 4m et de hauteur 4n avec des tétrominos[15],[16].

Si G est un graphe planaire, alors est égal à la somme des partitions eulériennes pondérées du graphe médial de G, où le poids d'une orientation est 2 à la puissance le nombre de points selle de cette orientation (c'est-à-dire le nombre de sommets où l'orientation des arêtes incidentes est cycliquement « in, out, in out »)[17].

Modèles de Potts et d'Ising[modifier | modifier le code]

Les fonctions de partition pour le modèle d'Ising et pour les modèles de Potts à 3 et 4 états, dans le plan de Tutte..

On considère l'hyperbole définie par :

.

Sur cette hyperbole, le polynôme de Tutte se spécialise en la fonction de partition du modèle d'Ising étudiée en physique statistique. Plus précisément, le long de l'hyperbole ces deux fonctions sont liées par l'équation[18] :

Le lien avec l'hyperbole vient de ce que

pour tout nombre complexe .

Plus généralement, si on définit, pour tout entier q, l'hyperbole :

,

le polynôme de Tutte se spécialise en la fonction de partition du modèle de Potts (en) à q états. Diverses quantités physiques analysées dans le contexte du modèle de Potts se transposent en des parties spécifiques de .

Correspondances entre le modèle de Potts et le plan de Tutte[19]
modèle de Potts polynôme de Tutte
Ferromagnétisme Branche positive de
Antiferromagnétisme Branche négative de avec
Haute température Asymptote de pour
Basse temperature ferromagnétique Branche positive de asymptotique en
Température nulle antiferromagnetique q-coloration du graphe, c'est-à-dire .

Polynôme de flot[modifier | modifier le code]

Le polynôme de flot dans le plan de Tutte.

Au point , le polynôme de Tutte se spécialise en un polynôme appelé polynôme de flot et étudié en combinatoire. Pour un graphe non orienté connexe G et un entier k, un k-flot qui ne s'annule pas (« nowhere-zero flow ») est une affectations de valeurs de « flot » aux arcs d'une orientation arbitraire de G telle que le flot total entrant et le flot total sortant de chaque sommet sont égaux modulo k. Le polynôme de flot compte le nombre de ces k-flots. Cette valeur est étroitement liée au polynôme chromatique : en fait, si G est un graphe planaire, le polynôme chromatique de G est équivalent au polynôme de flot de son graphe dual par le théorème de Tutte suivant :

.

La connexion avec le polynôme de Tutte est donnée par :

.

Polynôme de fiabilité[modifier | modifier le code]

Le polynôme de fiabilité dans le plan de Tutte.

Pour , le polynôme de Tutte se spécialise en le polynôme de fiabilité entre tous points étudié en théorie des réseaux. Pour un graphe connexe G, on affecte une probabilité p à la suppression d'une arête, pour modéliser un réseau sujet à des pannes aléatoires d'arêtes. Le polynôme de fiabilité est le polynôme en p qui donne la probabilité qu'un couple de sommets de G reste connecté après des pannes sur les arêtes. Le lien avec le polynôme de Tutte est donné par :

.

Polynôme dichromatique[modifier | modifier le code]

Tutte a aussi défini une généralisation en deux variables des polynômes chromatiques voisine de ces derniers ; il les appelle polynômes dichromatiques. Pour un graphe , c'est le polynôme

est le nombre de composantes connexes du graphe partiel . Ce polynôme est relié au « corank-nullity polynomial » par

Le polynôme dichromatique ne se généralise pas aux matroïdes parce n'est pas une propriété de matroïdes : des graphes avec même matroïde peuvent avoir un nombre différent de composantes connexes.


Polynôme de Martin[modifier | modifier le code]

Le polynôme de Martin d'un graphe orienté 4-régulier a été défini par Pierre Martin in 1977[20]. Il a prouvé que si est un graphe planaire et si est son graphe médial orienté, alors

Algorithmes[modifier | modifier le code]

Suppression–contraction[modifier | modifier le code]

L'algorithme de suppression–contraction applique au graphe diamant. Les arêtes rouges sont supprimées dans l'enfant gauche, contractées dans l'enfant droit. Le polynôme résultant est la somme des monômes des feuilles, . Basé sur (Welsh et Merino 2000).

La formule de récurrence de suppression–contraction pour le polynôme de Tutte s'écrit :

,

lorsque n'est ni une boucle ni un isthme. Ceci donne immédiatement un algorithme récursif pour le calcul du polynôme : choisir une arête quelconque e et appliquer la formule jusqu'à ce que toutes les arêtes restantes sont soit des boucles, soit des isthmes. Les cas restants sont alors faciles à évaluer.

À un facteur polynomial près, le temps de calcul t de cet algorithme peut être exprimé en fonction du nombre n de sommets et du nombre m d'arêtes du graphe,

 ;

cette relation est similaire à celle des nombres de Fibonacci avec pour solution[21] :

L'analyse peut être améliorée par un facteur polynomial en le nombre d'arbre couvrants du graphe donné[22]. Pour des graphes creux, avec ce temps de calcul est en . Pour les graphes réguliers de degré k, le nombre d'arbres couvrants peut être borné par

avec ,

de sorte que l'algorithme de suppression–contraction s'exécute avec un facteur polynomial de cette borne. On a par exemple[23] :

En pratique, un test d'isomorphisme de graphes est utilisé pour éviter des appels récursifs. Cette approche fonctionne bien pour des graphes que sont creux et qui présentent beaucoup de symétries ; l'efficacité dépend alors de l'heuristique utilisé pour choisir l'arête e[22],[24],[25]. Le choix de l'arête à supprimer s'avère primordial[26].

Élimination de Gauss-Jordan[modifier | modifier le code]

Dans certains cas particuliers, le polynôme de Tutte peut être calculé en temps polynomial parce que l'élimination de Gauss-Jordan permet un calcul efficace des opérations matricielles comme le déterminant et le pfaffien. Ces algorithmes constituent eux-mêmes des résultats importants en théorie algébrique des graphes et en mécanique statistique .

est égal au nombre d'arbres couvrants d'un graphe connexe. Cette valeur peut se calculer en temps polynomial en tant que déterminant de la sous-matrice principale maximale de la matrice laplacienne de G, un résultat ancien de la théorie algébrique des graphes connu sous le nom de théorème de Kirchhoff. De façon similaire, la dimension de l'espace à deux cycles de peut être calculé par élimination de Gauss en temps polynomial.

Pour les graphes planaires, la fonction de partition du modèle d'Ising, c'est-à-dire le polynôme de Tutte sur l'hyperbole , peut être exprimée comme un pfaffian et calculé de façon efficace au moyen de l'algorithme FKT. Cette idée a été développée par Fisher, Kasteleyn et Temperley pour calculer le nombre de pavages par des dominos pour paver un modèle en treillis plan.

Méthode de Monte-Carlo[modifier | modifier le code]

En utilisant une méthode de Monte-Carlo par chaînes de Markov, le polynôme de Tutte peut être approximé d'arbitrairement près sur la branche positive de , qui est aussi la fonction de partition du modèle d'Ising ferromagnétique. Elle exploite le lien étroit entre le modèle d'Ising et le problème du comptage des couplages dans un graphe. L'idée derrière ce célèbre résultat de Jerrum et Sinclair[27] est de définir une chaîne de Markov dont les états sont les couplages du graphe donné. Les transitions sont définies en choisissant aléatoirement des arêtes et en modifiant le couplage en conséquence. La chaîne de Markov obtenue est rapidement mélangeante et conduit à des couplages « suffisamment aléatoires » qui peuvent être utilisées pour récupérer la fonction de partition par un échantillonnage aléatoire. L'algorithme obtenu est un schéma d'approximation en temps polynomial (FPRAS).

Complexité algorithmique[modifier | modifier le code]

Plusieurs problèmes algorithmiques sont associés aux polynômes de Tutte. Le plus direct est le calcul des coefficients du polynôme de Tutte pour un graphe donné.

Ceci permet l’évaluation du polynôme de Tutte en tout point ; par exemple l'évaluation de équivaut au calcul du nombre de 3-coloriages de G. Ce problème est #P-complet, même restreint à la famille des graphes planaires, de sorte que le calcul des coefficients d'un polynôme de Tutte d'un graphe est #P-difficile même pour les graphes complets.

Un ensemble de problèmes appelés de manière abrégée « Tutte » a été étudié ; il s'agit, étant donné un graphe de calculer, pour un couple de nombres complexes , la valeur . La difficulté du problème varie avec les coordonnées .

Calcul exact[modifier | modifier le code]

Le plan de Tutte. À chaque point du plan réel est associée la complexité du calcul de  :
en un point rouge, le calcul est en temps polynomial ;
sur une courbe bleue, le problème est #P-difficile en général, mais polynomial pour les graphes planaires ;
● pour tout point des régions blanches, le problème est #P-difficile même pour les graphes planaires bipartis.

Si x et y sont tous deux des entiers positifs ou nuls, le problème du calcul de est dans la classe #P. Pour des entiers relatifs, le polynôme de Tutte contient des termes négatifs, ce qui place le problème dans la classe de complexité GapP (en), qui est la fermeture par soustraction de la classe #P. Pour prendre en compte des coordonnées rationnelles, on peut définir une classe correspondant à #P[28].

La complexité algorithmique du calcul exact de se partage en deux classes selon la valeur de . Le problème est #P-difficile sauf si est un point de l'hyperbole ou est l’un des points de l'ensemble suivant (avec ) :

,

auquel cas le calcul est en temps polynomial[29]. Pour les graphes planaires, le problème devient polynomial aussi pour les points sur l'hyperbole , mais reste #P-difficile pour les autres points, même pour les graphes planaires bipartis[30]. Dans la conclusion de son article sur la dichotomie pour les graphes planaires[31], Vertigan note que le même résultat reste valable même pour la restriction supplémentaire aux graphes de degré au plus trois, à l'exception du point , qui compte les flots « nowhere-zero » pour lequel le polynôme se calcul en temps polynomial.

Ces résultats admettent divers cas particuliers intéressants. Par exemple, le problème du calcul de la fonction de partition du modèle d'Ising est #P-difficile en général, même avec l'algorithme d'Onsager et Fisher de résolution dans les treillis planaires. De même, les polynômes de Jones sont #P-difficiles à évaluer. Enfin, le calcul du nombre de 4-coloriages d'un graphe planaire est #P-complet, alors que le problème de décision est trivial par le théorème des quatre couleurs. En revanche, compter le nombre de 3-coloriages d'un graphe planaire est #P-complet parce que le problème de décision est connu pour être NP-complet.

Approximation[modifier | modifier le code]

Les points qui possèdent un bon algorithme d'approximation sont peu connus. En dehors des points pour lequel le calcul exact peut être fait en temps polynomial, le seul algorithme d'approximation connu pour est l'algorithme FPRAS de Jerrum et Sinclair, qui est valable pour les points sur l'hyperbole d'Ising pour y > 0. Si le graphe donné est restreint aux instances denses, avec degré , alors il existe un FPRAS pour x ≥ 1, y ≥ 1[32]. Même si les résultats sont moins complets que pour le calcul exact, on sait que de larges portions du plan sont difficiles à approximer[28].

Annexes[modifier | modifier le code]

Notes[modifier | modifier le code]

  1. Courtiel 2014 présente une synthèse.
  2. Bollobás 1998, chapter 10.
  3. Biggs 1993, chapter 13.
  4. Godsil et Royle 2004, chap. 15.
  5. La définition de Tutte est comme suit. On suppose les arêtes du graphe totalement ordonnées. Une arête externe (interne) d’un arbre couvrant donné est active si elle est minimale dans son cycle (cocycle) fondamental. L'activité interne (ou externe) est le nombre d'arêtes internes ou externes actives. Les ensembles d'arêtes internes ou externes actives d’un arbre couvrant dépendent de l’ordre total choisi, mais la somme de leurs nombres, sur les arbres couvrants, est indépendante de l’ordre choisi.
  6. Sokal 2005.
  7. Sokal 2005, eq. (2.26).
  8. a et b Tutte 2004.
  9. Expression de Welsh.
  10. a et b Farr 2007.
  11. Fortuin et Kasteleyn 1972.
  12. a et b Welsh 1999.
  13. Lass 2001.
  14. Las Vergnas 1980.
  15. Korn et Pak 2004.
  16. Korn et Pak 2003 donnent des interprétations combintoires du polynôme de Tutte en de nombreux autres points.
  17. Las Vergnas 1988.
  18. Welsh 1993, p. 62.
  19. Welsh et Merino 2000.
  20. Martin 1977.
  21. Wilf 1986, p. 46.
  22. a et b Sekine, Imai et Tani 1995.
  23. Chung et Yau 1999, d'après Björklund et al. 2008.
  24. Haggard, Pearce et Royle 2010.
  25. Pearce, Haggard et Royle 2010.
  26. Monagan 2012.
  27. Jerrum et Sinclair 1993.
  28. a et b Goldberg et Jerrum 2008.
  29. Jaeger, Vertigan et Welsh 1990.
  30. Vertigan et Welsh 1992.
  31. Vertigan 2005.
  32. Pour x ≥ 1 et y = 1, voir Annan 1994. Pour x ≥ 1 et y > 1, voir Alon, Frieze et Welsh 1995.

Références[modifier | modifier le code]

Livres et monographies
Articles
  • (en) Michael Monagan, « Computing Tutte Polynomials », dans Hiro-Fumi Yamada etNantel Bergeron, 24th International Conference on Formal Power Series and Algebraic Combinatorics (FPSAC 2012)), Nagoya, Japon, coll. « Discrete Mathematics and Theoretical Computer Science (DMTCS), DMTCS Proceedings », (HAL hal-01283134, lire en ligne), p. 839-850.
  • (en) David J. Pearce, Gary Haggard et Gordon Royle, « Edge-selection heuristics for computing Tutte polynomials », Chicago Journal of Theoretical Computer Science,‎ , Article 6, 14 (MR 2659710, lire en ligne).
  • (en) Kyoko Sekine, Hiroshi Imai et Seiichiro Tani, « Computing the Tutte polynomial of a graph of moderate size », dans Algorithms and computations (Cairns, 1995), Springer, coll. « Lecture Notes in Computer Science » (no 1004), (DOI 10.1007/BFb0015427, MR 1400247), p. 224–233.
  • (en) Alan D. Sokal et Bridget S. Webb, « The multivariate Tutte polynomial (alias Potts model) for graphs and matroids », dans Surveys in Combinatorics, Cambridge University Press, coll. « London Mathematical Society Lecture Note Series » (no 327), (DOI 10.1017/CBO9780511734885.009, arXiv math/0503607), p. 173–226.

Articles liés[modifier | modifier le code]

Liens externes[modifier | modifier le code]