Introduction à l'algorithmique et ses bases
Introduction à l'algorithmique et ses bases
[Link]
I NTRODUCTION ÀÙ L ’A LGORITHMIQUE
« Un langage de programmation est une convention pour donner des ordres à un
ordinateur. Ce n’est pas censé être obscur, bizarre et plein de pièges subtils. Ça, ce sont les caractéristiques de la
magie. » - Dave Small
« C'est illogique, Capitaine » - Mr Spock
L’algorithmique est un terme d’origine arabe, comme algèbre, amiral ou zénith. Ce n’est pas une excuse pour
massacrer son orthographe, ou sa prononciation.
Ainsi, l’algo n’est pas « rythmique », à la différence du bon rock’n roll. L’algo n’est pas non plus « l’agglo ».
Alors, ne confondez pas l’algorithmique avec l’agglo rythmique, qui consiste à poser des parpaings en cadence.
1. Qu’est-ce que l’algomachin ?
Avez-vous déjà ouvert un livre de recettes de cuisine ? Avez-vous déjà déchiffré un mode
d’emploi traduit directement du coréen pour faire fonctionner un magnétoscope ou un répondeur téléphonique
réticent ? Si oui, sans le savoir, vous avez déjà exécuté des algorithmes.
Plus fort : avez-vous déjà indiqué un chemin à un touriste égaré ? Avez-vous fait chercher un
objet à quelqu’un par téléphone ? Ecrit une lettre anonyme stipulant comment procéder à une remise de rançon ? Si
oui, vous avez déjà fabriqué – et fait exécuter – des algorithmes.
Comme quoi, l’algorithmique n’est pas un savoir ésotérique réservé à quelques rares initiés
touchés par la grâce divine, mais une aptitude partagée par la totalité de l’humanité. Donc, pas d’excuses…
Un algorithme, c’est une suite d’instructions, qui une fois exécutée
correctement, conduit à un résultat donné. Si l’algorithme est juste, le résultat est le résultat voulu, et le
touriste se retrouve là où il voulait aller. Si l’algorithme est faux, le résultat est, disons, aléatoire, et
décidément, cette saloperie de répondeur ne veut rien savoir.
Complétons toutefois cette définition. Après tout, en effet, si l’algorithme,
comme on vient de le dire, n’est qu’une suite d’instructions menant celui qui l’exécute à résoudre un
problème, pourquoi ne pas donner comme instruction unique : « résous le problème », et laisser
l’interlocuteur se débrouiller avec ça ? A ce tarif, n’importe qui serait champion d’algorithmique sans faire
aucun effort. Pas de ça Lisette, ce serait trop facile.
Le malheur (ou le bonheur, tout dépend du point de vue) est que
justement, si le touriste vous demande son chemin, c’est qu’il ne le connaît pas. Donc, si on n’est pas un
goujat intégral, il ne sert à rien de lui dire de le trouver tout seul. De même les modes d’emploi
contiennent généralement (mais pas toujours) un peu plus d’informations que « débrouillez-vous pour que
ça marche ».
Pour fonctionner, un algorithme doit donc contenir uniquement des
instructions compréhensibles par celui qui devra l’exécuter. C’est d’ailleurs l’un des points délicats
pour les rédacteurs de modes d’emploi : les références culturelles, ou lexicales, des utilisateurs, étant
variables, un même mode d’emploi peut être très clair pour certains et parfaitement abscons pour d’autres.
En informatique, heureusement, il n’y a pas ce problème : les choses auxquelles on doit donner des
instructions sont les ordinateurs, et ceux-ci ont le bon goût d’être tous strictement aussi idiots les uns que
les autres.
il faut avoir une certaine intuition, car aucune recette ne permet de savoir a priori quelles
instructions permettront d’obtenir le résultat voulu. C’est là, si l’on y tient, qu’intervient la forme
« d’intelligence » requise pour l’algorithmique. Alors, c’est certain, il y a des gens qui possèdent au départ
davantage cette intuition que les autres. Cependant, et j’insiste sur ce point, les réflexes, cela s’acquiert. Et
ce qu’on appelle l’intuition n’est finalement que de l’expérience tellement répétée que le raisonnement, au
départ laborieux, finit par devenir « spontané ».
il faut être méthodique et rigoureux. En effet, chaque fois qu’on écrit une série d’instructions
qu’on croit justes, il faut systématiquement se mettre mentalement à la place de la machine qui va les
exécuter, armé d'un papier et d'un crayon, afin de vérifier si le résultat obtenu est bien celui que l’on
voulait. Cette opération ne requiert pas la moindre once d’intelligence. Mais elle reste néanmoins
indispensable, si l’on ne veut pas écrire à l’aveuglette.
Et petit à petit, à force de pratique, vous verrez que vous pourrez faire de
plus en plus souvent l’économie de cette dernière étape : l’expérience fera que vous « verrez » le résultat
produit par vos instructions, au fur et à mesure que vous les écrirez. Naturellement, cet apprentissage est
long, et demande des heures de travail patient. Aussi, dans un premier temps, évitez de sauter les étapes :
la vérification méthodique, pas à pas, de chacun de vos algorithmes
représente plus de la moitié du travail à accomplir... et le gage de vos progrès.
Retour Haut de Page
3. L’ADN, les Shadoks, et les ordinateurs
Quel rapport me direz-vous ? Eh bien le point commun est : quatre mots de
vocabulaire.
L’univers lexical Shadok, c’est bien connu, se limite aux termes
« Ga », « Bu », « Zo », et « Meu ». Ce qui leur a tout de même permis de formuler quelques fortes
maximes, telles que : « Mieux vaut pomper et qu’il ne se passe rien, plutôt qu’arrêter de pomper et risquer
qu’il se passe quelque chose de pire » (pour d’autres fortes maximes Shadok, n’hésitez pas à visiter
leur site Internet, il y en a toute une collection qui vaut le détour).
L’ADN, qui est en quelque sorte le programme génétique, l’algorithme à la
base de construction des êtres vivants, est une chaîne construite à partir de quatre éléments invariables. Ce
n’est que le nombre de ces éléments, ainsi que l’ordre dans lequel ils sont arrangés, qui vont déterminer si
on obtient une puce ou un éléphant. Et tous autant que nous sommes, splendides réussites de la Nature,
avons été construits par un « programme » constitué uniquement de ces quatre briques, ce qui devrait nous
inciter à la modestie.
Enfin, les ordinateurs, quels qu’ils soient, ne sont fondamentalement capables de comprendre que
quatre catégories d'ordres (en programmation, on n'emploiera pas le terme d'ordre, mais plutôt celui
d'instructions). Ces quatre familles d'instructions sont :
l’affectation de variables
la lecture / écriture
les tests
les boucles
4. Algorithmique et programmation
Pourquoi apprendre l’algorithmique pour apprendre à programmer ? En quoi
a-t-on besoin d’un langage spécial, distinct des langages de programmation compréhensibles par les
ordinateurs ?
Parce que l’algorithmique exprime les instructions résolvant un problème
donné indépendamment des particularités de tel ou tel langage. Pour prendre une image, si un
programme était une dissertation, l’algorithmique serait le plan, une fois mis de côté la rédaction et
l’orthographe. Or, vous savez qu’il vaut mieux faire d’abord le plan et rédiger ensuite que l’inverse…
Apprendre l’algorithmique, c’est apprendre à manier la structure
logique d’un programme informatique. Cette dimension est présente quelle que soit le langage de
programmation ; mais lorsqu’on programme dans un langage (en C, en Visual Basic, etc.) on doit en plus
se colleter les problèmes de syntaxe, ou de types d’instructions, propres à ce langage. Apprendre
l’algorithmique de manière séparée, c’est donc sérier les difficultés pour mieux les vaincre.
A cela, il faut ajouter que des générations de programmeurs, souvent
autodidactes (mais pas toujours, hélas !), ayant directement appris à programmer dans tel ou tel
langage, ne font pas mentalement clairement la différence entre ce qui relève de la structure logique
générale de toute programmation (les règles fondamentales de l’algorithmique) et ce qui relève du langage
particulier qu’ils ont appris. Ces programmeurs, non seulement ont beaucoup plus de mal à passer ensuite
à un langage différent, mais encore écrivent bien souvent des programmes qui même s’ils sont justes,
restent laborieux. Car on n’ignore pas impunément les règles fondamentales de l’algorithmique… Alors,
autant l’apprendre en tant que telle !
Bon, maintenant que j’ai bien fait l’article pour vendre ma marchandise, on
va presque pouvoir passer au vif du sujet…
Définition 1.1. Un algorithme est une procédure de calcul bien définie qui prend en entrée un ensemble de
valeurs et qui délivre en sortie un ensemble de valeurs.
Exemple 1.1
Problème : trier une suite de nombres entiers dans l'ordre croissant.
Entrée : suite de n nombres entiers ( a1,a2… ana1,a2… an)
Sortie : une permutation de la suite donnée en entrée (a′1,a′2… a′na1′,a2′… an′) telle que a′1≤a′2≤⋯≤a
′na1′≤a2′≤⋯≤an′. À partir de la suite (6,9,2,4), un algorithme de tri fournira le résultat (2,4,6,9).
Définition 1.2. Une valeur particulière de l'ensemble des valeurs données en entrée est appelée instance du
problème.
Définition 1.3. Un algorithme est correct, si pour toute instance du problème il se termine et produit une
sortie correcte.
Les algorithmes peuvent être spécifiés en langage humain ou tout langage informatique. Dans ce qui suit, nous
utiliserons un langage proche du langage naturel. Nous donnerons une implémentation en Python (voir
cours MISMI MIS 102).
Définition 1.4. Une heuristique est une procédure de calcul correcte pour certaines instances du problème
(c'est-à-dire se termine ou produit une sortie correcte).
Pour qu'un algorithme puisse être décrit et s'effectue, les données d'entrées doivent être organisées.
Définition 1.5. Une structure de données est un moyen de stocker et d'organiser des données pour faciliter
leur stockage, leur utilisation et leur modification.
bio-informatique ;
moteur de recherche sur Internet ;
commerce électronique ;
affectation de tâches.
Définition 1.6. L'efficacité d'un algorithme est mesurée par son coût (complexité) en temps et en mémoire.
Un problème NP-complet est un problème pour lequel on ne connaît pas d'algorithme correct efficace, c'est-à-
dire réalisable en temps et en mémoire. Le problème le plus célèbre est le problème du voyageur de
commerce.
L'ensemble des problèmes NP-complets ont les propriétés suivantes :
L'efficacité d'un algorithme est fondamentale pour résoudre effectivement des problèmes.
Exemple1.2.
Supposons que l'on dispose de deux ordinateurs. L'ordinateur A est capable d'effectuer 10 9 instructions par
seconde. L'ordinateur B est capable d'effectuer 10 7 instructions par seconde. Considérons un même problème
(de tri par exemple) dont la taille des données d'entrées est n. Pour l'ordinateur A, on utilise un algorithme qui
réalise 2n2 instructions. Pour l'ordinateur B, on utilise un algorithme qui réalise 50nlog(n)50nlog(n)instructions.
Pour traiter une entrée de taille 10 6 : l'ordinateur A prendra 2000 s et l'ordinateur B prendra 100 s. Ainsi, même
si la machine B est médiocre, elle résoudra le problème 20 fois plus vite que l'ordinateur A.
en temps, le nombre d'opérations élémentaires effectuées pour traiter une donnée de taille n ;
en mémoire, l'espace mémoire nécessaire pour traiter une donnée de taille n.
Dans ce cours, nous considérerons que la complexité des instructions élémentaires les plus courantes sur un
ordinateur ont un temps d'exécution que l'on considérera dans ce cours comme constant égal à 1. Les
instructions élémentaires sont : addition, multiplication, modulo et partie entière, affectation, instruction de
contrôle. Ce qui intéresse fondamentalement l'algorithmique, c'est l'ordre de grandeur (au voisinage de l'infini)
de la fonction qui exprime le nombre d'instructions. Les courbes de références sont ici.
Il est nécessaire de disposer d'un langage qui soit non lié à l'implémentation. Ceci permet une description plus
précise des structures de données ainsi qu'une rédaction de l'algorithme plus souple et plus « lisible ». Le
langage EXALGO est un exemple de ce qui peut être utilisé et qui sera utilisé dans ce cours. Il est composé de
chaînes de caractères alphanumériques, de signes opératoires (+, -, *, /, <, <=, >=, >, <>, ==, =, ou, non,
et), de mot-clés réservés, et de signes de ponctuation : ''=, ;, (, ), début, fin, //. Les
balises début et fin peuvent être remplacées par { et }.
Python n'utilise pas de marqueurs de fin. Le caractère « : » est le marqueur de début et quand l'indentation
cesse Python considère que c'est un marqueur de fin.
d'un nom ;
d'un ensemble de valeurs ;
d'un ensemble d'opérations définies sur ces valeurs.
On écrit en EXALGO
Sélectionnez
var NomDeVariable : Type;
Type est à prendre pour l'instant dans l'ensemble {entier, car, booléen, réel}.
Définition 2.3. Les Expressions sont constituées à l'aide de variables déjà déclarées, de valeurs, de
parenthèses et d'opérateurs du (des) type(s) de variables concernées.
Définition 2.4. L'affectation est l'instruction qui permet de stocker une valeur dans une variable.
On écrit :
Sélectionnez
NomDeVariable = ExressionDuTypeDeLaVariable;
Une variable de type booléen prend comme valeur VRAI ou FAUX. Les opérations usuelles sont ET, OU et NON
qui sont données dans les tables qui suivent :
Une variable de type entier peut prendre comme valeur l'ensemble des nombres entiers signés. Les opérations
associées sont les opérations usuelles +,-,*,/.
Une variable de type réel peut prendre comme valeur l'ensemble des nombres réels. Les opérations associées
sont les opérations usuelles +,-,*,/.
Une variable de type car peut prendre comme valeur l'ensemble des caractères imprimables. On notera les
valeurs entre guillemets. On considère souvent que les caractères sont ordonnés dans l'ordre alphabétique.
Les valeurs :
sont différentes et ne seront pas codées de la même manière dans la mémoire de la machine.
Les opérateurs <, ≤, ==, !=, >, ≥ permettent de comparer les valeurs de type entier, réel et caractère. Le
résultat de cette comparaison est une valeur booléenne.
Bloc d'instructions :
Sélectionnez
début
instruction1
instruction2
..........
fin
Alternative :
o Alternative simple :
Traduction Python
Sélectionnez
Sélectionnez
si ExpressionBooléenne alors
BlocInstructions1
sinon
BlocInstructions2
finsi;
Alternative multiple (traduction Python) :
Traduction Python
Sélectionnez
Sélectionnez
selon que
cas cas1 : BlocInstructions1
cas cas2 : BlocInstructions2
..........
autrement : BlocInstruction
finselonque
Répétition
L'instruction exit permet d'arrêter la répétition.
o Le bloc d'instructions peut ne pas être exécuté :
Traduction Python
Sélectionnez
Sélectionnez
tant que ExpressionBooléenne faire
BlocInstructions
fintantque;
Le bloc d'instructions peut ne pas être exécuté et il y a une variable indicatrice :
Traduction Python
Sélectionnez
Sélectionnez
pour VariableIndicatrice
allant de ValeurInitiale à ValeurFinale
par pas de ValeurPas faire
BlocInstructions
finpour;
Le bloc d'instructions est exécuté au moins une fois (ne se traduit pas directement en Python) :
Sélectionnez
répéter
BlocInstructions
jusqu'à ExpressionBooléenne finrépéter;
Une fonction est une section d'algorithme qui a un objectif bien défini et un nom. En général, elle communique
avec l'extérieur par le biais de paramètres typés. Elle possède des variables locales qui ne sont pas visibles à
l'extérieur de la fonction. Ces variables peuvent être des fonctions. Une fonction retourne une valeur par
l'instruction simple retourne(Expression). L'expression peut être :
vide, tout s'est bien passé, mais il n'y a pas de résultat à retourner : retourne() ;
sans résultat, il est impossible de retourner un résultat suite à un cas de figure de
l'instance : retourne(NUL).
Écriture de la fonction :
Sélectionnez
fonction NomDeFonction (ListeParamètres) :TypeRésultat;
// déclarations des variables ou fonctions locales autres que
les paramètres
début
// partie instruction qui contient l'appel à retourne
fin
finFonction
Liste des paramètres
Les paramètres sont passés :
o par référence ref, on écrit :
Sélectionnez
ref ListeVariable NomDeType
Sélectionnez
val ListeVariable:NomDeType
Sélectionnez
NomDeFonction(ListeInstanceParamètres)
Sélectionnez
fonction exemple(val n :entier;ref m : entier) :vide;
début
n = 5;
m = 7;
fin
finFonction
Sélectionnez
var p,q :entier;
début
p = 1;
q = 2;
exemple(p,q);
fin
EXALGO permet de fixer les quelques règles élémentaires permettant d'écrire des algorithmes en
s'affranchissant l'implémentation.
Le langage EXALGO est composé de chaînes de caractères alphanumériques, de signes opératoires, de mot-clés
réservés, et de signes de ponctuation : =, ;,(,), début, fin, //. Les marqueurs de fin, début et fin peuvent être
remplacés par { et } lorsqu'il y a encombrement.
Sélectionnez
var NomDeVariable: TypePrédéfini;
Constituées à l'aide de variables déjà déclarées, de parenthèses et d'opérateurs du (des) type(s) des variables
concernées.
affectation :
Sélectionnez
NomDeVariable = ExressionDuTypeDeLavariable;
Bloc d'instruction :
Bloc d'instruction
Sélectionnez
Instruction1
instruction2
.............
Alternative :
Sélectionnez
si ExpressionBooléenne alors
BlocInstruction1
sinon
BlocInstruction2
finsi;
Alternative multiple :
Sélectionnez
selon que
cas cas1 : BlocInstruction1
cas cas2 : BlocInstruction2
.............
autrement : BlocInstruction
finselonque
Répétition : exit permet d'arrêter la répétition
o le bloc d'instruction peut ne pas être exécuté
Sélectionnez
tant que ExpressionBooléenne faire
BlocInstruction
fintantque;
le bloc d'instruction peut ne pas être exécuté et il y a une variable indicatrice
Sélectionnez
pour VariableIndicatrice
allant de ValeurInitiale à ValeurFinale
par pas de ValeurPas faire
BlocInstruction
finpour;
Sélectionnez
répéter
BlocInstruction
jusqu'à ExpressionBooléenne finrépéter;
Une fonction retourne une valeur par l'instruction simple (retourne(Expression)). Une
fonction s'utilise dans le calcul d'une expression ou comme instruction simple.
Écriture de la fonction :
Sélectionnez
fonction NomDeFonction (ListeParamètres):TypeRésultat;
// déclarations des variables locales autres que les paramètres
début
// partie instruction qui contient l'appel à retourne()
fin
finFonction
Liste des paramètres .
Les paramètres sont passés :
o par référence ref, on écrit :
Sélectionnez
ref ListeVariable:NomDeType
par valeur val, on écrit :
Sélectionnez
val ListeVariable:NomDeType
Un type structuré est constitué à partir de types de base ou d'autres types déclarés.
Sélectionnez
type NomDeType: structure
champ1:NomDeType1
champ2:NomDeType2
…
finstructure
Après la déclaration :
Sélectionnez
var E:NomDeTypeEnregistrement
on accède aux différents champs par le nom de la variable suivi d'un point suivi du nom de champ (E.champ1).
Sélectionnez
var P:^NomDeType
alors on peut obtenir un objet accessible par allouer(P). Lorsqu'on n'utilise plus l'objet, il faut libérer l'espace qu'il
utilise par desallouer(P).
Définition 3.1. Une séquence sur un ensemble E est une suite d'éléments (e1,e2,…en) d'éléments de E.
Exemple 3.1 (3,5,8,2,12,6) : est une séquence d'éléments de N, ensemble des entiers naturels.
("a","z","T","A","a") est une séquence sur l'ensemble des caractères imprimables (char).
Il existe plusieurs variantes de séquences suivant les opérations de manipulation autorisées : accès par l'indice
de l'élément ou non, accès à la fin de la séquence ou non…
On utilisera en général des noms particuliers dépendant des caractéristiques de la séquence.
Exemple 3.2 : Un vecteur peut être défini par une séquence dans laquelle l'accès aux éléments se fait par son
indice et la taille de la séquence dépend de l'espace dans lequel on se trouve. On dit aussi qu'on a un accès
direct à l'élément. Dans la plupart des langages de programmation, le vecteur existe sous le nom d'array.
Sélectionnez
fonction fac(val n :entier) :entier;
début
si n <= 1 alors
retourner(1)
sinon
retourner(n * fac(n-1))
finsi
fin
finfonction
Programme Python
Sélectionnez
La séquence des valeurs de n au cours des appels récursifs doit être mémorisée. Supposons l'appel fac(4) alors :
Définition 3.2. Une structure sur F1×F2×⋯×FpF1×F2×⋯×Fp est une séquence (f1, f2,…, fk)(f1, f2,
…, fk) telle que
∀i∈[1..k],fi∈Fi∀i∈[1..k],fi∈Fi.
Les structures sont des cas particuliers de séquences. En algorithmique, chaque ensemble FiFi peut être un type
de base ou une structure. Ce mécanisme permet de définir de nouveaux types plus complexes que les types de
base. En EXALGO, on écrit :
Sélectionnez
nom_du_type = structure
nom_champs_1 :type1;
nom_champs_2 :type2;
…
nom_champs_k :typek;
finstructure
Cela signifie que lorsqu'une variable est déclarée de ce type, elle référence k variables en même temps. Soit V
une variable dont le type est une structure, on désigne un des champs par V. suivi du nom du champ.
Exemple 3.4 : Une date de naissance est un exemple de structure. On peut écrire :
Sélectionnez
dateDeNaissance = structure
jourDeNaissance: entier;
moisDeNaissance: entier;
annéeDeNaissance: entier;
finstructure
Soit la déclaration :
Sélectionnez
var I :individu
alors [Link] sera un booléen et [Link] sera un entier. Ainsi les instructions suivantes ont un sens :
Sélectionnez
[Link] = 12;
[Link] = faux;
Définition 3.3 : Soit F un ensemble. Une table d'association à clé unique est une séquence d'éléments
de N×FN×F (N est l'ensemble des entiers naturels), ((c1,f1), (c2,f2),…, (ck,fk))((c1,f1), (c2,f2),
…, (ck,fk)) telle que :
∀i,j∈[1..k],i≠j,ci≠cj∀i,j∈[1..k],i≠j,ci≠cj
Les tables d'association sont un cas particulier de séquences d'éléments structurés. La structure se décrit en
EXALGO :
Sélectionnez
association = structure
cle :entier;
valeur :type_prédéfini;
finstructure
Exemple 3.5 : Lors de l'activation du compte électronique, l'étudiant de l'Université Bordeaux 1 fournit un
numéro INE qui sera associé à un mot de passe. On a donc quelque part dans le système de gestion des comptes
une table d'association à index unique dont l'élément de séquence est :
Sélectionnez
Etudiant = structure
INE :entier;
motDePasse :typeMotDePasse;
finstructure
k,a>0k,a>0 tels que ∀x>a,|f(x)|≤k|g(x)|∀x>a,|f(x)|≤k|g(x)|.
Définition 4.1. (Notation de Landau). On dit que f=O(g)f=O(g) s'il existe deux nombres
réels
Exemple 4.1. Si le nombre d'instructions est égal à f(n)=an2+bn+cf(n)=an2+bn+c avec a,b,c des
constantes réelles, alors f(n)=O(n2)f(n)=O(n2).
Les figures permettent de comparer les fonctions usuelles utilisées pour décrire la complexité d'un algorithme
en fonction de la taille n des données d'entrées. Parmi les fonctions usuelles, le log à base 2
de nlog2(n)nlog2(n) joue un rôle important. Pour un algorithme A, notons CA(D)CA(D), le coût de l'algorithme
A pour une instance D.
Soit DnDn l'ensemble des instances de taille n. Si toutes les instances sont équiprobables, on a :
Parfois, il est nécessaire d'étudier la complexité en mémoire lorsque l'algorithme requiert de la mémoire
supplémentaire (donnée auxiliaire de même taille que l'instance en entrée par exemple).
opérations élémentaires +, -, *, / ;
test d'expression booléenne ;
appel de fonctions.
Source Python
Sélectionnez
On a :
C>suite(n)=C<suite(n)=Csuite(n)=O(n)Csuite>(n)=Csuite<(n)=Csuite(n)=O(n)
Entrée : un entier n
Sortie : « vrai » si on rencontre une pile, « faux » sinon.
La fonction suivante retourne « vrai » lorsque l'un des lancers est égal à 6 et « faux » sinon.
Sélectionnez
fonction jeuDePile(val n :integer) :booléen;
var i : entier;
début
pour i allant de 1 à n faire
f = résultat_lancer_pièce()
si (f==pile) alors
retourner(vrai)
finsi
finpour
retourner(faux)
fin
finFonction
Source Python
Sélectionnez
On écrit en EXALGO :
Sélectionnez
nom_tableau = tableau[min_indice..max_indice] de type_predefini;
La taille du tableau est donc max_indice - min_indice + 1. Pour accéder à un élément d'un tableau T d'indice I, on
écrit T[I]. La complexité de l'accès à un élément du tableau est O(1)O(1).
Soit min_indice< i<j <max_indice, on notera T[i..j] la séquence des éléments de T (T[i],T[i+1],…,T[j]).
Beaucoup d'algorithmes peuvent être décrits sans préciser un type particulier. Dans ce cas, on écrira à la place
de type_prédéfini le mot élémentet on précisera les valeurs possibles pour élément.
TC = tableau[1..10] de car; ;
TE = tableau[1..10] d'entiers;.
L'algorithme qui permet de trier TC et TE est le même. Seul diffère le type de l'élément manipulé. On écrira dans
ce cas un algorithme sur un tableau.
Sélectionnez
T = tableau[1..10] d'éléments;
Les paramètres tableaux doivent, sauf raison majeure, être passés en paramètre par référence afin d'éviter
la recopie.
Sélectionnez
fonction init(ref T :tableau[min_indice..max_indice] d'éléments;
val valeurInitiale :élément) :vide;
var i :entier;
début
pour i allant de min_indice à max_indice faire
T[i] = valeurInitiale
finpour
fin
finfonction
Complexité : O(n)O(n)
Sélectionnez
fonction taille(ref T :tableau[min_indice..max_indice] d'éléments) :entier;
début
retourner(max_indice - min_indice + 1)
fin
finfonction
Complexité : O(1)O(1).
Sélectionnez
fonction echange(ref T :tableau[min_indice..max_indice] d'éléments;
val indice1,indice2 : entier) :vide;
var e :élément;
début
e = T[indice1];
T[indice1] = T[indice2];
T[indice2] = e;
fin
finfonction
Complexité : O(1)O(1).
Sélectionnez
fonction copie(ref T1,T2 :tableau[min_indice..max_indice] d'élément;
val indiceT1_1,indiceT1_2,indiceT2 : entier) :booléen;
var i :entier;
début
si indiceT2+indiceT1_2-indiceT1_1>max_indice alors
retourner(faux)
sinon
pour i allant de indiceT1_1 à indiceT1_2 faire
T2[indiceT2] = T1[i];
indiceT2 = indiceT2 + 1;
finpour
retourner(vrai)
fin
finfonction
Complexité :
minimum : O(1)O(1)
maximum : O(n)O(n)
Source Python
Sélectionnez
Sélectionnez
fonction somme(ref T :tableau[min_indice..max_indice] d'entiers) :entier;
var s,i :entier;
début
s = 0;
pour i allant de min_indice à max_indice faire
s = s + T[i]
finpour
retourner(s)
fin
finfonction
Complexité : O(n)O(n)
Propriété 5.2. Soit i,j deux entiers, i<=j. Soit T un tableau d'éléments d'indice variant entre i et j. Pour tout
élément e, appartenant au tableau T, on a :
Sélectionnez
fonction cherche(ref T :tableau[min_indice..max_indice] d'éléments;val
e :élément) :entier;
var i :entier;
début
pour i allant de min_indice à max_indice faire
si T[i]==e alors
retourner(i)
finsi
finpour
retourner()
fin
finfonction
Complexité :
minimum : O(1)O(1)
maximum : O(n)O(n)
Source Python
Sélectionnez
Propriété 5.3. Soit i,j deux entiers, i<=j. Soit T un tableau d'entiers d'indice variant entre i et j. Soit m l'élément
minimum du tableau, on a :
Sélectionnez
fonction minimum(ref T :tableau[min_indice..max_indice] d'entier) :entier;
var i,sauv :entier;
début
sauv = min_indice;
pour i allant de min_indice+1 à max_indice faire
si T[i]<T[sauv] alors
sauv = i
finsi
finpour
retourner(sauv)
finfonction
Complexité : O(n)O(n)
Sélectionnez
var M :tableau[1..n] de tableau [1..m] d'éléments;
Sélectionnez
fonction initMatrice(ref M :tableau[1..n] de tableau [1..m] d'éléments;
val valeurInitiale :élément) :vide;
var i,j :entier;
début pour i allant de 1 à n faire
pour j allant de 1 à m faire
M[i][j] = valeurInitiale
finpour
finpour
retourner()
finfonction
Complexité : O(nm)O(nm)
Sélectionnez
fonction sommeMatrice(ref M1,M2 :tableau[1..n] de tableau [1..m] de réels) :
tableau[1..n] de tableau [1..m] de réels;
var i,j :entier;
var M :tableau[1..n] de tableau [1..m] de réels;
début
pour i allant de 1 à n faire
pour j allant de 1 à m faire
M[i][j] = M1[i][j] + M2[i][j];
finpour
finpour
retourner(M)
finfonction
Complexité : O(nm)O(nm)
On considérera dans tout ce chapitre que l'on manipule des entiers. L'objet du tri est d'ordonner une séquence
de N entiers. On considérera que ces entiers sont rangés dans un tableau.
Sélectionnez
var T :tableau[1..N] d'entiers;
Ce tri est basé sur l'algorithme de recherche du minimum On adapte cet algorithme pour pouvoir effectuer la
recherche dans un sous-tableau. On a le déroulement ici.
Sélectionnez
fonction minimumSoustableau(ref T :tableau[1..N] d'entiers, val
Imin,Imax :entier) :entier;
var sauv :entier;
début
sauv = Imin;
pour i allant de Imin+1 à Imax faire
si T[i]<T[sauv] alors
sauv = i;
finsi
finpour
retourner(sauv);
fin
finfonction
Propriété 6.1. La complexité de l'algorithme triSelection sur une instance de taille N est O(n2)O(n2)
Propriété 6.2. Soit T un tableau d'entiers trié d'indice variant entre i et j. Soit e un entier quelconque, alors on a
l'une des propriétés suivantes :
e ≤ T[i] ;
il existe un unique entier k dans [i..j-1] tel que T[k] < e ≤ T[k+1] ;
e > T[j].
Sélectionnez
fonction triInsertion(ref T :tableau[1..N] d'entiers) :vide;
var i,j,cle :entier;
début
pour i allant de 2 à N faire
cle = T[i];
j = i-1;
tant que j>0 et T[j]>cle faire
T[j+1] = T[j];
j = j-1;
fintantque
T[j+1] = cle;
finpour
fin
finfonction
On a le déroulement ici.
Propriété 6.3. La complexité de l'algorithme triInsertion sur une instance de taille N est :
au minimum en O(N)O(N) ;
au maximum et en moyenne en O(N2)O(N2).
Idée de la démonstration
La boucle « pour » s'effectue systématiquement et demandera O(N)O(N) opérations.
La boucle « tant que » effectue au minimum O(1)O(1) opération (cas où les nombres sont déjà triés) et au
maximum O(N)O(N).
La boucle « tant que » effectue en moyenne O(N/2)O(N/2) opérations.
Sélectionnez
fonction triBulle(ref T :tableau[1..N] d'entiers) :vide;
var i,j,cle :entier;
début
pour i allant de 1 à N-1 faire
pour j allant de N à i+1 par pas de -1 faire
si T[j] < T[j-1] alors
echange(T, j, j-1);
finsi
finpour
finpour
fin
finfonction
On a le déroulement ici.
Propriété 6.4. La complexité de l'algorithme triBulle sur une instance de taille N est O(N2)O(N2).
Lorsque deux tableaux T1 et T2 sont triés, il est aisé de construire un nouveau tableau contenant la séquence
triée regroupant les séquences correspondantes à T1 et T2.
Sélectionnez
PREMIERE VERSION
fonction fusion(ref T1 :tableau[1..N1] d'entier;
ref T2 :tableau[1..N2] d'entier
) :tableau[1..N1+N2] d'entier;
var I1,I2,i :entier;
var T :tableau[1..N1+N2] d'entier;
début
I1 = 1;
I2 = 1;
pour i allant de 1 à N1+N2 faire
si T1[I1]≤T2[I2] alors
T[i] = T1[I1];
I1 = I1+1;
sinon
T[i] = T2[I2];
I2 = I2+1;
finsi
finpour
retourner(T)
fin
finfonction
Complexité : O(n)O(n)
Cette version ne fonctionne pas toujours. Par exemple, si I1 a dépassé N1 et vaut par exemple N1+1, on
comparera T1[N1+1] à T2[I2] ce qui n'a pas de sens. Il faut donc utiliser un algorithme exact. On a le
déroulement ici.
Soit une séquence d'éléments de [0..k], il est alors possible de réaliser l'histogramme des valeurs. Par la suite le
tri des éléments de la séquence se fait en temps linéaire O(n)O(n).
Sélectionnez
fonction triHisto(ref T :tableau[1..N] d'entiers) :vide;
var H :tableau[0..maximum(T)] d'entier;
var i,j,k,max : entier;
début
init(H,0);
pour i allant de 1 à N faire
H[T[i]] = H[T[i]] + 1;
finpour
i = 1;
max = maximum(T);
pour j allant de 0 à max faire
pour k allant de 1 à H[j] faire
T[i] = j;
i = i+1;
finpour
finpour
fin
finfonction
On a le déroulement ici.
On considère une nouvelle fonction copie qui copie un tableau dans un autre même s'ils n'ont pas la même
définition. L'en-tête de la fonction est :
Sélectionnez
fonction copie(ref T1:tableau[1..N1] d'élément;
ref T2:tableau[1..N2] d'élément;
val indiceT1_1,indiceT1_2,indiceT2: entier):vide;
Sélectionnez
fonction fusion(ref T1:tableau[1..N1] d'entier;
ref T2:tableau[1..N2] d'entier
):tableau[1..N1+N2] d'entier;
var I1,I2,i:entier;
var T:tableau[1..N1+N2] d'entier;
début
Initialiser I1,I2,i
tant que I1 ≤ N1 et I2 ≤ N2 faire
// ... On compare tant qu'il reste des éléments dans T1 et T2
fintantque
si I1 ≤ N1 alors
// il n'y a plus d'éléments dans T2 :copier(..........)
sinon
// il n'y a plus d'éléments dans T1 :copier(..........)
finsi
retourner(T)
fin
finfonction
Sélectionnez
fonction fusion(ref T1:tableau[D1..N1] d'entier;
ref T2:tableau[D2..N2] d'entier
):tableau[1..N1+N2-D1-D2+2] d'entier;
var I1,I2,i: entier;
var T:tableau[1..N1+N2-D1-D2+2] d'entier;
début
i = 1;
I1 = D1;
I2 = D2;
tant que I1 ≤ N1 et I2 ≤ N2 faire
si T1[I1] ≤ T2[I2] alors
T[i] = T1[I1];
I1 = I1+1;
sinon
T[i] = T2[I2];
I2 = I2+1;
finsi
i = i+1
fintantque
si I1 ≤ N1 alors
copier(T1,T,I1,N1,i)
sinon
copier(T2,T,I2,N2,i)
finsi
retourner(T)
fin
finfonction
Sélectionnez
fonction fusion(ref T1:tableau[D1..N1] d'entier;
ref T2:tableau[D2..N2] d'entier;
val DT1,FT1,DT2,FT2:entier):
tableau[1..FT1+FT2-DT1-DT2+2] d'entier;
var I1,I2,i:entier;
var T:tableau[1..FT1+FT2-DT1-DT2+2] d'entier;
début
i = 1;
I1 = DT1;
I2 = DT2;
tant que I1 ≤ FT1 et I2 ≤ FT2 faire
si T1[I1] ≤ T2[I2] alors
T[i] = T1[I1];
I1 = I1+1;
sinon
T[i] = T2[I2];
I2 = I2+1;
finsi
i = i+1
fintantque
si I1 ≤ FT1 alors
copier(T1,T,I1,FT1,i)
sinon
copier(T2,T,I2,FT2,i)
finsi
Retourner(T)
fin
finfonction
Comme vu au chapitre Codage et structures de contrôle, on peut déclarer dans une fonction des variables et
des fonctions locales :
Sélectionnez
fonction NomDeFonction (ListeParamètres) :TypeRésultat;
// déclarations des variables ou fonctions locales
début
// partie instruction qui contient l'appel à retourner
fin
finFonction
La multi-imbrication possible des fonctions entraîne l'existence de problèmes de visibilité : entre les variables et
entre les fonctions.
Règle 1 : une variable V (locale ou non) est visible depuis sa déclaration jusqu'au
marqueur finFonction de la fonction F où elle a été déclarée.
Règle 2 : si une fonction G est locale à F et déclare une variable V déjà déclarée
dans F alors la variable originelle est momentanément cachée.
Soit la fonction P suivante :
Sélectionnez
fonction P (....) :....;
var x,y,z : entier ;
début
x = 1;
y = 2;
z = 3;
R() …
Q(z);
fin
finFonction
On a le déroulement ici.
Une fonction est visible depuis la fin de son entête jusqu'au finFonction de la fonction où elle a été déclarée.
Cependant comme pour les variables, elle peut momentanément être cachée par une autre fonction ayant le
même entête (surcharge).
La fonction P suivante est annotée pour préciser la visibilité des fonctions Q,R,T.
Sélectionnez
fonction P(....) :....;
.....
fonction Q(....) :.....;
.....
fonction R(...) :.....;
....
début
....// on peut utiliser P,Q,R
fin
finFonction ;
début
....// on peut utiliser P,Q,R
fin
finFonction
fonction T(...) :...;
début
....// on peut utiliser P,Q,T mais pas R
fin
finFonction ;
début
... //// on peut utiliser P,Q,T mais pas R
fin
finFonction
La récursivité consiste à remplacer une boucle par un appel à la fonction elle-même. Considérons la suite
factorielle, elle est définie par :
0!=10!=1 ;
n!=n(n−1)!n!=n(n−1)!.
Sélectionnez
fonction factorielle(val n :entier) :entier;
début
si (n == 0)
retourne(1)
sinon
retourne(factorielle(n-1) * n)
finsi
fin
finfonction;
On a le déroulement ici. On peut décrire sur le papier les changements et les appels sous la forme suivante :
Plusieurs appels à la fonction peuvent être exécutés dans son corps. Soit la suite dite de Fibonacci définie par :
u0=1u0=1 ;
u1=1u1=1 ;
un=un+1+un+2 pour n>2un=un+1+un+2 pour n>2.
Sélectionnez
fonction fibo(val n :entier) :entier;
début
si (n == 0) ou (n == 1) alors
retourne(1)
sinon
retourne(fibo(n-1) + fibo(n-2))
finsi
fin
finfonction;
On a le déroulement ici. On peut décrire sur le papier les changements et les appels sous la forme suivante :
Examinons la suite définie par :
u1=1u1=1 ;
un=un−1+n pour n>1un=un−1+n pour n>1.
Sélectionnez
fonction suite(val n :entier) :entier;
var i,s :entier;
début
s = 0;
pour i allant de 1 à n faire
s = s+i;
finpour;
retourner(s)
fin
finfonction;
Sélectionnez
fonction suiteR(val n :entier) :entier;
début
si n == 1 alors
retourne(1)
sinon
retourne(suiteR(n-1) + n)
finsi
fin
finfonction;
La complexité en nombre d'opérations de suite et suiteR est en O(n)O(n). On aurait donc tendance à
préférer suiteR pour sa lisibilité. Cependant, si on examine la complexité en mémoire, suite est
en O(1)O(1) alors que suiteR est en O(n)O(n).
La programmation non récursive est donc plus efficace. L'utilisation de la récursivité ne doit pas se faire
au détriment de l'efficacité.
Chaque fois que l'on désire programmer une fonction récursive, on doit répondre aux questions suivantes :
Sélectionnez
fonction cherche(ref T :tableau[min_indice..max_indice] d'éléments;
val e : élément) :entier;
début
si T[min_indice] == e alors
retourner(min_indice)
sinon
si min_indice == max_indice alors
retourner(NUL)
sinon
retourner(cherche(T[min_indice+1..max_indice], e))
finsi
finsi
fin
finfonction
Sélectionnez
fonction minimumTableau(ref T :tableau[1..N] d'entiers;
val Imin :entier) :entier;
var sauv :entier;
début
si Imin == N alors
retourner(T[N])
sinon
sauv = minimumTableau(T, Imin+1];
si T[Imin] < sauv alors
retourner(T[Imin])
sinon
retourner(sauv)
finsi
finsi
fin
finfonction
La dichotomie fait partie des méthodes dites « diviser pour régner ». Elle consiste pour un objet de taille N à
exécuter un algorithme de façon à réduire le problème à un objet de taille N/2. On répète alors l'algorithme de
réduction sur ce dernier objet. Ainsi, il suffit de connaître la résolution pour un problème de taille faible
(typiquement N=1 ou N=2) pour obtenir la totalité de la ré[Link] type d'algorithme est souvent implémenté
de manière récursive. Lorsque cette technique est utilisable, elle conduit à un algorithme très efficace et très
lisible.
Il est parfois nécessaire de traiter les données avant d'appeler la fonction récursive. La fonction récursive est
alors une fonction locale à la fonction d'appel.
Soit g une fonction croissante sur un intervalle [a,b] et telle que f(a)≤0 et f(b)≥0. L'algorithme ci-dessous permet
de trouver la valeur x de [a,b] telle que f(x)=0 avec une précision e.
Sélectionnez
fonction zero(ref g(val n :réel) :fonction;val a,b,e :réel) :réel;
var M :réel;
début
M = g((a+b)/2);
si |M| < e alors
retourne((a+b)/2)
sinon
si M > 0 alors
zero(g, a, (a+b)/2, e)
sinon
zero(g, (a+b)/2, b, e)
finsi
finsi
fin
finfonction
Nous avons déjà traité cet algorithme sous une autre forme au chapitre Tableaux.
Propriété 8.1. T un tableau d'entiers triés d'indice variant entre d et f. Posons m=⌊(d+f)/2⌋m=⌊(d+f)/2⌋.
Soit e un entier appartenant à la séquence contenue dans T. On a l'une des propriétés suivantes :
T[m] = e ;
e est dans la séquence contenue dans T[d..m-1] ;
e est dans la séquence contenue dans T[m+1..f].
Sélectionnez
fonction cherche(ref T :tableau[1..N] d'entiers; val e :entier) :entier;
var d,f :entier;
fonction chercheRec(ref T :tableau[1..N] d'entiers; val d,f,e :entier) :entier;
var m;
début
si f == d alors
si T[d] == e alors
retourner(f)
sinon
retourner(NUL)
finsi
sinon
m = partieEntiere((d+f)/2);
si T[m] < e alors
retourner(chercheRec(T,m+1,f,e))
sinon
retourner(chercheRec(T,d,m,e))
finsi
finsi
fin
finfonction
début
d = 1;
f = N;
retourner(chercheRec(T,d,f,e))
fin
finfonction
Idée de la preuve : la complexité de la fonction cherche est donnée par la complexité de chercheRec.
Soit f(n)f(n) le nombre de tests effectués par cette fonction.
On a :
f(n)f(1)=2+f(⌊n/2⌋)=2f(n)=2+f(⌊n/2⌋)f(1)=2
La complexité en temps de l'algorithme est donc déterminée par une équation de récurrence de la forme :
C(n)=aC(n/b)+d(n)C(n)=aC(n/b)+d(n)
qui après résolution permet de montrer que cette méthode conduit à des algorithmes plus efficaces en nombre
d'opérations. Cependant, cela ne doit pas occulter l'aspect mémoire. La complexité en mémoire doit rester d'un
ordre raisonnable. (cf. récursivité).
On considérera dans tout ce chapitre que l'on manipule des entiers. L'objet du tri est d'ordonner une séquence
de N entiers. On considérera que ces entiers sont rangés dans un tableau :
Sélectionnez
var T :tableau[1..N] d'entiers;
Cet algorithme consiste à diviser la séquence d'entiers en deux sous-séquences, à les trier de manière récursive,
puis à fusionner les deux sous-séquences triées. On utilise la fonction fusion vue au chapitre tris non
récursifs.
Sélectionnez
fonction triFusion(ref T :tableau[1..N] d'entiers) :vide;
var d,f :entier;
fonction fusionLocal(ref T :tableau[1..N] d'entier;
val d,m,f :entier) :vide;
var C :tableau[1..f-d+1] d'entier;
début
C = fusion(T,T,d,m,m+1,f);
copie(C,T,d,f,d);
fin
finfonction
fonction triFusionRec(ref T :tableau[1..N] d'entiers; val d,f :entier) :vide;
début
si d<f alors
m = partieEntiere((d+f)/2);
triFusionRec(T,d,m);
triFusionRec(T,m+1,f);
fusionLocal(T,d,m,f);
finsi
fin
finfonction
début
trifusionRec(T,1,N)
fin
finfonction
On a le déroulement ici.
Propriété 9.1. La complexité du tri fusion pour une séquence de n éléments est O(nlog2(n))O(nlog2(n)).
idée de la preuve : la complexité de la fonction triFusion est donnée par la complexité de triFusionRec.
Soit f(n)f(n) le nombre de tests effectués par cette fonction. On a :
f(n)f(1)=1+2f(⌊n/2⌋)+3n=0f(n)=1+2f(⌊n/2⌋)+3nf(1)=0
On en déduit :
f(n)=∑i=0p2i+3pnf(n)=∑i=0p2i+3pn
et
f(n)=2p+1+3pn−1f(n)=2p+1+3pn−1
Cet algorithme consiste à utiliser une valeur x de la séquence pour diviser la séquence d'entiers en deux sous-
séquences :
Sélectionnez
fonction triRapide(ref T :tableau[1..N] d'entiers) :vide;
fonction diviserSequence(ref T :tableau[1..N] d'entier;
val d,f :entier) :entier;
var x,i :entier;
début
x = T[f];
i = d-1;
pour j allant de d à f-1 faire
si T[j] ≤ x alors
i = i+1;
echanger(T,i,j);
finsi
finpour
echanger(T,i+1,f);
retourner(i+1);
fin
finfonction
fonction triRapideRec(ref T :tableau[1..N] d'entiers; val d,f :entier) :vide;
var p :entier;
début
si d<f alors
p = diviserSequence(T,d,f);
triRapideRec(T,d,p-1);
triRapideRec(T,p+1,f);
finsi
fin
finfonction
début
triRapideRec(T,1,N)
fin
finfonction
On a le déroulement ici.
Propriété 9.2. La complexité du tri rapide pour une séquence de n éléments est :
au maximum en O(n2)O(n2) ;
en moyenne et au minimum en O(nlog(n))O(nlog(n)).
Décrire une structure permettant de gérer des polynômes définis sur les réels. Écrire un ensemble de primitives
associées permettant les principales opérations.
Sélectionnez
constante Taille = 200;
type polynome = structure
coeff:tableau[0..200] de réels;
degre: entier;
finstructure
Sélectionnez
fonction ajoute(ref p1,p2:polynome):polynome;
var p:polynome;
var m,i:entier;
début
m = min([Link],[Link]);
pour i de 0 à m faire
[Link][i] = [Link][i] + [Link][i]
finpour;
si m < [Link] alors
pour i de m+1 to p1[degre] do
[Link]][i] = [Link]][i];
finpour:
[Link] = [Link];
sinon
pour i de m+1 to p2[degre] do
[Link][i] = [Link][i];
finpour:
[Link] = [Link];
finsi
retourner(p);
fin;
finfonction;
Complexité : O(n)O(n)
Sélectionnez
fonction multiplie(ref p1,p2:polynome):polynome;
var p:polynome;
var i,j:entier;
début
init(p);
[Link]=[Link] + [Link];
pour i de 0 à [Link];
pour j de 0 à [Link] faire
[Link][i+j] = [Link][i+j] + [Link][i] * [Link][j];
finpour
finpour
retourner(p);
fin
finfonction
Complexité : O(n2)O(n2)
Soit :
P(x)=∑i=0naixiP(x)=∑i=0naixi
Q(x)=∑i=0n−pixiQ(x)=∑i=0n−pixi
Sélectionnez
fonction moins(ref p:polynome):polynome;
var i:entier;
var m:polynome;
début
[Link]=[Link];
pour i de 0 à [Link] faire
[Link][i] =- [Link][i];
finpour
retourner(m)
fin
finfonction
Complexité : O(n)O(n)
Sélectionnez
fonction decale(ref p:polynome;val n:entier):polynome
var i:entier;
var d:polynome;
début
[Link] = [Link] + n;
si n>0 then
pour i de 0 à n-1 FAIRE
[Link][i] = 0;
finpour
pour i de 0 à [Link] faire
[Link][i+n] = [Link][i];
finpour
sinon
pour i de -n à [Link] faire
[Link][i+n] = [Link][i];
finpour
finsi
retourner(d)
fin
finfonction
Complexité : O(n)O(n)
Sélectionnez
fonction deriv(ref p1:polynome):polynome;
var p:polynome;
var i:entier;
début
si [Link] == 0 alors
[Link] = 0;
[Link][0] = 0;
sinon
[Link] = [Link]-1;
pour i de 1 à [Link] faire
[Link][i-1] = [Link][i]
finpour;
finsi;
retourner(p)
fin
finfonction
Complexité : O(n)O(n)
Sélectionnez
fonction valeur(ref p:polynome;val x: réel):réel;
var f,s,i:réel;
début
s = 0;
f = 1;
pour i allant de 0 à [Link] faire
s = s + f * [Link][i];
f = f * x
finpour
retourner(s)
fin
finfonction
Complexité : O(n)O(n)
Sélectionnez
fonction integraleDefinie(ref p1:polynome;val x,y: réel):réel;
var p:polynome;
var s:réel;
var i:entier;
début
[Link] = [Link]+1;
[Link][0] = 0;
pour i allant de 0 à [Link] faire
[Link][i+1] = [Link][i] / (i+1);
finpour;
retourner(valeur(p,y) - valeur(p,x))
fin
finfonction
Même si en première approche, la complexité ne prend en compte que le nombre d'opérations (+, *), en seconde
analyse, les multiplications sont beaucoup plus coûteuses que les additions. Cela est pris en compte dans les
algorithmes ci-dessous.
L'algorithme énoncé au paragraphe précédent effectue 2n multiplications. Le schéma d'Horner d'un polynôme
permet d'effectuer n multiplications seulement. Le schéma d'Horner repose sur la propriété suivante :
P(x)=∑i=0naixiP(x)=∑i=0naixi
P(x)=A(x)x+a0P(x)=A(x)x+a0
P(x)=((…((anx+an−1)x+an−2)…)x+a1)x+a0P(x)=((…((anx+an−1)x+an−2)…)x+a1)x+a0
Sélectionnez
fonction valeur(ref p:polynome;val x: réel):réel;
var s:réel;
début
s = [Link][[Link]];
pour i allant de [Link]-1 à 0 pas de -1 faire
s = s * x + [Link][i]
finpour
retourner(s)
fin
finfonction
Une méthode « diviser pour régner » permet d'améliorer cet algorithme. Elle est basée sur l'égalité suivante :
(ay+b)(cy+d)=acy2+((a+b)(c+d)−ac−bd)y+bd(ay+b)(cy+d)=acy2+((a+b)(c+d)−ac−bd)y+bd
Cette égalité signifie, entre autres, que si deux polynômes sont de degré 1, il suffit de trois multiplications de
réels pour obtenir leur produit. Les quantités a, b, c, d, y étant quelconque, celles-ci peuvent être elles-mêmes
des polynômes. Soit le polynôme :
P(x)=∑i=0npixiP(x)=∑i=0npixi
P(x)=A(x)x1+⌊n/2⌋+B(x)P(x)=A(x)x1+⌊n/2⌋+B(x)
avec :
A(x)=∑i=0n−1−⌊n/2⌋pixiA(x)=∑i=0n−1−⌊n/2⌋pixi
B(x)=∑i=0⌊n/2⌋pixiB(x)=∑i=0⌊n/2⌋pixi
De même, on a :
Q(x)=∑i=0nqixiQ(x)=∑i=0nqixi
Q(x)=C(x)x1+⌊n/2⌋+D(x)Q(x)=C(x)x1+⌊n/2⌋+D(x)
C(x)=∑i=0n−1−⌊n/2⌋qixiC(x)=∑i=0n−1−⌊n/2⌋qixi
D(x)=∑i=0⌊n/2⌋qixiD(x)=∑i=0⌊n/2⌋qixi
De plus, on est amené à calculer des produits de polynômes de degré au plus n/2. Il s'agit donc d'une méthode :
« diviser pour régner ». La complexité en nombre de multiplications est alors O(nlog23)O(nlog23). Comme on
notera ci-dessous, l'algorithme est plus complexe à écrire, mais il est bien plus efficace aussi.
Un prétraitement permet de considérer des polynômes de même degré. Il faut donc une fonction « chapeau ».
On définit de plus deux autres fonctions utiles pour le calcul :
etend : si le degré de P(x)P(x) est supérieur au degré de Q(x)Q(x), alors Q(x)Q(x) est modifié
de manière à ce que les coefficients manquants soient à zéro, et les degrés des deux
polynômes deviennent ainsi égaux ;
tronque : permet d'initialiser un polynôme avec les premiers termes d'un autre polynôme (calcul
de B(x)B(x) et D(x)D(x)).
Sélectionnez
fonction multiplie(ref p,q:polynome):polynome;
var p1,p2:polynome;
fonction etend(ref p:polynome;val n:entier):polynome;
var i:entier;
var r:polynome;
début
r = decale(p,0);
[Link] = n;
pour i de [Link]+1 à n faire
[Link][i] = 0;
finpour;
retourner(r);
fin
finfonction
fonction tronque(ref p:polynome;val n:entier)
var c:polynome;
var i:entier;
début
[Link] = n;
pour i de 0 à n faire
[Link][i] = [Link][i];
finpour;
retourner(c);
fin
finfonction
fonction multirec(ref p,q:polynome)
var a,b,c,d:polynome;
var c0,c1,c2:réel;
var C0,C1,C2:réel;
var m:entier;
début
selon que
cas [Link] == 0
retourner(polynome([[Link][0]*[Link][0]]))
cas [Link] == 1
c0 = [Link][0] * [Link][0];
c2 = [Link][1] * [Link][1];
c1 = ([Link][0] + [Link][1]) * ([Link][0] + [Link][1]);
c1 = c1 - c0 - c2;
retourner(polynome([c0,c1,c2]));
autrement
m = partieEntière([Link]/2);
a = decale(p, -(m+1));
b = tronque(p, m);
c = decale(q, -(m+1));
d = tronque(q, m);
C2 = multirec(a, c);
C0 = multirec(b, d);
C1 = multirec(ajout(a,b), ajout(c,d));
C1 = ajout(C1, ajout(moins(C0), moins(C2)));
C1 = decale(C1, 1+m);
C2 = decale(C2, 2+2*m);
C0 = ajout(C0, ajout(C1,C2));
retourner(C0);
finselonque:
fin
finfonction;
début
si [Link] > [Link] alors
p1 = p;
p2 = etend(q, [Link]);
sinon
p1 = q;
p2 = etend(p, [Link]);
finsi;
retourner(multirec(p1,p2));
fin
finfonction
Définition 6.1. Une liste est une table d'association à clé unique telle que :
La complexité de l'accès à un élément par son pointeur est O(1)O(1). Si p est un pointeur vers un élément
alors contenu(p) est l'élément lui-même. Un pointeur qui n'adresse aucun élément a pour valeur NIL. On écrit
en EXALGO pour déclarer un pointeur :
Sélectionnez
nom_pointeur=^type_predefini;
Sélectionnez
type_liste=liste de type_predefini;
La manipulation des éléments de la liste dépend des fonctions définies comme s'exécutant en temps O(1)O(1).
Définition 6.2. Une liste est dite simplement chainée si les opérations suivantes s'effectuent en O(1)O(1) :
accès :
Sélectionnez
fonction premier(val L :type_liste) :^type_predefini;
fonction suivant(val L :type_liste;
val P :^type_predefini) :^type_predefini;
fonction listeVide(val L :type_liste) :booléen;
modification :
Sélectionnez
fonction créer liste(ref L :type_liste) :vide;
fonction insérerAprès(val x :type_prédéfini;
ref L :type_liste;
P :^type_predefini) :vide;
fonction insérerEnTete(val x :type_prédéfini;
ref L :type_liste) :vide;
fonction supprimerAprès(ref L :type_liste;val P :^type_predefini) :vide;
fonction supprimerEnTete(ref L :type_liste) :vide;
On écrira en EXALGO listeSC pour préciser qu'il s'agit d'une liste simplement chaînée.
Sélectionnez
fonction estDernier(ref L :listeSC de type_prédéfini;
ref P :^type_prédéfini) :booléen;
début
retourner(suivant(L,P) == NIL)
fin
finfonction
Sélectionnez
fonction chercher(ref L :listeSC de type_prédéfini;
ref E :type_prédéfini) :^type_predefini;
début
var p :^type_prédéfini;
début
si listeVide(L) alors
retourner(NIL)
sinon
p = premier(L);
tant que non(estDernier(L,p)) et (contenu(p) != e) faire
p = suivant(L,p);
fintantque
si (contenu(p) != e) alors
retourner(NIL)
sinon
retourner(p)
finsi
finsi
fin
finfonction
Complexité :
minimum : O(1)O(1) ;
maximum : O(n)O(n).
Sélectionnez
fonction trouverDernier(ref L :listeSC de type_prédéfini) :^type_predefini;
var p :^type_prédéfini;
début
si listeVide(L) alors
retourner(NIL)
sinon
p = premier(L);
tant que non(estDernier(L,P)) faire
p = suivant(L,p);
fintantque
retourner(p)
finsi
fin
finfonction
Complexité : O(n)O(n).
Définition 6.3. Une liste doublement chaînée est une liste pour laquelle les opérations en
temps O(1)O(1) sont celles des listes simplement chaînées auxquelles on ajoute les fonctions d'accès.
Sélectionnez
fonction dernier(val L :type_liste) :^type_predefini;fonction précédent(val
L :type_liste;
val P :^type_predefini) :^type_predefini;
On écrira en EXALGO listeDC pour préciser qu'il s'agit d'une liste doublement chaînée.
Sélectionnez
fonction supprimer(ref L :listeDC de type_predefini;
val e : type_prédéfini) :booléen;
var p,prec,suiv :^type_prédéfini;
début
p = chercher(L,e);
si p == NIL alors
retourner(FAUX)
sinon
si estPremier(L,p) alors
supprimerEnTete(L)
sinon
prec = précédent(L,p);
supprimerApres(L,prec);
finsi
retourner(VRAI)
finsi
fin
finfonction
Complexité :
minimum : O(1)O(1) ;
maximum : O(n)O(n).
Sélectionnez
fonction taille(val L :type_liste) :entier;
var p :^type_prédéfini;
var t :entier;
début
si listeVide(L) alors
retourner(0)
sinon
retourner(1 + taille( suivant(L, premier(L)) ))
finsi
fin
finfonction
Complexité : O(n)O(n).
Sélectionnez
fonction insertionTrie(ref L :listeDC de type_prédéfini;
val e : type_prédéfini) :vide;
var p :^type_prédéfini;
début
si listeVide(L) alors
insererEnTete(e,L)
sinon
si contenu(premier(L))>e alors
insererEnTete(e,L)
sinon
insererTrie(suivant(L,premier(L)), e)
finsi
finsi
fin
finfonction
On considérera dans tout ce chapitre que l'on a des valeurs qui correspondent à un caractère.
Pour certaines structures de données, l'ensemble des langages de programmation proposent une traduction
immédiate. Pour d'autres, il n'existe pas de traduction immédiate. Il faut alors définir explicitement l'algorithme
de chacune des primitives.
Exemple - les listes. On doit définir le stockage de la liste, et en fonction de ce stockage comment s'effectue
par exemple l'adjonction.
L'implémentation doit respecter la complexité des primitives à part celle d'initialisation (celle-ci ne s'exécutera
qu'une fois).
Exemple - les listes. On utilise souvent les fonctions ajouter et supprimer, mais une seule fois creerListe.
Ici nous allons choisir de ranger les éléments dans un tableau « suffisamment grand ». Chaque élément du
tableau est une paire (valeurElement, pointeurSuivant). Un pointeur est la valeur d'un index du tableau ; ainsi l'accès
au suivant est en complexité O(1)O(1). La zone de stockage peut donc être décrite par :
Sélectionnez
elementListe = structure
valeur :car;
suivant :entier;
finstructure;
stockListe = tableau[1..tailleStock] d'elementListe;
La valeur du pointeur (champ suivant) est donc un entier compris entre 0 et tailleStock. La valeur 0 correspondant
à l'absence d'élément suivant. Le premier élément doit être accessible en O(1)O(1), il faut donc conserver son
index. Si la liste est vide, par convention, l'index du premier sera 0. On peut donc représenter une liste par la
structure suivante :
Sélectionnez
listeSC_Car = structure
tailleStock :entier;
premier :entier;
vListe :stockListe;
finstructure;
Le tableau de stockage étant grand, mais pas illimité, il faudra prévoir que l'espace de stockage puisse être
saturé.
Sélectionnez
fonction premier(val L :listeSC_Car) :entier;
début
retourner [Link];
fin;
finfonction
Pour ajouter un élément, il faut pouvoir trouver un élément « libre » dans le tableau. Une première solution
consiste à marquer les éléments libres du tableau (par exemple champ suivant de l'élément a pour valeur -1).
Dans ce cas, il faudra parcourir le tableau (complexité O(n/2)O(n/2) en moyenne). Par suite, la
primitive insérerAprès ne sera plus en complexité O(1)O(1) puisqu'il faudra d'abord trouver un élément libre. Une
solution compatible avec la complexité des primitives consiste à gérer cet espace de stockage en constituant la
liste des cellules libres. On modifie donc en conséquence la description de listeSC_Car :
Sélectionnez
listeSC_Car = structure
tailleStock :entier;
premier :entier;
premierLibre :entier;
vListe :stockListe;
finstructure;
Par convention, l'espace de stockage sera saturé lorsque l'index premierLibre vaut 0 (la liste des cellules libres est
vide). On définit donc la fonction de test :
Sélectionnez
fonction listeLibreVide(val L :listeSC_Car) :booléen;
début
retourner [Link] == 0;
fin
finfonction
Sélectionnez
fonction prendreCellule(ref L :listeSC_Car,ref tete :entier) :entier;
var nouv :entier;
début
nouv = tete;
tete = suivant(L,nouv);
retourner nouv;
fin
finfonction
Sélectionnez
fonction créer_liste(ref L :listeSC_Car;val tailleMax :entier) :vide;
var i :entier;
début
[Link] = tailleMax;
[Link] = 0;
[Link] = 1;
pour i allant de 1 à [Link]-1 faire
[Link][i].suivant = i+1;
finpour
[Link][tailleStock].suivant = 0;
fin
finfonction