Scalable Gaussian Processes on Discrete Domains