Voir la source de Générateur congruentiel linéaire

{{Sources|date=septembre 2014}}
Un '''générateur congruentiel linéaire''' est un [[générateur de nombres pseudo-aléatoires]] dont l'algorithme, introduit en 1948 par [[Derrick Lehmer]], sous une forme réduite, pour produire des nombres aléatoires, est basé sur des [[Congruence sur les entiers|congruences]] et une [[fonction affine]].

== Définition ==
Les nombres pseudo aléatoires forment une suite dont chaque terme dépend du précédent, selon la formule :
:<math>X_{n+1}=(a\,X_n+c)\mod m</math>
où ''{{math|a}}'' est le ''multiplicateur'', ''{{math|c}}'' l'''incrément'' et ''{{math|m}}'' le ''module''.

Le terme initial <math>X_0~</math> est appelé la graine (''seed'' en anglais). C'est elle qui va permettre de générer une suite apparemment aléatoire. Pour chaque graine, on aura une nouvelle suite. Cependant, il est possible que certaines graines permettent d'obtenir une suite plus aléatoire que d'autres.

Du fait de l'opération ''mod'' (reste dans la [[division euclidienne]] de <math>(a\,X_n+c)</math> par ''{{math|m}}''), les termes de cette suite sont compris entre 0 et <math>(m-1)</math>. De plus, comme chaque terme dépend entièrement du précédent, si un nombre apparaît une deuxième fois, toute la suite se reproduit à partir de ce nombre. Or le nombre de valeurs que le nombre peut prendre étant fini (égal à ''{{math|m}}''), la suite est amenée à se répéter au bout d'un certain temps. On dit qu'elle est [[fonction périodique|ultimement périodique]].

== Du bon choix des paramètres ''a'', ''c'' et ''m'' ==

Dans ce type d'algorithme, la qualité du générateur va entièrement dépendre du choix de ''a'', ''c'' et ''m'' car on ne peut pas se satisfaire d'un générateur qui produit des suites plus ou moins aléatoires, sans aucune raison apparente, selon le choix de X<sub>0</sub>.

=== Un mauvais exemple ===
Choisir les valeurs ''a'', ''c'' et ''m'' « au hasard » n'est pas une bonne idée. Prenons un exemple, soit le générateur congruentiel linéaire employant les valeurs ''a'' = 25, ''c'' = 16, ''m'' = 256, nous obtenons :
* avec ''X<sub>0</sub> = 125'', la suite : {{Vert|125, 69, 205, 21, 29, 229, 109, 181, 189, 133, 13, 85, 93, 37, 173, 245, 253, 197, 77, 149, ...}}
* avec ''X<sub>0</sub> = 96'', la suite : {{Vert|96, 112, 0, 16, 160, 176, 64, 80, 224, 240, 128, 144, 32, 48, 192, 208, '''96, 112, 0, 16, ...''', }}
* avec ''X<sub>0</sub> = 50'', la suite : {{Vert|50, 242, 178, 114, '''50, 242, 178, 114, 50, 242, ...'''}}
* avec ''X<sub>0</sub> = 10'', la suite : {{Vert|10, '''10, 10, 10, ...'''}}

Il est clair que ces suites ne peuvent être considérées comme aléatoires.<br />
On voit donc bien que l'on doit choisir avec précaution les paramètres du générateur si l'on espère obtenir des nombres qui s'approchent de l'aléa parfait. Il faut alors se demander comment choisir ''a'', ''c'' et ''m'' convenablement.

=== Choix du module ===

Les générateurs congruentiels font intervenir un calcul modulo m, et donc a priori une division euclidienne, ce qui peut avoir un coût de calcul important dans le cadre d'une utilisation fréquente du générateur. La solution la plus simple est d'utiliser un module de type m = {{2e}}. En effet, les ordinateurs calculant naturellement en base binaire, un tel choix est tout à fait transparent pour eux, ce qui rend inutile une division euclidienne.

Toutefois, ce choix présente une limite importante : les bits dits de poids faible (les bits les plus à droite) sont beaucoup moins aléatoires que ceux de la partie de poids fort (les bits les plus à gauche). En effet, si ''d'' est un diviseur de ''m'', alors la suite ''Y<sub>n</sub>'', telle que :
<center><math>Y_n = X_n \mod d</math></center>
satisfait à la suite congruentielle linéaire :
<center><math>Y_{n+1} = (a \cdot Y_n + c) \mod d</math></center>
en particulier, avec d = 2<sup>k</sup>, pour ''k'' fixé, compris entre 1 et ''e'', on voit que les ''k'' chiffres de poids faible, ont une période maximale de 2<sup>k</sup>, évidemment inférieure à m.<br />
Pour remédier à ce problème, on peut ne garder que les bits de poids fort, c'est-à-dire garder les bits les plus à gauche du nombre obtenu. Si l'on tronque les k derniers bits, on aura alors un générateur pseudo aléatoire de nombres compris entre 0 et <math>2^{e-k}-1</math>.

Une autre solution consiste à prendre un module du type <math>m = 2^e \pm 1</math>, ce qui permet encore d'éviter une division euclidienne par une astuce (cf. ''The Art Of Computer Programming'' de D. E. Knuth).

=== Choix du multiplicateur et de l'incrément ===

Afin de pouvoir choisir une graine <math>X_0</math> sans contraintes entre <math>0</math> et <math>m - 1</math>, il faut chercher à maximiser la période du générateur. Or il se trouve que les valeurs de <math>a</math> et <math>c</math> sont connues, ce qui permet d'obtenir une période maximale (égale à <math>m</math>).

Si <math> c \not= 0 </math> la période d’un générateur congruentiel linéaire est '''maximale''' si et seulement si :
:# <math>c</math> est premier avec <math>m</math>. <math>pgcd(c,m) = 1</math>.
:# Pour chaque nombre premier <math>p</math> divisant <math>m</math>, <math>(a-1)</math> est un multiple de <math>p</math>.
:# <math>(a-1)</math> est un multiple de 4 si <math>m</math> en est un.

On remarquera que l'on possède un théorème d'équivalence (si et seulement si) : la période est maximale pour toutes les valeurs qui possèdent les propriétés du théorème et seulement pour celles-là.

On possède, de plus, de conditions suffisantes dans le cas particulier où ''c'' est nul.
Ainsi, pour c nul la période d'un générateur congruentiel linéaire est '''maximale''' si{{Citation nécessaire}} :
:# <math> m </math> est premier.
:# <math> a^{m-1} - 1 </math> est un multiple de <math>m</math>.
:# <math> a^j - 1 </math> n'est pas divisible par <math>m</math>, pour <math> j = 1,...,m-2 </math>

Attention: ce n'est pas une équivalence !

== Le potentiel ==

Le potentiel est une notion qui permet d'écarter certains générateurs insuffisamment aléatoires. Il est toujours défini que pour une suite possédant une période maximale (et ce, d'après la deuxième propriété énoncée plus haut). On ne peut pas toujours le définir sur les autres générateurs, et il existe de bons générateurs sur lesquels le potentiel n'est pas défini.

Le potentiel <math>s</math> d’une suite possédant une période maximale est défini comme le plus petit entier tel que :
<center><math>(a-1)^s \equiv 0 \mod m</math></center>

Pour une suite de période maximale, on peut prendre <math>X_0 = 0</math>, et donc :
<center><math>X_n = \frac {a^n - 1} {a - 1} \cdot c \mod m</math></center>
ainsi, avec <math> a^n = ((a-1) + 1)^n </math>:
<center><math>X_n = c \cdot \left ( n + {n \choose 2} \cdot (a-1) + \ldots + {n \choose s} \cdot (a-1)^{s-1} \right ) \mod m</math></center>
et un potentiel inférieur ou égal à trois, apparaît immédiatement comme trop faible car la suite n’est pas suffisamment aléatoire. En fait un potentiel de 5 ou plus est conseillé.

Rappelons cependant que la notion de potentiel est uniquement là pour écarter quelques mauvais générateurs. Un générateur possédant un potentiel supérieur à 5 ne peut en aucun cas être d’emblée considéré comme bon, il doit d’abord subir quelques tests.

== Tests statistiques ==
[[Fichier:Lcg 3d.gif|vignette|Générateur congruentiel linéaire en 3D]]
Comme tous les générateurs de nombres pseudo-aléatoires, le générateur congruentiel linéaire doit être soumis à une série de tests avant d’être homologué.

Le plus connu d’entre eux, le test de fréquence n'est généralement pas un problème pour un générateur congruentiel linéaire de période maximale, mais il existe bien d'autres tests : Test des séries, Gap test, Run test, ... ; <br />
l'un des plus discriminant, car aucun générateur congruentiel linéaire démontré mauvais ne l’a réussi, est le [[Test spectral]].

== Exemples ==
{{À vérifier|thème=informatique|date=décembre 2016}}
{| class=wikitable
!Algorithme
!''a''
!''c''
!''m''  
!Remarques
|-
|[[RANDU]]
|65539
|0
|2<sup>31</sup>
|Biaisé et fortement déconseillé
|-{{ligne grise}}
|générateur de [[Robert Sedgewick]]
|31415821
|1
|10<sup>8</sup>
|Intéressant mais déconseillé (essayer avec <math>X_0 = 0~</math>)
|-
|[[Standard minimal]]
|16807
|0
|2<sup>31</sup>-1
|
|}
* RANDU défini par <math>X_{n+1} = 65539 X_n \mod 2^{31}</math>. Implanté sur les [[IBM]] System/370 et {{Référence souhaitée|réputé mauvais par tous ceux qui ont à l'utiliser sérieusement}}.
* Le générateur de [[Turbo Pascal]] défini par <math>X_{n+1} =(129 X_n + 907633385) \mod 2^{32}</math>. On sait dans 1 cas sur 129 que <math>X_{n+1} > X_n</math>, ce qui est un biais énorme pour une application scientifique. En plus comme 129 est de la forme 2^n + 1, les nombres produits ne sont pas si « aléatoires que ça » [cf. Knuth p.23].
* Le générateur d'UNIX (dont le comité ANSI C a heureusement recommandé de ne conserver que les 16 bits de poids fort) :
<math>X_{n+1} = (1103515245 X_n + 12345) \mod 2^{31}</math>. Même la documentation du Berkeley 4.2 admet que les bits de poids faible des nombres produits ne sont pas vraiment aléatoires.

=== Autres exemples ===
Knuth recense dans « Seminumerical Algorithms » un certain nombre de générateurs congruentiels linéaires de qualité. S'il n'est pas possible d'utiliser le Standard Minimal, on pourra essayer par exemple :
* <math>X_{n+1} = ( 137 X_n + 187 ) \mod 2^8</math> (une suite cyclique de 256 nombres, est-ce encore du hasard ?)
* <math>X_{n+1} = ( 1664525 X_n + 1013904223 ) \mod 2^{32}</math> : générateur de Knuth & Lewis
* <math>X_{n+1} = 69069 X_n \mod 2^{32}</math> : générateur de Marsaglia au multiplicande remarquable, utilisé sur le [[VAX]]
* <math>X_{n+1} = 31167285 X_n \mod 2^{48}</math> : générateur de Lavaux & Jenssens
* <math>X_{n+1} = 6 364 136 223 846 793 005 X_n \mod 2^{64}</math> : générateur de Haynes

Dans tous ces cas, comme le modulant est une [[Puissance de deux|puissance de 2]], les bits de poids faible des nombres produits ne sont pas eux très « aléatoires ». Il est par conséquent préférable de ne conserver que les poids forts (16 bits par exemple), comme le fait le générateur de [[Borland]] C++ par exemple :

* <math>X_{n+1} = (22 695 477 X_n + 1) \bmod 2^{32}</math> et on sort <math>(X_{n+1} \gg 16)\bmod 2^{15}</math>
Ce faisant, la période du générateur reste la même mais on ne va produire que des nombres compris entre 0 et 65535, chacun se retrouvant 32768 fois dans la suite.

Un autre exemple serait de mixer la sortie d'un générateur congruentiel linéaire et d'un registre à décalage à rétroaction linéaire<ref>{{lien web |langue=en |titre=Alexpukall/prng64 : pukall prng 64 : mix the output of an lfsr 64 and an  lcg 64 (c language or 6502 assembly) |url=https://github.com/AlexPukall/PRNG64 |site=GitHub |consulté le=12-11-2023}}.</ref>.

== Références ==
{{Références}}

== Voir aussi ==
=== Bibliographie ===
{{en}} [[Donald Knuth|Donald E. Knuth]], ''[[The Art of Computer Programming]]'', vol. 2 : ''Seminumerical Algorithms'', {{3e}} éd., [[Addison-Wesley]], Boston, 1998

=== Articles connexes ===
* [[Générateur de nombres aléatoires]]
* {{Lien|trad=Lehmer random number generator|Générateur de Lehmer}}
{{Palette|Générateur de nombres pseudo-aléatoires}}

{{Portail|informatique}}

[[Catégorie:Nombre aléatoire]]