Surrogate Models for Diffusion on Graphs: A High-Dimensional Polynomial Approach