Given the upper Shift Matrix, which for e.g. dimension $5$ is $$ {\bf E}_{\,{\bf 5}} = \left( {\matrix{ 0 & 1 & 0 & 0 & 0 \cr 0 & 0 & 1 & 0 & 0 \cr 0 & 0 & 0 & 1 & 0 \cr 0 & 0 & 0 & 0 & 1 \cr 0 & 0 & 0 & 0 & 0 \cr } } \right) $$ then its non-negative integral powers are just given by a shift of the non-null diagonal, till ${\bf E}_{\,{\bf 5}} ^{\,{\bf 5}} ={\bf 0}$.
I know that the Jordan decomposition of ${\bf E}_{\,h} ^{\,{\bf n}}$ is given by
$$
{\bf E}_{\,h} ^{\,{\bf n}} = {\bf P}_{\,h} (n)\;{\bf C}_{\,h} (n)\;{\bf P}_{\,h} (n)^{\, - \;{\bf 1}}
$$
where
- ${\bf P}_{\,h} (n)$ is a permutation matrix;
- ${\bf C}_{\,h} (n)$ is actually a ${\bf E}_{\,h}$ with $n-1$ ones missing in certain positions, i.e.
it is also expressible as a permutation.
For instance $$ \eqalign{ & {\bf E}_{\,{\bf 5}} ^{\,{\bf 2}} = \left( {\matrix{ 0 & 0 & 1 & 0 & 0 \cr 0 & 0 & 0 & 1 & 0 \cr 0 & 0 & 0 & 0 & 1 \cr 0 & 0 & 0 & 0 & 0 \cr 0 & 0 & 0 & 0 & 0 \cr } } \right) = {\bf P}_{\,{\bf 5}} (2)\;{\bf C}_{\,{\bf 5}} (2)\;{\bf P}_{\,{\bf 5}} (2)^{\, - \;{\bf 1}} = \cr & = \left( {\matrix{ 1 & 0 & 0 & 0 & 0 \cr 0 & 0 & 0 & 1 & 0 \cr 0 & 1 & 0 & 0 & 0 \cr 0 & 0 & 0 & 0 & 1 \cr 0 & 0 & 1 & 0 & 0 \cr } } \right)\left( {\matrix{ 0 & 1 & 0 & 0 & 0 \cr 0 & 0 & 1 & 0 & 0 \cr 0 & 0 & 0 & 0 & 0 \cr 0 & 0 & 0 & 0 & 1 \cr 0 & 0 & 0 & 0 & 0 \cr } } \right)\left( {\matrix{ 1 & 0 & 0 & 0 & 0 \cr 0 & 0 & 1 & 0 & 0 \cr 0 & 0 & 0 & 0 & 1 \cr 0 & 1 & 0 & 0 & 0 \cr 0 & 0 & 0 & 1 & 0 \cr } } \right) \cr} $$
But, after various trials, I could not yet find an effective way to express the type of P and C permutations wrt $h$ and $n$.
Most probably the subject has been already studied or it is easily assessable with the appropriate approach, and I am asking
for hints in this respect.