Пояснение относительно того, откуда взялась эллиптическая кривая
и какая от неё польза.
Пусть
, (
натуральное число) - рациональная диофантова тройка. Попытаемся расширить её до диофантовой четверки.
Для этого нужно найти рациональное число
, так, чтобы три числа
были квадратами.
Тогда
, где
рациональное число.
- уравнение эллиптической кривой.
Заменой
кривая
превращается в эквивалентную ей
.
Эта кривая имеет 8 рациональных точек конечного порядка
.
Рациональные точки бесконечного порядка могут на ней быть, могут не быть. Это зависит от ранга кривой
.
Если ранг положительный, таких точек на кривой бесконечно много, если равен нулю, таких точек нет.
В первой сотне натуральных чисел
62 кривых
имеют ранг больше нуля и наименьшее значение
. (Это число и появилось в условии задачи).
Легко посчитать, что для точек конечного порядка числа
не являются квадратами.
Значит, из рассмотрения выпадают кривые с нулевым рангом. Так, если взять
, то система уравнений
не имеет решений в рациональных числах.
На каждой кривой
с ненулевым рангом бесконечно много рациональных точек, для которых числа
являются квадратами.
Эти точки получаются суммированием всех удвоенных рациональных точек с одной из двух точек
, причем X-координаты этих точек не равны
.
Возможен следующий порядок вычисления
:
1. Найти рациональную точку
бесконечного порядка на кривой
.
2. Если точка не является удвоенной, то удвоить её и найти точку
.
3. Сложить точки
и
и получить точку
.
(Точки
на
соответствуют точкам
на
).
4. Взять x-координату точки
и вычислить
.
Если
,(равенство совершенно невероятно, учитывая теорему Пуанкаре-Гурвица о плотности рациональных точек на связной компоненте эллиптической кривой при наличии на ней хотя бы одной рациональной точки бесконечного порядка) то
- искомое число.
Если имеет место равенство, выбрать другую рациональную точку на
и и повторить п.п.2-4.
Замечу, что все эти действия можно совершить с помощью PARI-GP.
Во всяком случае, приведенное мной выше решение для
найдено именно так.
Код:
1.gp > E=ellinit([0,81+1/81,0,1,0]);
gp > ellgenerators(E)
%4 = [[-289/81, 7616/243]]
2.gp > elladd(ellinit([0,1/81+81,0,1,0]),[-289/81,7616/243],[-289/81,7616/243])
%5 = [5784025/165173904, -778241937005/2122815014208]
3.gp > elladd(ellinit([0,1/81+81,0,1,0]),[5784025/165173904,-778241937005/2122815014208],[-1,80/9])
%6 = [-48846121/598340521, -89118069146480/131724067357629]
(Если заменить в п.3
на
, то получим решение
hippie)
Далее выполняется п.4 и находится
.
Относительно расширения рациональной диофантовой тройки
до диофантовой четверки, из вышеизложенного следует, что это возможно тогда и только тогда, когда ранг кривой
больше нуля.