Let $x$ and $y$ be given real numbers. We may suppose that $2\leqslant x \leqslant y$ and that $u:= \log(y)/\log(x)$ remains bounded in a compact set away from $1$ as $x,y\to\infty$. An integer $n$ is said to be $x$-rough if all its prime factors are at least $x$. Let us denote the set of all $x$ rough numbers as $R_x$.
I am interested in obtaining asymptotic formula for the sum $$ \sum_{n\leqslant y,\ n\in R_x} \frac{1}{n}. $$
In fact, for my purpose, an upper bound of the correct order of magnitude should suffice. A naive upper bound maybe obtained as follows; From Buchstab's theorem (see Chapter 7 of Montgomery-Vaughan) states that there are $\sim w(u)y/\log(x)$ integers $n\leqslant y$ such that $n\in R_x$. Here $w(u)$ is a weight function depending only on $u$. Therefore
$$ \sum_{n\leqslant y,\ n\in R_x} \frac{1}{n}\leqslant 1+ \sum_{n=x}^{x + w(u) y/\log(x)} \frac{1}{n} \ll (u-1)\log(x) - \log\log(x) + \mathcal{O}_u(1). $$
I am inclined to believe that this is not the correct order of magnitude. For example if $u\leqslant 2$, then the sum is over primes and the growth is $\sim \log\log(x)$.
The analogous question regarding smooth numbers seems to have received a lot of attention (see this question and the references therein).
Any help in this regard is greatly appreciated.