Bonjour tout le monde;

Je suis entrain de coder ma formule sur la décomposition d'un nombre en produit de facteurs premiers;

Pour l'instant le code marche avec les petites valeurs sauf à ce niveau j'ai un petit souci pour la lecture clavier qui ne répond pas; donc je suis toujours obligé d'initialiser f à décomposer:
Code : Sélectionner tout - Visualiser dans une fenêtre à part
f: int = int(input("Veuillez saisir le nombre entier à décomposer : "))
Bon; il y a aussi le reste de code qui ne fonctionne pas si le nombre à décomposer n'est pas un carré parfait;
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
if __name__ == '__main__': 
 
   tabpremier=[]
   tabcomposer=[]
# On  demande à l'utilisateur à saisir le nombre  qu'on convertit immédiatement en entier s'il saisie un chaine
 
   #f: int = int(input("Veuillez saisir le nombre entier à decomposer : "))
 
   f=21
 
   if miller_rabin(f, 50)==True:
         print("ce nombre est premier")
   else:
 
           print("La decomposition commence!")
           k1,k2=decompose(f)
           f1=2*k1+1
           f2=2*k2+1
           print(f"f={f1}*{f2}")
           if miller_rabin(f1,50)==True:              
              tabpremier.append(f1)
           else:              
              tabcomposer.append(f1)
 
           if miller_rabin(f2,50)==True:
              tabpremier.append(f2)
           else:
              tabcomposer.append(f2)
 
   while len(tabcomposer) != 0:
         for x in tabcomposer:
            k1,k2=decompose(x)
            f1=2*k1+1
            f2=2*k2+1
            tabcomposer.remove(x)
 
            if miller_rabin(f1,50)==True:
                tabpremier.append(f1)               
            else:
                tabcomposer.append(f1)  
 
            if miller_rabin(f2,50)==True:
               tabpremier.append(f2)              
            else:
                 tabcomposer.append(f2)
 
print("Resultat:")
 
print(tabpremier)
l'erreur est au niveau de la fonction miller-rabbin:
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
def miller_rabin(n, k):
 
    # Implementation uses the Miller-Rabin Primality Test
    # The optimal number of rounds for this test is 40
    # See http://stackoverflow.com/questions/6325576/how-many-iterations-of-rabin-miller-should-i-use-for-cryptographic-safe-primes
    # for justification
 
    # If number is even, it's a composite number
 
    if n == 2:
        return True
 
    if n % 2 == 0:
        return False
 
    r, s = 0, n - 1
    while s % 2 == 0:
        r += 1
        s //= 2
    for _ in range(k):
        a = random.randrange(2, n - 1)
        x = pow(a, s, n)
        if x == 1 or x == n - 1:
            continue
        for _ in range(r - 1):
            x = pow(x, 2, n)
            if x == n - 1:
                break
        else:
            return False
    return True
Erreur: f=21

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
La decomposition commence!
delta: 64
racine= 8.0
racin delta is not decimal; k1 et k2: 1 3
f=3*7
Traceback (most recent call last):
  File "D:\Users\DAFFE\AppData\Local\Programs\Python\Python36\Doc\decomposition", line 147, in <module>
    if miller_rabin(f1,50)==True:              
  File "D:\Users\DAFFE\AppData\Local\Programs\Python\Python36\Doc\decomposition", line 116, in miller_rabin
    a = random.randrange(2, n - 1)
  File "D:\Users\DAFFE\AppData\Local\Programs\Python\Python36\lib\random.py", line 198, in randrange
    raise ValueError("empty range for randrange() (%d,%d, %d)" % (istart, istop, width))
ValueError: empty range for randrange() (2,2, 0)
[Finished in 731ms]
avec f=25:

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
 
La decomposition commence!
delta: 0
f=5*5
Resultat:
[5, 5]
[]
[Finished in 1.4s]
le dernier problème est la conversion du grand nombre;
avec f=22112825529529666435281085255026230927612089502470015394413748319128822941402001986512729726569746599085900330031400051170742204560859276357953757185954298838958709229238491006703034124620545784566413664540684214361293017694020846391065875914794251435144458199

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
La decomposition commence!
delta: 4026587744844483487536318615131830859270168820589264937687931746069624082190901751073970800735693961268773532158456941497347383268550394184237379432857564428082918697714124782611156420627755895559074161507012059201937483627100670423971972472547282768135772469857150537429931278057182940716774200363536136141887673102130914062979710654116235678493238795584989353758044570014034679663088386688017502671153712917960234419471976554033363387117536910299216717624957837239475202332338138290772
Traceback (most recent call last):
  File "D:\Users\DAFFE\AppData\Local\Programs\Python\Python36\Doc\decomposition", line 143, in <module>
    k1,k2=decompose(f)
  File "D:\Users\DAFFE\AppData\Local\Programs\Python\Python36\Doc\decomposition", line 58, in decompose
    racine_delta=math.sqrt(delta)       
OverflowError: int too large to convert to float
[Finished in 1.4s]
c'est vraiment une case tête, trop mesquin
merci d'avance pour votre aide