Esta es sólo una respuesta parcial, pero espero que proporciona el inicio de un camino posible para resolver el problema. Como ya se señaló en la pregunta, la $2$ polinomio factores deben ser $7$ veces un cuadrado. En particular,
$$x^3 - x + 1 = 7w^2 \tag{1}\label{eq1}$$
$$x^2 + x + 1 = 7y^2 \tag{2}\label{eq2}$$
para algunos enteros positivos $w$ e $y$. En el caso de que $w = y$ ha sido discutido en la respuesta por Martin Hansen. Nota \eqref{eq2} puede ser hecho en la forma genérica $ax^2 + bxy + cy^2 + dx + ey + f = 0$ moviendo el $7y^2$ a la izquierda. Conectar el adecuado coeficientes en Genérico de dos enteros de la variable de la ecuación de solver da la base de la solución positiva de $x = 2$ e $y = 1$, con la recursividad de las fórmulas de
$$x_{n+1} = 127 x_n + 336 y_n + 63 \tag{3}\label{eq3}$$
$$y_{n+1} = 48 x_n + 127 y_n + 24 \tag{4}\label{eq4}$$
y también:
$$x_{n+1} = 127 x_n - 336 y_n + 63 \tag{5}\label{eq5}$$
$$y_{n+1} = -48 x_n + 127 y_n - 24 \tag{6}\label{eq6}$$
Si usted no está seguro de que se trate o que estas son todas las posibles soluciones, luego otro enfoque es tener en cuenta que \eqref{eq2} también puede ser tratada como una ecuación de segundo grado en $x$, por lo que el determinante debe ser un cuadrado perfecto. En particular, esto significa que
$$1 - 4\left(1 - 7y^2\right) = a^2 \; \Rightarrow \; a^2 - 28y^2 = -3 \tag{7}\label{eq7}$$
para algún entero positivo $a$. De acuerdo a AoPS en Línea de la ecuación de Pell, esta es una Pell-como la ecuación porque es de la forma $x^2 - Dy^2 = k$. Con respecto a la solución, incluyendo la determinación de cuáles son todas las posibles soluciones, esto es a pedido y con varias buenas respuestas dadas en el MSE del que Hace la Pell-como la ecuación de $X^2-dY^2=k$ tienen un simple recursividad como $X^2-dY^2=1$?. Además, el uso Genérico de dos enteros de la variable de la ecuación de solver de nuevo ahora da la base de la solución positiva de $x = 5$ e $y = 1$, con la recursividad de las fórmulas de
$$x_{n+1} = 127 x_n + 672 y_n \tag{8}\label{eq8}$$
$$y_{n+1} = 24 x_n + 127 y_n \tag{9}\label{eq9}$$
y también:
$$x_{n+1} = 127 x_n - 672 y_n \tag{10}\label{eq10}$$
$$y_{n+1} = -24 x_n + 127 y_n \tag{11}\label{eq11}$$
El uso de cualquier método, usted ahora tiene un conjunto relativamente pequeño de posibles soluciones para comprobar mediante el uso de los valores resultantes de $x$ en \eqref{eq1} para determinar si el resultado $w$ es un número entero. Por desgracia, no sé si hay alguna analítica para comprobar si hay alguna solución, mucho menos la menos uno mayor que $w = 1$ si alguno no existen.
Nota también podría considerar \eqref{eq1} ser una ecuación cúbica en $x$ y, a continuación, determinar qué condiciones se requieren para que un integrante de la raíz, y las restricciones que éste impone a $w$. Sin embargo, yo no lo he probado todavía, así que no sé de improviso si te ayuda o no.