pub fn div_mod(n: &BigInt, m: &BigInt, p: &BigInt) -> Result<BigInt, MathError>
Finds a nonnegative integer x < p such that (m * x) % p == n.