Let $P_{ij}$ be variables, and let $A \in \mathbb{R}^{n\times n}$ be the matrix defined by $$ A_{ij} = \begin{cases} -P_{ij} & i \neq j,\\ P_{i1} + P_{i2} + \dots + P_{in} & i=j. \end{cases} $$ For instance for $n=3$ $$ A = \begin{bmatrix} P_{1,1}+P_{1,2}+P_{1,3} & -P_{1,2} & -P_{1,3}\\ -P_{2,1} & P_{2,1}+P_{2,2}+P_{2,3} & -P_{2,3}\\ -P_{3,1} & -P_{3,2} & P_{3,1}+P_{3,2}+P_{3,3} \end{bmatrix}. $$ Then, computationally I observe that $\det(A)$ is the sum of $(n+1)^{n-1}$ distinct terms, each of which is the product of $n$ terms $P_{ij}$. For instance for $n=3$ we get the 16 terms \begin{align*} \det(A) &= P_{1,1}\,P_{2,1}\,P_{3,1}+P_{1,1}\,P_{2,1}\,P_{3,2}+P_{1,1}\,P_{2,2}\,P_{3,1}+P_{1,1}\,P_{2,1}\,P_{3,3} \\ &+ P_{1,1}\,P_{2,2}\,P_{3,2}+P_{1,1}\,P_{2,3}\,P_{3,1}+P_{1,2}\,P_{2,2}\,P_{3,1}+P_{1,1}\,P_{2,2}\,P_{3,3} \\ &+P_{1,2}\,P_{2,2}\,P_{3,2}+P_{1,1}\,P_{2,3}\,P_{3,3}+P_{1,2}\,P_{2,2}\,P_{3,3}+P_{1,3}\,P_{2,1}\,P_{3,3} \\& +P_{1,3}\,P_{2,2}\,P_{3,2}+P_{1,2}\,P_{2,3}\,P_{3,3}+P_{1,3}\,P_{2,2}\,P_{3,3}+P_{1,3}\,P_{2,3}\,P_{3,3}. \end{align*} Is this result known? It looks like it might have a nice combinatorial interpretation, since the number of terms matches Cayley's formula for the number of trees on $n+1$ labeled nodes (https://oeis.org/A000272).