Sobre el determinante de a $A$, es un caso particular de algo que se llama Cauchy doble alternante.
En general, dados cualesquiera $2n$ números de $x_1, \ldots, x_n, y_1, \ldots, y_n$. Si uno
la construcción de un $n \times n$ matriz con la entrada $\frac{1}{x_i + y_j}$ $i^{th}$ fila $j^{th}$ columna, uno tiene
$$\det\left[ \frac{1}{x_i + y_j}\right] = \frac{\prod_{1\le i < j\le n}(x_i - x_j)(y_i - y_j)}{\prod_{1\le i, j \le n}(x_i + y_j)}\tag{*1}$$
En particular, esto significa que si $x_i$ $-y_j$ $2n$ distintos números, el determinante no es cero y por lo tanto la matriz es invertible.
No es difícil demostrar $(*1)$ nosotros mismos, considere la posibilidad de la expresión
$$\prod_{1\le i, j \le n}(x_i + y_j)\times\det\left[ \frac{1}{x_i + y_j}\right]$$
Si uno de ampliar el determinante, es fácil ver que esta expresión es un polinomio en a $x_i, y_j$. Desde el determinante se desvanece cuando $x_i$ es igual a algunos $x_j$ o $y_i$ es igual a algunos $y_j$. Este polinomio contiene $\prod_{1\le i < j\le n}(x_i - x_j)(y_i - y_j)$ como un factor. Por coincidencia el grado del polinomio en ambos lados, nos encontramos con
$$\det\left[ \frac{1}{x_i + y_j}\right] = \lambda_n \frac{\prod_{1\le i < j\le n}(x_i - x_j)(y_i - y_j)}{\prod_{1\le i, j \le n}(x_i + y_j)}$$
para algunas constantes $\lambda_n$ sólo depende de $n$.
Si se sustituye $x_1$ por $\epsilon x_1$, $y_1$ por $\epsilon y_1$, envíe $\epsilon$ $0$y un vistazo a la limitación de los comportamientos de ambos lados, se puede concluir $\lambda_n = \lambda_{n-1}$. Uno puede comprobar $\lambda_1 = \lambda_2 = 1$ a mano y, por tanto, $\lambda_n = 1$ todos los $n$.
Aplicar esto a la de Hilbert como matriz $A$ a mano, obtenemos
$$\det\left[\frac{1}{i+j}\right] = \frac{(c_n c_{n+1})^2}{c_{2n+1}} \ne 0
\quad\text{ donde }\quad c_n = \prod_{k=1}^{n-1} k!
$$
y, por tanto, $A$ es invertible.