I want to achieve this kind of simplification: e+ac+ad+bc+bd = e+(a+b)(c+d) . None of SymPy simplification functions worked this way. Is there any other method
e+ac+ad+bc+bd = e+(a+b)(c+d)
You can use collect(expr, e, func=factor).
collect(expr, e, func=factor)
In [5]: expr = e + a*c + a*d + b*c + b*d In [6]: collect(expr, e, func=factor) Out[6]: e + (a + b)⋅(c + d)