In a field $\mathbb{F}_4 = \{0, 1, a, a+1\}$ where $1 + 1 = 0$ and $a^2 = a + 1$,
I am to find the characteristic polynomial of $A = \begin{pmatrix} a+1 & 0 & 0 & a \\ 0 & 1 & 0 & 0 \\ a+1 & a & a & 1 \\ a & a & a+1 & a+1 \end{pmatrix}$.
I first swap row 1 with row 2, and take the determinant of $(A-\lambda I)$ which yields
$(A-\lambda I) = -(1-\lambda)\begin{pmatrix} a+1-\lambda & 0 & a \\ a+1 & a-\lambda & 1 \\ a & a+1 & a+1-\lambda \end{pmatrix}$,
but I just get stuck in very long expressions where I can't easily determine its eigenvalues (because I will later on write the Jordan normal form of $A$).
I arrive at, after using the fact that $1 + 1 = 0$ and $a^2 = a + 1$
$-(1-\lambda) \left( ((a+1-\lambda)(a-\lambda)(a+1-\lambda)-(a+1)) + a((a+1)^2-a(a-\lambda)) \right) = 0$
Which is the determinant of the above $3\times 3$ matrix,
and the first term I can see means one eigenvalue is $\lambda=1$, but then I get an expression I cannot make sense of inside the parentheses.
What am I doing wrong?