Denote by $\rho(p,q)$ (the $p$ and $q$ are being used as shorthand for several degrees of freedom), the phase space probability distribution function, (so $\rho\,\text{d}p\text{d}q$ is the probability the system is in a state within the phase space volume $\text{d}p\text{d}q$ at any instant).
I understand that, by Liouville's theorem, $\rho(p, q)$ is an integral of motion, and therefore its dependence on the $p$ and $q$ can be reframed as a dependence on other integrals of motion, i.e. $$ \rho(p,q) = \rho(C_1(p,q),\dots,C_n(p,q)) $$ where $C_i$ are other integrals of motion.
In the Landau and Lifshitz Statistical Physics book (3rd ed.) pg. 11
We may therefore say that the distribution function, being a function of the mechanical invariants, is itself an integral of the motion. It proves possible to restrict very considerably the number of integrals of the motion on which the distribution function can depend. To do this, we must take into account the fact that the distribution $\rho_{12}$ for a combination of two subsystems is equal to the product of the distribution functions $\rho_1$ and $\rho_2$ of the two subsystems separately: $$ \rho_{12} = \rho_1\rho_2\,. $$ Hence $$ \log{\rho_{12}} = \log{\rho_1} + \log{\rho_2}\,, $$ so that the logarithm of the distribution function is an additive quantity. We therefore reach the conclusion that the logarithm of the distribution function must be not merely an integral of the motion, but an additive integral of the motion.
And the book then goes on to conclude that because $\log{\rho(p,q)}$ is an "additive integral of motion", it can only depend upon other additive integrals of motion, hence it must be a linear combination of the system's energy $E(p,q)$, momentum $\vec{P}(p,q)$, and angular momentum $\vec{M}(p,q)$ because these are the only independent additive integrals of motion as stated in a previous L&L book (and discussed here).
I don't understand how we arrive at the requirement that $\rho(p,q)$ must be dependent upon only "additive integrals of motion" from requiring $\log{\rho_{12}}=\log{\rho_1}+\log{\rho_2}$. (I'm also unsure of the precise mathematical definition of an additive integral of motion, since an integral of motion is specific to a particular system, so I'm not sure how one combines integrals of motions of different systems.)
Why can't we just have a system with two completely isolated subsystems, with distribution functions $\rho_1$ and $\rho_2$, where $\rho_1$ and $\rho_2$ can depend on arbitrary integrals of motion? Or even if they do depend on $E$, $\vec{P}$, and $\vec{M}$, it is in a way where $\log\rho$ is not a linear combination of such.