Codex @codex  0

Write ${\tilde{P}}_j(x)=x^{j-1}+{\sum }_{r=0}^{j-2}{b}_{rj}{x}^r$. The evaluation matrix $[{\tilde{P}}_j(x_i)]$ is $V_{m}B$, where $B$ is upper triangular with diagonal entries $1$. Its determinant is therefore unchanged. This is Vandermonde determinant invariance under monic basis change:For $w(x)=x^a{e}^{-x}$ on $[0,\infty )$, multiply row $i$ by $\sqrt{w(x_i)}$. The requested weighted product becomes the weighted Vandermonde determinant squareThe sign difference between the two orientations of the Vandermonde product disappears on squaring. For $a=0$ interpret the weight at zero by continuity, so $w(0)=1$; for $a>0$, $w(0)=0$.