Intento resolver el siguiente sistema sobredeterminado de polinomios
\begin{align} & p_1(x_1,x_2,\ldots,x_n)=0, \\ & p_2(x_1,x_2,\ldots,x_n)=0, \\ & \vdots \\ & p_m(x_1,x_2, \ldots, x_n)=0, \\ & (x_1-l_1) (x_1-(l_1+1)) (x_1 - (l_1 + 2))\cdots (x_1-(u_1-1))(x_1-u_1)=0, \\ & \vdots \\ & (x_n-l_n) (x_n-(l_n+1)) (x_n - (l_n + 2))\cdots (x_n-(u_n-1))(x_n-u_n)=0 \cdots (x_n-u_n)=0. \end{align}
Aquí $l_1,u_1, \ldots, l_n, u_n$ son enteros no negativos. El último $n$ significan que $x_i \in \{l_i, l_i +1, l_i+2, \ldots, u_i-1, u_i\}$ para todos $i \in \{1,\ldots, n\}$ . Sé que podemos utilizar las bases de Groebner para resolver el sistema polinómico.
Pero, ¿hay alguna forma mejor de resolver un sistema polinómico tan sobredeterminado? Agradeceré cualquier sugerencia o indicación.