Lagrangian time-discretization of the Hunter–Saxton equation