Optimal monomial quadratization for ODE systems