Solution (source code)

= Solution

For
$$
|\Psi\rangle_{AR}
=\sum_i\sqrt{\lambda_i}|a_i\rangle_A|b_i\rangle_R,
$$
orthonormality of the $|b_i\rangle$ gives
$$
\operatorname{Tr}_R|\Psi\rangle\langle\Psi|
=\sum_i\lambda_i|a_i\rangle\langle a_i|
=\rho_A.
$$
Normalization follows from $\sum_i\lambda_i=1$, so this is a <purification of a density operator>.