In classical statistical mechanics, the Liouville's theorem tells that for a system in equilibrium the Poisson bracket $$\{\rho(q,p),H(q,p)\}=0.\tag{1}$$
1. Does it necessarily imply $\rho(q,p)=\rho(H(q,p))$? In other words, I'm not asking whether $\rho(H(q,p))$ satisfies (1). It of course does. But is it not possible to have other functional forms $\rho(q,p)$ which are nor necessarily of the form $\rho(H(q,p))$ but still satisfying (1)?
2. In the answer to question 1 is in negative, is it possible to mathematically derive or restrict the possible functional forms of $\rho(H)$? If not mathematically, can we rule out physically some functional forms of $\rho(H)$ (if not all)? After all, apparently any function $\rho(H(q,p))$ will satisfy (1).
A coherent, simplistic answer addressing all the questions above will be highly appreciated.