Irreducibility and Splitting Fields Show that over any field $F$, the polynomial $x^3-3x+1$ is either irreducible or splits into linear factors.

Edited: 
This is my attempt: Let $f(x)=x^3-3x+1$. Let $a_1,a_2,a_3$ be the roots of $f$. Suppose char $F\neq 2,3$. Suppose also that $f$  is neither irreducible nor splits in $F$. Then $f$ is reducible which implies that $a_1 \in F$. i.e. $f(x)=(x-a_1)g(x)$, where $g(x)\in F[x]$ is irreducible with deg $g=2$. Let $K$ be the splitting field of $g$. The $K$ is Galois over $F$.
So if $\sigma \in $ Aut($K/F$), then $\sigma (a_1)=a_1$ since $\sigma $ fixes $F$ and $a_1 \in F$. Since $\sigma$ permutes the roots of $f$, WLOG suppose $\sigma (a_2)=a_3$. Then
$\sigma(\triangle) = \sigma((a_1-a_2)(a_1-a_3)(a_2-a_3))=-\triangle$.  
But $\triangle^2=D(f)=81$, so $\triangle = \pm 9 \in F$, so $\triangle \in F$. Therefore,
$9=-9$ $\implies 1=-1 \implies$ char $F =2$ which is a contradiction. So $f$ is either irreducible or splits in $F$. 
Next suppose char $F=2$. Well, I'm not exactly sure what I can say about $f$.  
I would like to know if my approach is correct and also what to do in the second case. Thanks.  
ADDED:
If char $F=2$, then $f=x^3+x+1$. Suppose $b$ is a root of $f$. Then $b^2 $ is also a root, since $f(b^2)=(b^2)^3+b^2+1=(b+1)^2+b^2+1=2b^2+2=0$. Is it enough to conclude that $f$ splits?  

Here is one way to do it using almost no theory, just playing with algebra.
Suppose the polynomial has a root $a$ in $F$.  If you divide $x^3 - 3x + 1$ by $x-a$ the quotient is the polynomial
$$
x^2 + ax + a^2 - 3
$$
in $F(a)[x]$.  From the quadratic formula (assuming characteristic $\neq 2$ here, for the moment, I guess) you can see that this will have its roots in $F(a)$ if and only if $12-3a^2$ is a square in $F(a)$.
You may know that any element in $F(a)$ can be written in the form $pa^2 + qa + r$ for some $p,q,r$ in $F$.  The idea is to square this symbolic expression, rewrite it as a polynomial in $a$ of degree at most $2$, and then see if values $p, q, r$ in $F$ can be found making this equal to $-3a^2 + 0a + 12$.
Using the identity $a^3 = 3a - 1$ (from the fact that $a$ is a root of the given polynomial) and $a^4 = 3a^2 - a$ (from multiplying the previous identity by $a$) one gets
$$
(pa^2 + qa + r)^2 = (3p^2 + 2pr + q^2) a^2 + (-p^2 + 2qr + 6pq) a + (r^2 - 2pq).
$$
So the goal is to find $p,q,r$ satisfying $3p^2 + 2pr + q^2 = -3$, $-p^2 + 2qr + 6pq = 0$, and $r^2 - 2pq = 12$.  Since we don't know anything about $F$, we might optimistically look for these $p, q, r$ in $\mathbb{Z}$.  This system of $3$ polynomial equations in $3$ integer unknowns can be fed to software, and one finds that e.g. $p = 2$, $q = 1$, and $r = -4$ give a solution (no matter what $F$ is!).  So the roots of $x^3 - 3x + 1$ are all in $F(a)$, and we can actually write formulas for them: $a$, $\frac{-a + (2a^2 + a -4)}{2}$, and $\frac{-a - (2a^2 + a - 4)}{2}$, which simplify to $a$, $a^2 - 2$, $-a^2 - a + 2$.  
Although we assumed characteristic $\neq 2$ to use the quadratic formula, we can immediately check that the single identity $a^3 - 3a + 1 = 0$ is indeed enough to ensure that
$$
(x-a)(x - (a^2 - 2)) (x - (-a^2 - a + 2)) = x^3 - 3x + 1.
$$
[In detail: expanding, the coefficient of $x^3$ is $1$ on the nose, the coefficient of $x^2$ is $0$ on the nose, and the coefficient of $x$ is a polynomial in $a$ that, when divided by $a^3 - 3a + 1$, has a remainder of $-3$; similarly the constant term is a polynomial in $a$ that, when divided by $a^3 - 3a + 1$, has remainder of $1$.]
So if $x^3 -3x + 1$ has one root in $F$, it splits in $F$ for the reason that we can explicitly write the other two roots as polynomials in the one that we already have.
