\[ \int (c x)^{-1-n-n p} \left (a+b x^n\right )^p \, dx \]
Optimal antiderivative \[ -\frac {\left (a +b \,x^{n}\right )^{1+p} \left (c x \right )^{-n \left (1+p \right )}}{a c n \left (1+p \right )} \]
command
integrate((c*x)**(-n*p-n-1)*(a+b*x**n)**p,x)
Sympy 1.10.1 under Python 3.10.4 output
\[ \frac {b^{p} c^{- n} c^{- n p} x^{- n} \left (\frac {a x^{- n}}{b} + 1\right )^{p} \Gamma \left (- p - 1\right )}{c n \Gamma \left (- p\right )} + \frac {b b^{p} c^{- n} c^{- n p} \left (\frac {a x^{- n}}{b} + 1\right )^{p} \Gamma \left (- p - 1\right )}{a c n \Gamma \left (- p\right )} \]
Sympy 1.8 under Python 3.8.8 output
\[ \text {Timed out} \]________________________________________________________________________________________