diophantine y^2+1=1250z