16
$\begingroup$

The question was stimulated by this one. Here it comes:

When you look at the sum $\sum\limits_{k=1}^N k!$ for $N\geq 10$, you'll always find $3$ and $11$ among the prime factors, due to the fact that $$ \sum\limits_{k=1}^{10}k!=3^2\times 11\times 40787. $$ Increasing $N$ will give rise to factors $3$ resp. $11$.

Are $3$ and $11$ the only common prime factors in $\sum\limits_{k=1}^N k!$ for $N\geq 10$?

I think, one has to show, that $\sum\limits_{k=1}^{N}k!$ has a factor of $N+1$, because the upcoming sum will always share the $N+1$ factor as well. This happens for $$ \underbrace{1!+2!}_{\color{blue}{3}}+\color{blue}{3}! \text{ and } \underbrace{1!+2!+\cdots+10!}_{3^2\times \color{red}{11}\times 40787}+\color{red}{11}! $$

  • 7
    Let $a_p = \sum_{k=1}^p k!$. If we assume that $a_p \equiv 0 \mod p$ with "probability" $1/p$, then the expected number of times this happens for $p is $\log \log N$. This quantity is not larger then $3$ until $N$ is about 500,000,000. So the fact that you have only found two examples with a short search does not convince me that there are no more.2012-03-08
  • 2
    By the way, I just checked up to the 500th prime, 3571, without finding another example.2012-03-08
  • 4
    Checked for other examples for primes up to 15000000; none found.2012-12-29
  • 0
    The search up to $3571$ is slightly more surprising: $\sum_{11 < p \le 3571} 1/p$ is just under $1.1$, while $\sum_{3571 < p \le 15000000} 1/p$ is closer to $0.7$. The (marginal) odds only get slimmer from here on out.2013-01-13
  • 0
    While fighting with the factorial sum, I found a nice closed formula at Mathworld (["Factorial sum"](http://mathworld.wolfram.com/FactorialSums.html)): $\frac{-e+Ei(1)+\Re[E_{n+2}(-1)]\Gamma(n+2)}e$. Attempts to prove it by induction led to $e+\Re(E_n(-1))=n\Re(E_{n+1}(1))$, which boils down to the recurrence relation (12) from Mathworld's [$E_n$ function](http://mathworld.wolfram.com/En-Function.html)...2018-02-21
  • 0
    The formula works fine. Put `table round(N[(Ei(1)+Real(ExpIntegralE[(n+2), -1])*Gamma(n+2))/e-1]) n=1 to 15` at http://www.wolframalpha.com/2018-02-21
  • 3
    There is really little known about the distribution of factorials modulo a prime $p$. For instance, it is conjectured that the sequence $(n!)$ attains $\approx (1-\frac{1}{e})p$ residues modulo $p$ and it is known unconditionally that it is at least $$\frac{p\log\log p}{\log\log\log p}.$$ However, this does not help for your question. If I should guess an answer, I'd bet (given the "random" behaviour of the primes) that there are infinitely many primes $p$ such that $p\mid 1!+2!+\cdots+(p-1)!$.2018-02-26
  • 1
    It might be a good idea to make the problem statement unambiguous by explicitly pointing out you're looking for prime factors $p$ shared by all the sums with $N\geq (p-1)$, rather than using $N\geq 10$ (which would make the question trivial).2018-02-26
  • 1
    What do you mean by "common prime factors?" Do you mean in common among the various sums?2018-02-26
  • 0
    @Chickenmancer yes common for all sums like 3 and 11 for N>102018-02-26

1 Answers 1