Discussion :

[laurent_ott]

Chers membres du club,

J'ai le plaisir de vous présenter ce tutoriel :

Comprendre la méthode de factorisation du Crible Quadratique
Une invention de Carl Pomerance

Cet article vous permet de comprendre la méthode de factorisation du crible quadratique. Vous trouverez dans le fichier joint les codes source en VBA du crible quadratique ainsi que d'autres fonctions utilisées pour la factorisation : le test de primalité Miller-Rabin, le crible d'Ératosthène, la factorisation RhoPollard, l'algorithme Tonelli-Shanks, mais aussi les algorithmes pour les opérations sur les grands nombres (additions, soustractions, multiplications, puissances, divisions, modulos, racines carrées…).

Bonne lecture

Retrouvez les meilleurs cours et tutoriels pour apprendre l'algorithmique.

[valentin03]

A quoi ça sert concrètement ?

[laurent_ott]

Bonjour.
Si la question est à quoi sert le crible quadratique, la réponse est : à factoriser "rapidement" des grands nombres.
Si la question est à quoi sert cette documentation, la réponse est : à expliquer "simplement" cette méthode de factorisation qui est une référence dans ce domaine.
Si la question est à quoi ça sert de savoir comment marche le crible quadratique, alors on s'éloigne peut-être de l'algorithmique et l'on s'approche de la philosophie.

Plus sérieusement, la question mérite en effet d'être posée. J'avoue ne pas avoir la réponse.
Et pour celles et ceux qui (se) demanderont à Qui ça sert, je peux répondre tout de suite : "je ne sais pas, en tout cas pas à moi."
Bonne lecture pour les plus curieux d'entre vous.

[valentin03]

Qui factorise et pourquoi ?
Tu dis ne pas avoir la réponse mais peut-être que quelqu'un ici l'a. Je suis curieux pathologique.

[aghrsa]

Bonjour,

Je n'ai pas compris dans l'exemple du début vous vous arrêtiez une fois arrivé à 51, ou dans le deuxième exemple, pourquoi à 60261, pourriez-vous clarifier ce détail ?

Merci beaucoup

[laurent_ott]

Le premier exemple est très simpliste et ne tient pas compte des améliorations possibles. La règle du pivot de Gauss dit que si l'on a au moins autant de lignes que de colonne une solution peut être trouvée.
60261 est la 7eme ligne sur 7 colonnes (-1 est exclu car non utilisé ici) donc on cherche une solution. Si pas trouvée on ajoutera une ligne. Ainsi de suite.
La recherche étant chronophage, c'est inutile de la lancer avant d'avoir le minimum de lignes requises, sauf à avoir de la chance.
Cordialement

[aghrsa]

Ah d'accord, on détermine les facteurs à l'aide de la friabilité et de Legendre, puis on arrête de chercher des équations une fois qu'on en a obtenu autant que l'on a de facteurs ?

[laurent_ott]

C'est comme ça que je l'ai compris, et programmé.

[jefflg]

Bonjour et merci pour votre article "Comprendre la méthode de factorisation du Crible Quadratique. Il répond d’une façon claire à des questions que je me posais.

Par contre, je me suis surpris à trouver une simplification, étonnamment simple(équation du second degré), à cette méthode : Il suffit de remplacer le calcul des X²- N par la fonction :

F(n) = n²-2(X-1) n+ ((X²-N) -2X+1) avec n entier naturel
On remarque que si on calcule le discriminant (simplifié car b pair)= b'²-ac = ((X-1)²-((X²-N) -2X+1) = N
Normal, sinon on pourrait calculer directement les facteurs premiers.

Exemple pour N = 48206621

X=6944 (racine de N arrondi à l’entier supérieur)

F(n) = n² + 2(6944-1)n+(12515-2*6944+1) = n²+13886n –1372

F(1) = 1+13886-1372= 12515
F(2)=4+2*13886-1372= 26404

Etc....

Cette fonction facilite le calcul du crible quadratique :

Il suffit de calculer la fonction, modulo p (avec p=7,11,19,23,..)- On notera au passage que le test de Legendre sert à éliminer les nombres premiers qui n'admettent pas de racine carrée au discriminant dans p


exemple p= 11 f(n)mod11 = n²+ 4n +3 d'où b'²-ac = 4-3 =1 d'où les deux racines
n1 = -b'+1/a= -2+1= -1= 10 mod11
et n2= -2-1=-3 = 8 mod 11

p= 7 f(n) = n²+ 5n racine n1 = 0 et n= -5 = 2 mod 7

P=19 f(n) = n²+16n+15 d'où b'²-ac = 64-15 = 49 n1= -8-7= -15 = 4 mod 19 et n2= -1 = 18.

Calculer le crible quadratique revient donc à résoudre l'équation du second degré de la f(n) mod(p)=0.

Peut-être que l'utilisation de cette fonction peut permettre un gain de calcul conséquent.

Cordialement

[jefflg]

P.S. Lire F(n) = n²+2(X-1) n+ ((X²-N) -2X+1) au lieu de F(n) = n²-2(X-1) n+ ((X²-N) -2X+1)

[jefflg]

Exemple de la cible quadratique obtenu avec N=48206621 en pdf

[vincent Boulou]

Bonsoir.
J'ai lu attentivement les 2 articles de MR OTT concernant la décomposition de grands nombres en facteurs premiers.
J'ai donc commencé par l'algorithme ECM (plus facile) que j'ai porté en langage assembleur en suivant la démarche proposée par Laurent OTT.
Tout va bien pour des petits nombres mais mon programme n'arrive pas à factoriser le nombre 3000004222000029407 alors que l'algorithme rho y arrive en quelques secondes.
Si quelqu'un à implanté cet algorithme dans un autre langage ou si Mr OTT peut tester en visual basic pour me dire si la factorisation réussi et en combien de passages. Car je pense à une erreur dans mon programme à partir d'une certaine taille de nombre.
Merci d'avance.

[laurent_ott]

Bonjour.
Je confirme, votre chiffre peut être décomposé avec l'algorithme ECM, en tout cas ça marche chez moi.
Il faut faire plusieurs passages. Difficile à dire combien car à chaque passage les variables de départ sont aléatoires. Donc non reproductible.
D'après mes rapides essais, une cinquantaine de passages.
Mais c'est parfois moins, parfois plus.
Il faudrait peut être programmer votre recherche en lui indiquant le temps que vous souhaitez consacrer à l'algorithme, et donc le laisser tourner autant de fois que ce délai n'est pas dépassé. Ou comme j'ai fait, lui indiquer un nombre maximum de passages avant de laisser tomber.
Ce qui revient au même.
Bonne programmation.
Laurent OTT.

[vincent Boulou]

Merci de votre réponse.
Je vais donc chercher l'erreur dans mon programme.

[tbc92]

Ton programme plante pour 3000004222000029407
Est ce qu'il donne des résultats corrects pour des nombres un peu plus grands ?
Dans beaucoup de langages, le traitement des grands nombres est très problématique. Les langages sont bâtis en disant : je sais traiter tous les nombres jusqu'à xxx, mais pas au-delà. Et quand on travaille sur des nombres supérieurs, on a des comportements très étranges. (prévisibles pour les spécialistes, mais étranges pour les novices)

[vincent Boulou]

Bonjour.
En effet c'est bien un problème de dépassement de capacité.
Je n'avais pas fait attention qu'un nombre modulo N peut être très grand quand N devient très grand.
Merci de vos messages.