Prove $$\color{blue}{ \Delta=\begin{vmatrix} (y+z)^2&xy&zx\\ xy&(x+z)^2&yz\\ xz&yz&(x+y)^2 \end{vmatrix}=2xyz(x+y+z)^3} $$ using elementary operations and the properties of the determinants without expanding.
My Attempt $$ \Delta\stackrel{C_1\rightarrow C_1+C_2+C_3}{=}\begin{vmatrix} xy+y^2+yz+zx+yz+z^2&xy&zx\\ x^2+xy+xz+xz+yz+z^2&(x+z)^2&yz\\ x^2+xy+xz+xy+y^2+yz&yz&(x+y)^2\\ \end{vmatrix}\\ =\begin{vmatrix} y(x+y+z)+z(x+y+z)&xy&zx\\ x(x+y+z)+z(x+y+z)&(x+z)^2&yz\\ x(x+y+z)+y(x+y+z)&yz&(x+y)^2\\ \end{vmatrix}\\ =(x+y+z)\begin{vmatrix} y+z&xy&zx\\ x+z&(x+z)^2&yz\\ x+y&yz&(x+y)^2\\ \end{vmatrix}\\ \stackrel{R_1\rightarrow R_1+R_2+R_3}{=} \\(x+y+z)\begin{vmatrix} 2(x+y+z)&x^2+xy+xz+xz+yz+z^2&x^2+xy+zx+xy+y^2+yz\\ x+z&(x+z)^2&yz\\ x+y&yz&(x+y)^2\\ \end{vmatrix}\\ =(x+y+z)^2\begin{vmatrix} 2&x+z&x+y\\ x+z&(x+z)^2&yz\\ x+y&yz&(x+y)^2\\ \end{vmatrix}\\ =2(x+y+z)^2\begin{vmatrix} 1&\frac{x+z}{2}&\frac{x+y}{2}\\ x+z&(x+z)^2&yz\\ x+y&yz&(x+y)^2 \end{vmatrix} \stackrel{R_2\rightarrow R_2-(x+z)R_1|R_2\rightarrow R_3-(x+y)R_1}{=} 2(x+y+z)^2\begin{vmatrix} 1&\frac{x+z}{2}&\frac{x+y}{2}\\ 0&\frac{(x+z)^2}{2}&\frac{yz-x^2-xy-xz}{2}\\ 0&\frac{yz-x^2-xy-xz}{2}&\frac{(x+y)^2}{2} \end{vmatrix} $$ How do I extract $x+y+z$ or $xyz$ from the determinant to find the solution ?
Note: In a similar problem How to solve this determinant there seems to be solutions talking about factor theorem and polynomials. I am basically looking for extracting the terms $(x+y+z)$ or $xyz$ from the given determinant to solve it only using the basic properties of determinants.
$$ \Delta=\begin{vmatrix} (y+z)^2&xy&zx\\ xy&(x+z)^2&yz\\ xz&yz&(x+y)^2 \end{vmatrix}=xyz\begin{vmatrix} \frac{(y+z)^2}{x}&x&x\\ y&\frac{(x+z)^2}{y}&y\\ z&z&\frac{(x+y)^2}{z} \end{vmatrix}= xyz\begin{vmatrix} \frac{(y+z)^2}{x}-x&x&0\\ y-\frac{(x+z)^2}{y}&\frac{(x+z)^2}{y}&y-\frac{(x+z)^2}{y}\\ 0&z&\frac{(x+y)^2}{z}-z \end{vmatrix}= xyz\begin{vmatrix} \frac{(y+z)^2-x^2}{x}&x&0\\ \frac{y^2-(x+z)^2}{y}&\frac{(x+z)^2}{y}&\frac{y^2-(x+z)^2}{y}\\ 0&z&\frac{(x+y)^2}{z} \end{vmatrix}= xyz\begin{vmatrix} \frac{[(y+z)-x][(y+z)+x]}{x}&x&0\\ \frac{[y-(x+z)][y+(x+z)]}{y}&\frac{(x+z)^2}{y}&\frac{[y-(x+z)][y+(x+z)]}{y}\\ 0&z&\frac{[(x+y)-z][(x+y)+z]}{z} \end{vmatrix}= xyz(x+y+z)^2\begin{vmatrix} \frac{[(y+z)-x]}{x}&x&0\\ \frac{[y-(x+z)]}{y}&\frac{(x+z)^2}{y}&\frac{[y-(x+z)]}{y}\\ 0&z&\frac{[(x+y)-z]}{z} \end{vmatrix}= (x+y+z)^2\begin{vmatrix} {[(y+z)-x]}&x^2&0\\ {[y-(x+z)]}&{(x+z)^2}&{[y-(x+z)]}\\ 0&z^2&{[(x+y)-z]} \end{vmatrix} =(x+y+z)^2\begin{vmatrix} {[(y+z)-x]}&x^2&0\\ {[y-(x+z)]}&{(x+z)^2}&{[y-(x+z)]}\\ 0&z^2&{[(x+y)-z]} \end{vmatrix}= (x+y+z)^2\begin{vmatrix} {[(y+z)-x]}&x^2&0\\ -2z&2xz&-2x\\ 0&z^2&{[(x+y)-z]} \end{vmatrix}= (x+y+z)^2\begin{vmatrix} {(y+z)}&x^2&\frac{x^2}{z}\\ 0&2xz&0\\ \frac{z^2}{x}&z^2&{(x+y)} \end{vmatrix}\\ =(x+y+z)^2.2xz.\Big[(x+y)(y+z)-\frac{x^2.z^2}{xz}\Big]\\ =(x+y+z)^2(2xz)(xy+xz+y^2+xz-xz)\\ =(x+y+z)^2(2xz)(xy+y^2+xz)\\=(x+y+z)^2.2xz.y(x+y+z)=2xyz(x+y+z)^2 $$