Clearly $ax + b \equiv 0 \pmod{p}$ is solvable for almost all primes $p$. In fact, if $p$ does not divide $a$, then $x \equiv -ba^{-1}$ is a solution. So for any linear polynomial $f(x)$ we get that for $100\%$ of the primes the equation $f(x) \equiv 0 \pmod{p}$ is solvable. With the use of quadratic reciprocity one can often relatively easy establish that for a given irreducible polynomial $f$ of degree two, for half of all primes the equation $f(x) \equiv 0 \pmod{p}$ has a solution, and I believe that this is true in general (?).
However, I now have the following polynomial of degree $d$:
$$f_d(x) = \sum_{i=0}^{d} \prod_{j = 0, j\neq i}^{d} (x+j)$$
And I would like to know the proportion of primes for which $f_d$ has a root. Is this something that can be calculated?
For some examples, with $d = 1, \ldots, 5$ we get:
$f_1(x) = 2x + 1$
$f_2(x) = 3x^2 + 6x + 2$
$f_3(x) = 4x^3 + 18x^2 + 22x + 6$
$f_4(x) = 5x^4 + 40x^3 + 105x^2 + 100x + 24$
$f_5(x) = 6x^5 + 75x^4 + 340x^3 + 675x^2 + 548x + 120$
For $f_1$ we have that $f_1(x) \equiv 0 \pmod{p}$ is solvable for all odd primes and it's a nice exercise to show that $f_2(x) \equiv 0 \pmod{p}$ is solvable whenever $p \equiv \pm 1 \pmod{12}$. Already for $f_3$ I have no idea what to do, let alone in general. From some searching around on the internet, it looks like Chebotarevs density theorem might be useful, but as of now I don't even understand its statement. Any help would be much appreciated!
Edit: as Will Jagy pointed out in the comments, for odd $d$ it looks like $f_d(x)$ is divisible by $2x + d$ or, alternatively, if $d = 2m+1$, then $f_d(x-m)$ is divisible by $2x+1$. Perhaps interestingly and at least for odd $d$ smaller than $10$, we have that the constant and linear terms of $f_d(x)/(2x + d)$ are equal to those of $f_{d-1}(x)$.
For even $d$ we can write $f_{d}(x-d/2) = \displaystyle \sum_{i=-d/2}^{d/2} \prod_{j = -d/2, j\neq i}^{d/2} (x+j)$ and then it is not hard to see that for $i_1 = -i_2$ the odd powers in the products for $i = i_1$ and $i = i_2$ cancel each other out, so $f_{d}(x+d/2)$ is the sum of even powers of $x$. For example:
$f_2(x) = 3(x+1)^2 - 1$
$f_4(x) = 5(x+2)^4 - 15(x+2)^2 + 4$
$f_6(x) = 7(x+3)^6 - 70(x+3)^4 + 147(x+3)^2 - 36$
$f_8(x) = 9(x+4)^8 - 210(x+4)^6 + 1365(x+4)^4 - 2460(x+4)^2 + 576$
$f_{10}(x) = 11(x+5)^{10} - 495(x+5)^8 + 7161(x+5)^6 - 38225(x+5)^4 + 63228(x+5)^2 - 14400$
Here one can note that these coefficients all seem divisible by small primes.
With the help of a computer I have checked that $f_d(x)$ is irreducible for even $d \le 200$ and for odd $d < 200$, $f_d(x)$ can be written as $(2x+d)$ times an irreducible.
– Woett Sep 20 '18 at 17:59Confirmation that the Galois group of $f_4(x)$ equals $D_4$ and that this implies density $3/8$ for $d = 4$.
Finding the density for $d = 6,8,10$.
(Mainly for theoretical and esthetic reasons) Proving that $2x + d$ divides $f_d(x)$ for all odd $d$l.