Function mul_div_ceil

pub fn mul_div_ceil(x: Number, y: Number, denom: Number) -> Number

Multiplies and divides with ceiling: ceil(x * y / denom). Uses 512-bit intermediate (U512) to avoid overflow; returns 0 if denom is 0. Implemented via quotient/remainder to avoid relying on addition overflow behavior.