-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathShor's Algorithm.py
More file actions
76 lines (50 loc) · 1.35 KB
/
Shor's Algorithm.py
File metadata and controls
76 lines (50 loc) · 1.35 KB
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
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
from random import randrange
from math import floor
def euclid_algo(x, y):
if x < y:
return euclid_algo(y, x)
while y != 0:
(x, y) = (y, x % y)
return x
def verificaIgualdade(num, lista):
for number in lista:
if number == num:
return True
return False
def contaFrequencia(number, lista):
count = 0
for num in lista:
count += 1
if num == number:
return count
def descobrirFrequencia(rand, num, grandeza = 15):
restosLista = []
for x in range(grandeza):
resto = (rand**x)%num
if not verificaIgualdade(resto, restosLista):
restosLista.append(resto)
else:
return contaFrequencia(resto, restosLista)
return 0
def decodificar(num, grandeza):
freq = 0
factors = []
while freq == 0 or not (freq%2 == 0):
rand = randrange(2, num)
freq = descobrirFrequencia(rand, num)
guess = rand**(freq/2)
factors.append(euclid_algo(guess+1,num))
factors.append(euclid_algo(guess-1,num))
if factors[0] == 1:
factors[0] = num/factors[1]
if factors[1] == 1:
factors[1] = num/factors[0]
if factors[0] == num or factors[1] == num:
factors = decodificar(num, grandeza)
return factors
def main():
num = int(input("The chosen number: "))
grandeza = int(input("Number of trys: "))
print(decodificar(num, grandeza))
if __name__ == '__main__':
main()