Numbers whose prime factors all have odd exponent Let $S$ denote the set of natural numbers $m$ with the property that for all prime powers $p^k || m$ we have $k \equiv 1 \pmod{2}$. 
What is the asymptotic density of $S$?
Note that $S$ contains all prime numbers and more generally, all square-free numbers, so that $\liminf_{X \rightarrow \infty} \frac{\# (S \cap [1,X])}{X} \geq \frac{6}{\pi^2}$. 
 A: The Dirichlet series $L (s)=\sum_{n \in S} n^{-s} $ has the Euler product expression
$$L (s) = \prod_p \frac{1+p^{-s}-p^{-2s}}{1-p^{-2s}}.$$
So $S $ has analytic density
$$\mu (S) = \operatorname{Res}_{s=1} L(s) = \lim_{s \to 1} \frac{L(s)}{\zeta(s)} = \prod_p \frac{1+1/p-1/p^2}{1+1/p}$$
as can also be seen by an informal argument. We can also rewrite this as
$$ \mu(S) = \prod_p \bigl(1-\frac{1}{p(p+1)}\bigr) = \prod_p \frac{1-2/p^2+1/p^3}{1-1/p^2}$$
I don't expect this Euler product to have a nice expression since the roots of the Euler factor $1+X-X^2$ in the numerator are not roots of unity, so the product cannot be expressed, at least in an obvious way, in terms of values of the Riemann zeta function.
A: As I see it, the desired set is the disjoint union of (the rad(n) multiples of) squarefree numbers times n^2 as n ranges over all positive integers, which is like (zeta(3) + 1/64 + 3/512 + 2/729 + 1/1728 + 7/4096 + ...)/zeta(2) for the density. 
Gerhard "Doing This On Minimal Caffeine" Paseman, 2018.10.04.
