Continuing the curiosity of my last couple questions: Is it the case that for every set of primes $F$, the asymptotic density of the integers coprime to all of $F$ is $\displaystyle \prod\_{p \in F} (1 - 1/p)$?


(The case of finite $F$ is obvious. Furthermore, as the product is an upper bound on the density, this holds whenever the product goes to zero; in particular, by the divergence of the harmonic series, this holds for all co-finite $F$. But that leaves a lot of ground untouched. Ideally, there'd be some way to extend the reasoning from the finite case to the general case, but it's not clear to me how to do so given the general lack of countable additivity of density.)



Relationship of Euler product to coprimality densities for arbitrary sets of primes
GH from MO is answering a more difficult question where one asks for uniform estimates for the number of integers up to $x$ not divisible by primes in a set $F$ (which may also depend on $x$). This is delicate, and the answer is indeed no. 


The problem, as I understand it, keeps $F$ fixed, and asks for the density as $x\to \infty$. This is easier, and the answer is in fact yes. As the OP observes, if the set $F$ is finite the result is true; and also if $\prod\_{p\in F} (1-1/p)=0$ then the result is true (take the restriction of $F$ to primes below some point $N$, 
and note that the density of the numbers coprime to $F\_N$ is an upper bound for 
what we want; and then take $N\to \infty$). 


Now consider a set $F$ with $\prod\_{p\in F}(1-1/p) >0$. Therefore if $G$ denotes the primes not in $F$ then $\prod\_{p\in G} (1-1/p)=0$. Put $\mu\_F(n)$ to be a multiplicative function defined by $\mu\_F(p)=-1$ if $p\in F$ and $\mu\_F(p) =0$ otherwise; put $\mu\_F(p^k)=0$ for all prime powers with $k\ge 2$. So this is a small generalization of the Mobius function. Then 
$$ 
\sum\_{\substack{{n\le x} \\ {(n,F)=1}} } 1 = \sum\_{n\le x} \sum\_{d|n} \mu\_F(d) 
= \sum\_{d\le x} \mu\_F(d) \Big(\frac xd +O(1)\Big).
$$
By assumption $\sum\_{p\in F} 1/p <\infty$ and so the Euler product $\prod\_{p\in F}(1-1/p)$ converges absolutely (to a non-zero real number) and equivalently the sum $\sum\_{n=1}^{\infty} \mu\_F(n)/n$ converges. Thus the main term above equals 
$$
x\sum\_{d=1}^{\infty} \frac{\mu\_F(d)}d + o(x)= x \prod\_{p\in F}\Big(1-\frac{1}{p}\Big) + o(x). 
$$
As for the remainder term, since $\mu\_F(d) =0$ unless $d$ is divisible only by primes in $F$, the remainder term is at most the number of integers up to $x$ that are coprime to all elements in $G$; but we know that the density of such integers tends to zero for large $x$. That completes the proof. 


More generally, Erdos asked whether for any real valued multiplicative function $f$ with $-1\le f(n)\le 1$ for all $n$ we have 
$$ 
\lim\_{x\to \infty} \frac{1}{x} \sum\_{n\le x} f(n) = \prod\_{p} \Big(1+\frac{f(p)}{p} + \frac{f(p^2)}{p^2}+\ldots \Big) \Big(1-\frac 1p\Big).
$$
Your question is a special easy case of this problem; Erdos's problem also includes the prime number theorem as the special case when $f(n)=\mu(n)$. Erdos's problem was fully solved by Wirsing, who showed that the limit does equal the expected Euler product. The extension to complex valued multiplicative functions was made by Gabor Halasz (another GH, not from MO!).


