The surrogate gradient algorithm for Lagrangian relaxation method