Otra solución más. Se basa en el hecho de que $$ f = x^4+x^3+x^2+x+1\in\Bbb F_7[x] $$ es recíproco, y trata de utilizar la idea de la OP. (No debería haber división en dos polinomios de grado dos).
Supongamos que existe una factorización $$ f=(x^2+ax+b)(x^2+cx+d)\ ,\qquad a,b,c,d\in\Bbb F_7\ .$$ Dejemos que $u\ne 0$ sea una raíz del primer factor en alguna extensión (de grado dos). Entonces $1/u$ también es una raíz. Distinguimos dos casos.
Primer caso : $1/u$ es también una raíz del primer factor. Entonces $b=1$ entonces $d=1$ , $c+a=1$ el producto $$(x^2+ax+1)(x^2+(1-a)x+1)$$ es recíproco, y tenemos que buscar una coincidencia del coeficiente en $x^2$ para algunos valores de $a$ (modulando la acción $a\leftrightarrow (1-a)$ ). No existe tal $a$ con $1+a(1-a)+1=1$ .
Segundo caso : $1/u$ es no una raíz del primer factor. Sea $u$ , $v$ sean las dos raíces de $x^2+ax+b$ . Entonces $1/u$ , $1/v$ son las raíces del segundo factor, que es el polinomio recíproco, hecho mónico, por lo que esperamos una igualdad de polinomios recíprocos $$ (x^2+ax+b)(1+ax+bx^2) = b(x^4+x^3+x^2+x+1)\ . $$ Las dos ecuaciones obtenidas, $ab+a=b$ y $1+a^2+b^2=b$ no tienen solución en $\Bbb F_7$ . (La sustitución de $a=b/(b+1)$ en la segunda ecuación conduce a... $(b^4+b^3+b^2+b+1)/(b+1)^2=0$ .)
Nota : Esto es posiblemente más corto que tomar todos los polinomios posibles $x^2+ax+b$ , $b\ne 0$ asociando $c=1-a$ , $d=1/b$ y comprobar si las demás condiciones entre $a,b,c,d$ están satisfechos.
Nota : Los programas de álgebra computacional muestran que tenemos un polinomio irreducible en nuestras manos. Por ejemplo, utilizando sage :
sage: R.<x> = PolynomialRing(GF(7))
sage: f = x^4 + x^3 + x^2 + x + 1
sage: f.is_irreducible()
True
Se puede utilizar entonces la potencia del ordenador también como sigue para concluir. (Con un esfuerzo un poco mayor que para el tecleo, también se puede concluir humanamente). Supongamos que $f$ se divide en dos (o más factores). Se puede comprobar que no hay ninguna raíz en $\Bbb F_7$ . Entonces hay una raíz en $\Bbb F_{7^2}=\Bbb F_{49}$ por lo que los polinomios $f$ y $x^{49}-x$ tienen un divisor común. El gcd es pero... escribimos
sage: gcd( x^49 - x, f )
1
sage: gcd( x^(7^4) - x, f )
x^4 + x^3 + x^2 + x + 1
Para ello, también se puede calcular fácilmente $(x^{49}-x,f)$ como humano al observar que $f$ divide $x^5-1$ Así que $(x^{49}-x,f)=(x^{49}-x^4+x^4-x,f)=(x^4-x,f)=(x^3-1,f)=(x^2+x+1,f)=(x^2+x+1,x+1)=(1,x+1)=1$ .