Function moderate_u128_divrem

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.