The book on Kinetic theory I'm reading derives the BBGKY hierarchy after introducing the reduced distribution functions
$f_s(q^1,p_1,q^2,p_2,\dots,q^s,p_s):=\int\ \rho\ \ \mathrm d q^{s+1} \mathrm d p_{s+1}\cdots \mathrm d q^N \mathrm d p_N,$
where $\rho$ is the phase space probability function. In fact, in the derivation, to get rid of most terms, he reasonably states that we must make the assumption that the $f_s$ are symmetric, i.e.
$\forall k.\ f_s(\dots,q^k,p_k,q^{k+1},p_{k+1},\dots)=f_s(\dots,q^{k+1},p_{k+1},q^k,p_k,\dots).$
The step from going from the particle trajectory $\pi$ (which solves the Hamiltonian equations of motions) to the distribution $\rho$ (which solves the continouty equation with the vector field, generated by multiple $\pi$-streamlines, as flow) can be a bit confusing. The fact that the author concerns himself two times with this and doesn't just say "we can clearly assume that.." suggests to me that it might not be so simple to put asside. He offers a symmetrization method $f_s\mapsto \frac{1}{s!}\sum_{\pi}f_s$, where the sum permutes the arguments.
Now firstly, I wonder how this works with multi-particle distribution functions. What do the equations look like for an $f_{s,t}$ which describe $s$ particle of one sort and $t$ particles of another sort?
Also, as this implies even the Boltzmann equation rests on this assumptions which therefore apparently isn't initially part of the Hamiltonian description - what information is lost in the process of symmetrization?