[leethargo]

The model could be extended by adding a penalty term to the objective. I guess the penalty would need to be nonlinear, as the current objective already minimizes the x_ij (for i != j) implictly? One could also introduce a fixed-cost offset for any nonzero redistribution amount, but that would change the problem type from linear program (LP) to mixed-integer program (MIP).