Optimisation linéaire en nombres entiers
L'optimisation linéaire en nombres entiers (OLNE) (ou programmation linéaire en nombres entiers (PLNE) ou Modèle:Langue (IP) ou Integer Linear Programming (ILP)) est un domaine des mathématiques et de l'informatique théorique dans lequel on considère des problèmes d'optimisation d'une forme particulière. Ces problèmes sont décrits par une fonction de coût et des contraintes linéaires, et par des variables entières.
La contrainte d'intégralité sur les variables, qui différencie l'OLNE de l'optimisation linéaire classique est nécessaire pour modéliser certains problèmes, en particulier des problèmes algorithmiques. Mais cette contrainte supplémentaire rend le problème plus complexe et demande des techniques particulières.
Définitions
Définitions générales
Un problème d'optimisation est un problème mathématique où, étant donné un ensemble de variables et des contraintes sur ces variables, on doit trouver une assignation qui maximise (ou minimise) une certaine fonction de coût. On parle de problème linéaire lorsque les contraintes et la fonction de coût sont combinaisons linéaires des variables et le problème est à nombres entiers si ces variables ne sont autorisées qu'à prendre des valeurs dans l'ensemble des entiers.
La contrainte qui force les variables à prendre des valeurs entières est appelée contrainte d'intégralité. Lorsque l'on gomme cette contrainte, on parle d'un problème relaxé ou de relaxation continue, et l'on a alors affaire à un problème d'optimisation linéaire. Le rapport entre l'optimal dans la version relaxée et dans la version entière est souvent appelé integrality gap.
Formes canonique et standard
Un problème d'OLNE peut être mis sous deux formes classiques : la forme canonique et la forme standard. La forme canonique pour une maximisation est :
- <math> \begin{align}
& \text{maximiser} && \mathbf{c}^\mathrm{T} \mathbf{x}\\ & \text{tel que} && A \mathbf{x} \le \mathbf{b}, \\ & && \mathbf{x} \ge \mathbf{0}, \\ & \text{et} && \mathbf{x} \in \mathbb{Z}, \end{align} </math>,
et la forme standard est :
- <math> \begin{align}
& \text{maximiser} && \mathbf{c}^\mathrm{T} \mathbf{x}\\ & \text{tel que} && A \mathbf{x} + \mathbf{s} = \mathbf{b}, \\ & && \mathbf{s} \ge \mathbf{0}, \\ & && \mathbf{x} \ge \mathbf{0}, \\ & \text{et} && \mathbf{x} \in \mathbb{Z}, \end{align} </math>
où <math>\mathbf{c}, \mathbf{b}</math> sont des vecteurs et <math>A</math> est une matrice ayant des valeurs entières.
Exemple
On donne un exemple de problème d'OLNE, illustré par l'image ci-contre.
- <math>
\begin{align}
\max & \text{ } y \\
-x +y & \leq 1 \\
3x +2y & \leq 12 \\
2x +3y & \leq 12 \\
x,y & \ge 0 \\
x,y & \in \mathbb{Z}
\end{align}
</math>
Il y a deux variables, les solutions sont donc des couples d'entiers. Les points rouges sont les couples qui vérifient les contraintes et les pointillés rouges montrent l'enveloppe convexe de ces points. Les solutions optimales de ce problèmes sont (1,2) et (2,2).
Les lignes bleues et l'axe des abscisse délimitent les couples de réels qui satisfont toutes les contraintes sauf la contrainte d'intégralité. L'optimum est meilleur dans cette version relaxée.
Modélisation par OLNE
De nombreux problèmes de recherche opérationnelle et d'algorithmique peuvent être traduits en problème d'OLNE. Par exemple le problème de couverture par ensembles est le suivant. Étant donné un ensemble A, on dit qu'un élément e est couvert par A si e appartient à A ; pour un ensemble U et une famille S de sous-ensembles de U, le problème consiste à couvrir tous les éléments U avec une sous-famille de S la plus petite possible. Ce problème se traduit naturellement sous la forme suivante<ref>Voir par exemple Modèle:Lien web</ref> :
| minimiser : | <math>\sum_{S \in \mathcal S}x_S</math> | (Minimiser le nombre de sous-ensembles) | |
| tel que : | <math>\sum_{S\colon e \in S} x_S \geqslant 1 </math> | <math>\forall e\in \mathcal U</math> | (Tous les éléments sont couverts) |
| <math>x_S \in \{0,1\}</math> | <math>\forall S\in \mathcal S</math>. | (Chaque sous-ensemble est, ou bien dans la couverture, ou bien pas) |
Résolution
Complexité
Au sens de la théorie de la complexité, l'optimisation linéaire en nombres entiers est considérée comme difficile car c'est un problème NP-difficile<ref>Il fait d'ailleurs partie de la liste des 21 problèmes NP-complets de Karp sous le nom de 0-1 INTEGER PROGRAMMING. Voir : Modèle:Reducibility Karp 1972.</ref>. Cette complexité est facilement déduite de la NP-difficulté du problème de couverture par ensembles. Une des conséquences pratiques est que, pour des problèmes de grande taille, le temps de calcul peut être très grand.
Cependant, la complexité est polynomiale quand le nombre de variables est fixé, comme montré par Lenstra en 1983 <ref>Modèle:Article</ref>.
Cas particuliers
Lorsque les contraintes prennent la forme d'une matrice totalement unimodulaire et entière, une résolution en temps polynomial est possible car les solutions du problème relaxé sont entières.
Si le nombre de variables est fixé, alors le problème est aussi polynomial<ref> Modèle:Article.</ref>.
Méthodes de résolution exacte
Parmi les méthodes classiques de résolution, on peut citer la méthode des plans sécants (notamment à l'aide de coupes de Gomory) et le principe de séparation et évaluation (Modèle:Langue). Depuis les années 1990, l'inclusion de coupes de Gomory a permis de grandement accélérer l'algorithme séparation et évaluation. Cela a donné naissance à une nouvelle classe d'algorithmes nommés [[Branch and cut|Modèle:Langue]] <ref>Modèle:Article</ref>.
Approximation
La théorie des algorithmes d'approximation fait souvent appel à une formulation OLNE des problèmes et essaie de borner l'integrality gap pour obtenir une solution approchée en temps polynomialModèle:Refnec.
Applications
Ce type d'optimisation est très utilisée en recherche opérationnelle. C'est aussi devenu une approche classique en bio-informatique<ref>Modèle:Site web.</ref>,<ref>Modèle:Lien web.</ref>.
Notes et références
Modèle:Traduction/Référence <references />
