div64_u64_safe
n = div64_u64_safe(pdo * puo * (puv - uv),
n = div64_u64_safe(pdo * puo * uv, pdo * (puv - uv) - puo * uv);