Integer linear programming formulation for the unified duplication-loss-coalescence model