I asked this question here yesterday and it was pointed out in the comments that if $f(x)=x^4+x+4$ were reducible, i.e. $f(x)=g(x)h(x)$ for some degree $2$ polynomials in $\mathbb{F}_5[x]$, then as $g$ has degree $2$, it has a root in degree $2$ extension of $\mathbb{F}_5$. Since degree $2$ extension of $\mathbb{F}_5$ are isomorphic to $ \mathbb{F}_{25}$, this gives us that $g$ and in turn $f$ has a root in $\mathbb{F}_{25}$.
Now how do we check if $f$ has any root in $\mathbb{F}_{25}$? $\mathbb{F}_{25}$ is only defined up to isomorphism. I suppose I could take it to be $K=\frac{\mathbb{F}_{5}[x]}{(x^2+x+1)}$, and explicitly write down the elements and compute $f(z)=z^4+z+4$ for every $z \in K$ using the relation $x^2+x+1=0$.
Is this the way to go forward or is there any way to avoid doing $24$ different computations?
I'm assuming you're talking about $x^4+x+4$ as in the linked question, and not $x^4+x+1$, which has the root $x=3$ in $\mathbb F_5$.
One way is to evaluate $\gcd(x^4+x+4,x^{5^2}-x)$. This can be done by first evaluating $x^{25}$ modulo $x^4+x+4$ using exponentiation by squaring, which is annoying but probably doable by hand. Then, you should have that $$\gcd(x^4+x+4,x^{5^2}-x)=\gcd(x^4+x+4,f(x))$$ for some $f\in\mathbb F_5[x]$ of degree at most $3$, from which the standard polynomial $\gcd$ algorithm will allow you to compute the $\gcd$ quickly.
If $x^4+x+4$ is the square of a quadratic, then this $\gcd$ will be that quadratic. If it's the product of two distinct quadratics, the $\gcd$ will be $x^4+x+4$ itself. If it's irreducible, the $\gcd$ will be $1$.
Note: This and its variants are actually used in practice to factor and find roots of polynomials over finite fields.