Make your matrix
mm[n_] := Table[If[i == 1, x^j, a[i - 1, i]^j], {i, 1, n}, {j, 0, n - 1}]
mm[3] // MatrixForm
and calculate its determinant
det[m_] := Product[
(m[[j, 2]] - m[[i, 2]]), {i, 1, Length[m] - 1}, {j, i + 1, Length[m]}]
Example
det[mm[5]] - Det[mm[5]] // Expand
The det-method is (of course) considerably faster
Timing[det[mm[10]];]
Timing[Det[mm[10]];]