Fast evaluation of discrete integral operators by Chebyshev and Leja polynomial approximation