Your sample proportion of having a "pair" is $\hat{p} = 0.019$, which is the best estimate for the population proportion (assuming random sampling). So we set $\hat{p} = p$, and therefore the number of people in the population that have a "pair" is $40,000,000 \times 0.019 = 760,000$.
Now assuming independent observations (your sample is about 3.75% of the population, so OK) and a sufficient number of both "failures" and "successes" (both have more than 10 observations, so OK), then:
$$ \hat{p} \sim N(p, \frac{p(1-p)}{N})$$
The standard error is defined as the standard deviation of the sampling distribution of $\hat{p}$, so in this case:
$$ SE = \sqrt{\frac{p(1-p)}{N}} = \sqrt{\frac{0.019(1-0.019)}{1,500,000}} = 0.0001$$
Note that this is so small, because N is so large! Now, in order to construct a 95% confidence interval:
$$\hat{p} \pm 1.96\times SE = 0.019 \pm 0.000196 = [0.018804, 0.019196]$$
Translating this back to the population, the "true" number of people in the population that have a "pair" is in the interval $[752160, 767840]$. (A lot more goes into the interpretation of the confidence interval, but I'm not going to elaborate here.)