Запоздалое послесловие.
Так или иначе, вопрос быстро сводится к уравнению

, которое имеет лишь конечное число решений, а если

- целый квадрат, вроде бы только одно:

.
Что касается непосредственно уравнения

, то целые решения у него (отличные от

) при полном квадрате

могут и быть.
Например, при

это

. Пожалуй, лучше тут обойтись без нахождения фундаментальной единицы кубического поля и алгорифма Вороного, что естественно напрашивается при решении уравнения

.
Для изначально рассматриваемой задачи, на мой вкус, проще сделать так:
Подробно:
1.Написать уравнение

(это сделал
Shadow),
2. Сделать замену

(а не

как предлагает
Andrey A) и получить уравнение

,
3. Сделать еще одну замену

и получить уравнение

4. Убедиться, что ранг

равен нулю, а именно:
Код:
gp > E=ellinit([0,0,0,0,194481]);
gp > ellanalyticrank(E)[1]
%13 = 0
5. Вычислить рациональные точки кручения на кривой с уравнением

Код:
E := EllipticCurve([0,0,0,0,194481]);
E;
G, h := TorsionSubgroup(E);
torsion_pts_E := [ h(g) : g in G ];
torsion_pts_E;
Результаты:
Код:
Elliptic Curve defined by y^2 = x^3 + 194481 over Rational Field
[ (0 : 1 : 0), (0 : 441 : 1), (0 : -441 : 1) ]
Точки кручения дают для уравнения

решение

.
Из п.4 следует, что других рациональных решений для

и, следовательно, для

нет.
Т.о. задача решена.