OFFSET
1,2
COMMENTS
This sequence is finite by Siegel's theorem on integral points. The next term, if it exists, is greater than 10^35. - David Radcliffe, Jan 01 2024
MAPLE
select(t -> issqr(t+1), [seq(i*(i+1)*(i+2)/6, i=0..10^6)]); # Robert Israel, Jan 02 2024
MATHEMATICA
t = {}; Do[tet = n (n + 1) (n + 2)/6; If[IntegerQ[Sqrt[tet + 1]], AppendTo[t, tet]], {n, 0, 100000}]; t (* T. D. Noe, Mar 18 2013 *)
PROG
(Python)
import math
for i in range(1<<33):
t = i*(i+1)*(i+2)/6 + 1
sr = math.isqrt(t)
if sr*sr == t:
print (t-1, sep=' ')
(PARI)
A000292(n) = n*(n+1)*(n+2)\6;
for(n=0, 10^9, t=A000292(n); if (issquare(t+1), print1(t, ", ") ) );
/* Joerg Arndt, Mar 16 2013 */
KEYWORD
nonn,more,fini
AUTHOR
Alex Ratushnyak, Mar 15 2013
STATUS
approved