How to evaluate $$\det\begin{bmatrix}1+x_1y_1&\cdots & 1+x_1y_n\\ \vdots & \ddots & \vdots \\ 1+x_ny_1 & \cdots & 1+x_ny_n\end{bmatrix}$$?
I tried factorisation but I have no idea. I could not factor anything out... I need a hint to kick off.
Edit:
A conjecture due to multilinearity: if we let $v_j=(x_1y_j,\cdots,x_ny_j)^T$, then my guess would be
$$\det\begin{bmatrix}1+x_1y_1&\cdots & 1+x_1y_n\\ \vdots & \ddots & \vdots \\ 1+x_ny_1 & \cdots & 1+x_ny_n\end{bmatrix}=\det(1,v_2,\cdots,v_n)+\det(v_1,1,\cdots,v_n)+\cdots+\det(v_1,v_2,\cdots,1)=\color{red}{0}$$
where $1$ denote $(1,\cdots,1)^T\in \Bbb K^n$,$\Bbb K$ the underlying number field.
Another way of seeing that the determinant is $0$ for $n>2$ is as follows:
The determinant is a polynomial in $x_1$. If $n\geq 3$ then $x_2$ and $x_3$ are roots of this polynomial. (Since if $x_1=x_i$ then we have two rows which are equal)
Thus $(x_1-x_2)(x_1-x_3)$ is a factor of the determinant.
But we can also see that the determinant should be linear in $x_1$ since $x_1$ only occurs in the first row.
This implies that the determinant must be $0$. (Otherwise we would have a linear polynomial with a quadratic factor)