cp's OEIS Frontend

A374833 Elliptic net associated to y^2 + y = x^3 + x^2 - 2*x, based on the non-torsion generator points P = [0, 0] and Q = [1, 0].

%I A374833 #19 May 30 2025 08:01:34
%S A374833 0,1,1,-1,1,-1,-3,2,3,-5,-11,-5,1,8,31,38,7,-13,-19,53,94,-249,89,36,
%T A374833 -41,33,479,-4335,-2357,-149,181,-151,-350,919,5959,18041,-8767,-4544,
%U A374833 1535,989,-493,-2591,-12016,182879,3085709,496035,-48259,-11811,-1466,6627,13751,-55287,-201383,5002782,-124991065
%N A374833 Elliptic net associated to y^2 + y = x^3 + x^2 - 2*x, based on the non-torsion generator points P = [0, 0] and Q = [1, 0].
%C A374833 The curve y^2 + y = x^3 + x^2 - 2*x is one of the rank-2 elliptic curves with smallest conductor.
%C A374833 The signs are defined by the Weierstrass sigma function. In the literature are other variants of sign assignment for this particular net presented.
%H A374833 Katherine E. Stange, <a href="https://eprint.iacr.org/2006/392.pdf">The Tate Pairing via Elliptic Nets</a>, Proc. of Pairing Conference 2007, T. Takagi et al. eds., LNCS, Vol. 4575, pp. 329-348, Springer-Verlag, Berlin, 2007.
%F A374833 A(n, k) = ws(z1 + z2)/(ws(z1)^(n^2 - k*n)*ws(z1 + z2)^(k*n)*ws(z2)^(k^2 - k*n)), where ws is the Weierstrass sigma function using the lattice parameters of y^2 + y = x^3 + x^2 - 2*x, z1 is the lattice point corresponding to P = [0, 0] and z2 corresponds to Q = [1, 0].
%F A374833 A(n*c1, n*c2) divides A((n*k)*c1, (n*k)*c2), where c1, c2 are some integer constants not equaling zero simultaneously and k >= 1.
%F A374833 A(n, k) = (A(n-1, k)*A(n-3, k)*A(2, 0)^2 - A(1, 0)*A(3, 0)*A(n-2, k)^2)/A(n-4, k), for n > 4.
%F A374833 A(n, k) = (A(n, k-1)*A(n, k-3)*A(0, 2)^2 - A(0, 1)*A(0, 3)*A(n, k-2)^2)/A(n, k-4), for k > 4.
%F A374833 A(n, n) = (A(n-1, n-1)*A(n-3, n-3)*A(2, 2)^2 - A(1, 1)*A(3, 3)*A(n-2, n-2)^2)/A(n-4, n-4), for n > 4.
%F A374833 |A(n, 0)| = |A178624(n)|.
%e A374833 A(n, k) is a square array read by ascending antidiagonals:
%e A374833 .
%e A374833       --> k*Q
%e A374833  n*P |  0,  1,  -1,   -5,   31,    94
%e A374833   |  |  1,  1,   3,    8,   53,   479
%e A374833   |  | -1,  2,   1,  -19,   33,   919
%e A374833  \ / | -3, -5, -13,  -41, -350, -2591
%e A374833      |-11,  7,  36, -151, -493, 13751
%e A374833      | 38, 89, 181,  989, 6627, 68428
%e A374833 .
%e A374833 n*P means elliptic point multiplication here. A(n, k)^2 is the denominator of the x coordinate from n*P + k*Q with point multiplication and addition under the elliptic group law for rational numbers.
%o A374833 (PARI)
%o A374833 T(n, k) = { local(E, z1); local(E, z2); E=ellinit([0,1,1,-2,0]); z1=ellpointtoz(E,[0,0]); z2=ellpointtoz(E, [1,0]);round(ellsigma(E, n*z1+k*z2)/(ellsigma(E, z1)^(n^2-k*n)*ellsigma(E, z1+z2)^(k*n)*ellsigma(E, z2)^(k^2-k*n) )) }
%o A374833 A(size) = { my(si = max(0, size-5)); M = matconcat([matrix(5,5,m,k,T(m-1,k-1)),matrix(5, si);matrix(si, 5),matrix(si, si)]);
%o A374833 for(k = 1, 5, for(n = 6, size, M[n, k] = (M[n-1, k]*M[n-3, k]*M[3, 1]^2 - M[2, 1]*M[4, 1]*M[n-2, k]^2)/M[n-4, k]));
%o A374833 for(k = 1, 5, for(n = 6, size, M[k, n] = (M[k, n-1]*M[k, n-3]*M[1,3]^2 - M[1, 2]*M[1, 4]*M[k, n-2]^2)/M[k, n-4]));
%o A374833 for(k = 6, size, for(n = 6, size, M[n, k] = (M[n-1, k]*M[n-3, k]*M[3, 1]^2 - M[2, 1]*M[4, 1]*M[n-2, k]^2)/M[n-4, k])) }; M;
%o A374833 sd1(P)=sqrtint(denominator(P[1]));
%o A374833 Pnm(n, m, E, P1, P2) = elladd(E, ellmul(E, P1, n), ellmul(E, P2, m));
%o A374833 Aunsigned(size) = my(E=ellinit([0,1,1,-2,0]), P=[0,0], Q=[1,0]); matrix(size, size, m, n, sd1(Pnm(m-1, n-1, E, P, Q)));
%Y A374833 Cf. A178624.
%K A374833 sign,tabl
%O A374833 0,7
%A A374833 _Thomas Scheuerle_, Jul 21 2024