Sympy:如何简化总和的幂?

Sympy : How is it possible to simplify power of sum?

考虑这种形式的表达式:

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)

这是与

相同的一般问题
>>> 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)