sympy reach the specific polinomial expression from roots

Question

I have the polynomial : x^4+x^3+x^2+x+1. We can factorize the expression below like: (x^2+ax+1)(x^2+bx+1). I found 4 roots by: solve((x**4)+(x**3)+(x**2)+(x)+1)