Разложить квадраты на кубы - это круто.
... разделим почленно уравнение на

Соответствующее уравнение в рациональных числах имеет полное

-параметрическое решение:

С другой стороны, если для произвольного рационального аргумента

находятся

такие, что

видим:

. Домножая последнее на

, получаем в правой части рациональный квадрат, который не может не быть целым числом, ведь в левой части у нас разность целых квадратов. Отсюда определено

, и задача полностью сводится к уравнению

, которое можно записать еще так:

Некоторые его решения следуют из разложения

: если дробь

соответствует ур. Пелля, то дробь

возвращает нужные пары

. Решений как минимум две бесконечные серии, но могут быть и другие, хотя бы потому что пара

не обязана быть взаимно простой. Подробней об этом было
здесь. Пост написан давно и несколько сумбурен. Но уж как есть.