$x^{18}-x^3=x^3(x^{15}+1)=x^3(x^{5\cdot3}+1)=x^3(x^5+1)(x^{10}+x^5+1)$
$x^5+1$ has the root $1$ in $\mathbb{F}_2$ so using the factor theorem I got: $x^5+1=(x+1)(x^4+x^3+x^2+x+1)$ (how to prove that the second factor is irreducible over $\mathbb{F}_2$ ?).
Noticing that the multiplicative order is $4$, so the initial polynomial can be factorized only with polynomials of degree $1$, $2$ and $4$ (the divisors of $4$). Since $0$ and $1$ are not roots of $x^{10}+x^5+1$, its factors has to be of degree $2$ and $4$, the only irreducible polynomial deg $2$ over $\mathbb{F}_2$ is $x^2+x+1$, so dividing $x^{10}+x^5+1$ by it, I got $x^8+x^7+x^5+x^4+x^3+x+1$, which it has to decompose in two degree 4 irreducible polynomials. But this method needs too much time.
Another way is to compute the table of $\mathbb{F}_{16}$ (splitting field of the initial polynomial over $\mathbb{F}_2$), the cyclotomic cosets, etc.
Is there any faster method ?