Optimisation linéaire en nombres entiers — Wikipédia

Pour les articles homonymes, voir ILP.

L'optimisation linéaire en nombres entiers (OLNE) (ou programmation linéaire en nombres entiers (PLNE) ou integer programming (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.

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.

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 :

${\displaystyle {\begin{aligned}&{\text{maximiser}}&&\mathbf {c} ^{\mathrm {T} }\mathbf {x} \\&{\text{tel que}}&&A\mathbf {x} \leq \mathbf {b} ,\\&&&\mathbf {x} \geq \mathbf {0} ,\\&{\text{et}}&&\mathbf {x} \in \mathbb {Z} ,\end{aligned}}}$,

et la forme standard est :

${\displaystyle {\begin{aligned}&{\text{maximiser}}&&\mathbf {c} ^{\mathrm {T} }\mathbf {x} \\&{\text{tel que}}&&A\mathbf {x} +\mathbf {s} =\mathbf {b} ,\\&&&\mathbf {s} \geq \mathbf {0} ,\\&&&\mathbf {x} \geq \mathbf {0} ,\\&{\text{et}}&&\mathbf {x} \in \mathbb {Z} ,\end{aligned}}}$

où ${\displaystyle \mathbf {c} ,\mathbf {b} }$ sont des vecteurs et ${\displaystyle A}$ est une matrice ayant des valeurs entières.

On donne un exemple de problème d'OLNE, illustré par l'image ci-contre.

${\displaystyle {\begin{aligned}\max &{\text{ }}y\\-x+y&\leq 1\\3x+2y&\leq 12\\2x+3y&\leq 12\\x,y&\geq 0\\x,y&\in \mathbb {Z} \end{aligned}}}$

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ème sont (1,2) et (2,2).

Les lignes bleues et les axes des abscisses et des ordonnées 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.

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^[1] :

minimiser :	${\displaystyle \sum _{S\in {\mathcal {S}}}x_{S}}$		(Minimiser le nombre de sous-ensembles)
tel que :	${\displaystyle \sum _{S\colon e\in S}x_{S}\geqslant 1}$	${\displaystyle \forall e\in {\mathcal {U}}}$	(Tous les éléments sont couverts)
	${\displaystyle x_{S}\in \{0,1\}}$	${\displaystyle \forall S\in {\mathcal {S}}}$.	(Chaque sous-ensemble est, ou bien dans la couverture, ou bien pas)

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^[2]. 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 ^[3].

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^[4].

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 (branch and bound). 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 ^[5].

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 polynomial^{[réf. nécessaire]}.

Ce type d'optimisation est très utilisée en recherche opérationnelle. C'est aussi devenu une approche classique en bio-informatique^[6],[7].

• (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Integer programming » (voir la liste des auteurs).

• (en) Bradley, Hax, and Magnanti, « Integer programming », dans Applied Mathematical Programming, Addison-Wesley, 1977

• Leo Liberti et Ruslan Sadykov, « Présentation : Introduction à la Programmation Linéaire en Nombres Entiers », sur École polytechnique, 2006

v · m Optimisation: théorie et algorithmes
Non linéaire	Méthode du nombre d'or Recherche linéaire Méthode de Nelder-Mead Critères de Wolfe Méthode de Broyden-Fletcher-Goldfarb-Shanno Algorithme à régions de confiance Pénalisation Algorithme du gradient Algorithme du gradient stochastique Méthode de Newton Algorithme de Gauss-Newton Algorithme de Levenberg-Marquardt Algorithme du lagrangien augmenté
Convexe	Optimisation complètement positive Optimisation copositive Optimisation SDP Méthode des plans sécants Algorithme de Frank-Wolfe Méthode de l'ellipsoïde Linéaire Optimisation conique Algorithme du simplexe Méthodes de points intérieurs Décomposition de Benders Génération de colonnes Quadratique Optimisation quadratique successive
Combinatoire	Algorithme d'approximation Programmation dynamique Algorithme glouton Optimisation linéaire en nombres entiers
Métaheuristique	Algorithme évolutionniste Algorithme génétique Essaims Forêts aléatoires Boosting
Liste

• Portail des mathématiques
• Portail de l'informatique théorique