Sympy：如何简化求和的幂？答案

【问题标题】：Sympy : How is it possible to simplify power of sum?Sympy：如何简化求和的幂？
【发布时间】：2021-11-22 22:54:57
【问题描述】：

考虑这种形式的表达式：

x,y,n=sp.symbols("x y n",positive=True,real=True)
sp.Pow(x+y+x**2,n+1)*sp.Pow(x+2*y+4*y**3,-n-1)

如何将其简化为具有共同的权力？（即sp.Pow((x+y+x**2)/(x+2*y+y**3),n+1)）

【问题讨论】：

github.com/sympy/sympy/issues/22546
@OscarBenjamin。你认为修改 powsimp.py 中的第 393 行 (exp_c, exp_t = e.as_coeff_Mul(rational=True) to exp_c, exp_t = factor_terms(e).as_coeff_Mul(rational=True) 是解决方案？（评论也添加到您的问题中）。它适用于我的问题，但需要进行大量测试，我不知道该怎么做。

标签： sympy

【解决方案1】：

这是和here一样的普遍问题

>>> var('z', positiv=True)
z
>>> expr = sp.Pow(x+y+x**2,n+1)*sp.Pow(x+2*y+4*y**3,-n-1)
>>> powsimp(expr.subs(n + 1, var('z',positive=1))).subs(z, n + 1)
((x**2 + x + y)/(x + 4*y**3 + 2*y))**(n + 1)

【讨论】：