Yo, I couldn't solve this exercise after thinking for a while.
For every $A \in GL_{3} (\mathbb{C})$ and $n$, there's a $B \in Mat_{3, 3}(\mathbb{C})$ such that $B^n = A$
The previous exercise was that for every nilpotent $N \in Mat_{3, 3} (\mathbb{C})$ and every $n$, $C = 1 + \frac{1}{n}N + \frac{1-n}{2n^2}N^2$ satisfies $C^n = 1 + N$, so I suppose there's a trick using this result.
I tried to play a little with the splitting of of $A$ as a nilpotent plus a semisimple, however I couldn't get anything useful.
Thanks in advance.
If you write $A = D + N$ as semisimple + nilpotent (where $D$ and $N$ commute), then $D$ is invertible and $$A = D(I + D^{-1}N),$$ where $D^{-1}N$ is nilpotent (because $D$ and $N$ commute). Now $D$ has a $n$th root (because we are in $\mathbb{C}$ so it's diagonalizable), and so does $I + D^{-1}N$ by the previous exercise. The product of these two $n$th roots is your desired $B$ (because they commute).
More generally, for matrices of any size, you can put $B = \exp(\tfrac1n \log A)$. Here, $\exp$ is defined by the usual power series, and $\log A$ is any matrix such that $\exp(\log A) = A$. If $A$ is invertible then this exists. Indeed, with $A = D(I + D^{-1}N)$ as above, then $\log D$ exists (clearly we can take the logarithm of any invertible diagonal matrix), and for the other factor we can use the power series $$\log (I+X) = \sum_{k=1}^\infty \frac{(-1)^{k-1}}{k} X^k$$ for $X = D^{-1}N$. The power series converges whenever the spectral radius of $X$ is $< 1$; in particular, it converges (after finitely many terms) when $X$ is nilpotent. Then $$\log A = \log D + \sum_{k=1}^\infty \frac{(-1)^{k-1}}{k} (D^{-1}N)^k.$$ It's worth showing that the $C$ from your previous exercise is just $C = \exp(\tfrac1n \log(I+N))$.
Actually, it's even easier to derive $C$ using the binomial series $$ (I+X)^\alpha = \sum_{k=0}^\infty \binom{\alpha}{k} X^k. $$