Solution:
If x=0 we have
z2−22y=5t⟺(z+2y)(z−2y)=5t
Putting z+2y=5a and z−2y=5b with a+b=t we get 5a−5b=2y+1. This gives us b=0 and now we have 5t−1=2y+1. If y≥2 then consideration by modulo 8 gives 2∣t. Putting t=2s we get (5s−1)(5s+1)=2y+1. This means 5s−1=2c and 5s+1=2d with c+d=y+1. Subtracting we get 2=2d−2c. Then we have c=1,d=2, but the equation 5s−1=2 has no solutions over nonnegative integers. Therefore so y≥2 in this case gives us no solutions. If y=0 we get again 5t−1=2 which again has no solutions in nonnegative integers. If y=1 we get t=1 and z=3 which gives us the solution (t,x,y,z)=(1,0,1,3).
Now if x≥1 then by modulo 3 we have 2∣t. Putting t=2s we get
3x4y=z2−52s⟺3x4y=(z+5s)(z−5s)
Now we have z+5s=3m2k and z−5s=3n2l, with k+l=2y and m+n=x≥1. Subtracting we get
2⋅5s=3m2k−3n2l
Here we get that min{m,n}=0. We now have a couple of cases.
Case 1. k=l=0. Now we have n=0 and we get the equation 2⋅5s=3m−1. From modulo 4 we get that m is odd. If s≥1 we get modulo 5 that 4∣m, a contradiction. So s=0 and we get m=1. This gives us t=0,x=1,y=0,z=2.
Case 2. min{k,l}=1. Now we deal with two subcases:
Case 2a. l>k=1. We get 5s=3m−3n2l−1. Since min{m,n}=0, we get that n=0. Now the equation becomes 5s=3m−2l−1. Note that l−1=2y−2 is even. By modulo 3 we get that s is odd and this means s≥1. Now by modulo 5 we get 3m≡22y−2≡1,−1(mod5). Here we get that m is even as well, so we write m=2q. Now we get 5s=(3q−2y−1)(3q+2y−1).
Therefore 3q−2y−1=5v and 3q+2y−1=5u with u+v=s. Then 2y=5u−5v, whence v=0 and we have 3q−2y−1=1. Plugging in y=1,2 we get the solution y=2,q=1. This gives us m=2,s=1,n=0,x=2,t=2 and therefore z=13. Thus we have the solution (t,x,y,z)=(2,2,2,13). If y≥3 we get modulo 4 that q,q=2r. Then (3r−1)(3r+1)=2y−1. Putting 3r−1=2e and 3r+1=2f with e+f=y−1 and subtracting these two and dividing by 2 we get 2f−1−2e−1=1, whence e=1,f=2. Therefore r=1,q=2,y=4. Now since 24=5u−1 does not have a solution, it follows that there are no more solutions in this case.
Case 2b. k>l=1. We now get 5s=3m2k−1−3n. By modulo 4 (which we can use since 0<k−1=2y−2) we get 3n≡−1(mod4) and therefore n is odd. Now since min{m,n}=0 we get that m=0,0+n=m+n=x≥1. The equation becomes 5s=22y−2−3x. By modulo 3 we see that s is even. We now put s=2g and obtain (2y−1−5g)(5g+2y−1)=3x. Putting 2y−1−5g=3h,2y−1+5g=3i, where i+h=x, and subtracting the equations we get 3i−3h=2y. This gives us h=0 and now we are solving the equation 3x+1=2y.
The solution x=0,y=1 gives 1−5g=1 without solution. If x≥1 then by modulo 3 we get that y is even. Putting y=2y1 we obtain 3x=(2y1−1)(2y1+1). Putting 2y1−1=3x1 and 2y1+1=3x2 and subtracting we get 3x2−3x1=2. This equation gives us x1=0,x2=1. Then y1=1,x=1,y=2 is the only solution to 3x+1=2y with x≥1. Now from 2−5g=1 we get g=0. This gives us t=0. Now this gives us the solution 1+3⋅16=49 and (t,x,y,z)=(0,1,2,7).
This completes all the cases and thus the solutions are (t,x,y,z)=(1,0,1,3),(0,1,0,2),(2,2,2,13), and (0,1,2,7).