P(x) = x4 + x³ + x² + x + 1 = (x5 - 1)/(x - 1) since the LHS is a geometric series
To find the zeros of P(x), we need find the zeroes of (x5 - 1) for x =/= 1 (i.e. non real solutions of x)
So solving x5 - 1 = 0 for x=/=1
x = cis 2π/5, cis 4π/5, cis 6π/5, cis 8π/5 (ignore cis 0 = 1 as we require non-real solutions)
= cis 2π/5, cis 4π/5, cis -2π/5, cis -4π/5 since co-efficients are real, we have complex conjugate roots
Hence P(x) = (x - cis 2π/5)(x - cis -2π/5)(x - cis 4π/5)(x - cis -4π/5)
However, we require the factorisation to be over the real field, so note that:
*cis θ + cis (-θ) = cos θ + isin θ + cos (-θ) + isin (-θ)
= cos θ + isin θ + cos θ - isin θ (as cos (-θ) = cos θ and sin (-θ) = -sin θ)
= 2cos θ
AND
*cis θ.cis (-θ) = cis 0 = 1
Thus,
P(x) = (x² - 2xcos 2π/5 + 1)(x² - 2xcos 4π/5 + 1)