pub fn canon_egcd(a: i64, b: i64, c: i64) -> Option<(i64, i64, i64)>
Assuming a != 0, finds smallest coef_b >= 0 such that a * coef_a + b * coef_b = c.
Panics if a == 0.