code: def condensed(network, state):
      
      result = []
      covered_nodes = set()
  
      for c in reversed(sorted(complexes(network, state))):
          if not any(n in covered_nodes for n in c.subsystem.node_indices):
              result.append(c)
              covered_nodes = covered_nodes | set(c.subsystem.node_indices)
  
      return result
docs: Return a list of maximal non-overlapping complexes.
  
      Args:
          network (Network): The |Network| of interest.
          state (tuple[int]): The state of the network (a binary tuple).
  
      Returns:
          list[SystemIrreducibilityAnalysis]: A list of |SIA| for non-overlapping
          complexes with maximal |big_phi| values.