For each $m\in\mathbb{N}$ and fixed $a>0,\theta\in\mathbb{R}$, I want to factorizate the polynomial $p_m(x) = x^{2m} - 2a^m\cos (m\theta)x^m + a^{2m}$ into $m$ polynomials of second order. Using basic algebra I get for $m=2$ the following decomposition
$$ p_2 (x)= x^4 - 2a^2\cos(2\theta)x^2 + a^4 = (x^2+2a\cos(\theta)x +a^2)(x^2 - 2a\cos(\theta)x +a^2) $$ I used the formula $\cos (2\ \cdot) = \cos^2 (\cdot) - \sin^2(\cdot)$. I want to do the general case $m\in\mathbb{N}$ and I am afraid that complex numbers are a good tool here, but I can't see how to manage this problem using complex numbers. Please, any help or comment is welcome. Thank you in advance.