Let the n-th Polynomial $P_n(x)$ be $$P_n(x) = n + \prod_{k=1}^{n} (x - k)= (x - 1)(x - 2)(x - 3)\cdots(x - n) + n$$
How to prove that for every $n$, $P_n(x)$ is irreducible over integers, or more generally rational field?
I tried using Wolfram Alpha to check some values, and it results in irreducible. As $n$ gets bigger and bigger, it is hard to check if it is irreducible or not. And, I tried to use the fact that $1$ to $n$ is a root of $P_n(x) - n$, but got stuck.
The following is only a partial answer, but I believe it is relevant to the question.
For $\,n \ge 13\,$ the result follows from a far more general, albeit lesser known, irreducibility criterion of Pólya from Verschiedene Bemerkungen zur Zahlentheorie (1919). Below is the English translation of Prasolov's Russian translation in his book Polynomials.
With $\,a_i=i\,$, the criterion implies that all polynomials $\,(x-1)(x-2)\ldots(x-n) \pm k\,$ are irreducible for $\,k = 1, 2, \ldots, \left\lfloor\dfrac{m!}{2^m}\right\rfloor\,$ where $\,m = \left\lfloor\dfrac{n+1}{2}\right\rfloor\,$. For $\,n \ge 13\,$ it is straightforward to show that $\,n \lt \left\lfloor\dfrac{m!}{2^m}\right\rfloor\,$, so the conclusion in OP's question follows.
I can't read German, but the very specific example $\,x(x-1)\ldots(x-n+1) + \lambda n\,$ appears to be mentioned in Pólya's paper as a concrete solution $\lambda=\pm 1$ to Hilbert's existential irreducibility theorem in this case.
An English translation of Pólya's proof (adapted to the present case):
Assume contrariwise that $P_n(x)=g(x) h(x)$ for some $g,h\in\Bbb{Q}[x]$. By Gauss' Lemma we can w.l.o.g. assume that $g,h\in\Bbb{Z}[x]$. One of the factors, say $g(x)$, has degree at least $m\ge\lfloor\dfrac{n+1}2\rfloor$. Observe that $m+1=\deg g(x)+1\le n$ as otherwise the other factor $h(x)$ is a constant.
By Lagrange interpolation formula we can write $$ g(x)=\sum_{i=1}^{m+1} g(i) \Delta_i(x), $$ where for all $i, 0< i\le m+1$, $$\Delta_i(x)=\frac{\prod_{0< j\le m+1, j\neq i}(x-j)}{\prod_{0< j\le m+1, j\neq i}(i-j)}.$$ The leading coefficient $c_i$ of $\Delta_i(x)$ has absolute value $$|c_i|=\frac1{i!(m-i)!}=\frac1{m!}\binom m i.$$ Furthermore, $g(i)$ must be a factor of $P_n(i)=n$, so $|g(i)|\le n$.
Putting all these pieces together, the triangle inequality says that the leading coefficient $c=1$ of $g(x)$ has an upper bound $$ |c|\le\sum_{i=1}^{m+1} |g(i)|\cdot |c_i|\le \frac{n}{m!}\sum_{i=0}^m\binom m i=\frac{n\cdot 2^m}{m!}. $$ As $c=1$, this is a contradiction whenever $n\cdot 2^m<m!$.