Let's rename the linear operator as $T$ and the non-zero vectors corresponding to the cyclic decomposition as $x^1,x^2,\ldots,x^r$(force of habit :)).
We'll take $W_0 = \{0\}$(with this, we will always have $p_1$(the $T$-annihilator of $x^1$)$=p$(the minimal polynomial), which greatly helps our calculation).
You've determined that the characteristic polynomial, $f=(t-1)^3$(I prefer to keep the leading coefficient as 1). Verify that the minimal polynomial is $p=(t-1)^2=p_1$.
We have $\deg(p_1)=\dim(Z(x^1;T))=2$, and hence $\dim(Z(x^2;T))=1$. We have no more cyclic factors of $x^i$ in our cyclic decomposition.
Since $\dim(Z(x^2;T))=1$, $x^2$ must be an eigenvector of $T$.
Verify that the only two linearly independent eigenvectors of $T$(which correspond to the eigenvalue $1$) are $v^1=\begin{pmatrix} 2 \\ 0 \\ 1 \end{pmatrix}$ and $v^2=\begin{pmatrix} 2 \\ 1 \\ 0 \end{pmatrix}$. Both of these are permitted vectors for $x^2$. We'll take $\boxed{x^2=\begin{pmatrix} 2 \\ 0 \\ 1 \end{pmatrix}}$
Let's look at how we can determine $x^1$. The following must be true regarding $x^1$:
(i). It must be linearly independent of $x^2$
(ii). It must generate a cyclic subspace of dimension precisely $2$, i.e. $Z(x^1;T)=span(\{x^1,T(x^1)\})$, in particular, $T^2(x^1)$ must be a linear combination of $x^1$ and $T(x^1)$(and consequently higher composition orders of $T$ as well).
(iii). We have $p_1=p=(t-1)^2$, i.e. $(T-2I)^2$ must be the $T$-annihilator of $x^1$.
Verify that $\boxed{x^1=\begin{pmatrix} 1 \\ 0 \\ 0 \end{pmatrix}}$ is allowed.
Note: $x^1,x^2$ are certainly not unique, there can be many(even infinitely many) possible candidates for them. Only the number of cyclic factors and $T$-annihilators are unique(and hence the later is also referred to as invariant factors).