We presented some parallel NewtonKrylovSchwarz (NKS) algorithm to solving the complementarity problems. Using semismooth function, the solution of the optimization problem can be obtained by solving a large sparse nonlinear system of algebraic equations. Numerical results show that the efficiency can be achieved by the proposed method.