A continuación es bastante brutal enfoque : cada vez que veo los inversos de los coeficientes binomiales, trato de usar la siguiente relación entre la Beta y Gamma funciones (me puede proporcionar un enlace si es necesario) :
$$\mbox{With } 0 < m \le n,\quad \quad \ \frac{1}{\binom{m+n}{m}} = \frac{mn}{m+n} \cdot \displaystyle{\int_0^1} t^{m-1}(1-t)^{n-1} dt$$
Su suma se convierte en :
$$S_n = 1 + \sum \limits_{k=1}^{n^2} (-1)^k \binom{n^2}{k} \frac{kn}{k+n} \displaystyle{\int_0^1} t^{k-1}(1-t)^{n-1} dt$$
$$S_n = 1 + n \displaystyle{\int_0^1} (1-t)^{n-1} \cdot \Big( \sum \limits_{k=1}^{n^2} \frac{(-1)^k k}{k+n} \binom{n^2}{k} t^{k-1} \Big) dt$$
Breve interludio : podemos simplificar esta suma. Denotar $f : t \mapsto t^{n+1} \sum \limits_{k=1}^{n^2} \frac{(-1)^k k}{k+n} \binom{n^2}{k} t^{k-1}$. A continuación,$f'(t) = \sum \limits_{k=1}^{n^2} (-1)^k k \binom{n^2}{k} t^{k+n-1} = t^n \frac{d}{dt} \big( \sum \limits_{k=0}^{n^2} \binom{n^2}{k} (-t)^k \Big) = t^n \frac{d}{dt} \Big( (1-t)^{n^2} \Big) = -n^2 t^n (1-t)^{n^2-1}$.
Por lo tanto $$S_n = 1 + n \displaystyle{\int_0^1} (1-t)^{n-1} \cdot \frac{1}{t^{n+1}} \displaystyle{\int_0^{t}}-n^2 s^n (1-s)^{n^2-1} ds dt$$
y por el intercambio de las integrales :
$$S_n = 1 - n^3 \displaystyle{\int_0^1} s^n (1-s)^{n^2-1} \cdot \displaystyle{\int_s^1} \frac{(1-t)^{n-1}}{t^{n+1}} dtds$$
$$S_n = 1 - n^3 \displaystyle{\int_0^1} s^n (1-s)^{n^2-1} \cdot \frac{(1-s)^n}{n s^n}ds$$
$$S_n = 1 -n^2 \displaystyle{\int_0^1} (1-s)^{n^2+n-1} = 1 - \frac{n^2}{n^2+n}$$
Puede (debe) ser una forma más elegante, pero con esto podemos concluir que el $S_n = \frac{1}{n+1}$.