Abstract
Fermi-Dirac and Bose-Einstein functions arise as quantum statistical distributions. The Riemann zeta function and its extension, the polylogarithm function, arise in the theory of numbers. Though it might not have been expected, these two sets of functions belong to a wider class of functions whose members have operator representations. In particular, we show that the Fermi-Dirac and Bose-Einstein integral functions are expressible as operator representations in terms of themselves. Simpler derivations of previously known results of these functions are obtained by their operator representations.