For reasons of readability let us write An=(n2+11n−4)⋅n!+33⋅13n+4. First, consider the value of An modulo 8 for n≥4. We have 8∣1⋅2⋅3⋅4⋅⋯⋅n=n!, so
An=(n2+11n−4)⋅n!+33⋅13n+4≡0+1⋅5n+4(mod8).
Since 52≡1(mod8), we have 52k≡1(mod8) and 52k−1≡5(mod8) for all k∈N. For even n we get An≡5(mod8) and for odd n we have An≡1(mod8). On the other hand a perfect square can only give the remainder of 0,1 or 4 when divided by 8, so all even n≥4 are out of consideration.
Now, consider the value of An modulo 7 for n≥7. We have 7∣n!, so
An=(n2+11n−4)⋅n!+33⋅13n+4≡0+5⋅(−1)n+4(mod7).
Since (−1)2=1(mod7), we have (−1)2k=1(mod7) and (−1)2k−1=−1(mod7) for all k∈N. For even n we get An≡2(mod7), for odd n we get An≡6(mod7). A perfect square can only give the remainders of 0,1,2 or 4 when divided by 7, so we can exclude all odd n≥7.
The remaining options are n=1,2,3 and 5.
If n=3 we have
A3≡3⋅1+3⋅33+4≡3(mod5).
A perfect square can only give the remainder of 0,1 or 4 when divided by 5, so A3 is not a perfect square.
A similar reasoning helps us see that n=5 does not work, since
A5≡(−4)⋅0+3⋅35+4≡3(mod5).
If n=1, then A1=441=212. If n=2, then A2=5625=752. Thus, n=1 and n=2 are the only solutions.