AAnother proof by induction (on the dimension $n$ of the vector-space $X\simeq K^n$ or $\mathbb{C}^n$).
The basis step, $n=1$, is trivial.
Suppose we have proceeded up to and including dimension $n$ and suppose $X$ is an $(n+1)$$n$-dimensional $K$-vectorspace and $T:\DeclareMathOperator\End{End}\End(X) \to \End(X)$ preserves the determinant. Let $\{x_j\}_{0\leq j \leq n}$$\{x_j\}_{1\leq j \leq n}$ be a basis of $X$ and $\{x_j^*\}_{0\leq j \leq n}$ a$\{x_j^*\}_{1\leq j \leq n}$ its dual basis.
$T$ is injective and preserves rank. For rank 1 matrices, one of the implications is that $$\forall j \in \{0,\ldots,n\}:T(x_0x_j^*) = y_0 \varphi_j^* \text{ and }T(x_jx_0^*)=y_j\varphi_0^*,\\ \text{or}, \forall j \in \{0,\ldots,n\}:T(x_0x_j^*) = y_j \varphi_0^*\text{ and }T(x_jx_0^*)=y_0\varphi_j^* \qquad(1)$$$$\forall j \in \{1,\ldots,n\}:T(x_1x_j^*) = y_1 \varphi_j^* \text{ and }T(x_jx_1^*)=y_j\varphi_1^*,\\ \text{or}, \forall j \in \{1,\ldots,n\}:T(x_1x_j^*) = y_j \varphi_1^*\text{ and }T(x_jx_1^*)=y_1\varphi_j^* \qquad(1)$$ where $\{y_j\}_{0\leq j \leq n}$$\{y_j\}_{1\leq j \leq n}$ is a basis of $X$ ($\{y_j^*\}_j$ the corresponding dual basis) and $\{\varphi_j^*\}_{0\leq j \leq n}$$\{\varphi_j^*\}_{1\leq j \leq n}$ is a basis of $X^*$ ($\{\varphi_j\}_j$ the corresponding dual basis, keeping $X^{**}\simeq X$ in mind). We can restrict ourselves to the former scenario outlined in (1) by replacing $T$ in the latter scenario with $T\circ t$, where $t$ is a "coordinate-transpose", i.e. $t(x_jx_k^*):=x_kx_j^*$ and $t$ linear. If we then define the invertible matrices $U_1=\sum_{j=0}^n x_j y_j^*,\,V_1=\sum_{j=0}^n \varphi_jx_j^*$$U_1=\sum_{j=1}^n x_j y_j^*,\,V_1=\sum_{j=1}^n \varphi_jx_j^*$, then $T_1(.):=U_1T(.)V_1$ fixes both $\{x_0x_j^*\}_j$$\{x_1x_j^*\}_j$ and $\{x_jx_0^*\}_j$$\{x_jx_1^*\}_j$, i.e. $T_1$ fixes matrices whose non-zero entries occur only in the zeroth1st row and/or zerothfirst column. (Note that $T_1$ multiplies the determinant of its input by an uncertain factor $\lambda^2:=\det(U_1 V_1)\neq 0$$\det(U_1 V_1)\neq 0$, but later on the proof will indirectly show that this factor must be 1 anyway)
Let $W=\text{span}(x_1,\ldots,x_n)\subseteq X$, $\pi:X\to W: c_0x_0+c_1x_1+\ldots +c_nx_n\mapsto c_1x_1+\ldots +c_nx_n$ and $i: W \to X: w \mapsto w$. Let $\tilde{T}:\End(W) \to \End(W): A \mapsto \pi \circ T_1(i\circ A\circ \pi) \circ i$. Using the Laplace expansion for the determinant on the zeroth row or column, and using the fact that $T_1$ sends singularmaps rank 1 matrices to singularrank 1 matrices. At this stage, we seethis implies that $\forall A \in \End(W)$ $$\lambda\det(A) = \lambda\det(x_0x_0^*+i\circ A \circ \pi)=\det(T_1(x_0x_0^*+i\circ A \circ \pi))=\det(x_0x_0^*+T_1(i\circ A \circ \pi))=\det(\tilde{T}(A)).$$ So according to the induction hypothesis there exist $U_2$,$\forall j,k \in \{2,\ldots,n\}:\,\exists a_{jk}\in \mathbb{C}:\,T_1(x_jx_k^*)=a_{jk}x_jx_k^*$ $V_2$ with(test $\det(U_2)=\det(V_2)=\lambda$$T_1$ on ($\lambda$ is the square root of$(x_1+x_j)x_k^*$ and on $\lambda^2$ closest$x_j(x_1^*+x_k^*)$ to the $\mathbb{R}^+$ axisarrive at this conclusion) s, i.te. $\tilde{T}(.) \equiv U_2(.)V_2$ or $\tilde{T}(.) \equiv U_2(t(.))V_2$.
It is now easy to check$T_1$ performs a 'pointwise' multiplication on the matrix entries of its input (using onlyviewed in the property$\{x_j\}_j$-basis). To see that the below$a$-definedcoefficients all have to be equal to 1, test $T_2$ preserves$T_1$ once more on the rank 1 matrices)matrix $(x_1+x_j)(x_1^*+x_k^*)$. So we conclude that $$T_2(.):=(x_0x_0^* + i\circ U_2^{-1}\circ \pi) *T_1(.)*(x_0x_0^* +i\circ V_2^{-1}\circ \pi)$$ is $T_1$ is the identity map on $\End(X)$. (The latter scenario that was hypothesized at the end of step 4) is found to be impossible$\End{X}$ and, if I'm not mistaken retracing our steps, you should findwe see that our adjustments in step 3) have already forced $U_2=I=V_2$ a posteriori)$T(.)\equiv U_1^{-1}(.)V_1^{-1}$ or $T(.)\equiv U_1^{-1}(t(.))V_1^{-1}$.