Voir la source de Équilibres et jeux matriciels

En [[théorie des jeux]], '''un équilibre de Nash''' définit une situation dans laquelle aucun joueur ne peut modifier seul sa stratégie sans affaiblir sa position personnelle<ref group="Note">Défini en détail plus loin.</ref>. La notion d'équilibre a pris un essor considérable au {{s-|XX|e}} et se trouve désormais au centre de la théorie des jeux. Depuis 1928, lorsque [[John von Neumann]] initie la théorie des jeux avec son [[Théorème_du_minimax_de_von_Neumann|théorème du minimax]], les équilibres deviennent un outil essentiel d'analyse dans des domaines variés tels que l'économie, les sciences sociales et les sciences politiques. En 1958, [[John Forbes Nash]] décrit un théorème qui assure que tout jeu à nombre de joueurs et stratégies finis possède un [[Équilibre_de_Nash|équilibre]] pour les stratégies mixtes. La théorie des jeux se consacre, entre autres, à calculer si ces états peuvent être atteints, efficacement en pratique, avec des agents rationnels. En 1994, [[Christos Papadimitriou]] crée la [[PPAD (complexité)|classe de complexité PPAD]] dans laquelle sont repris les calculs des équilibres de Nash pour les jeux à plusieurs joueurs.

== Notion d'équilibre == 
Dans la [[théorie des jeux]], l'équilibre ou l''''équilibre de Nash''', nommé d'après [[John Forbes Nash]], est un [[concept de solution]] dans lequel l'ensemble des choix faits par plusieurs joueurs, connaissant leurs stratégies réciproques, est devenu stable du fait qu'aucun ne peut modifier seul sa stratégie sans affaiblir sa position personnelle.

Cela signifie que n'importe quel jeu pourrait, en principe, atteindre un état de calme où personne n'aurait l'intention de changer son comportement. Cette situation soulève immédiatement plusieurs questions : est-ce possible dans la réalité de notre monde économique ou politique ? Est-il possible de calculer cet équilibre efficacement ?
Il est important de résoudre cette deuxième question. Car ce principe d'équilibre est au-delà de tout outil conceptuel, c'est une prédiction de comportements stratégiques rationnels par des agents en situation de conflit, un contexte pourtant dénué de logique et de calculs. La crédibilité de cette solution envers les économistes dépend directement de l'efficacité avec laquelle la théorie des jeux peut la proposer.

== Matrices de gain == 
La [[matrice des gains]], matrice de coût ou encore matrice des paiements (anglicismes issus de ''payoff matrix'') est un moyen de représenter un jeu sous forme stratégique par un tableau indiquant les gains (ou paiements) associés à chaque action en fonction des actions de l'autre joueur.

Typiquement, une matrice de gain pour '''les jeux à somme nulle''' se représente comme suit :

{| class="wikitable" border="1" cellspacing="2" style="margin: 1em auto 1em auto; text-align:center; width:25%"
|-
|
| colspan="3" |j2
|-
| rowspan="3" |j1
| 
| stratégie 1
| stratégie 2
|-
| stratégie 1
| -2
| 3
|-
| stratégie 2
| 3
| -4
|}

Lorsque le joueur 1 choisit sa stratégie 2 et le joueur 2 choisit sa stratégie 1 : le joueur 1 gagne 3 et le joueur 2 perd 3.

Lorsque le joueur 1 choisit sa stratégie 1 et le joueur 2 choisit sa stratégie 1 : le joueur 1 perd 2 et le joueur 2 gagne 2.

Une matrice de gain pour '''les jeux à somme non nulle''' se représente comme suit :

{| class="wikitable" border="1" cellspacing="2" style="margin: 1em auto 1em auto; text-align:center; width:25%"
|-
|
| colspan="3" |j2
|-
| rowspan="3" |j1
| 
| stratégie 1
| stratégie 2
|-
| stratégie 1
| (5 ; 3)
| (5 ; -10)
|-
| stratégie 2
| (-10 ; 5)
| (-5 ; -5)
|}


Lorsque le joueur 1 choisit sa stratégie 2 et le joueur 2 choisit sa stratégie 1 : le joueur 1 perd 10 et le joueur 2 gagne 5.

Nous prendrons soin de toujours placer le joueur 1 à gauche et le joueur 2 en haut.

== Équilibre dans les jeux deux joueurs à somme nulle ==
=== Pour les stratégies pures ===
Soit un joueur qui dispose d'un ensemble de choix fini <math>X</math>; une stratégie pure est un choix <math>x \in X </math>.

