Today I was solving an exercise and one of the things I tried (which later turned out to be useless) involved considering the following:
Is there a simple way to describe in terms of $n$ the range of the inverse harmonic mean of two integers?
Let $m, k \in \mathbb N$, $m,k \le n$. Consider
$$\frac 1m + \frac 1k = \frac 2{\mathcal H(m,k)}$$
For which $f: \mathbb N \mapsto \mathbb Q$ can I find $m,k$ such that
$$\frac 1m + \frac 1k = \frac 2{\mathcal H(m,k)} = f(n)$$?
So you are looking for a function $f:\mathbb{N}\rightarrow\mathcal{P}(\mathbb{Q})$,where $\mathcal{P}(\mathbb{Q})$ is the power set of the rationals, rather than a function $f:\mathbb{N}\rightarrow\mathbb{Q}$. The way you are posing the problem causes some confusion.
$\mathcal H(m,k)$ obtains its min value for $m=k=1$ and then $\mathcal H(1,1)=1$ and its max value for $m=k=n$ and then $\mathcal H(n,n)=n$. So, generally, $$f(n)=\left\{\frac{2mk}{m+k}|k,m=1,2,...n\right\}\subset\mathbb{Q}$$