This is an application of Poonen’s Bertini theorem
[106] by Bucur–David–Feigon–Lalín
[17]. Here’s where this is coming from. The quantity
\(q^2+q+1\) is the number of points on
\(\PP^2(\FF_q)\text{.}\) We can think of each point as a variable, and ask if that point is a rational point of a given plane curve
\(X\text{.}\) At the point 0, we can locally expand out our curve as being cut out by the equation
\(a+bx+cy+\dots\text{.}\) So 0 is on the curve exactly when
\(a=0\text{,}\) and it’s a smooth point if
\(a=0\) and
\(bc\neq 0\text{.}\) If we exclude the case where
\(a=b=c=0\) (by sieving), all other possibilities are equally likely. Now we have
\(q^3-1\) total possibilities, of which
\(q^2-1\) are good. So the probability that a given point is on a random curve is
\((q+1)/(q^2+q+1)\text{,}\) and one uses Poonen’s theorem to ensure that each point contributes independently to the count.