Quelques variables sont nécessaires :
- Le tableau des nombres premiers que l'on a trouvé, que l'on initialise avec les nombres premiers inférieurs à soixante.
- Trois tableaux contenant les nombres qu'il reste à tester, et qui ont certains restes dans la division 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 (notons que les nombres avec certains restes sont ignorés) :
- 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 solutions de 3x - y = n.
- Si le compte est impair, 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 restant dans la liste des nombres premiers.