x^2 (1 +x) == y (-1 + 3 y) in integers