A scalable multigrid method for solving indefinite Helmholtz equations with constant wave numbers