I might venture into some math which is slightly above what you need, but I wanted to show you the reason behind the correct equation to use. To start with, the pressure due to radiation, in the most general sense, is given by
$$P_{\mathrm{absorption}} = \frac{\langle S\rangle}{c} \cos(\theta)$$
$$P_{\mathrm{reflection}} = 2\frac{\langle S\rangle}{c} \cos^2(\theta)$$
The top equation represents the case where some light hits your object and the light gets absorbed, and the bottom equation represents the case where light hits and bounces off your object (really it's re-emitted). Your case is the bottom equation so I'll work with that from now on and drop the reflection subscript. Without needing to explain exactly what it is, just accept that $\langle S\rangle$ is the average energy of the photon being absorbed or emitted. Now, that's for a single photon. If we have a lot of photons, generally the energy distribution of all the photons together is said to be described by some function, such as $I_\lambda$ (the subscript refers to wavelength). We can define the full pressure then by not only integrating over all wavelengths of all the photons, but importantly all the directions that the photons may hit or be emitted.
$$P = \frac{2}{c}\int_0^\infty\int_0^{2\pi}\int_0^{\pi/2} I_\lambda(\lambda) \cos^2(\theta)\sin(\theta)\:d\theta\:d\phi\:d\lambda$$
Note now that we've replaced $\langle S\rangle$ with $I_\lambda(\lambda)$. Also, since we're assuming this is pressure on a flat surface, we're not integrating over every angle, but only half of them $-$ any photon, in order to hit the solar sail, has to come from the direction the sail is facing, and can only reflect in the direction it came from, i.e., we're discluding photons "behind" the sail. For that reason, the $d\theta$ integral only goes to $\pi/2$.
If you don't understand the math at this point, that's fine but here's the important takeaway: If we only integrated over all wavelengths, you'd get the equation without the three that you listed. What that means is that the equation without the three doesn't take into account all the possible angles and orientations of the incoming photons. It simply assumes all your photons hit your surface perpendicularly, and then bounce off perpendicularly. This is of course not physically accurate as the photons can come in from a variety of angles.
If you do the angle integrals out, you find that you get the following:
$$P = \frac{2}{c}\int_0^{2\pi}\:d\phi \int_0^{\pi/2} \cos^2(\theta)\sin(\theta)\:d\theta\int_0^\infty I_\lambda(\lambda)\:d\lambda$$
$$P = \frac{2}{c}\frac{2\pi}{3}\int_0^\infty I_\lambda(\lambda)\:d\lambda$$
You can see there in the last line that the three popped out of those integrals over all the possible angles (specifically the $d\theta$ integral). Now, you get the rest of this by assuming some sort of function for $I_\lambda(\lambda)$ which describes the energy distribution of the photons as a function of wavelength. Since your light comes from a star, you know that's just going to be the Planck distribution (specifically the equation for $B_\lambda(\lambda,T)$). This integral is non-trivial, but if you again run through the math, you find
$$P = \frac{2}{c}\frac{2\pi}{3}\int_0^\infty \frac{2hc^2}{\lambda^5}\frac{1}{e^{hc/\lambda k_BT}-1}\:d\lambda = \frac{2}{c}\frac{2\pi}{3}\frac{\sigma T^4}{\pi}$$
$$\boxed{P = \frac{4\sigma T^4}{3c}}$$
This is the true equation for radiation pressure on the solar sail, assuming reflection is occurring. This includes the three because we need to account for the photons hitting and bouncing off at all possible angles (for the solar sail). Note one important differences here is that my equation already has the extra "2" due to reflection factored in. Your last equation differs from mine in that it throws in that extra "2" for reflection. So either it erroneously included the factor of "2" twice, or (more likely) its account for reflection at every angle, not just all possible angles. Remember above, we only integrated the $d\theta$ angle up to $\pi/2$ for the reason explained there. If we integrated up to $\pi$ to include every angle, then we'd get an extra "2" as well and my equation would match your last one. Your solar sails though, don't receive or re-emit photons at every angle, and so my equation differs by that factor of "2".