Consider a set of fractions $\left\{1, \frac{1}{2}, \frac{1}{3}, \ldots, \frac{1}{n}\right\}$. How many subsets of this set have sum at most 1? I'm interested in the asymptotics of this number.


Clearly, any subset of $\left\{\frac{1}{\lceil n/2 \rceil}, \ldots, \frac{1}{n}\right\}$ works, hence the answer is $\Omega(2^{n / 2})$. Can we show $\Omega(2^{\beta n})$ for $\beta > \frac{1}{2}$? Can we determine $\beta$ exactly? From numeric estimates of [OEIS sequence](https://oeis.org/A212657) $\beta$ seems to at least $0.88$ (link to the sequence and correction of the estimate due to Max Alekseyev).


The question arose while I was thinking about upper bounds for [this question](https://mathoverflow.net/questions/84207/asymptotic-growth-of-antichains-in-divisibility-posets?noredirect=1&lq=1). Clearly, every divisibility antichain $I$ of $[n]$ must satisfy $\sum\_{x \in I} \lfloor\frac{n}{x}\rfloor \leq n$, which is a very similar condition.


**Post-mortem**: while I accepted Lucia's answer (simply because it was the first to contain the correct answer and some reasoning to why it is correct), the whole discussion here is very valuable. Be sure to also check out js21's answer with an approach based on large deviations method, and RaphaelB4's answer for a more off-the-ground explanation of the method. In a comment Jay Pantone shared a link to a paper on series analysis, in particular, the differential approximation method allows to obtain the same answer with high precision and is, without doubt, a great practical tool. Kudos to all of you guys! What a great day to learn.



Sets of unit fractions with sum $\leq 1$
Let $n\_0$ be the smallest number such that the sum of the reciprocals of the integers from $n\_0+1$ to $n$ is $<2$. It is easy to see that $n\_0 \approx n/e^2$, since $\sum\_{j>n/e^2}^{n} 1/j \approx \log n - \log (n/e^2) 
=2$. Now for any subset $A$ of $\{n\_0 +1, \ldots, n\}$ either the sum of the reciprocals of elements in $A$ or the sum of the reciprocals of its complement must be $<1$. Therefore there are at least 
$$ 
\frac 12 2^{n-n\_0} \asymp 2^{n(1-1/e^2)}
$$ 
possible sets. My guess is that this exponent $1-1/e^2$ is correct -- note that $1-1/e^2 = 0.86466\ldots$. 


Maybe my first guess is not right! Here's an upper bound, which gives an exponent around $0.91\ldots$ (my numerical calculations are pretty rough). For any positive $x$, an upper bound on the quantity we want is 
$$ 
e^x \prod\_{j=1}^n (1+e^{-x/j}). 
$$ 
To see this, just expand out the product and terms with sum of reciprocals less than $1$ will contribute at least $1$, and the rest are positive. Now choose $x$ so as to minimize the above (a standard idea, known in analytic number theory as Rankin's trick). 


Calculus shows that one must choose $x$ so that 
$$
1= \sum\_{j=1}^{n} \frac 1j \frac{1}{1+e^{x/j}}. 
$$
It is natural to guess that $x$ is of the shape $\alpha n$ for a 
constant $\alpha$, and then for large $n$ the condition on $\alpha$ becomes 
$$ 
1= \int\_0^1 \frac{1}{1+e^{\alpha/y}} \frac{dy}{y} = \int\_1^\infty \frac{1}{1+e^{\alpha y}} \frac{dy}{y}. 
$$ 
If I calculated right, this gives $\alpha \approx 0.1273$. For this choice of $\alpha$ (and so $x$), one obtains the bound (approximately) 
$$ 
\exp\Big(n\Big( \alpha + \int\_0^1 \log (1+e^{-\alpha/y}) dy\Big)\Big), 
$$
which seems to be about 
$$ 
\exp(-.631n) \approx 2^{0.911n}.
$$ 
(I won't swear to the numerics -- someone should check.) 


My second guess is that the upper bound is tight (and I think this could be proved with some effort). The idea is to choose $j$ to be in your set with probability $1/(1+\exp(x/j))$ with the same $x$ as in the upper bound. The expected value of $1/j$ with this distribution is $1$, by the choice of $x$. An entropy calculation for this distribution then gives the exponent. (More generally, in all the situations I know, the Rankin upper bound is pretty close to optimal.)


