code: def add_cut(problem, indicators, bound, Constraint):
      
      cut = Constraint(sympy.Add(*indicators), ub=bound)
      problem.add(cut)
      return cut
docs: Add an integer cut to the problem.
  
      Ensure that the same solution involving these indicator variables cannot be
      found by enforcing their sum to be less than before.
  
      Parameters
      ----------
      problem : optlang.Model
          Specific optlang interface Model instance.
      indicators : iterable
          Binary indicator `optlang.Variable`s.
      bound : int
          Should be one less than the sum of indicators. Corresponds to P - 1 in
          equation (14) in [1]_.
      Constraint : optlang.Constraint
          Constraint class for a specific optlang interface.
  
      References
      ----------
      .. [1] Gevorgyan, A., M. G Poolman, and D. A Fell.
             "Detection of Stoichiometric Inconsistencies in Biomolecular
             Models."
             Bioinformatics 24, no. 19 (2008): 2245.