С помощью PARI/GP можно решить уравнение Туэ, но при заданных
, для получения всех решений
придется перебрать все целые
. (Кстати, решения для всех целых
у меня пока нет, хотя для
такие решения отсутствуют, но есть пара сомниnельных
например, 978 с переполнением стека, некогда ими заниматься).
Приведу пример получения решений итогового уравнения
с кодами Pari
, а
меняется от
до
Код:
(15:19) gp > {A=[];
a=115;
b=116;
c=117;
P=(b-a)^2;
Q=4*a*(b-a);
M=6*a^2-2*a*b-2*a*c-2*b*c;
U=4*a*(c-a);
V=(c-a)^2;
for(N=1,1000,
r=thue(thueinit(P*t^4-Q*t^3+M*t^2-U*t+V),N^2);if(r!=A,print1("S=",N," ");print(r)))}
S=1 [[-1, 0], [-1, 1], [1, -1], [1, 0]]
S=2 [[-2, 1], [0, -1], [0, 1], [2, -1]]
S=4 [[-2, 0], [-2, 2], [2, -2], [2, 0]]
S=8 [[-4, 2], [0, -2], [0, 2], [4, -2]]
S=9 [[-3, 0], [-3, 3], [3, -3], [3, 0]]
S=16 [[-4, 0], [-4, 4], [4, -4], [4, 0]]
S=18 [[-6, 3], [0, -3], [0, 3], [6, -3]]
S=25 [[-5, 0], [-5, 5], [5, -5], [5, 0]]
S=32 [[-8, 4], [0, -4], [0, 4], [8, -4]]
S=36 [[-6, 0], [-6, 6], [6, -6], [6, 0]]
S=49 [[-7, 0], [-7, 7], [7, -7], [7, 0]]
S=50 [[-10, 5], [0, -5], [0, 5], [10, -5]]
S=53 [[-3, 1], [-1, 2], [1, -2], [3, -1]]
S=64 [[-8, 0], [-8, 8], [8, -8], [8, 0]]
S=72 [[-12, 6], [0, -6], [0, 6], [12, -6]]
S=81 [[-9, 0], [-9, 9], [9, -9], [9, 0]]
S=98 [[-14, 7], [0, -7], [0, 7], [14, -7]]
S=100 [[-10, 0], [-10, 10], [10, -10], [10, 0]]
S=106 [[-4, 1], [-2, 3], [2, -3], [4, -1]]
S=121 [[-11, 0], [-11, 11], [11, -11], [11, 0]]
S=128 [[-16, 8], [0, -8], [0, 8], [16, -8]]
S=144 [[-12, 0], [-12, 12], [12, -12], [12, 0]]
S=162 [[-18, 9], [0, -9], [0, 9], [18, -9]]
S=169 [[-13, 0], [-13, 13], [13, -13], [13, 0]]
S=196 [[-14, 0], [-14, 14], [14, -14], [14, 0]]
S=199 [[-7, 3], [-1, 4], [1, -4], [7, -3]]
S=200 [[-20, 10], [0, -10], [0, 10], [20, -10]]
S=212 [[-6, 2], [-2, 4], [2, -4], [6, -2]]
S=225 [[-15, 0], [-15, 15], [15, -15], [15, 0]]
S=242 [[-22, 11], [0, -11], [0, 11], [22, -11]]
S=256 [[-16, 0], [-16, 16], [16, -16], [16, 0]]
S=288 [[-24, 12], [0, -12], [0, 12], [24, -12]]
S=289 [[-17, 0], [-17, 17], [17, -17], [17, 0]]
S=324 [[-18, 0], [-18, 18], [18, -18], [18, 0]]
S=338 [[-26, 13], [0, -13], [0, 13], [26, -13]]
S=361 [[-19, 0], [-19, 19], [19, -19], [19, 0]]
S=392 [[-28, 14], [0, -14], [0, 14], [28, -14]]
S=398 [[-8, 1], [-6, 7], [6, -7], [8, -1]]
S=400 [[-20, 0], [-20, 20], [20, -20], [20, 0]]
S=424 [[-8, 2], [-4, 6], [4, -6], [8, -2]]
S=441 [[-21, 0], [-21, 21], [21, -21], [21, 0]]
S=450 [[-30, 15], [0, -15], [0, 15], [30, -15]]
S=477 [[-9, 3], [-3, 6], [3, -6], [9, -3]]
S=484 [[-22, 0], [-22, 22], [22, -22], [22, 0]]
S=512 [[-32, 16], [0, -16], [0, 16], [32, -16]]
S=529 [[-23, 0], [-23, 23], [23, -23], [23, 0]]
S=576 [[-24, 0], [-24, 24], [24, -24], [24, 0]]
S=578 [[-34, 17], [0, -17], [0, 17], [34, -17]]
S=625 [[-25, 0], [-25, 25], [25, -25], [25, 0]]
S=648 [[-36, 18], [0, -18], [0, 18], [36, -18]]
S=676 [[-26, 0], [-26, 26], [26, -26], [26, 0]]
S=722 [[-38, 19], [0, -19], [0, 19], [38, -19]]
S=729 [[-27, 0], [-27, 27], [27, -27], [27, 0]]
S=784 [[-28, 0], [-28, 28], [28, -28], [28, 0]]
S=796 [[-14, 6], [-2, 8], [2, -8], [14, -6]]
S=800 [[-40, 20], [0, -20], [0, 20], [40, -20]]
S=841 [[-29, 0], [-29, 29], [29, -29], [29, 0]]
S=848 [[-12, 4], [-4, 8], [4, -8], [12, -4]]
S=882 [[-42, 21], [0, -21], [0, 21], [42, -21]]
S=900 [[-30, 0], [-30, 30], [30, -30], [30, 0]]
S=954 [[-12, 3], [-6, 9], [6, -9], [12, -3]]
S=961 [[-31, 0], [-31, 31], [31, -31], [31, 0]]
S=968 [[-44, 22], [0, -22], [0, 22], [44, -22]]
Ну, а далее разбираться на любителя.
Хотя вот ещё строки с кодами из Maple для вычисления
при
Код:
solve({F = (x^2 - 1)/a, G = (y^2 - 1)/b, H = (z^2 - 1)/c, a = 115, b = 116, c = 117, k = -7/3, l = 9/199, m = 2, x = l*(k^2 - m), y = l*(k^2 + 2*k + m), z = l*(k^2 + (2*k + 1)*m)}, {F, G, H, a, b, c, k, l, m, x, y, z})
{F = -336/39601, G = -336/39601, H = -336/39601, a = 115, b = 116, c = 117, k = -7/3, l = 9/199, m = 2, x = 31/199, y = 25/199, z = -17/199}