In this paper, we propose a splitting least-squares mixed finite element method for
the approximation of elliptic optimal control problem with the control constrained by pointwise
inequality. By selecting a properly least-squares minimization functional, we derive equivalent two
independent, symmetric and positive definite weak formulation for the primal state variable and
its flux. Then, using the first order necessary and also sufficient optimality condition, we deduce
another two corresponding adjoint state equations, which are both independent, symmetric and
positive definite. Also, a variational inequality for the control variable is involved. For the
discretization of the state and adjoint state equations, either RT mixed finite element or standard
C^0 finite element can be used, which is not necessary subject to the Ladyzhenkaya-Babuska-Brezzi
condition. Optimal a priori error estimates in corresponding norms are derived for the control,
the states and adjoint states, respectively. Finally, we use some numerical examples to validate
the theoretical analysis.