Solution (source code)

= Solution

Each term is a <Hermitian> rank-one <projection> and $p_j$ is real, so $\rho^\dagger=\rho$. Using the rank-one <trace>,
$$
\operatorname{Tr}\rho=\sum_jp_j\langle\psi_j|\psi_j\rangle
=\sum_jp_j=\boxed{1}.
$$
The remaining admissibility condition follows from $\langle v|\rho|v\rangle=\sum_jp_j|\langle\psi_j|v\rangle|^2\ge0$: \b[$\rho$ is <positive semidefinite>].