Вы можете разбить уравнение на два неравенства и применить метод big-M:
ui + qj <= uj + M(1-xij)
ui + qj >= uj - M(1-xij)
Модели с константами big-M обычно бывают слабыми и численно нестабильными, поэтому я предлагаю выбрать константукак можно меньше (т. е. сделать так, чтобы M зависело от ij, если это возможно).Чтобы узнать больше об этом, взгляните на Опасности "Big M" .