Background
I'm reading some notes in multivariate data analysis, in particular factor analysis.
A data vector $X_{p\times 1}$, with $E(X) = \mu$
A vector $F_{m \times 1} $ of factors,
A matrix $L_{p\times m}$ of loadings, and
A vector $\varepsilon_{p \times 1}$ of $p$ errors with a diagonal covariance matrix $Var(\varepsilon) = \Psi$
which gives the model
$$ X-\mu = LF + \varepsilon $$
The model has $mp + p$ parameters after considering $L$ and $\Psi$.
The solutions are not unique, and are only determined up to orthogonal rotations $L^* = LT$ where $T$ is an orthogonal matrix. This is used profitably to rotate the factors in a way that provides better interpretation.
Question
Then, the notes say:
after rotating with $T$, there are $\frac{m(m-1)}{2}$ fewer parameters.
I just can't figure out where that comes from. Can someone please explain?