Not very elegant of a solution there, here is my own alternatively:
From the manipulation at the start and observation, we can see that x=i, x=-i are not solutions. Hence when we factorise over R we just add the conjugates and express as quadratic factors: