OFFSET
1,1
COMMENTS
y-x = 1 is not included because (x+1)!/x! + 1 = x+2.
y-x = 4 is not included because (x+4)!/x! + 1 = (x^2+5*x+5)^2 = A028387(x+1)^2 + 1.
y-x = 2 is easily seen not to be possible.
y-x = 3 corresponds to the elliptic curve (x + 1)*(x + 2)*(x + 3) + 1 = k^2.
According to Sage, the only integral points on that curve are (x,k) = (-3,+-1),(-2,+-1),(-1,+-1),(1,+-5),(3,+-11),(54,+-419).
Brocard's problem is the case x=1.
If a(6) exists, it is greater than 10^16.
LINKS
EXAMPLE
4!/1! + 1 = 5^2.
6!/3! + 1 = 11^2.
7!/1! + 1 = 71^2.
57!/54! + 1 = 419^2.
14!/7! + 1 = 4159^2.
MAPLE
N:= 10^15: # to get all terms <= N
L:= N^2-1:
Res:= {5, 11, 419}:
for x from 1 do
t:= (x+1)*(x+2)*(x+3)*(x+4)*(x+5);
if t > L then break fi;
if issqr(t+1) then
Res:= Res union {sqrt(t+1)};
fi;
for y from x+6 do
t:= t*y;
if t > L then break fi;
if issqr(t+1) then
Res:= Res union {sqrt(t+1)};
fi;
od;
od:
sort(convert(Res, list));
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Robert Israel, Oct 15 2019
STATUS
approved