$\phi_n(x,y,z)\in FO$, we can use only three variables, but we can requantify them. Only available relation is edge relation $E$. For each finite directed cycle: $$(C, x : a, y : b, z : c) \models \phi_n(x,y,z)\leftrightarrow \text{$C$ has $3n$ edges, directed distances: dist(a,b) = dist(b,c) = dist(c,a) = n}$$
So, lets look at my trial:
$\phi_n(x,y,z) \equiv \psi_{n}(x,y)\wedge \psi_{n}(y,z) \wedge \psi_{n}(z,x)$
$\psi_0(x,y) \equiv x = y$
$\psi_n(x,y) \equiv \exists!_z(E(x,z) \wedge \psi_{n-1}(z,y))$
$\exists!$ means that there exists exactly one. Can someone check my solution ?