Python中算术表达式的因式分解

Python中是否有一种直接的方法来构建函数 factorize 例如,给定一个简单的算术表达式作为输入,只使用乘法和加法(例如。 'a*b*d+a*b*c*d+e' ), 返回给定因子列表的因式分解表达式?

例如

x = 'a*b*d+a*b*c*d+e'
y = ['a','b']
print(factorize(x,y))

会打印出类似的内容吗 '(a*b)*(d+cd)+e'

是否有一些内置功能可以帮助此类操作?