Cet article concerne le polynôme de Tutte d'un graphe. Pour le polynôme de Tutte d'un matroïde, voir Matroïde.
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é.
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].
Définitions et exemples
Définitions équivalentes
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 :
On peut donner une formulation légèrement différente de la même définition en posant
,
où 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
où 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
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
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 :
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
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.
Pour , le polynôme de Tutte se spécialises en le polynôme chromatique :
où 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 :
si G a n sommets et pas d'arêtes, alors .
si G contient une boucle, alors .
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.
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
(2,1)
compte le nombre de forêts, c'est-à-dire le nombre des sous-ensembles d'arêtes sans cycle.
(1,1)
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)
compte le nombre de sous-graphes couvrants (ensemble d'arêtes qui ont même nombre de composantes connexes que G).
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].
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]
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 :
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
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
où 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
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
Suppression–contraction
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,
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
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
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
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
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
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].
↑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.
Julien Courtiel, Combinatoire du polynôme de Tutte et des cartes planaires (thèse de doctorat en informatique), Université de Bordeaux I, Labri, (arXiv1411.0737, lire en ligne).
(en) N. Alon, A. Frieze et D. J. A. Welsh, « Polynomial time randomized approximation schemes for Tutte-Gröthendieck invariants: The dense case », Random Structures and Algorithms, vol. 6, no 4, , p. 459–478 (DOI10.1002/rsa.3240060409).
(en) J. D. Annan, « A Randomised Approximation Algorithm for Counting the Number of Forests in Dense Graphs », Combinatorics, Probability and Computing, vol. 3, no 3, , p. 273–283 (DOI10.1017/S0963548300001188).
(en) Jordan Awan et Olivier Bernardi, « Tutte polynomials for directed graphs », Journal of Combinatorial Theory, Series B, vol. 140, , p. 192–247 (DOI10.1016/j.jctb.2019.05.006, arXiv1610.01839).
(en) Andreas Björklund, Thore Husfeldt, Petteri Kaski et Mikko Koivisto, « Computing the Tutte polynomial in vertex-exponential time », dans Proc. of the 47th annual IEEE Symposium on Foundations of Computer Science (FOCS 2008), (ISBN978-0-7695-3436-7, DOI10.1109/FOCS.2008.40), p. 677–686.
(en) Graham E. Farr, « Tutte-Whitney polynomials: some history and generalizations », dans Geoffrey Grimmett et Colin McDiarmid, Combinatorics, complexity, and chance. A tribute to Dominic Welsh, Oxford University Press, coll. « Oxford Lecture Series in Mathematics and its Applications » (no 34), (ISBN0-19-857127-5, zbMATH1124.05020), p. 28–52.
(en) Gary Haggard, David J. Pearce et Gordon Royle, « Computing Tutte polynomials », ACM Transactions on Mathematical Software, vol. 37, no 3, , Art. 24, 17 (DOI10.1145/1824801.1824802, MR2738228).
(en) F. Jaeger, D. L. Vertigan et Dominic J. A. Welsh, « On the computational complexity of the Jones and Tutte polynomials », Mathematical Proceedings of the Cambridge Philosophical Society, vol. 108, , p. 35–53 (DOI10.1017/S0305004100068936).
Bodo Lass, « Orientations Acycliques et le Polynôme Chromatique », European Journal of Combinatorics, vol. 22, no 8, , p. 1101–1123 (ISSN0195-6698, DOI10.1006/eujc.2001.0537).
(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 », (HALhal-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 (MR2659710, 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), (DOI10.1007/BFb0015427, MR1400247), 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), (DOI10.1017/CBO9780511734885.009, arXivmath/0503607), p. 173–226.
(en) D. L. Vertigan et D. J. A. Welsh, « The Computational Complexity of the Tutte Plane: the Bipartite Case », Combinatorics, Probability and Computing, vol. 1, no 2, , p. 181–187 (DOI10.1017/S0963548300000195, lire en ligne).