Si ce bandeau n'est plus pertinent, retirez-le. Cliquez ici pour en savoir plus.
Les informations figurant dans cet article ou cette section doivent être reliées aux sources mentionnées dans les sections «Bibliographie», «Sources» ou «Liens externes» ().
Le crible d'Ératosthène consiste à trouver les valeurs pouvant se réduire à la forme quadratique binaire réduite x⋅y (produit de deux entiers strictement supérieurs à 1). Le crible d'Atkin consiste lui à dénombrer les valeurs d'une forme quadratique binaire non-réduite. Par exemple, un entier sans facteur carrép qui est congru à 1 modulo 4 est premier si et seulement si l'équation
4x2 + y2 = p
admet un nombre impair de solutions positives (x, y).
Algorithme
Quelques variables sont nécessaires:
le tableau des nombres premiers que l'on a trouvés, que l'on initialise avec les nombres premiers inférieurs à soixante;
trois tableaux contenant les nombres qu'il reste à tester, selon le reste dans la division euclidienne par soixante; ils sont initialisés avec tous les entiers entre soixante et le dernier nombre qui nous intéresse, tel que le reste modulo soixante est dans l'une de ces trois listes (on ne considère que les restes premiers avec 60):
reste valant 1, 13, 17, 29, 37, 41, 49, ou 53,
reste valant 7, 19, 31, ou 43,
reste valant 11, 23, 47, ou 59.
Après, pour chaque nombre n restant:
dans la première classe, on compte le nombre de couples x > 0 et y > 0 solutions de 4x² + y² = n;
dans la deuxième classe, on compte le nombre de couples x > 0 et y > 0 solutions de 3x² + y² = n;
dans la troisième classe, on compte le nombre de couples x > 0 et y > 0, avec x > y, solutions de 3x² – y² = n;
si le compte est pair, on supprime n de la liste.
Ensuite, on applique un crible d'Ératosthène modifié: pour chaque nombre premier p supérieur à 7, on supprime tous les multiples de p².
Enfin, on place les nombres restants dans la liste des nombres premiers.
Complexité
Avec cet algorithme, la recherche des nombres premiers inférieurs à N a une complexité en temps en o(N) et une complexité en mémoire en O(Nα) avec α < 1.