Because $x^2 + ax + b$ has all real coefficients, we know the roots must be complex conjugates. So let $\lambda e^{i \theta}$ be a root, so we know that $(x - \lambda e^{i \theta})(x - \lambda e^{-i \theta}) = x^2 - \lambda x (e^{i \theta} + e^{-i\theta}) + \lambda^2$ is our polynomial.
We also know that $(\lambda e^{i \theta})^3 = 7^3 \implies \lambda e^{i \theta} = 7$, and since $|\lambda e^{i \theta}| = |\lambda| = 7,$ we know that $\lambda = 7$.
This means that $e^{i\theta} = 1$ and so $\theta = 0$. But clearly this doesn't give the right answer, since our polynomial isn't supposed to have real roots. Where is the flaw in my reasoning?
$(\lambda e^{i \theta})^3 = 7^3 \implies e^{3i \theta} = (\frac{7}{\lambda})^3e^{0}$
Using De Moivre's theorem
$\implies e^{i \theta} = (\frac{7}{\lambda}) e^{\frac{2k\pi}{3}i}, \text{k = 0,1,2}$
$\implies \lambda=7$
But as given in the question, roots are not real, so $k$ can not be $0$, therefore,
Roots are $7e^{\frac{2\pi}{3}i}$ and $7e^{\frac{4\pi}{3}i}$
$\implies 7(-\frac{1}{2} + \frac{\sqrt3}{2}i) \text{ and } 7(-\frac{1}{2} - \frac{\sqrt3}{2}i)$