There are obvious counterexamples [hint: $\,f(0)=p\Rightarrow p\mid f(pn)\,$] and this idea leads to a general proof that no such prime-producing polynomial exists.
However, you probably wonder why the first $32$ values are all prime. To explain that, notice that if you shift x by $16$ then it becomes $\rm\:f(x) = x^2 + x + 17.\:$ Now one may apply the following theorem on prime-producing polynomials (generalizing the famous example $\rm\:x^2+x+41\:$ noticed by Euler), noting that the ring of integers of $\rm\:\mathbb Q(\sqrt{-67})\:$ is a UFD.
Theorem $\ $ The polynomial $\rm\ f(x)\ =\ (x-\alpha)\:(x-\alpha')\ =\ x^2 + x + k\ $ assumes only prime values for $\rm\ 0\ \le\ x\ \le\ k-2 \ \iff\ \mathbb Z[\alpha]\ $ is a PID.
Hint $\: (\Rightarrow)\ $ Show all primes $\rm\ p\le\sqrt{n},\, n=1\!-\!4k\,$ have $\rm \left[\frac{n}p\right]\! = -1,\,$ i.e. $\rm\! \bmod p\!:\rm\, n\,$ is nonsquare, so no primes split/ramify.
For proofs, see e.g. Cohn, Advanced Number Theory, pp. 155-156, or Ribenboim, My numbers, my friends, 5.7 p.108. Note: both proofs employ the bound $\rm\ p < \sqrt{n}\ $ without explicitly mentiioning that this is a consequence of the Minkowski bound - presumably assuming that is obvious to the reader based upon earlier results. Thus you'll need to read the prior sections on the Minkowski bound. Compare Stewart and Tall, Algebraic number theory and FLT, 3ed, Theorem 10.4 p.176 where the use of the Minkowski bound is mentioned explicitly.
Alternatively see the self-contained paper [1] which proceeds a bit more simply, employing Dirichlet approximation to obtain a generalization of the Euclidean algorithm (the Dedekind-Rabinowitsch-Hasse criterion for a PID). If memory serves correct, this is close to the approach originally employed by Rabinowitsch when he first published this theorem in 1913.
[1] Daniel Fendel, Prime-Producing Polynomials and Principal Ideal Domains,
Mathematics Magazine, Vol. 58, 4, 1985, 204-210