# Cryptomath Module # https://www.nostarch.com/crackingcodes (BSD Licensed) from typing import Optional def gcd(a: int, b: int) -> int: """ Euclid's Algorithm Return the Greatest Common Divisor of a and b """ while a != 0: a, b = b % a, a return b def findModInverse(a: int, m: int) -> Optional[int]: # Return the modular inverse of a % m, which is # the number x such that a*x % m = 1 if gcd(a, m) != 1: # Ah dynamic typing! return None # No single mod inverse exists if a & m aren't relatively prime. # Calculate using the Extended Euclidean Algorithm: u1, u2, u3 = 1, 0, a v1, v2, v3 = 0, 1, m while v3 != 0: # Note that // is the integer division operator q = u3 // v3 v1, v2, v3, u1, u2, u3 = (u1 - q * v1), (u2 - q * v2), (u3 - q * v3), v1, v2, v3 return u1 % m