I would take a different approach altogether and show that for each $n\in\Bbb N$ there is a unique $\langle i,j\rangle\in\Bbb N\times\Bbb N$ such that $f(i,j)=n$ by actually calculating $i$ and $j$ from $n$. Here’s a suggestion for how this can be done.
First note that for each $k\ge 2$ there are $k-1$ pairs $\langle i,j\rangle$ of positive integers such that $i+j=k$. Then observe that
$$\frac{(i+j-2)(i+j-1)}2=\sum_{k=1}^{i+j-2}k=\sum_{k=2}^{i+j-1}(k-1)$$
is the number of pairs $\langle u,v\rangle$ of positive integers such that $u+v<i+j$. Now we consider the $i+j-1$ pairs $\langle u,v\rangle$ such that $u+v=i+j$: they are
$$\langle i+j-1,1\rangle,\langle i+j-2,2\rangle,\ldots,\langle 1,i+j-1\rangle\;.$$
The pair $\langle i,j\rangle$ is the $j$-th pair in that list. In other words, $f(i,j)$ is the number of pairs $\langle u,v\rangle$ of positive integers such that $u+v<i+j$, or $u+v=i+j$ and $v\le j$. (It’s the position of $\langle i,j\rangle$ in a diagonal enumeration of $\Bbb N\times\Bbb N$.)
Now work backwards. Given a positive integer $n$, show that there is a unique $m$ such that
$$\frac{m(m+1)}2=\sum_{k=1}^mk<n\le\sum_{k=1}^{m+1}k=\frac{(m+1)(m+2)}2\;.$$
Then show that if $f(i,j)=n$, then necessarily $i+j=m+2$ and $j=n-\frac{m(m+1)}2$.
If you get stuck, take a look at Wikipedia’s discussion of the pairing function; its version differs slightly from yours, because its version is a pairing function for the non-negative integers rather than for the positive integers – the $\Bbb N$ there includes $0$, while yours does not – but the principle is exactly the same, and the sketch should be helpful.