On the other hand, evaluating (1) for every pair (i, j) is not feasible as it would make the complexity of the model quadratic.
we let J i contain all indices j for which x ij > 0, as well as a random sample of indices for which x ij = 0.