An optimal control problem governed by an elliptic variational inequality is studied. The feasible set of the problem is relaxed, and a pathfollowing-type method is used to regularize the constraint on the state variable. First order optimality conditions for the relaxed-regularized subproblems are derived, and convergence of stationary points with respect to the relaxation and regularization parameters is shown. In particular, C- and strong stationarity as well as variants thereof are studied. The subproblems are solved by using semismooth Newton methods. The overall algorithmic concept is provided, and its performance is discussed by means of examples, including problems with bilateral constraints and a nonsymmetric operator.