In the kinetic theory of dense fluids the many-particle collision bracket integral is given in terms of a classical collision operator defined in the phase space. To find an algorithm to compute the collision bracket integrals, we revisit the eigenvalue problem of the Liouville operator and re-examine the method previously reported [Chem. Phys. 1977, 20, 93]. Then we apply the notion and concept of the eigenfunctions of the Liouville operator and knowledge acquired in the study of the eigenfunctions to cast collision bracket integrals into more convenient and suitable forms for numerical simulations. One of the alternative forms is given in the form of time correlation function. This form, on a further manipulation, assumes a form reminiscent of the Chapman- Enskog collision bracket integrals, but for dense gases and liquids as well as solids. In the dilute gas limit it would give rise precisely to the Chapman-Enskog collision bracket integrals for two-particle collision. The alternative forms obtained are more readily amenable to numerical simulation methods than the collision bracket integrals expressed in terms of a classical collision operator, which requires solution of classical Lippmann-Schwinger integral equations. This way, the aforementioned kinetic theory of dense fluids is made fully accessible by numerical computation/simulation methods, and the transport coefficients thereof are made computationally as accessible as those in the linear response theory.