Dado un polinomio en $n$ variables de la forma
$$P(x_1,x_2,\dots,x_n)=\left(\sum_{i,j}a_{ij}x_ix_j+\sum_{i}b_{i}x_i+c\right)^2$$
¿hay alguna manera de encontrar un polinomio también en $n$ variables de grado 2 como máximo
$$Q(y_1,y_2,\dots,y_n)=\sum_{i,j}d_{ij}y_iy_j+\sum_{i}e_{i}y_i+f$$
con las mismas raíces y que es positivo para todos los valores del dominio? Puedes suponer que todas las variables sólo toman los valores $0$ o $1$ y todos los coeficientes son reales.
Si no es posible encontrar un polinomio en $m$ variables con $m>n$
$$R(y_1,y_2,\dots,y_m)=\sum_{i,j}g_{ij}y_iy_j+\sum_{i}h_{i}y_i+k$$
tal que si $(z_1,z_2,\dots,z_n)$ es una raíz de $P$ entonces $R(z_1,z_2,\dots,z_n,z_{n+1},\dots,z_m)=0$ para algún valor de $z_{n+1},\dots,z_m$ y R es positivo en el dominio?
Motivación: Tengo un montón de ecuaciones que van a alimentar a un ordenador cuántico, pero éste sólo puede manejar interacciones de 2 qubits, por lo que necesito reducir las expresiones a polinomios de grado no superior a 2. Por el momento simplemente uso Mathematica para encontrar una instancia de un polinomio que satisfaga las restricciones anteriores, pero falla cuando $n$ es grande. ¿Existe un procedimiento general para encontrar $R$ ? ¿En qué condiciones $Q$ ¿dejar de existir?