cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-1 of 1 results.

A364248 For n >= 3, a(n) is the least r >= 0 such that the elliptic equation y^2 = n^3 + n^2 + 2*r*n + r^2 has an integer solution.

Original entry on oeis.org

0, 2, 5, 9, 14, 0, 27, 5, 44, 12, 65, 21, 0, 32, 119, 9, 152, 35, 21, 77, 230, 0, 275, 117, 54, 14, 377, 41, 434, 32, 55, 221, 0, 27, 629, 285, 52, 20, 779, 49, 860, 11, 21, 437, 1034, 0, 1127, 75, 34, 182, 1325, 27, 110, 159, 19, 725, 1652, 10, 1769, 837, 0, 320, 195, 99, 2144, 374
Offset: 3

Views

Author

Ctibor O. Zizka, Sep 01 2023

Keywords

Comments

The equation y^2 = n^3 + A*n^2 + B*n + C, where A = 1, B = 2*r, C = r^2 is a minimal model of an elliptic curve with integral coefficients, for details see the Links section.
For n >= 4, n*(n - 3)/2 >= a(n) >= 2.

Examples

			n = 6: y^2 = 6^3 + (6 + r)^2 is valid for r = 9, 19, 47. Least r is 9 thus a(6) = 9 and [y, n, n+r] is [21, 6, 15].
n = 7: y^2 = 7^3 + (7 + r)^2 is valid for r = 14, 164. Least r is 14 thus a(7) = 14 and [y, n, n+r] is [28, 7, 21].
		

Crossrefs

Programs

  • PARI
    a(n)=vecmin(((select((x)->x[1]>=0&&x[2]>=n,thue(thueinit(x^2-1,1),n^3)))~[1]))-n \\ Thomas Scheuerle, Sep 03 2023
    
  • Sage
    def a(n):
        for d in Integer(n^3).divisors():
            if ((d-n^3/d)%2 == 0) and ((d-n^3/d)/2 >= n):
                return (d-n^3/d)/2 - n  # Robin Visser, Sep 30 2023

Formula

a(p) = p*(p - 3)/2, p >= 5 a prime number.
a(n) = 0 for n from A005563.
Showing 1-1 of 1 results.