I think it's easier to focus on the roots of the polynomial, which should be fine since that tells you pretty much everything you need to know about the polynomial anyway. The characteristic polynomial of A is det(A-xI), correct? Say the roots of that are x=x
1, x
2, and x
3. Now look at the characteristic polynomial of A
2:
\det (A^2 - \lambda I) = \det \left( \left(A-\sqrt{\lambda}\right) \cdot \left(A+\sqrt{\lambda}\right) \right) = \det \left(A-\sqrt{\lambda}\right) \cdot \det \left(A+\sqrt{\lambda}\right)
But, you already know the roots of det(A-xI), so the roots of the above equation just come from the equation
\pm \sqrt{\lambda}=x \Longrightarrow \lambda = x^2
So in this case, the roots of the characteristic polynomial of A
2 would be x
12, x
22, and x
32
So if the characteristic polynomial for A looked like this:
p(x) = (x-x_1)(x-x_2)(x-x_3)
then the characteristic polynomial for A
2 would look like:
p_2(x) = (x-x_1^2)(x-x_2^2)(x-x_3^2)
In fact, the characteristic polynomial for A
n is:
p_n(x) = (x-x_1^n)(x-x_2^n)(x-x_3^n)
Now, using eigenvalues rather than determinants, this property is actually easier to derive. You know that if x is a root of the characteristic polynomial of A, then x is an eigenvalue, so there exists a non-zero vector v such that:
A\vec{v}=x\vec{v}
Left multiply both sides by A to get
A^2\vec{v}=x(A\vec{v})
But you know that Av=xv, so:
A^2\vec{v} = x(x\vec{v}) = x^2\vec{v}
That tells you that x
2 is an eigenvalue of A
2, and you can easily show by induction (just left multiply by A, n times) that x
n is an eigenvalue of A
n. Since an eigenvalue is a root of the characteristic polynomial, x
n must be a root of the characteristic polynomial of A
n. So once you know the characteristic polynomial of A, you know that the roots of the characteristic of A
n are the roots of the characteristic of A raised to the n
th power.
That may or may not translate to something nice in the p(x) = ax
3+bx
2+cx+d form (I'll let you try that...

), but it tells you something very nice about the roots which is actually true for any mxm matrix, not just a 3x3 matrix. If you can turn the characteristic of an mxm A into the form:
p(x) = (x-x_1)\cdot (x-x_2) \cdot (x-x_3) \cdots (x-x_m)
Then you immediately know that the characteristic of A
n is:
p(x) = (x-x_1^n)\cdot (x-x_2^n) \cdot (x-x_3^n) \cdots (x-x_m^n)
Hope that helps answer your question, even if I didn't answer it directly.