Gauss' law for two point charges?
The \$4\pi\$ is the number of steradians on the surface of a sphere.

If you had two points charges within the surface of integration, you could do the sums for each charge individually, and the use superposition to combine them. This of course generalises to any number of charges, and any arbitrary charge distribution, as long as it's entirely within the relevant surface.