Después de haber visto la respuesta dada por joriki,creo que es el momento de compartir mi propia respuesta(agradezco su respuesta,y necesito algo de tiempo para digerir...).
Vamos
\begin{equation}
f(x)=
\begin{vmatrix}
x^{2n+1}&x^{2n}&\cdots&x&1\\
x_1^{2n+1}&x_1^{2n}&\cdots&x_1&1\\
\vdots&\vdots&\cdots&\vdots\\
x_n^{2n+1}&x_n^{2n}&\cdots&x_n&1\\
(2n+1)x^{2n}&2nx^{2n-1}&\cdots&1&0\\
\vdots&\vdots&\cdots&\vdots&\cdots\\
(2n+1)x_{n}^{2n}&2nx_n^{2n-1}&\cdots&1&0\\
\end{vmatrix}
\end{equation}
Es fácil comprobar que $f(x)$ es un polinomio de grado $4n$.Y
\begin{equation}
f(x_1)=\cdots =f(x_n)=0
\end{equation}
Por lo $(x-x_1)(x-x_2)\cdots (x-x_n)|f(x)$.Y
\begin{align*}
f'(x)= \begin{vmatrix}
x^{2n+1}&x^{2n}&\cdots&x&1\\
x_1^{2n+1}&x_1^{2n}&\cdots&x_1&1\\
\vdots&\vdots&\cdots&\vdots\\
x_n^{2n+1}&x_n^{2n}&\cdots&x_n&1\\
2n(2n+1)x^{2n-1}&(2n-1)2nx^{2n-2}&\cdots&0&0\\
\vdots&\vdots&\cdots&\vdots&\cdots\\
(2n+1)x_{n}^{2n}&2nx_n^{2n-1}&\cdots&1&0\\
\end{vmatrix}
\end{align*}
Por lo $f'(x_1)=\cdots f'(x_n)=0$.Por lo que $(x-x_1)^2(x-x_2)^2\cdots
(x-x_n)^2|f(x)$.Y
\begin{align*}
f''(x)= \begin{vmatrix}
(2n+1)x^{2n}&2nx^{2n-1}&\cdots&1&0\\
x_1^{2n+1}&x_1^{2n}&\cdots&x_1&1\\
\vdots&\vdots&\cdots&\vdots\\
x_n^{2n+1}&x_n^{2n}&\cdots&x_n&1\\
2n(2n+1)x^{2n-1}&(2n-1)2nx^{2n-2}&\cdots&0&0\\
\vdots&\vdots&\cdots&\vdots&\cdots\\
(2n+1)x_{n}^{2n}&2nx_n^{2n-1}&\cdots&1&0\\
\end{vmatrix}+
\begin{vmatrix}
x^{2n+1}&x^{2n}&\cdots&x&1\\
x_1^{2n+1}&x_1^{2n}&\cdots&x_1&1\\
\vdots&\vdots&\cdots&\vdots\\
x_n^{2n+1}&x_n^{2n}&\cdots&x_n&1\\
(2n-1)2n(2n+1)x^{2n-2}&(2n-2)(2n-1)2nx^{2n-3}&\cdots&0&0\\
\vdots&\vdots&\cdots&\vdots&\cdots\\
(2n+1)x_{n}^{2n}&2nx_n^{2n-1}&\cdots&1&0\\
\end{vmatrix}
\end{align*}
Es fácil comprobar que $f''(x_1)=\cdots f''(x_n)=0$.Así
\begin{equation}
(x-x_1)^3(x-x_2)^3\cdots (x-x_n)^3|f(x)
\end{equation}
Y también es fácil de averiguar $f'''(x)$,por lo que es fácil comprobar que
\begin{equation}
f'''(x_1)=\cdots =f'''(x_n)=0
\end{equation}
Así
\begin{equation}
(x-x_1)^4(x-x_2)^4\cdots (x-x_n)^4|f(x)
\end{equation}
Debido a $f(x)$ es un polinomio de grado $4n$,lo $f(x)=a(x-x_1)^4(x-x_2)^4\cdots (x-x_n)^4$.De acuerdo a la simetría ammong $x_0,x_1,\cdots,x_n$ en este determinante,es fácil comprobar que
\begin{equation}
f(x_0)=c\prod_{n\geq i>j\geq 0}(x_i-x_j)^4
\end{equation}
Ahora podemos determinar la constante de $c$.Lo hacemos por inducción.Es fácil comprobar que Cuando se $n=1$,
\begin{equation} \det \begin{pmatrix} x_0^3&x_0^2&x_0&1\\ x_1^3&x_1^2&x_1&1\\ 3x_0^2&2x_0&1&0\\ 3x_1^2&2x_1&1&0\\ \end{pmatrix}=-(x_0-x_1)^4 \end{equation}Entonces cuando $n=2$,Vamos a ver el determinante
\begin{equation} \det\begin{pmatrix} x_0^5&x_0^4&x_0^3&x_0^2&x_0&1\\ x_1^5&x_1^4&*x_1^3&*x_1^2&*x_1&*1\\ x_2^5&x_2^4&*x_2^3&*x_2^2&*x_2&*1\\ 5x_0^4&4x_0^3&3x_0^2&2x_0&1&0\\ 5x_1^4&4x_1^3&*3x_1^2&*2x_1&*1&*0\\ 5x_2^4&4x_2^3&*3x_2^2&*2x_2&*1&*0\\ \end{pmatrix} \end{equation}
El elemento marcado * también en forma de un determinante,sabemos que la constante de este determinante es -1,por lo que el término constante de la determinante \begin{equation} \det\begin{pmatrix} x_0^5&x_0^4&x_0^3&x_0^2&x_0&1\\ x_1^5&x_1^4&*x_1^3&*x_1^2&*x_1&*1\\ x_2^5&x_2^4&*x_2^3&*x_2^2&*x_2&*1\\ 5x_0^4&4x_0^3&3x_0^2&2x_0&1&0\\ 5x_1^4&4x_1^3&*3x_1^2&*2x_1&*1&*0\\ 5x_2^4&4x_2^3&*3x_2^2&*2x_2&*1&*0\\ \end{pmatrix} \end{equation} es $-1\times (4-5)=1$
......
Así que por inducción,sabemos que $c=(-1)^n$.
En mi respuesta yo hacer uso de la siguiente resultado
Si \begin{equation} f(x)=\begin{vmatrix} f_{11}(x)&f_{12}(x)&\cdots&f_{1n}(x)\\ f_{21}(x)&f_{22}(x)&\cdots&f_{2n}(x)\\ \vdots&\vdots& &\vdots\\ f_{n1}(x)&f_{n2}(x)&\cdots&f_{nn}(x)\\ \end{vmatrix} \end{equation}
A continuación,\begin{equation} f'(x)=\sum_{i=1}^n \begin{vmatrix} f_{11}(x)&f_{12}(x)&\cdots&f_{1n}(x)\\ \vdots&\vdots & &\vdots\\ f'_{i1}(x)&f'_{i2}(x)&\cdots&f'_{in}(x)\\ \vdots&\vdots&&\vdots\\ f_{n1}(x)&f_{n2}(x)&\cdots&f_{nn}(x)\\ \end {vmatrix} \end{equation}