\section{Question 2}
\label{sec:q2}

L'objectif de ce problème est de concevoir un horaire de travail avec minimisation de la perte de profit en utilisant le solveur Choco 4 \cite{chocoSolver}.

\subsection{Constantes}
\label{sec:q2constantes}

Plusieurs constantes sont définies par l'énoncé du problème. On définit les périodes $t=1,\ldots,16$ de la journée comme étant des périodes de 30 minutes débutant à 9h00 et se terminant à 16h59.

\begin{itemize}
\item $N$ = Nombre d'employés disponibles
\item $MIN_H$ = Nombre minimal de périodes travaillées par l'employé, incluant la pause
\item $MAX_H$ = Nombre maximal de périodes travaillées par l'employé, incluant la pause
\item $MIN_{P}$ = Nombre minimal de période dans le bloc de travail précédant et suivant la pause
\item $e_j^{req}$ = Nombre d'employés requis à la période $j \in \lbrace 1,16 \rbrace$.
\item $e_j^{souh}$ = Nombre d'employés souhaités à la période $j \in \lbrace 1,16 \rbrace$.
\item $N_P$ = Nombre de périodes dans l'horaire. Constante initialisée à 16.
\item $V_{PERTE}$ = Valeur de la perte lié à un écart d'un employé par rapport au nombre souhaité. Constante initialisée à 20.
\end{itemize}

On définit aussi ces constantes supplémentaires pour la création de contraintes, car Choco 4 ne permet pas d'effectuer d'opérations arithmétiques dans une contrainte.

\begin{itemize}
\item $MAX_P = MAX_H-MIN_P-1$ = Nombre maximal de périodes dans le bloc de travail précédant et suivant la pause.
\item $MAX_{PR} = MAX_P-MIN_H-1$ = Nombre maximal de périodes dans le bloc de repos débutant ou terminant la journée.
\item $MIN_{HT} =  MIN_H-1$ = Nombre minimal de périodes travaillées par l'employé, excluant la pause
\item $MAX_{HT} =  MAX_H-1$ = Nombre maximal de périodes travaillées par l'employé, excluant la pause
\end{itemize}

\subsection{Variables}
\label{sec:q2variables}

On définit une matrice de variables $E_{N \times N_P}$ où $e_{ij}$ correspondant à une valeur booléenne prenant la valeur $1$ si l'employé $i$ travaille à la période $j$ et $0$ sinon.

\begin{align}
  \label{eq:v1}
  &e_{ij} & \forall 1 \leq i \leq N, 1 \leq j \leq N_P\\
  &dom(e_{ij}) = \lbrace 0,1 \rbrace & \forall 1 \leq i \leq N, 1 \leq j \leq N_P\nonumber
\end{align}

Afin de résoudre le problème d'optimisation, on doit définir des variables supplémentaires, car les composantes de la fonction d'optimisation sont aussi exprimées sous forme de contraintes.

$e_{j}^{eff}$, le nombre d'employés effectif à la période $j$, et $e_{j}^{dist}$ la différence absolue entre le nombre d'employés effectif et souhaité à la période $j$.
\begin{align}
  \label{eq:v2}
  e_{j}^{eff} &= \sum_{i=1}^{N} e_{ij} & \forall 1 \leq j \leq N_P\\
  \label{eq:v3}
  e_{j}^{dist} &= \lvert e_{j}^{eff} - e_j^{souh} \rvert & \forall 1 \leq j \leq N_P
\end{align}

Enfin, nous créons une variable $N_{PERTE}$ représentant le total des unités de perte, équivalant
à la somme des différences absolues. Le nombre total d'unités de perte prend une valeur dans l'intervalle $[0,N \times N_P]$.
\begin{align}
  N_{PERTE} = \sum_{j=1}^{N_P} e_j^{dist}
\end{align}


\subsubsection{Nombre total de variables}
\label{sec:q2nbtotvar}

Le nombre total de variables pour les définitions \eqref{eq:v1}, \eqref{eq:v2},
\eqref{eq:v3} :

\begin{align}
  \label{eq:vartot}
  N \times N_P + N_P + N_P &= (N+2) \times N_p \in \theta(N)
\end{align}

\subsection{Contraintes}
\label{sec:q2contraintes}

