Eisenstein's Criterion arises from exploiting factorization structure in a "simpler" image ring. $ $ Here the structure exploited is simply that the image reduces to a $\rm\color{#0a0}{prime\ power}$, and prime products always factor $\rm\color{#0a0}{uniquely}$ in domains. Below we explore this viewpoint.
The first part of the hypothesis says: $\bmod p\!:\ f \equiv a x^n\not\equiv 0\,$ is ($a$ times) a $\rm\color{#0a0}{prime\ power}$ $\:\!x^n$
so a $\rm\color{#c00}{proper}$ factorization has form $\, gh\equiv (b x^i) (c x^j)\,$ for $\,\color{#c00}{i,j \ge 1},\,$ by $\rm\color{#0a0}{uniqueness}$ ($\rm\color{#90f}{Lemma}$)
But $\,\color{#c00}{i,j \ge 1}\,\Rightarrow\, p\mid g(0),h(0)\,\Rightarrow\, p^2\mid g(0)h(0)\!=\!f(0),\,$ contra second part of hypothesis.
Said in ideal language it is simply $\,(p,x)^2\equiv (p^2)\,\pmod{\! x}$
$\rm\color{#90f}{Lemma}$ $\:x$ prime, $\,ax^n = fg\Rightarrow f = bx^i,\, g = cx^j,\ i\!+\!j=n,\ a=bc,\,$ in any domain.
Proof $ $ Induct on $\,n.$ $\,n\!=\!0\,$ clear. Else prime $\,x\:\!|\:\! fg\,\Rightarrow\, x\:\!|\:\! f$ or (wlog) $\,x\:\!|\:\! g,\,$ so $\,ax^{n-1}\! = f(g/x),\,$ so by induction $f = bx^i,\ g/x = cx^j,\, i\!+\!j=n\!-\!1,\, a=bc,\,$ so $\,g = cx^{j+1},\, i+(j\!+\!1)=n$.
Recall $x$ is prime in $D[x]\!\iff\! D[x]/x\cong D\,$ is a domain.
Remarks $ $ With this view in mind one can easily generalize Eisenstein's criterion to other polynomials of similar shape, e.g. replace $x$ by a prime polynomial $q$ with $\,f\equiv aq^n\!\pmod{\! p},\,$ where $\,\deg f = \deg q^n,\,$ and $\,p^2\nmid f\bmod q,\,$ i.e. use $\,(p,q)^2\equiv (p)^2\pmod{\! q}\,$ similar to above.
This view also leads to the discriminant-based test for finding shifts $x\mapsto x\!+\!c$ that are Eisenstein, e.g. see this answer.
However, there is a limit to how far one can go since the criterion corresponds to totally ramified primes, but there are number fields without such primes (see e.g. Keith Conrad's note Totally ramified primes and eisenstein polynomials, and note that we can also apply Eisenstein to reversed (reciprocal) polynomials, as here)
In order to find fruitful generalizations one may apply valuation-theoretic ideas. The Eisenstein and related irreducibility tests are all special cases of much more general techniques exploiting Newton polygons. They yield the master theorem behind all these related results. An excellent comprehensive introduction can be found in Filaseta's lecture notes - see the links below. See also this MathOverflow question.
[1] http://www.math.sc.edu/~filaseta/gradcourses/Math788F/latexbook/
[2] http://www.math.sc.edu/~filaseta/gradcourses/Math788F/NewtonPolygonsTalk.pdf
[3] Newton Polygon Applet
http://www.math.sc.edu/~filaseta/newton/newton.html
[4] Abhyankar, Shreeram S.
Historical ramblings in algebraic geometry and related algebra.
Amer. Math. Monthly 83 (1976), no. 6, 409-448.
http://mathdl.maa.org/mathDL/22/?pa=content&sa=viewDocument&nodeId=2964&pf=1