The gravitational wave (GW) detection of a binary neutron star inspiral made by the Advanced LIGO and Advanced Virgo paves the unprecedented way for multimessenger observations. The propagation speed of this GW can be scrutinized by comparing the arrival times between GW and neutrinos or photons. It provides the constraint on the mass of the graviton. f(R) gravity theories have the habitual non-zero mass gravitons in addition to usual massless ones. Previously, we show that the model independent f(R) gravity theories can be constructed from the both background evolution and the matter growth with one undetermined parameter. We show that this parameter can be constrained from the graviton mass bound obtained from GW detection. Thus, the GW detection provides the invaluable constraint on the validity of f(R) gravity theories.The action of the general f(R) gravity theories are given bywhere κ 2 = 8π G/c 4 , f (R) is a general function of the Ricci scalar R, and L (m) is the matter Lagrangian. The so-called metric formalism gravitational field equation is obtained from the variation of action, Eq. (1) with respect to the metricwhereμν is an energymomentum tensor of the matter. The trace of the field equation is given byIn order to invoke the gravitational waves, one needs to investigate the linearized theory of f(R) in vacuum, T (m) = 0. The a e-mail: skylee@kias.re.kr linear perturbations on the metric h μν is written aswhereḡ μν is the background metric. One expands the Ricci tensor and scalar curvature up to the first order