Problem
Let $k$ be an odd integer. Let $A$ be a matrix such that $A^{k} = A+ I$. Prove that $\det(A) > 0$.
My attempt:
We know that $\det(A) = \lambda_1\dots\lambda_n$ where $\lambda_1,\dots,\lambda_n$ are the complex eigenvalues of $A$. We know that $\lambda_1,\dots,\lambda_n$ are roots of $X^{2k+1} - X -1$. I can see graphically that the only real root of $X^{2k+1} -X -1$ is positive by plotting some graphs, and since the product of two conjugate complex numbers is positive, I think it has something to do with it.
The issue is that if $\lambda$ is an eigenvalue of $A$, I'm not sure that $\bar{\lambda}$ is an eigenvalue as well since $A$ possibly has complex coefficients (it's not specified in the statement).
Let $n$ be such that $A$ is an $n×n$ matrix. Let us also assume that $k$ is positive.
Then $$\det(A) =\prod_i \lambda_i^{c_i},$$ where the $\lambda_i$s are the algebraic eigenvalues of $A$ and $c_i$, for each such integer $i$, is the algebraic multiplicity of $\lambda_i$ as an algebraic eigenvalue of $A$. Then each $\lambda_i$ is the root of the polynomial $X^{2k+1}-X-1$. The $c_i$s as above are positive integers that sum to $n$.
The challenge of the problem, however, is this: We do not know a priori precisely which subset of the roots of the equation $X^{2k+1}-X-1$, are the algebraic eignevalues of $A$, and furthermore, with which multiplicity does it appear as an algebraic eigenvalue of $A$. For example, the real root of the polynomial $X^{2k+1}-X-1$ may be the one algebraic eignvalue of $A$, with multiplicity $n$.
However, if $A$ is real, then one can look at the polynomial $\det(A-\lambda I)$ [a real polynomial in $\lambda$] and note the following:
If $\lambda_i$ is a complex non real algebraic eigenvalue of $A$ with multiplicity $c_i$ for some positive intrger $c_i$, then $\bar{\lambda}_i$ is also an algebraic eigenvalue of $A$, with the same multiplicity $c_i$. However, for each such $c_i$ note that $\lambda_i^{c_i}\bar{\lambda}_i^{c_i}$ is $|\lambda_i|^{2c_i}$ which is necessarily positive.
So from 1. above it follows that, assuming $A$ is real, that $\det(A)$ must be of an integral power of the real root $y$ of the polynomial $X^{2k+1}-X-1$, times a positive number [make sure you see why]. More precisely, $\det(A)$ must be of the form
$$\det(A) = y^c \prod_i \lambda_i^{c_i}\bar{\lambda}_i^{c_i}$$ $$= \ y^c \prod_i |\lambda_i|^{2c_i}.$$ But as already observed in the OP, $y$ is positive. So from this it follows that $\det(A)$ must also be positive, if $A$ is real.
If $A$ is allowed to not be real then the claimed result does not always hold. Indeed, let $A$ be a $1×1$ matrix where the one entry of $A$ is a complex non-real root of $X^{2k+1}-X-1$. Then $\det(A)$ is simply that nonreal root, which is not positive.