Der Ansatz
(2·x + 1)^2 + (70·x + 2)·8 = y^2
ist doch sehr gut. Ich hätte jetzt eine Wertetabelle gemacht für
y = √(4·x^2 + 564·x + 17)
und diese auf ganzzahlige Lösungen untersucht.
Eine andere Möglichkeit ist Wolframalpha für sich schuften zu lassen.
https://www.wolframalpha.com/input/?i=%282*x%2B1%29%5E2%2B%2870*x%2B…
Wolframalpha findet die Lösungen
{{x == -2554, y == ±4965},
{x == -1313, y == ±2481},
{x == -268, y == ±369},
{x == -179, y == ±165},
{x == 38, y == ±165},
{x == 127, y == ±369},
{x == 1172, y == ±2481},
{x == 2413, y == ±4965}}