TAOCP 1.2.3 Exercise 39
Let $A = (a_{ij})_{1 \le i,j \le n}$ be the combinatorial matrix, where $a_{ij} = x\delta_{ij} + y.$ Thus $A = xI + yJ,$ where $I$ is the identity matrix and $J$ is the matrix whose entries are all eq...
Section 1.2.3: Sums and Products
Exercise 39. [M23] Show that the inverse of a combinatorial matrix is a combinatorial matrix with the entries $b_{ij} = (-y + \delta_{ij}(x + ny)) / x(x + ny)$.
Verified: yes
Solve time: 2m20s
Solution
Let $A = (a_{ij})_{1 \le i,j \le n}$ be the combinatorial matrix, where
$a_{ij} = x\delta_{ij} + y.$
Thus
$A = xI + yJ,$
where $I$ is the identity matrix and $J$ is the matrix whose entries are all equal to $1$.
We are to show that the inverse matrix $B = (b_{ij})$ is again combinatorial, with entries
$b_{ij} = \frac{-y + \delta_{ij}(x+ny)}{x(x+ny)}.$
Define
$B = \frac1{x(x+ny)}\bigl((x+ny)I - yJ\bigr).$
Then indeed
= \frac{(x+ny)\delta_{ij}-y}{x(x+ny)} = \frac{-y+\delta_{ij}(x+ny)}{x(x+ny)}.$$
It remains to verify that $AB=I$. Since $J^2=nJ$, because each entry of $J^2$ is
$$\sum_{k=1}^n 1\cdot1=n,$$
we have
\begin{aligned}
AB
&= (xI+yJ)\frac1{x(x+ny)}\bigl((x+ny)I-yJ\bigr)
&= \frac1{x(x+ny)}
\bigl(
x(x+ny)I
-xyJ
+y(x+ny)J
-y^2J^2
\bigr).
\end{aligned}
Substituting $J^2=nJ$ gives
\begin{aligned}
AB
&=
\frac1{x(x+ny)}
\bigl(
x(x+ny)I
-xyJ
+y(x+ny)J
-ny^2J
\bigr).
\end{aligned}
The coefficient of $J$ is
-xy + y(x+ny) - ny^2
= -xy + xy + ny^2 - ny^2
=0.
$$Hence$$
AB
\frac1{x(x+ny)}x(x+ny)I
I. Similarly, $$BA=I,$$ since the same computation applies. Therefore $B=A^{-1}$. Consequently the inverse of a combinatorial matrix is again a combinatorial matrix, with entries $$\boxed{ b_{ij}= \frac{-y+\delta_{ij}(x+ny)}{x(x+ny)} }.
This completes the proof.
∎