solve (2k+1)(2k^2+1) = y^2 over the integers