==== Point-selle du Minimax ====
{{Théorème|nom=Proposition<ref>{{harvsp|id=BINMORE2007|Binmore|2007|p=220-221}}</ref>|Une matrice'' <math>(n,k)</math> de réels <math>A=(a_{ij})</math> ''possède un point d'équilibre (ou point-selle) si et seulement si''

<center><math>\max_{1\leq j\leq k}(\min_{1\leq i\leq n} a_{ij})=\min_{1\leq i\leq n}(\max_{1\leq j\leq k} a_{ij})</math></center>

''En tout point d'équilibre de la matrice, l'entrée correspondante du tableau est égale à ce minimax''.|style=background-color:#ffeedd;border-color:#ff8822;}}

'''Jeu 1'''
<hr/>

Soit ''k'' le nombre de stratégies pures disponibles pour le joueur 1 et ''n'' le nombre de stratégies pures disponibles pour le joueur 2. On numérote de 1 à ''k'' les stratégies à la disposition du J1 et de 1 à ''n'' celles du J2. Prenons un exemple facile de jeu où ''n'' = ''k'' = 3.

Voici la matrice de gain :

{| class="wikitable" border="1" cellspacing="2" style="margin: 1em auto 1em auto; text-align:center"
|- 
| J1/J2
| joue sa stratégie 1
| joue sa stratégie 2
| joue sa stratégie 3
|-
| joue  sa stratégie 1
| -1
| 1,5
| 0
|-
| joue sa stratégie 2
| 2
| 2
| 1
|-
| joue sa stratégie 3
| 2
| -3
| 0
|}

Le joueur 1 se dit : « Si je joue 1, je risque de perdre 1 unité, et si je joue 3 d'en perdre 3; mais si je joue 2, je gagne une unité dans le pire des cas. »

Le joueur 2 se dit : « Si je joue 1, je risque de perdre 2 unités, et si je joue 2 d'en perdre 2; mais si je joue 3 mes pertes sont limitées à une unité dans le pire des cas.


Le joueur 1, reconstituant le raisonnement du joueur 2 à choisir sa stratégie 3 pour limiter ses pertes, est conforté dans son choix.

Le joueur 2, reconstituant le raisonnement du joueur 1 à choisir sa stratégie 2, est conforté dans son choix (il perdrait 2 s'il choisissait autrement).

Aucun des joueurs n'est prêt à changer sa stratégie car ce sont les choix les plus rationnels. 

Par exemple, si le joueur 1 devient gourmand et change sa stratégie pour gagner 1,5 ou 3 unités, sachant que le joueur 2 joue intelligemment et donc sa stratégie 3, le joueur 1 ne gagnerait plus que 0 unité, au lieu de 1. Il n'a donc aucun intérêt à changer sa stratégie. Ce raisonnement est réciproque pour le joueur 2.

En effet, selon la proposition, nous trouvons bien un équilibre :

<math>\max_{1\leq j\leq k}(\min_{1\leq i\leq n} a_{ij})=\max\{-1,1,-3\}=1</math> tandis que <math>\min_{1\leq i\leq n}(\max_{1\leq j\leq k} a_{ij})=\min\{2,1.5, 1\}=1.</math>

Le minimax et le maximin sont égaux; la ligne où le minimax est atteint et la colonne où le maximin est atteint indiquent un point d'équilibre.


'''Jeu 2'''
<hr/>

{| class="wikitable" border="1" cellspacing="2" style="margin: 1em auto 1em auto; text-align:center"
|-
| J1/J2
| joue sa stratégie 1
| joue sa stratégie 2
| joue sa stratégie 3
|-
| joue sa stratégie 1
| 0
| -1,03
| 1,02
|-
| joue sa stratégie 2
| 1,05
| 0
| -1,01
|-
| joue sa stratégie 3
| -1,07
| 1,04
| 0
|}

Similairement au jeu 1, le joueur 1 va choisir sa stratégie 2 pour minimiser les pertes, ici 1,01 unité. Le joueur 2 choisira sa stratégie 3 où il perdra au pire 1,02 unité. 

Le joueur 1 anticipe la décision du joueur 2 et se rend compte qu'il va bel et bien perdre 1,01 unité et choisira alors sa stratégie 1 pour gagner 1,02 unité. Le joueur 2 va lui aussi anticiper cette décision et choisir sa stratégie 2 pour gagner 1,03 unité, etc. Aucune stratégie pure de la part des deux joueurs ne parviendra à s'imposer.

Toujours par rapport à la proposition, nous pouvons remarquer qu'il n'y a aucun équilibre : 

<math>\max_{1\leq j\leq k}(\min_{1\leq i\leq n} a_{ij})=-1,01</math> tandis que <math>\min_{1\leq i\leq n}(\max_{1\leq j\leq k} a_{ij})=1,02</math>

=== Pour les stratégies mixtes ===
Soit un joueur qui dispose d'un ensemble de choix fini <math>X</math>; une stratégie mixte est l'attribution d'une certaine distribution de probabilité sur l'ensemble des stratégies pures <math>x \in X</math>. Bien entendu, chaque stratégie pure <math>x \in X</math> est une stratégie mixte qui choisit <math>x</math> avec une probabilité de 1.

==== Théorème du Minimax ====
En 1928,  [[John von Neumann]] démontre son [[Théorème_du_minimax_de_von_Neumann|théorème du minimax]]. Il assure que, pour un jeu non coopératif opposant deux joueurs, à nombre fini de stratégies pures et à somme nulle, il existe au moins une situation d'interaction stable, à savoir une situation dans laquelle aucun des deux joueurs n'a intérêt à changer sa stratégie mixte si l'autre ne la change pas.
Ce théorème est un cas particulier du [[Équilibre de Nash#Théorème|théorème fondamental de la théorie des jeux à n joueurs]] de [[John Forbes Nash]], démontré en 1950.

À la  différence du cas des points-selle vu précédemment, il existe toujours un équilibre dans un jeu à somme nulle et à deux joueurs aux stratégies mixtes.

'''Explications :'''

Soit la stratégie mixte du joueur 1 un vecteur colonne <math>X =(p_1\,\ldots\,p_m)^T</math> et le joueur 2 choisit la colonne j, alors le gain moyen du J1 est <math>\sum_{1\leq i\leq m}p_ia_{ij}</math>. Similairement, soit la stratégie mixte du J2 représentée par le vecteur colonne <math>Y = (q_1\,\ldots\,q_n)^T</math> et le J1 choisit la ligne i, alors le gain moyen du J1 est <math>\sum_{1\leq j\leq n}a_{ij}q_j</math>.
De manière plus générale, si J1 utilise sa stratégie mixte <math>X</math> et J2 utilise sa stratégie mixte <math>Y</math>, le gain moyen pour le J1 est de : 

<center><math>Y^TAX = \sum_{1\leq i\leq m\atop 1\leq j\leq n}p_ia_{ij}qj</math></center>

Désormais, un couple de stratégies mixtes <math>(X_0,Y_0)\in\Delta_k\times\Delta_n</math> est un équilibre de Nash pour <math>A</math> lorsque :

*Pour tout <math>Y\in\Delta_n</math>, <math>Y_0^TAX_0\leq Y^TAX_0</math> (si J2 change de stratégie tandis que J1 conserve la stratégie <math>X_0</math>, J2 ne peut qu'y perdre)
*Pour tout <math>X\in\Delta_k</math>, <math>Y_0^TAX\leq Y_0^TAX_0</math>  (si J1 change de stratégie tandis que J2 conserve la stratégie <math>Y_0</math>, J1 ne peut qu'y perdre).

{{Théorème|nom=Proposition<ref>{{harvsp|id=BINMORE2007|Binmore|2007|p=228}}</ref>|''Soit'' <math>A=(a_{ij})</math> ''une matrice'' <math>(n,k)</math> ''de réels. Il y a un équilibre de Nash pour'' <math>A</math> ''(au sens défini ci-dessus) si et seulement si''

<center><math>\max_{X\in\Delta_k}\min_{Y\in\Delta_n}Y^TAX=\min_{Y\in\Delta_n}\max_{X\in\Delta_k}Y^TAX.</math></center>

''En tout équilibre de Nash'' <math>(X_0,Y_0)</math> ''pour'' <math>A</math>, ''la valeur de'' <math>Y_0^TAX_0</math> ''est alors égale à ce minimax''.|style=background-color:#ffeedd;border-color:#ff8822;}}

'''Exemple :'''

Soit la matrice de gain <math>A</math> :

{| class="wikitable" border="1" cellspacing="2" style="margin: 1em auto 1em auto; text-align:center; width:25%"
|-
|
| colspan="2" |j2
|-
| rowspan="2" |j1
| -2
| 3
|-
| 3
| -4
|}

Analysons premièrement ce jeu du point de vue du joueur 1. Supposons qu'il joue sa stratégie "une" 3/5 du temps et "deux" 2/5 du temps de manière aléatoire (<math>X</math> vaut ici <math>(3/5\,\,\,2/5)^T</math>). Dans ce cas :

*Si J2 choisit sa stratégie 1, J1 perd 2 3/5 du temps et gagne 3 2/5 du temps. En moyenne, il gagne -2(3/5) + 3(2/5) = 0
*Si J2 choisit sa stratégie 2, J1 gagne 3 3/5 du temps et perd 4 2/5 du temps. En moyenne, il gagne 3(3/5) - 4(2/5) = 1/5

Sur la durée, le J1 est assuré d'être rentable. Mais peut-il améliorer cette stratégie pour qu'il ait un gain positif quoi que J2 choisisse ?

Définissons <math>p</math> comme la quantité de fois que J1 choisit sa stratégie 1. Commençons par trouver p tel que J1 gagne le même montant qu'il choisisse sa stratégie 1 ou sa stratégie 2 :

<center><math>-2p + 3(1 - p) = 3p - 4(1 - p)</math>

<math>3 - 5p = 7p - 4</math>

<math>12p = 7</math>

<math>p = \frac7{12}</math>
</center>

Donc <math>X</math> = <math>(7/12\,\,\,5/12)^T</math>. En moyenne, J1 gagne −2(7/12) + 3(5/12) = 1/12 quoi que J2 choisisse.

Peut-il faire mieux ? Pas si J2 joue intelligemment. En fait, J2 peut utiliser la même procédure et il se retrouvera avec <math>Y</math> = <math>(7/12\,\,\,5/12)^T</math> qui lui garantit une perte moyenne d'au plus 1/12.

On dit alors que la valeur du jeu <math>v</math> vaut <math>\frac1{12}</math>.

==== Résolution des matrices de gain 2x2 ====
Soit <math>A=\begin{pmatrix}
a & b\\
d & c\\
\end{pmatrix}</math>

Nous supposons qu'il n'y a aucun point-selle.

Si <math>a \ge b</math> alors <math>b<c</math> sinon b serait un point-selle. Si <math>b<c</math>, alors <math>c>d</math> sinon c serait un point-selle.

Donc, <math>d<a</math> et <math>a>b</math>.

En résumé, si <math>a>b</math>, <math>a>b<c>d<a</math>.

De manière symétrique, si <math>a \le b</math>, alors <math>a<b>c<d>a</math>.

{{Théorème|nom=Proposition<ref name="fergusonGT">Thomas S. Ferguson ([[Université de Californie à Los Angeles]]), [http://www.math.ucla.edu/~tom/Game_Theory/mat.pdf Game Theory], {{p.|10,11}}.</ref>|S'il n'y a aucun point-selle, alors soit <math>a>b</math>, <math>b<c</math> ; <math>c>d</math> et <math>d<a</math> ; soit <math>a<b</math>, <math>b>c</math> ; <math>c<d </math> et <math>d>a</math>.|style=background-color:#ffeedd;border-color:#ff8822;}}


Comme vu précédemment, reprenons la formule du gain moyen pour le j1 avec <math>X</math> = <math>(p\,\,\,1-p)^T</math> lorsque le j2 utilise sa colonne 1 et 2 :

<center><math>ap+d(1-p) = bp+c(1-p)</math></center>

Lorsqu'on résout <math>p</math>, on trouve : 

<center><math> p = \frac{c-d}{(a-b)+(c-d)} </math></center>

Vu qu'il n'y a aucun point-selle, <math>(a-b)</math> et <math>(c-d)</math> sont soit tous deux strictement positifs soit tous deux strictement négatifs, donc <math>0<p<1</math>. Le gain moyen de j1 pour utiliser sa stratégie vaut :

<center><math>v=ap+d(1-p) = \frac{ac-bd}{a-b+c-d}</math></center>

Le calcul de la perte moyenne de j2 est similaire et nous obtenons :

<center><math>v=aq+b(1-q)= \frac{ac-bd}{a-b+c-d}</math></center>

Cette similarité prouve que le jeu a une valeur et que les joueurs ont leur stratégie optimale (propriété du théorème minimax pour les jeux à stratégies finies).

'''Exemple 1 : '''

{| cellpadding="2" style="background-color:white;margin: 1em auto 1em auto;"
|-
| <math>A=\begin{pmatrix}
-2 & 3\\
3 & -4\\
\end{pmatrix}</math>&nbsp;&nbsp; || <math>p = \frac{-4-3}{-2-3-4-3} = 7/12</math><br/>
<math>q = p</math><br/>
<math>v = \frac{-8-9}{-2-3-4-3} = 1/12</math>
|}
'''Exemple 2 : '''
{| cellpadding="2" style="background-color:white;margin: 1em auto 1em auto;"
|-
| <math>A=\begin{pmatrix}
0 & -10\\
1 & 2\\
\end{pmatrix}</math>&nbsp;&nbsp; || <math>p = \frac{2-1}{0+10+2-1} = 1/11</math><br/>
<math>q = \frac{2+10}{0+10+2-1} = 12/11</math>
|}

Dans ce dernier exemple, q devrait être entre 0 et 1. Il ne l'est pas car il existe un point-selle. Dans la matrice A, "1" est un point-selle donc p vaut 0, q vaut 1 et v vaut 1.

==== Résolution par suppression de la stratégie dominante ====
Parfois, de grandes matrices peuvent être réduites en taille jusqu'à, dans le meilleur des cas, une matrice 2x2. Pour y parvenir, il faut supprimer les lignes et les colonnes qui sont mauvaises (inutiles) pour les deux joueurs. 

{{Théorème|nom=Définition<ref name ="fergusonGT" />|La i-ième ligne d'une matrice <math> A = (a_{ij})</math> domine la k-ième ligne si <math>a_{ij} \geq a_{kj}</math> pour tout j.  La i-ième ligne de A domine strictement la k-ième ligne si <math>a_{ij}>a_{kj}</math> pour tout j. Similairement, la j-ième colonne de A domine (domine strictement) la k-ième colonne si <math>a_{ij} \leq a_{ik}</math> (<math>a_{ij} < a_{ik}</math>) pour tout i.|style=background-color:#ffeedd;border-color:#ff8822;}}

'''Tout ce qui peut être atteint avec une ligne ou une colonne dominante peut être atteint au moins aussi bien avec la ligne ou la colonne dominée.'''
Donc, la suppression d'une colonne ou d'une ligne dominante ne change pas la valeur du jeu. Cependant, il existe peut-être une stratégie optimale qui utilisait cette ligne et cette colonne dominante et disparaîtra avec sa suppression, mais il restera toujours au moins une stratégie optimale autre part.

Dans le cas d'une suppression d'une colonne ou d'une ligne strictement dominante, l'ensemble des stratégies optimales reste inchangé.

'''Exemple : '''

Nous allons itérer la procédure successivement sur les lignes et colonnes sur la matrice A.

{| cellpadding="2" style="background-color:white;"
|-
| <math>A=\begin{pmatrix}
2 & 0 & 4\\
1 & 2 & 3\\
4 & 1 & 2\\
\end{pmatrix}</math>&nbsp;&nbsp; || La dernière colonne est strictement dominée celle du milieu.
|}
<br/>
{| cellpadding="2" style="background-color:white;"
|-
| <math>A=\begin{pmatrix}
2 & 0\\
1 & 2\\
4 & 1\\
\end{pmatrix}</math>&nbsp;&nbsp; || La ligne une est dominée par la ligne trois.
|}
<br/>
{| cellpadding="2" style="background-color:white;"
|-
| <math>A=\begin{pmatrix}
1 & 2\\
4 & 1\\
\end{pmatrix}</math>&nbsp;&nbsp; || Cette matrice n'a pas de point-selle donc :<br/>
<math>p=3/4</math><br/>
<math>q=1/4</math><br/>
<math>v=7/4</math><br/>
<math>X</math> = <math>(0\,\,\,3/4\,\,\,1/4)^T</math> et <math>Y</math> = <math>(1/4\,\,\,3/4\,\,\,0)^T</math>
|}

==== Généralisation aux matrices NxN ====
'''Connaissances supposées''' : la [[Optimisation_linéaire|programmation linéaire]].

De manière informelle, nous voulons trouver des variables qui sont des distributions de probabilité sur différents choix <math>p_1\,\ldots\,p_n</math>. Ces probabilités doivent suivre deux inégalités linéaires : <math>p_1+\ldots+p_n = 1</math> et <math>p_i\geq0</math> pour tout <math>i</math>.

Notre but est de maximiser le pire cas (le minimum) de notre gain, sur toutes les colonnes que l'adversaire peut jouer. Plus formellement, on assignera une variable <math>v</math> qui représente le minimum et on donnera comme contrainte que notre gain soit au moins <math>v</math> pour chaque colonne et que notre objectif soit de maximiser <math>v</math>. En résumé nous avons :

'''Variables : ''' <math>p_1\,\ldots,p_n</math> et <math>v</math>.

'''Objectifs : ''' Maximiser <math>v</math>.

'''Contraintes :'''
* <math>p_i\geq0</math> pour tout <math>i</math> et <math>\sum_{i}p_i=1</math>
* pour toutes les colonnes <math>j</math>, nous avons <math>\sum_{i}{p_ia_{ij}}\geq v</math>

Ce problème est un programme linéaire et peut donc être résolu. 

===== Méthode du simplexe =====
'''Connaissances supposées''' : [[Algorithme_du_simplexe|l'algorithme du simplexe]].

L'algorithme le plus courant est la méthode du simplexe, introduite par [[George Dantzig]] à partir de 1947. La plupart des exemples sont résolus en temps polynomial mais ce n'est pas garanti. Certains cas peuvent dépasser cette complexité. Nous poursuivrons avec un exemple :

Déclinons <math>v</math> pour les deux joueurs : <math>\lambda</math> que j1 cherche à maximiser et <math>\mu</math> que j2 cherche à minimiser.
Nous avons donc :
* <math>\lambda - \sum_{1\leq i\leq n}{x_ia_{ij}}\leq 0</math> pour le j1
* <math>\mu - \sum_{1\leq j\leq m}{y_ja_{ij}}\geq 0</math> pour le j2

On remarque que le problème du j2 est le "dual" du problème du j1. Ces deux problèmes peuvent donc être résolus ensemble grâce à la [[Algorithme_du_simplexe#L.27algorithme_du_simplexe_dual|méthode du simplexe dual]].

Notez que si <math>A</math> est la matrice de gain, c'est la négation de la transposée de <math>A</math> qui est placée dans le tableau du simplexe :

{| "border="0" cellspacing="0" style="text-align:center;margin: 1em auto 1em auto;"
|-
|style="border-bottom:1px solid black;"|  ||style="border-left:1px solid black;border-bottom:1px solid black;"| <math>x_1</math> ||style="border-bottom:1px solid black;"| <math>x_2</math> || style="border-bottom:1px solid black;"|<math>\ldots</math> || style="border-bottom:1px solid black;"|<math>x_n</math> || style="border-bottom:1px solid black;"|<math>\lambda</math> || style="border-left:1px solid black;border-bottom:1px solid black;"|
|-
| <math>y_1</math> || style="border-left:1px solid black;" | <math>-a_{11}</math> || <math>-a_{21}</math> || <math>\ldots</math> || <math>-a_{n1}</math> || <math>1</math> || style="border-left:1px solid black;"|<math>0</math>
|-
| <math>\vdots</math> || style="border-left:1px solid black;"|<math>\vdots</math>||  ||  ||  || <math>\vdots</math> || style="border-left:1px solid black;"|<math>\vdots</math>
|-
| <math>y_m</math> || style="border-left:1px solid black;"|<math>-a_{1m}</math> || <math>-a_{2m}</math> || <math>\ldots</math> || <math>-a_{nm}</math> || <math>1</math> || style="border-left:1px solid black;"|<math>0</math>
|-
|<math>\mu</math> || style="border-left:1px solid black;"|<math>1</math> || <math>1</math> || <math>\ldots</math> || <math>1</math> || <math>0</math> || style="border-left:1px solid black;"|<math>1</math>
|-
|style="border-top:1px solid black;"| || style="border-left:1px solid black;border-top:1px solid black;"|<math>0</math> || style="border-top:1px solid black;" |<math>0</math> || style="border-top:1px solid black;" |<math>\ldots</math> || style="border-top:1px solid black;" |<math>0</math> || style="border-top:1px solid black;" |<math>-1</math> || style="border-left:1px solid black;border-top:1px solid black;"|&nbsp;&nbsp;&nbsp;<math>0</math>&nbsp;&nbsp;&nbsp;
|}


Soit la matrice :
<center>
<math>A=\begin{pmatrix}
1 & 4 & 0 & -1\\
2 & -1 & 3 & 1\\
-1 & 1 & 0 & 3\\
\end{pmatrix}</math>
</center>

Le tableau devient :

{| "border="0" cellspacing="0" style="text-align:center;margin: 1em auto 1em auto;"
|-
|style="border-bottom:1px solid black;"|  ||style="border-left:1px solid black;border-bottom:1px solid black;"| <math>x_1</math> ||style="border-bottom:1px solid black;"| <math>x_2</math> || style="border-bottom:1px solid black;"|<math>x_3</math> || style="border-bottom:1px solid black;"|<math>\lambda</math> || style="border-left:1px solid black;border-bottom:1px solid black;"|
|-
| <math>y_1</math> || style="border-left:1px solid black;" | <math>-1</math> || <math>-2</math> || <math>1</math> || <math>1</math> ||  style="border-left:1px solid black;"|<math>0</math>
|-
| <math>y_2</math> || style="border-left:1px solid black;"|<math>-4</math> || <math>1</math> || <math>-1</math> || <math>1</math> || style="border-left:1px solid black;"|<math>0</math>
|-
| <math>y_3</math> || style="border-left:1px solid black;"|<math>0</math> || <math>-3</math> || <math>0</math> || '''<span style="font-size:22px;">&#9312;</span>''' || style="border-left:1px solid black;"|<math>0</math>
|-
|<math>y_4</math> || style="border-left:1px solid black;"|<math>1</math> || <math>-1</math> || <math>-3</math> || <math>1</math> ||  style="border-left:1px solid black;"|<math>0</math>
|-
|<math>\mu</math> || style="border-left:1px solid black;"|<math>1</math> || <math>1</math> || <math>1</math> || <math>0</math> ||  style="border-left:1px solid black;"|<math>1</math>
|-
|style="border-top:1px solid black;"| || style="border-left:1px solid black;border-top:1px solid black;"|<math>0</math> || style="border-top:1px solid black;" |<math>0</math> || style="border-top:1px solid black;" |<math>0</math> || style="border-top:1px solid black;" |<math>-1</math> || style="border-left:1px solid black;border-top:1px solid black;"|<math>0</math>
|-
| style="padding: 0px 20px;"|  ||  style="padding: 0px 20px;"| ||  style="padding: 0px 20px;"| || style="padding: 0px 20px;"|  || style="padding: 0px 20px;"| || style="padding: 0px 20px;"|
|}

Nous désirerions pivoter <math>\lambda</math> et <math>\mu</math>. Malheureusement, comme toujours pour les matrices de gain pour les jeux, la valeur du pivot vaut 0. Donc on doit pivoter en deux fois : premièrement pour déplacer <math>\mu</math> en haut et ensuite pour déplacer <math>\lambda</math> à gauche. On inter-change alors <math>y_3</math> et <math>\lambda</math> et supprimons la ligne <math>\lambda</math>, puis nous inter-changeons <math>x_3</math> et <math>\mu</math> et supprimons la colonne <math>\mu</math>.



{| "border="0" cellspacing="0" style="text-align:center;margin: 1em auto 1em auto;"
|-
| valign="top"|

{| "border="0" cellspacing="0" style="text-align:center;"
|-
|style="border-bottom:1px solid black;"|  ||style="border-left:1px solid black;border-bottom:1px solid black;"| <math>x_1</math> ||style="border-bottom:1px solid black;"| <math>x_2</math> || style="border-bottom:1px solid black;"|<math>x_3</math> || style="border-bottom:1px solid black;"|<math>y_3</math> || style="border-left:1px solid black;border-bottom:1px solid black;"|
|-
| <math>y_1</math> || style="border-left:1px solid black;" | <math>-1</math> || <math>1</math> || <math>1</math> || <math>-1</math> ||  style="border-left:1px solid black;"|<math>0</math>
|-
| <math>y_2</math> || style="border-left:1px solid black;"|<math>-4</math> || <math>4</math> || <math>-1</math> || <math>1</math> || style="border-left:1px solid black;"|<math>0</math>
|-
| <math>y_4</math> || style="border-left:1px solid black;"|<math>1</math> || <math>2</math> || <math>-3</math> || <math>-1</math> || style="border-left:1px solid black;"|<math>0</math>
|-
|<math>\mu</math> || style="border-left:1px solid black;"|<math>1</math> || <math>1</math> || '''<span style="font-size:22px;">&#9312;</span>''' || <math>0</math> ||  style="border-left:1px solid black;"|<math>1</math>
|-
|style="border-top:1px solid black;"| || style="border-left:1px solid black;border-top:1px solid black;"|<math>0</math> || style="border-top:1px solid black;" |<math>-3</math> || style="border-top:1px solid black;" |<math>0</math> || style="border-top:1px solid black;" |<math>1</math> || style="border-left:1px solid black;border-top:1px solid black;"|<math>0</math>
|-
| style="padding: 0px 20px;"|  ||  style="padding: 0px 20px;"| ||  style="padding: 0px 20px;"| || style="padding: 0px 20px;"|  || style="padding: 0px 20px;"| || style="padding: 0px 20px;"|
|}
|| 
&nbsp;&nbsp;&nbsp;<math>\longrightarrow</math>&nbsp;&nbsp;&nbsp;
| valign="top" |
{| "border="0" cellspacing="0" style="text-align:center;"
|-
|style="border-bottom:1px solid black;"|  ||style="border-left:1px solid black;border-bottom:1px solid black;"| <math>x_1</math> ||style="border-bottom:1px solid black;"| <math>x_2</math> || style="border-bottom:1px solid black;"|<math>y_3</math> || style="border-left:1px solid black;border-bottom:1px solid black;"|
|-
| <math>y_1</math> || style="border-left:1px solid black;" | <math>-2</math> || <math>0</math> || <math>-1</math> || style="border-left:1px solid black;"|<math>-1</math>
|-
| <math>y_2</math> || style="border-left:1px solid black;"|<math>-3</math> || <math>5</math> || <math>-1</math> || style="border-left:1px solid black;"|<math>1</math>
|-
| <math>y_4</math> || style="border-left:1px solid black;"|<math>4</math> || <math>5</math> || <math>-1</math> || style="border-left:1px solid black;"|<math>3</math>
|-
|<math>x_3</math> || style="border-left:1px solid black;"|<math>1</math> || <math>1</math> || <math>0</math> ||  style="border-left:1px solid black;"|<math>1</math>
|-
|style="border-top:1px solid black;"| || style="border-left:1px solid black;border-top:1px solid black;"|<math>0</math> || style="border-top:1px solid black;" |<math>-3</math> ||style="border-top:1px solid black;" |<math>1</math> || style="border-left:1px solid black;border-top:1px solid black;"|<math>0</math>
|-
| style="padding: 0px 20px;"|  ||  style="padding: 0px 20px;"| ||  style="padding: 0px 20px;"| || style="padding: 0px 20px;"| || style="padding: 0px 20px;"|
|}
|}


Maintenant, nous utilisons la règle des pivots pour la méthode du simplexe jusqu'à obtenir, par exemple :

{| "border="0" cellspacing="0" style="text-align:center;margin: 1em auto 1em auto;"
|-
|style="border-bottom:1px solid black;"|  ||style="border-left:1px solid black;border-bottom:1px solid black;"| <math>y_4</math> ||style="border-bottom:1px solid black;"| <math>y_2</math> || style="border-bottom:1px solid black;"|<math>y_1</math> || style="border-left:1px solid black;border-bottom:1px solid black;"|
|-
| <math>y_3</math> || style="border-left:1px solid black;" |  ||  || || style="border-left:1px solid black;"|<math>3/7</math>
|-
| <math>x_2</math> || style="border-left:1px solid black;"| ||  ||  || style="border-left:1px solid black;"|<math>16/35</math>
|-
| <math>x_1</math> || style="border-left:1px solid black;"| ||  ||  || style="border-left:1px solid black;"|<math>2/7</math>
|-
|<math>x_3</math> || style="border-left:1px solid black;"| || ||  ||  style="border-left:1px solid black;"|<math>9/35</math>
|-
|style="border-top:1px solid black;"| || style="border-left:1px solid black;border-top:1px solid black;"|<math>13/35</math> || style="border-top:1px solid black;" |<math>8/35</math> ||style="border-top:1px solid black;" |<math>14/35</math> || style="border-left:1px solid black;border-top:1px solid black;"|<math>33/35</math>
|-
| style="padding: 0px 30px;"|  ||  style="padding: 0px 30px;"| ||  style="padding: 0px 30px;"| || style="padding: 0px 30px;"| || style="padding: 0px 30px;"|
|}


Donc, <math>v = 33/35</math>. L'unique stratégie optimale pour le J1 vaut <math>(2/7\,\,\,16/35\,\,\,9/35)^T</math>. L'unique stratégie optimale pour le J2 vaut <math>(14/35\,\,\,8/35\,\,\,0\,\,\,13/35)^T</math>.

===== Méthode de l'ellipsoïde =====
Cette méthode a été inventée par [[Leonid Khatchian]] dans les années 1980 en Russie. L'algorithme ne permet de résoudre que la "faisabilité" du problème mais peut se poursuivre avec une recherche binaire sur la fonction objectif pour résoudre le problème d'optimisation. 

Le principe est de commencer avec une ellipse (appelée ellipsoïde pour les grandes dimensions) qui contient à coup sûr la région de "faisabilité". Ensuite, on essaye le centre de l'ellipse et on observe si des contraintes sont enfreintes. Dans la négative, on a terminé. Si oui, on regarde les contraintes enfreintes. Désormais, nous savons que la solution (si elle existe) se trouve dans au plus la moitié restante de l'ellipse. On trouve alors une nouvelle ellipse plus petite qui contient cette moitié et on répète l'opération.

On estime qu'on diminue à chaque étape le volume de l'ellipse d'au moins un facteur de <math>(1-1/n)</math> par rapport à l'ellipse originale. Donc, pour  chaque n-étape, le volume diminue d'un facteur de <math>1/e</math>. 

C'était le premier algorithme qui garantissait la résolution d'un programme linéaire en temps polynomial mais reste assez lent en pratique.

===== Algorithme de Karmarkar =====
L'algorithme de Karmarka est un algorithme introduit par [[Narendra Karmarkar]] en 1984 et décrit une complexité polynomiale beaucoup plus rapide que la méthode de l'ellipsoïde en pratique. 

L’idée de base de la méthode est d’utiliser des fonctions barrières pour décrire l’ensemble des solutions qui est convexe par définition du problème. À l’opposé de l’algorithme du simplexe, cette méthode atteint l’optimum du problème en passant par l’intérieur de l’ensemble des solutions réalisables. Il utilise la [[Méthodes_de_points_intérieurs|méthode des points intérieurs]].

== Équilibre pour les jeux à n joueurs ==
Prenons le célèbre exemple du [[dilemme du prisonnier]] pour introduire les jeux à somme non nulle. Nous prenons encore un jeu à 2 joueurs par simple souci de clarté dans les matrices de gains. Les concepts de ce chapitre s'appliquent pour n joueurs.

'''Dilemme du prisonnier'''
<hr/>
Deux prisonniers sont accusés d'un crime et chacun est confronté au choix de trahir (T) ou de rester silencieux (S). Leur choix va définir leur sentence qui se traduit par des années passées en prison. Voici la matrice de coût associée à ces choix que les prisonniers veulent minimiser :

{| class="wikitable" border="1" cellspacing="2" text-align="left"
|-
| J1\J2
| (T)
| (S)
|-
| (T)
| (4 ; 4)
| (1 ; 5)
|-
| (S)
| (5 ; 1)
| (2 ; 2)
|}

La seule solution stable est celle où les deux prisonniers trahissent. Dans les trois autres possibilités, un prisonnier peut passer du silence à trahir et améliorer son gain. Les prisonniers n'oseront donc jamais jouer le payement 2 et 1 par peur de récolter la sentence maximale.

=== Solution par stratégie dominante ===
Dans le jeu du prisonnier, chaque joueur détient une stratégie optimale unique, indépendante de la stratégie des autres. Un jeu qui a cette propriété a une solution par stratégie dominante.

{{Théorème|nom=Proposition <ref name="nissanGT">N. Nisan, T. Roughgarden, E. Tardos, V. V. Vazirani (2007) ([[Université de Cambridge]]), [http://www.cambridge.org/journals/nisan/downloads/Nisan_Non-printable.pdf  Algorithmic Game Theory], {{p.|10,11,29-53}}.</ref> | Soit <br/>
- un vecteur de stratégie <math>s \in \mathbf{S}</math> où <math>s_i</math> est la stratégie du joueur <math>i</math> <br/>
- <math>s_{-i}</math> pour décrire le vecteur de dimension <math>(n - 1)</math> des stratégies des autres joueurs<br/>
- <math>u_i(s)</math> une valeur représentant l'issue de la situation, typiquement le gain ou le coût.<br/>

Un vecteur de stratégie <math>s \in \mathbf{S}</math> est ''' solution par stratégie dominante ''' si, pour chaque joueur <math>i</math>, et chaque vecteur de stratégie alternatif <math>s'\in \mathbf{S}</math>, nous avons :

<math>u_i(s_i,s'_{-i}) \geqslant u_i(s'_i,s'_{-i})</math>.|style=background-color:#ffeedd;border-color:#ff8822;}}

Il est important de noter qu'une solution à stratégie dominante ne donne pas forcément le meilleur gain. Nous l'avons vu, par exemple, dans le jeu du prisonnier.

Par ailleurs, ce cas particulier où tous les joueurs  ont leur propre stratégie dominante est assez rare. Peu de jeux satisfont cette propriété. La cause est que l'évaluation de la valeur <math>u_i(s)</math> peut être très complexe, jusqu'à devenir exponentielle.

=== Pour les stratégies pures ===
Parce que des situations de stratégies dominantes sont rares, il est nécessaire de chercher un concept de solutions plus largement applicable. Une solution acceptable selon la théorie des jeux est l'idée qu'aucun joueur, poursuivant toujours ses intérêts, ne peut améliorer individuellement sa situation sans s'écarter de sa stratégie. C'est l'idée que capture l'équilibre de Nash.

{{Théorème|nom=Proposition <ref name="nissanGT" /> | Soit un vecteur de stratégie <math>s \in \mathbf{S}</math> est ''' un équilibre de Nash ''' si, pour chaque joueur <math>i</math>, et chaque stratégie alternative <math>s'_i\in \mathbf{S_i}</math> nous avons : 

<math>u_i(s_i,s_{-i}) \geqslant u_i(s'_i,s_{-i})</math>|style=background-color:#ffeedd;border-color:#ff8822;}}

En d'autres termes, aucun joueur ne peut changer sa stratégie <math>s_i</math> en <math>s'_i</math> et améliorer son payement, en supposant que chacun des autres joueurs garde sa stratégie choisie dans s. 

Nous remarquons également qu'une solution par stratégie dominante est un équilibre de Nash.

=== Pour les stratégies mixtes ===
Désormais, les joueurs ne jouent plus leur stratégie de manière déterministe. Le problème se pose alors de comprendre comment les joueurs évaluent leur payement d'un choix probabiliste. Un joueur préférait-il un choix qui rapporte peu avec une forte probabilité, un choix qui rapporte beaucoup avec une faible probabilité ou encore un choix qui risque de faire perdre peu avec une grande probabilité ? Dans le cas des équilibres de Nash pour les stratégies mixtes, nous supposons qu'ils sont neutres par rapport au risque et qu'ils essayent uniquement de maximiser leur gain.

Pour définir ce choix de stratégie aléatoire formellement, on dit que chaque joueur peut définir une distribution de probabilité sur sa palette de choix, et que chacun choisira de manière indépendante un choix selon cette distribution de probabilité.

En 1951, Nash prouve que, sous ces conditions, tout jeu à nombre de joueurs et de stratégies mixtes fini possède un équilibre de Nash.

=== La classe PPAD ===
{{loupe|PPAD (complexité)}}
Calculer un équilibre de Nash peut s'avérer, pour certains jeux, extrêmement long. Traditionnellement, les problèmes se retrouvent dans deux catégories de complexité : la classe [[P (complexité)|P]] pour les problèmes qui ont un algorithme en temps polynomial et la classe [[Problème_NP-complet|NP-complet]]. Cependant, le concept de NP-complet ne peut être appliqué pour les rares problèmes "dont chaque instance possède une solution", comme les équilibres pour les jeux à nombre de joueurs et de stratégies mixtes fini. 

On démontre alors que trouver un équilibre de Nash est complet pour la [[PPAD (complexité)|classe de problème PPAD]] grâce au [[Théorème du point fixe de Brouwer]]. Chaque problème de cette classe partage cette même propriété que chacune de leurs instances possède une solution. PPAD est donc une sous-classe de NP qui contient les problèmes durs, les équilibres de Nash y compris<ref name="simplified">[[Constantinos Daskalakis]], Paul W. Goldberg, Christos H. Papadimitriou ([[Massachusetts Institute of Technology]]), [http://people.csail.mit.edu/costis/simplified.pdf The Complexity of Computing a Nash Equilibrium], {{p.|3 et suivantes }}.</ref>.

Plusieurs algorithmes de résolution ont été proposés depuis les années 1970, mais la plupart ont des complexités inconnues. C'est pourquoi désormais, les jeux sont séparés en "classes" de jeux pour étudier les différents cas et les différentes complexités. Principalement, les recherches espèrent établir l'existence d'un algorithme sous-exponentiel  - <math>O(n^{log n/ \epsilon^2})</math> - où n est la taille du jeu.

=== Équilibres de Nash approximés ===
Pour les jeux à deux joueurs, les mesures de l'équilibre de Nash sont des nombres rationnels (la distribution de leur stratégie mixte) et sa solution est donc claire. Mais, selon Nash, lorsque le nombre de joueurs augmente, il ne peut y avoir que des solutions irrationnelles. Les calculs deviennent alors numériquement imprécis, c'est pourquoi on introduit le concept d'équilibre de Nash approximé.

Si un équilibre de Nash est une situation dans laquelle personne ne souhaite dévier de sa stratégie, un ϵ-équilibre de Nash est un ensemble de stratégies mixtes dans laquelle les joueurs peuvent espérer augmenter leurs gains d'au plus ϵ en changeant de stratégie. Un ϵ-équilibre de Nash n'est pertinent que si ϵ est petit.

La question se pose de savoir s'il existe un [[schéma d'approximation en temps polynomial]] pour les équilibres de Nash approximés, c'est-à-dire un algorithme qui tourne en temps polynomial en rapport à la taille du jeu, en acceptant des dépendances sur temps sur 1/ϵ.

== Annexes ==
=== Notes et références ===
; Notes
{{Références|colonnes=2|groupe="Note"}}
; Références
{{Références}}

=== Liens externes ===
* {{en}} {{pdf}} T. S. Ferguson [http://www.math.ucla.edu/~tom/LP.pdf Linear Programming A Concise Introduction].
* {{en}} {{pdf}} T. S. Ferguson [http://www.math.ucla.edu/~tom/LP.pdf Game Theory Part II. Two-Person Zero-Sum Games].
* {{en}} {{pdf}} A. Blum [http://www.cs.cmu.edu/~avrim/451f11/lectures/lect1101.pdf Linear Programming]

=== Bibliographie ===
* {{Ouvrage
 | langue        = en
 | auteur1       = [[Kenneth Binmore|Ken Binmore]]
 | titre         = Playing for Real
 | sous-titre    = A Text on Game Theory
 | éditeur       = Oxford University Press US
 | année         = 2007
 | pages totales = 639
 | isbn          = 978-0-19-530057-4
 | lire en ligne = https://books.google.com/books?id=h5QRDAAAQBAJ&printsec=frontcover
 | id            = BINMORE2007
}}{{Commentaire biblio|Un cours d'initiation à la théorie des jeux et à son formalisme mathématique, par un des spécialistes mondiaux du sujet.}}
* {{en}} N. Nisan, T. Roughgarden, E. Tardos, V. V. Vazirani (2007). ''Algorithmic Game Theory''

{{Portail|économie|mathématiques|sociologie}}

[[Catégorie:Théorie des jeux]]