Function fast_u128_divrem

pub fn fast_u128_divrem(
    n: u128,
    d: u64,
    fast: u128,
    fast_shr: u32,
    factor: u128,
    factor_shr: u32,
) -> (u128, u64)

Fast division/remainder algorithm for u128, without a fast native approximation.