pub fn mod_log(a: u64, b: u64, m: u64) -> Option<u64>
aˣ = b (mod m)を満たすxを求める。
Time complexity $O(\sqrt{m})$