Other way, more elementary in my opinion: because they are polynomials it is enough to decompose the polynomials in it factors and, due to the fact that the determinant is a multiplicative function, then we can achieve the same result, that is
\( \displaystyle \det(A^2/2+A+I)=\frac1{2^3}\det(A^2+2A+2I)=\frac1{8}\det((A-\lambda I)(A-\bar\lambda I))=\frac18\det(A-\lambda I)\det(A-\bar\lambda I)=\frac18\lambda^3\bar\lambda^3=\frac{|\lambda|^6}8=1 \)
for \( \lambda \) and \( \bar\lambda \) the two roots of the polynomial \( x^2+2x+2 \), that is \( \lambda=-1+i \) and \( \bar\lambda =-1-i \), then clearly \( |\lambda|^2=\lambda\cdot\bar\lambda=2 \) and the result follows.
The important step is to note that \( \det(A-\lambda I)=-\lambda ^3 \) because a nilpotent matrix can be transformed, under a change of basis, in an upper-triangular matrix where the diagonal is full of zeros. A change of basis doesnt change the value of the determinant so we can see \( A \) in the form that we please.