For the integral
$$\int \frac{x^m}{x^n + 1} \,dx ,$$
the substitution $u = x^k$, $k \in \{2, 3, \ldots\}$, will only really be useful for simplifying this integral (more precisely, will not introduce any fractional powers) if $k \mid n$ and $m \equiv -1 \pmod k$.
If $n$ is prime, as in the example $m = 3, n = 5$, the only option is $k = n$, which means a substitution requires $m \equiv -1 \pmod n$, which is not the case for that example.
Once we've reached the point where no more substitutions are available, we can use that the denominator factors over $\Bbb R$ as
$$x^n + 1 = (x + 1) \prod_{k = 1}^{(n - 1) / 2} \left(x^2 - 2 \cos \frac{(2 k + 1)\pi }{n} \cdot x+ 1\right)$$ if $n$ is odd and as $$x^n + 1 = \prod_{k = 1}^{n / 2} \left(x^2 - 2 \cos \frac{(2 k + 1) \pi}{n} \cdot x+ 1\right)$$
if $n$ is even.
In our running example, $n = 5$, so
$$x^5 + 1 = (x + 1) \left(x^2 - 2 \cos \frac\pi5 \cdot x + 1\right)\left(x^2 - 2 \cos \frac{3\pi}5 \cdot x + 1\right) = (x + 1) (x^2 - \phi x + 1) (x^2 - \bar\phi x + 1) ,$$
where $\phi := \frac12(1 + \sqrt 5)$ is the Golden Ratio and $\bar\phi := \frac12(1 - \sqrt 5)$ is its conjugate.
So, our partial fraction decomposition has the form
$$\frac{x^3}{x^5 + 1} = \frac{A}{x + 1} + \frac{B x + C}{x^2 - \phi x + 1} + \frac{D x + E}{x^2 - \bar\phi x + 1} .$$
Solving the resulting $5 \times 5$ linear system gives
$$A = -\frac15, \quad B = E = \frac15 \phi, \quad C = D = \frac15 \bar\phi .$$