Efficient Computation of Counterfactual Bounds