Solve in positive integers: $n!=m(m+1)$

I'm pretty sure this is open. As suggested from Brocard's problem, it is interesting to investigate the Diophantine equations $$n!=P(m)$$ for polynomials $P$. You can see the paper "On polynomial-factorial diophantine equations", by D. Berend, J.E. Harmse where they make some advances on the problem and prove that this equation has finitely many solutions for many classes of polynomials (irreducible, with an irreducible factor of large degree or with an irreducible factor to a large power).

So by their results it is known that the equation $n!=m^r(m+1)$ has finitely many solutions if $r\geq 4$. But for $r\in \{1,2,3\}$ the problem is open.


I thought I should share with you the results from my computations from the past two days. All these computations could be carried out to significantly higher limits with better code and a bit more time. However, my motivation was just to play with them.

Please check the code before taking these computations as facts! I did not check the code too carefully.


Computation 1

The equation $$n! = m(m+1)$$ does not have integer solutions for $3 < n < 10^9$.

I used Mark Sapir's observation that if $n$ is a solution then $4n!+1$ is a square. I wrote a small program which calculates Legendre symbols $$\left(\frac{4n!+1}{p}\right),$$ with a set of primes $p>n$, to check if there might be a solution for $n$. Confirming that for $3 < n < 10^9$ there are no solutions took about one hour on my laptop using 40 primes. The code is available at http://users.jyu.fi/~tamaraja/temp/brogen.c for anyone interested to see it.


Computation 2

I wanted to test the slightly stronger guess in Mark Sapir's post that for $p>5$, $p^2$ does not divide any integer of the form $4n!+1$. This turned out not to be true. I checked that for $5 < p <10^6$ there is one exception, namely $$761671^2 ~|~ 4\cdot446142!+1.$$ Here I used only brute force and the trivial observation that $$p ~|~ 4n!+1 \Rightarrow n < p.$$ Testing this range for $p$ took a couple of hours. The code is available at http://users.jyu.fi/~tamaraja/temp/sqrdiv.c


Computation 3

Byron Schmuland pointed out the similarity to Brocard's problem in his comment. Because the previous computations I found for the equation $$n! + 1 = m^2$$ were carried out over a decade ago, I decided to extend them a bit.

With the same approach which was used by Berndt and Galway in 2000 (and in Computation 1 above) I confirmed that the equation has no solutions for $7 < n < 10^{10}$. This took about one day. The code was the same as in Computation 1.


Computation 4

After the previous 3 computations I did a little bit of the "missing" fourth. The equation $$n!+1 \equiv 0 \pmod{p^2}$$ has no solutions $(p,n)$ with $p$ prime and $613 < p < 10^6$.


Let us look for a probabilistic hint. Given $n\ge3$, define $N:=[\sqrt{n!}]$, the integer part of the square root $n!$ . Then $n!\in[N^2+1,\ldots,N^2+2N]$. The answer to the question is positive if and only if $n!=N^2+N$, because $m$ has to be $N$. At first glance, the probability of this event is $1/2N$. However, we know a priori that both $n!$ and $N(N+1)$ are even. Therefore this probability is $1/N\sim(n!)^{-1/2}$.

Since the series $$\sum_{n=2}^{\infty}\frac{1}{\sqrt{n!}}$$ converges, I expect that the number of solutions to this problem be finite. Actually, I checked that the answer is No for $4\le n\le 10$. Then the number of solutions with $n\ge4$ can be estimated by the series $$\sum_{n=11}^{\infty}\frac{1}{\sqrt{n!}}$$ Because this number is very small (not greater than $10^{-3}$), I bet that there does not exist a solution $n\ge4$.

This is the same kind of reasoning that is used to guess that there does not exist a prime number among Fermat numbers $F_m$ with $m\ge5$.