OFFSET
1,2
COMMENTS
Corresponding tetrahedral numbers are in A216268.
The curve 6*(x^2-1)-y*(y+1)*(y+2)=0 is elliptic, and has finitely many integral points by Siegel's theorem. - Robert Israel, Apr 22 2021
The above equation can be transformed to Y^2 = X^3 + 18*X^2 + 72*X + 1296 by setting X = 6*y and Y = 36*x, and then we can use for example SageMath to compute the integral points on this elliptic curve. - Jianing Song, Jan 09 2026
LINKS
MATHEMATICA
t = {}; Do[tet = n (n + 1) (n + 2)/6; If[IntegerQ[s = Sqrt[tet + 1]], AppendTo[t, s]], {n, 0, 100000}]; t (* T. D. Noe, Mar 18 2013 *)
PROG
(Python)
import math
for i in range(1<<30):
t = i*(i+1)*(i+2)//6 + 1
sr = math.isqrt(t)
if sr*sr == t:
print(sr)
(SageMath)
E = EllipticCurve([0, 18, 0, 72, 1296])
pts = E.integral_points()
set(-P[1]//36 for P in pts if (P[1] != 0) & (P[1] % 36 == 0)) # Jianing Song, Jan 09 2026
CROSSREFS
KEYWORD
nonn,fini,full
AUTHOR
Alex Ratushnyak, Mar 15 2013
STATUS
approved
