Guest User

Untitled

a guest
Jan 16th, 2018
85
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.72 KB | None | 0 0
  1. import numpy
  2. import math
  3. import random
  4.  
  5. def chk_cht(x):
  6. if x % 2 != 0:
  7. return x + 1
  8. return x
  9.  
  10. def Step1(q, t):
  11. ksi = random.normalvariate(0, 1)
  12. N = (2 ** (t-1)) / (q) + (2 ** (t-1) * ksi) / q
  13. N = chk_cht(N)
  14. u = 0
  15. p = Step3(q, t, N, u)
  16.  
  17. def Step3(q, t, N, u):
  18. p = (N + u) * q + 1
  19. if p > 2 ** t:
  20. p = Step1(q, t)
  21. else:
  22. p = Step5(q, t, N, u, p)
  23.  
  24. def Step5(q, t, N, u, p):
  25. temp1 = 2 ** (N - 1)
  26. temp2 = 2 ** (N + u)
  27. if (temp1 == numpy.mod(1, N)) and (temp2 != numpy.mod(1, N)):
  28. return p
  29. else:
  30. u = u + 2
  31. p = Step3(q, t, N, u)
  32.  
  33. t = 4
  34. q = 3
  35. p = 0
  36. u = 0
  37.  
  38. print(Step1(q, t))
  39.  
  40. RecursionError: maximum recursion depth exceeded in comparison
Add Comment
Please, Sign In to add comment