0, 1, 2, 5, 12, 13, 34, 70, 75, 89, 179, 233, 408, 507, 610, 1120, 1597, 2378, 2673, 2923, 3468, 4181, 6089, 10946, 13860, 15571, 16725, 19760, 23763, 28657, 39916, 51709, 80782, 75025, 113922, 162867, 206855, 196418, 249755, 353702
Offset: 1
n = 5: a(5) = k(5) = 12 because m(5) = A002559(5) = 29 with the triple MT(5) = (2, 5, 29). Whence 2*k1(5) - 5 == 0 (mod 29) for k1(5) = 17 < 29, and 5*k2(5) - 2 == 0 (mod 29) leads to k2(5) = 12. The smaller value is k2(5) = k(5) = 12. This leads to the form coefficients MF(5) = [29, 63, -31].
The forms MF(n) = [m(n), 3*m(n) - k(n), l(n) - 3*k(n)] with l(n) := (k(n)^2 + 1)/m(n) begin: [1, 3, 1], [2, 4, -2], [5, 11, -5], [13, 29, -13], [29, 63, -31], [34, 76, -34], [89, 199, -89], [169, 367, -181], [194, 432, -196], [233, 521, -233], [433, 941, -463], [610, 1364, -610], [985, 2139, -1055], [1325, 2961, -1327], [1597, 3571, -1597], [2897, 6451, -2927], [4181, 9349, -4181], [5741, 12467, -6149], [6466, 14052, -6914], [7561, 16837, -7639] ... .
The quadratic irrationals xi(n) = ((2*k(n) - 3*m(n)) + sqrt(D(n)))/(2*m(n)) begin: (-3 + sqrt(5))/2, -1 + sqrt(2), (-11 + sqrt(221))/10, (-29 + sqrt(1517))/26, (-63 + sqrt(7565))/58, (-19 + 5*sqrt(26))/17, (-199 + sqrt(71285))/178, (-367 + sqrt(257045))/338, (-108 + sqrt(21170))/97, (-521 + sqrt(488597))/466, (-941 + sqrt(1687397))/866, (-341 + sqrt(209306))/305, (-2139 + sqrt(8732021))/1970, (-2961 + sqrt(15800621))/2650, (-3571 + sqrt(22953677))/3194, (-6451 + sqrt(75533477))/5794, (-9349 + sqrt(157326845))/8362, (-12467 + 5*sqrt(11865269))/11482, (-3513 + 5*sqrt(940706))/3233, (-16837 + sqrt(514518485))/15122, ... .
The invariant M(xi(n)) = Lagrange(n) numbers begin with n >=1: sqrt(5), 2*sqrt(2), (1/5)*sqrt(221), (1/13)*sqrt(1517), (1/29)*sqrt(7565), (10/17)*sqrt(26), (1/89)*sqrt(71285), (1/169)*sqrt(257045), (2/97)*sqrt(21170), (1/233)*sqrt(488597), (1/433)*sqrt(1687397), (2/305)*sqrt(209306), (1/985)*sqrt(8732021), (1/1325)*sqrt(15800621), (1/1597)*sqrt(22953677), (1/2897)*sqrt(75533477), (1/4181)*sqrt(157326845), (5/5741)*sqrt(11865269), (10/3233)*sqrt(940706), (1/7561)*sqrt(514518485), ... .
Comments