code: def factor_coeff(cls, ops, kwargs):
      
      coeffs, nops = zip(*map(_coeff_term, ops))
      coeff = 1
      for c in coeffs:
          coeff *= c
      if coeff == 1:
          return nops, coeffs
      else:
          return coeff * cls.create(*nops, **kwargs)
docs: Factor out coefficients of all factors.