pub fn moderate_u128_divrem( n: u128, d: u64, factor: u128, factor_shr: u32) -> (u128, u64)
Fast division/remainder algorithm for u128, without a fast native approximation.