A Four-State HLL Riemann Solver for Numerical Simulation of Magneto-Hydrodynamics Based on the Least Squares Solution for the Middle Wave