\subsubsection{Nombre d'employés requis}
\label{sec:nbemployesreq}

La contrainte du nombre d'employés requis s'énonce comme suit: pour chaque période $i$, le nombre d'employés travaillant doit être supérieur ou égal au nombre d'employés minimum requis. Il s'agit d'une contrainte de type \textsc{Sum}.

\begin{align}
  \label{eq:q2c1}
  &\sum_{i=1}^{N} e_{ij} \geq e_j^{req} & \forall 1 \leq j \leq N_P
\end{align}

\subsubsection{Motif de la journée de travail}
\label{sec:q2automatefini}

Afin de représenter le motif correspondant à une journée de travail alternant les périodes de travail et celles de repos, on utilisera une contrainte de type \textsc{Regular} qui est paramétrée par un automate fini, c'est à dire une chaîne de caractère formant une expression régulière qui valide plusieurs contraintes définies par l'énoncé du problème.

Dans Choco 4, la contrainte \texttt{model.regular} prend en paramètre une séquence de variables entières à valider, et un automate fini défini par la classe \texttt{FiniteAutomaton}. Pour construire une instance de cette classe, on doit construire l'expression régulière en n'utilisant que des nombres comme valeurs. Comme nous avons choisi de représenter les différents états d'un employé par une variable boolenne, cette condition est satisfaite.

Voici comment nous avons construit l'expression:
\begin{itemize}
\item On débute par un bloc de repos d'au moins 0 et au plus $MAX_{PR}$ périodes.
\item On a ensuite un bloc de travail d'au moins $MIN_P$ et au plus $MAX_P$ périodes.
\item On a ensuite un bloc de repos d'une durée d'une période.
\item On a ensuite un bloc de travail d'au moins $MIN_P$ et au plus $MAX_P$ périodes.
\item On termine par un bloc de repos d'au moins 0 et au plus $MAX_{PR}$ périodes.
\end{itemize}

On obtient ainsi un automate fini défini par l'expression régulière suivante:
\begin{align}
  \label{eq:q2fa}
  {FA} = &0\lbrace 0,MAX_{PR} \rbrace\\
  &1 \lbrace MIN_P,MAX_P \rbrace\nonumber\\
  &0\nonumber\\
  &1 \lbrace MIN_P,MAX_P \rbrace\nonumber\\
  &0\lbrace 0,MAX_{PR} \rbrace\nonumber
\end{align}

L'ensemble des contraintes pour le motif d'horaire peut donc s'écrire
\begin{align}
  \label{eq:q2c2}
  &\textsc{Regular}(\lbrace e_{i1},\ldots,e_{ij} \rbrace, FA) &\forall 1 \leq i \leq N
\end{align}

\subsubsection{Nombre minimum et maximum de périodes travaillées par employé}
\label{sec:q2minmax}

On utilise une contrainte \textsc{Sum} pour appliquer le nombre de périodes minimum et maximum devant être travaillées, pour chaque employé.

\begin{align}
  \label{eq:q2c3}
   &MIN_{HT} \leq \sum_{j=1}^{N_P} e_{ij} \leq MAX_{HT} & \forall 1 \leq i \leq N
\end{align}

\subsubsection{Nombre total de contraintes}
\label{sec:q2totcontr}

Le nombre total de contraintes pour les définitions \eqref{eq:q2c1}, \eqref{eq:q2c2} et \eqref{eq:q2c3} est:
\begin{align}
  \label{eq:conttot}
  N_P + N + N &= N_P + 2N \in \theta(N)
\end{align}

\subsection{Optimisation}
\label{sec:q2optimisation}

Dans Choco, les contraintes d'optimisation sont crées de la même façon que les contraintes de satisfaction. Cependant, une des variables ne sera pas une constante, mais un intervalle dans lequel la fonction d'optimisation va effectuer sa recherche.

Nous implémenterons le nombre de travailleurs effectifs $e_j^{eff}$  à la période $j$ en utilisant une contrainte \textsc{Sum}:
\begin{align}
  &\textsc{Sum}(e_{1j},\ldots,e_{Nj},e_j^{eff}) & \forall 1 \leq j \leq N_P
\end{align}

Nous implémenterons ensuite la différence entre le nombre de travailleurs effectifs et souhaités $e_j^{dist}$ à la période $j$ en utilisant une contrainte \textsc{Distance}:
\begin{align}
  &\textsc{Distance}(e_j^{souh},e_j^{eff},e_j^{dist}) & \forall 1 \leq j \leq N_P
\end{align}

Enfin, la variable $N_{PERTE}$ est implémentée par une contrainte \textsc{Sum}:
\begin{align}
  \textsc{Sum}(e_1^{dist},\ldots,e_{N_P}^{dist},N_{PERTE}) 
\end{align}

\subsection{Paramètres}
\label{sec:q2parametres}

\subsection{Resultats}
\label{sec:q2resultats}


%%% Local Variables:
%%% mode: latex
%%% TeX-master: "rapport"
%%% End: