Abstract In this paper we consider harmonic functions on gradient shrinking Ricci solitons with constant scalar curvature. A Liouville theorem is proved without using gradient estimate: any bounded harmonic function is constant on gradient shrinking Ricci solitons with constant scalar curvature. As an application, we show that the space of harmonic functions with polynomial growth has finite dimension.