A generalized perturbation theory is derived, and a variational principle is presented, for the product of an arbitrary number of ratios of linear flux functionals and bilinear functionals in nuclear systems which are described by the homogeneous or the inhomogeneous Boltzmann equation.