MP ÉCOLE POLYTECHNIQUE FILIÈRE OPTION INFORMATIQUE CONCOURS D'ADMISSION 2008 COMPOSITION D'INFORMATIQUE (Durée : 4 heures) L'utilisation des calculatrices n'est pas autorisée pour cette épreuve. Le langage de programmation choisi par le candidat doit être spécifié en tête de la copie. Structure de corde La majorité des langages de programmation fournissent une notion primitive de chaînes de caractères. Si ces chaînes s'avèrent adaptées à la manipulation de mots ou de textes relativement courts, elles deviennent généralement inutilisables sur de très grands textes. L'une des raisons de cette inefficacité est la duplication d'un trop grand nombre de caractères lors des opérations de concaténation ou d'extraction d'une sous-chaîne. Or il existe des domaines où la manipulation efficace de grandes chaînes de caractères est essentielle (représentation du génôme en bio-informatique, éditeurs de texte, etc.). Ce problème aborde une alternative à la notion usuelle de chaîne de caractères connue sous le nom de corde. Une corde est tout simplement un arbre binaire, dont les feuilles sont des chaînes de caractères usuelles et dont les noeuds internes représentent des concaténations. Ainsi la corde représente le mot GATTACCATAGATTAC, obtenu par concaténation des cinq mots GATTAC, CAT, AG, ATT et AC. L'intérêt des cordes est d'offrir une concaténation immédiate et un partage possible de caractères entre plusieurs chaînes, au prix d'un accès aux caractères un peu plus coûteux. La partie I traite de fonctions préliminaires sur les mots. La partie II définit les principales opérations sur les cordes. Enfin, les parties III et IV étudient le problème de l'équilibrage des cordes, selon deux algorithmes différents dont le second, l'algorithme de Garsia-Wachs, est optimal. Les parties peuvent être traitées indépendamment. La partie IV utilise les notations de la partie III. Partie I. Préliminaires sur les mots On considère les mots dans construits sur un alphabet . Le mot x = a0 a1 · · · an-1 de longueur n est représenté dans ce problème par la liste des entiers codant ses caractères a0 , a1 , . . . an-1 . Ainsi le mot ATT est représenté par une liste h65, 84, 84i. Pour ne pas dépendre du codage des caractères, on identifie les caractères et leurs codes. Le type des mots est défini par : 1 (* Caml *) type mot == int list ; ; { Pascal } type mot = ^cellule ; cellule = record lettre : integer ; suite : mot ; end ; En Pascal la liste vide est nil et on pourra utiliser la fonction suivante pour construire des mots : function nouveauMot(a:integer; x:mot) : mot; var r : mot; begin new(r); r^.lettre := a; r^.suite := x; nouveauMot := r; end; Question 1 Écrire la fonction longueurMot qui calcule la longueur d'un mot. (* Caml *) longueurMot : mot -> int { Pascal } function longueurMot(x : mot) : integer Question 2 Écrire la fonction iemeCar qui prend en arguments un entier i et un mot x = a0 a1 . . . an-1 , et qui renvoie le caractère ai . On supposera 0 i < n. (* Caml *) iemeCar : int -> mot -> int { Pascal } function iemeCar(i : integer ; x : mot) : integer Question 3 Écrire la fonction prefixe qui prend en arguments un entier k et un mot x = a0 a1 . . . an-1 , et qui renvoie le mot a0 a1 . . . ak-1 c'est-à-dire le mot constitué des k premiers caractères de x. On supposera 0 k n. (* Caml *) prefixe : int -> mot -> mot { Pascal } function prefixe(k : integer ; x : mot) : mot Question 4 Écrire la fonction suffixe qui prend en arguments un entier k et un mot x = a0 a1 . . . an-1 , et qui renvoie le mot ak ak+1 . . . an-1 c'est-à-dire le mot obtenu en supprimant les k premiers caractères de x. On supposera 0 k n. (* Caml *) suffixe : int -> mot -> mot { Pascal } function suffixe(k : integer ; x : mot) : mot Partie II. Opérations sur les cordes Comme expliqué dans l'introduction, une corde est un arbre binaire dont les feuilles sont des mots. Plus précisément, une corde est soit vide, soit constituée d'un unique mot (une feuille), soit un noeud constitué de deux cordes et représentant leur concaténation. Pour des raisons d'efficacité, on conserve dans les feuilles aussi bien que dans les noeuds la longueur de la corde correspondante. On définit donc le type corde suivant : (* Caml *) type corde = | Vide | Feuille of int*mot | Noeud of int*corde*corde ; ; { Pascal } type corde = ^arbre ; nature = (Feuille, Noeud) ; arbre = record case indicateur : nature of Feuille : (n :integer ; x :mot) ; Noeud : (n :integer ; g, d :corde) ; end ; 2 En Pascal la corde vide est représentée par nil. Dans la suite, on garantira l'invariant suivant sur les cordes : dans une corde de la forme Feuille(n, x), on a n = longueurMot(x) et n > 0 ; dans une corde de la forme Noeud(n, c1 , c2 ), on a c1 6= Vide, c2 6= Vide et n est la longueur totale de la corde, c'est-à-dire la somme des longueurs de c1 et c2 . On notera en particulier que la corde de longueur 0 est nécessairement représentée par Vide. Question 5 Écrire la fonction longueur qui renvoie la longueur d'une corde. (* Caml *) longueur : corde -> int { Pascal } function longueur(c : corde) : integer Question 6 Écrire la fonction nouvelleCorde qui construit une corde à partir d'un mot. (* Caml *) nouvelleCorde : mot -> corde { Pascal } function nouvelleCorde(m : mot) : corde Question 7 Écrire la fonction concat qui construit la concaténation de deux cordes. (* Caml *) concat : corde -> corde -> corde { Pascal } function concat(c1 : corde ; c2 : corde) : corde Question 8 Écrire la fonction caractere qui prend en arguments un entier i et une corde c représentant le mot a0 a1 · · · an-1 , et qui renvoie le caractère ai . On supposera 0 i < n. (* Caml *) caractere : int -> corde -> int { Pascal } function caractere(i : integer ; c : corde) : integer Question 9 Écrire la fonction sousCorde qui prend en arguments un entier i, un entier m et une corde c représentant le mot a0 a1 · · · an-1 , et qui renvoie une corde représentant le mot ai ai+1 · · · ai+m-1 c'est-à-dire la sous-corde de c débutant au caractère i et de longueur m. On supposera 0 i < i + m n. (* Caml *) sousCorde : int -> int -> corde -> corde { Pascal } function sousCorde(i : integer ; m : integer ; c : corde) : corde On s'attachera à réutiliser dans la corde résultat autant de sous-arbres de la corde c que possible. Partie III. Équilibrage Le hasard des concaténations peut amener une corde à se retrouver déséquilibrée, c'est-à-dire à avoir certaines de ses feuilles très éloignées de la racine et donc d'accès plus coûteux. Le but de cette partie est d'étudier une stratégie de rééquilibrage a posteriori. Considérons une corde c composée de k + 1 feuilles, et donc de k noeuds internes. Notons ces k + 1 feuilles x0 , x1 , . . . xk lorsqu'on les considère de la gauche vers la droite, si bien que c représente 3 le mot x0 x1 . . . xk . La profondeur de la feuille xi dans c est notée prof(xi ) et est définie comme la distance de xi à la racine de c. Voici un exemple de corde pour k = 5 où la profondeur de chaque feuille est indiquée entre parenthèses : (1) (3) (3) (4) (3) (4) Le coût de l'accès à un caractère de la feuille xi est défini comme la profondeur de cette feuille dans c, soit prof(xi ) (on ne considère donc pas le coût de l'accès dans le mot xi lui-même). Le coût total d'une corde est alors défini comme la somme des coûts d'accès à tous ses caractères, et vaut donc Coût(c) = k X longueurMot(xi ) × prof(xi ) i=0 Un rééquilibrage consiste à construire une corde différente, dont les feuilles sont x0 , x1 . . . xk dans le même ordre (le mot représenté ne doit pas changer) et dont le coût est, éventuellement, meilleur. L'algorithme proposé est le suivant. On considère un tableau file de cordes dans lequel les feuilles de c vont être successivement insérées dans le sens des indices croissants. Les cases d'indices 0 et 1 ne sont pas utilisées. La case d'indice i contient soit la corde vide (Vide), soit une corde de hauteur inférieure ou égale à i - 2 et dont la longueur est comprise dans l'intervalle [Fi , Fi+1 [ où Fi désigne le i-ème terme de la suite de Fibonacci. La hauteur d'une corde c, notée hauteur(c), est la profondeur maximale de ses feuilles, c'est-à-dire : hauteur(Vide) = 0 hauteur(Feuille(n, x)) = 0 hauteur(Noeud(n, c1 , c2 )) = 1 + max{hauteur(c1 ), hauteur(c2 )} Pour équilibrer la corde c dont les feuilles sont les mots x0 , x1 , . . . xk , dans cet ordre, on procède ainsi : 1. On insère successivement chaque feuille xj dans le tableau file à partir de la case 2. L'insertion d'une feuille, et plus généralement d'une corde, à partir de la case d'indice i se fait ainsi : (a) La corde à insérer est concaténée à droite de la corde se trouvant dans la case i ; soit c le résultat. Si la longueur de c est comprise dans l'intervalle [Fi , Fi+1 [, on affecte c à la case i et on a terminé l'insertion de cette corde. (b) Sinon, on affecte Vide à la case i et on retourne à l'étape (a) pour effectuer l'insertion de c à partir de la case d'indice i + 1. On garantit l'invariant suivant : après l'insertion de la feuille xj , la concaténation successive de toutes les cordes contenues dans les cases de file, considérées dans le sens des indices décroissants, est égale à une corde représentant le mot x0 x1 . . . xj . 2. Le résultat est alors la corde résultant de la concaténation successive de toutes les cordes de file, considérées dans le sens des indices décroissants. 4 Question 10 Calculer le résultat de cet algorithme sur la corde de l'exemple précédent. Question 11 On rappelle que la suite de Fibonacci (Fn ) est définie par F0 = 0, F1 = 1, F n+2 = Fn+1 + Fn pour n 0. Afin d'éviter tout débordement arithmétique en calculant Fn , on limite la taille de file à 44 cases indexées de 0 à 43 (les cases 0 et 1 n'étant pas utilisées). On introduit la constante tailleMax = 44 et on calcule les valeurs de Fn pour 0 n tailleMax une fois pour toutes dans un tableau fib ainsi déclaré : { Pascal } (* Caml *) let tailleMax = 44 ; ; const tailleMax = 44 ; let fib = make_vect (tailleMax+1) 0 ; ; var fib :array[0..tailleMax] of integer ; Écrire la fonction initialiserFib qui initialise le tableau fib. Question 12 Le tableau file utilisé par l'algorithme est déclaré comme un tableau global contenant des cordes : { Pascal } (* Caml *) let file = make_vect tailleMax Vide ; ; var file :array[0..tailleMax-1] of corde ; Écrire la fonction inserer qui prend en arguments une corde c et un entier i, tels que 2 i < tailleMax, hauteur(c) i - 2 et longueur(c) Fi , et réalise l'insertion de c dans le tableau file à partir de la case d'indice i. (* Caml *) inserer : corde -> int -> unit { Pascal } procedure inserer(c : corde ; i : integer) Question 13 Montrer que l'invariant hauteur(ci ) i - 2 et longueur(ci ) Fi est préservé par cette fonction pour toutes les valeurs ci non vides des cases du tableau file (2 i < tailleMax). Question 14 Écrire la fonction equilibrer qui réalise l'équilibrage d'une corde par l'algorithme ci-dessus. (* Caml *) equilibrer : corde -> corde { Pascal } function equilibrer(c : corde) : corde Question 15 Soit c une corde non vide renvoyée par la fonction equilibrer ci-dessus. Soit n sa longueur et h sa hauteur. Montrer que l'on a n Fh+1 En déduire qu'il existe une constante K (indépendante de n) telle que Coût(c) n (log (n) + K) où est le nombre d'or (1 + 5)/2 et log désigne le logarithme à base . On admettra que l'on a Fi+1 i / 5 pour tout i 0. 5 Partie IV. Équilibrage optimal Bien que satisfaisant en pratique, l'équilibrage étudié dans la partie précédente n'est pas optimal. Le but de cette partie est d'étudier une stratégie optimale de rééquilibrage (algorithme de GarsiaWachs). Les notations sont celles de la partie III. L'algorithme proposé procède en deux temps : il commence par construire une corde de coût minimal, ayant les mêmes feuilles que c mais pas nécessairement dans le même ordre ; puis, dans un deuxième temps, il transforme cette corde en une autre de même coût où les feuilles sont maintenant dans le même ordre que dans c. La première partie de l'algorithme opère sur une liste de cordes q = hq0 , q1 , . . . qm i et procède de la manière suivante : 1. Initialement, la liste q est la liste hx0 , x1 , . . . xk i des k + 1 feuilles de c. 2. Tant que la liste q contient au moins deux éléments, on effectue l'opération suivante : (a) Déterminer le plus petit indice i tel que longueur(qi-1 ) longueur(qi+1 ) le cas échéant, et poser i = m sinon. (b) Ôter qi-1 et qi de la liste q et former leur concaténation ; soit c la corde obtenue. (c) Déterminer le plus grand indice j < i tel que longueur(qj-1 ) longueur(c ) le cas échéant, et poser j = 0 sinon. (d) Insérer c dans la liste q juste après qj-1 (et donc au début de la liste q si j = 0). 3. Le résultat est l'unique élément restant dans la liste q. Il est clair que le résultat de cet algorithme est une corde ayant les mêmes feuilles que c mais que ces feuilles ne sont pas nécessairement dans le bon ordre. On admettra le résultat suivant : l'arbre obtenu est de coût minimal. Question 16 Pour simplifier le codage, on suppose que le nombre k de feuilles est inférieur à une certaine valeur (ici maxf = 1000) et que la liste q est représentée dans un tableau global q : (* Caml *) let maxf = 1000 ; ; let q = make_vect maxf Vide ; ; { Pascal } const maxf :integer = 1000 ; var q :array[0..maxf-1] of corde ; Écrire la fonction initialiserQ qui prend en argument une corde c, remplit les k + 1 premiers éléments de q avec les feuilles x0 , . . . , xk de c (c'est-à-dire des cordes de la forme Feuille), et renvoie la valeur de k. On supposera c 6= Vide. (* Caml *) initialiserQ : corde -> int { Pascal } function initialiserQ(c : corde) : integer On admettra avoir effectué le reste de l'algorithme ci-dessus et avoir donc écrit une fonction phase1 qui prend en argument une corde c et renvoie la corde obtenue par l'algorithme ci-dessus. (* Caml *) phase1 : corde -> corde { Pascal } function phase1(c : corde) : corde 6 La deuxième étape de l'algorithme procède ainsi. Soit c1 la corde obtenue à l'issue de la première étape de l'algorithme. Chaque feuille xi de c se trouve dans c1 à une certaine profondeur ; notons pi cette profondeur. On admet la propriété suivante : il existe une corde c2 dont les feuilles sont exactement x0 , x1 , . . . xk dans cet ordre et où la profondeur de chaque xi est exactement pi . On peut alors construire c2 en ne connaissant que les pi , et c2 est dès lors un rééquilibrage optimal de c. Question 17 Les profondeurs pi seront stockées dans un tableau global prof : (* Caml *) let prof = make_vect maxf 0 ; ; { Pascal } var prof :array[0..maxf-1] of integer ; Écrire la fonction initialiserProf qui prend en argument les cordes c et c1 et range dans le tableau prof, à l'indice i, la profondeur de la feuille xi (de c) dans c1 pour 0 i k. On pourra avantageusement réutiliser le tableau q et la fonction initialiserQ. (* Caml *) initialiserProf : corde -> corde -> unit { Pascal } procedure initialiserProf(c : corde ; c1 : corde) Indication : on admettra que, pour comparer les feuilles de c et c1 , on peut utiliser l'égalité fournie par le langage, i.e. le symbole =. Question 18 Écrire la fonction reconstruire qui construit c2 à partir de la seule donnée des tableaux q et prof. Attention : pour des profondeurs pi quelconques, il n'existe pas nécessairement de corde où chaque xi a la profondeur pi . On demande ici un algorithme qui fonctionne uniquement sous l'hypothèse qu'une telle corde existe (ce qui est le cas ici). (* Caml *) reconstruire : unit -> corde { Pascal } function reconstruire : corde Question 19 Combiner les fonctions ci-dessus pour obtenir une fonction de rééquilibrage optimal. (* Caml *) equilibrerOpt : corde -> corde { Pascal } function equilibrerOpt(c : corde) : corde 7