Roundoff error is clearly not the problem. The long way is only an option because 3 is a small integer if we had
$$(-3+4\imath)^{31763}$$
or
$$(-3+4\imath)^\pi$$
we would not have "Guaranteed exact accuracy, no fractions, no decimals, no roundoff errors!"
I was surprised when I used my sine table in pi/48 increments it was not accurate enough
I needed interpolation to establish arcsin(4/5)~85pi/288
hence
(-3+4i)^3~125(cos(255pi/288)+isin(255pi/288))
Anyway in the original question we need to use quadrant knowledge as mentioned to find that
z^3 = 5^3(-cos(3⋅0.9273) + isin(3⋅0.9273))
or
z^3 = 5^3(cos(3⋅2.2143) + isin(3⋅2.2143))
One way to see this is to use the function arctan2 available on many apps and calculators
arctan2(5,-3)=pi-Arctan(4/3)~2.2143
it is a version of arctan that takes care of quadrants automatically
https://en.wikipedia.org/wiki/Atan2