pub fn gcd(a: u64, b: u64) -> u64
Greatest common divisor (binary/Euclid). gcd(0, n) = n.
gcd(0, n) = n