proof about commutative operators and T-cyclic vectors

With some general theory this can be done quite painlessly. By the structure theorem for finitely generated modules over a principal ideal domain (with $F[X]$ as PID, and $V$ considered as $F[X]$-module with $X$ acting as $T$), $V$ decomposes as a finite direct sum of cyclic modules; the decomposition is not unique, but the number and types of the modules are unique (and so $V$ has a cyclic vector if and only if there is (at most) one cyclic module in such a decomposition). The cyclic factors are either isomorphic to $F[X]$ as module over itself (but this cannot happen for $V$ because it is finite dimensional over $F$), or of the form $F[X]/P$ for a non-constant monic polynomial $P\in F[X]$; moreover these polynomials $P_i$ in the successive cyclic modules (called the invariant factors of the $F[X]$ module $V$) divide each other: $P_i\mid P_{i+1}$ whenever both polynomials exist.

Now the question is to show that if all linear operators commuting with $T$ are polynomials in $T$ (in other words they are just scalar multiplications in the $F[X]$-module, where "scalar" means element of $F[X]$) then there cannot be two or more cyclic modules in the decomposition (the converse is also true, see this question). So assume a decomposition given by the theorem has at least two factors $F[X]/P_1$, $F[X]/P_2$, with $1\neq P_1\mid P_2$. Since the factors in the decomposition are submodules, every scalar multiplication stabilises them; it will suffice to find a $F[X]$-module morphism $V\to V$ that does not stabilise all factors (the module morphism property says that it is $F$-linear and commutes with $T$). But we have a non-zero morphism $\pi:F[X]/P_2\to F[X]/P_1$ of reducing modulo $P_1$, which is defined because $P_2$ is a multiple of $P_1$. To make this into a $F[X]$-module endomorphism $V$, it suffices to first project $V\to F[X]/P_2$ parallel to all other cyclic factors, then apply$~\pi$, and finally inject the factor $F[X]/P_1$ back into $V$. The resulting endomorphism is easily seen to not stabilise the cyclic factor $F[X]/P_2$, as desired. QED

Here is another answer, less general that the one I gave five years ago as I am going to assume $F$ is algebraically closed (for instance $F=\Bbb C$), but also using less advanced theory. I will show that if every linear operator that commutes with $T$ is a polynomial in $T$, then $T$ has a cyclic vector.

Since $F$ is algebraically closed, $V$ decomposes into a direct sum of generalised eigenspaces, which are of the form $V_i=\ker((T-\lambda_i I)^{m_i})$, taken for the distinct eigenvalues $\lambda_i$ of$~T$ (the exponents $m_i$ just need to be big enough; taking the algebraic multiplicity of $\lambda_i$ or even just the fixed value $n=\dim(V)$ for them will do).

A first reduction is to show that if the restrictions to $T$ to each $V_i$ all admit a cyclic vector, then so does $T$; this will allow us to then reason for each of these restrictions separately. To prove this point, view $V=V_1\oplus V_2\oplus\cdots\oplus V_k$ as a Cartesian product, and taking a cyclic vector $v_i$ in each summand, form $v=(v_1,v_2,\ldots,v_k)$ which is going to be our cyclic vector. By hypothesis, the vectors $T^j(v_i)$ for $0\leq j<\dim(V_i)$ together span $V_i$. It suffices to show that for every $i,j$ there exists a polynomial $P_{i,j}$ such that $P_{i,j}[T](v)=(0,\ldots,0,T^j(v_i),0,\ldots,0)$ (with the nonzero component in position $i$), as then polynomials in $T$ applied to $v$ span everything. Now since the annihilator polynomial $(X-\lambda_i)^{m_i}$ of $V_i$ is relatively prime to the product $\prod_{l\neq i}(X-\lambda_l)^{m_l}$ of the annihilator polynomials of the other $V_l$, there exists by the Chinese remainder theorem a polynomial multiple $Q$ of that product that satisfies $Q\equiv 1\pmod{(X-\lambda_i)^{m_i}}$. Then $Q[T](v_l)=0$ for $l\neq i$ while $Q[T](v_i)=v_i$. But then taking $P_{i,j}=X^jQ$ does the job.

With that proven, focus on the restriction to a single $V_i$, or equivalently assume there is exactly only one eigenvalue$~\lambda$. We claim that our hypothesis of "commutation implies polynomial in $T$" forces the the Jordan decomposition to have only one block (equivalently, the eigenspace for$~\lambda$ must have dimension$~1$). For if there were more than one block (for$~\lambda$), take the first block to be at least as large as the second, and consider the linear operator that on the Jordan basis for this Jordan normal form vanishes on all basis vectors except those for the second block, which it sends to the corresponding basis vectors for the first block: this operator commutes with $T$, but is no polynomial in$~T$ (since those stabilise the individual blocks).

We are done: if $T$ has a single Jordan block, than it has a cyclic vector (the final basis vector of the Jordan basis).