period_finding

for i in [1..20]:
    p = next_prime(randint(2^1023,2^1024))
    q = next_prime(randint(2^1023, 2^1024))
    l = lcm(p-1, q-1)
    phi = (p-1) * (q-1)
    print(phi / l)