diophantine equation x^2-31y^2 = -3