For two lines to be parallel, the vectors defining their slopes have to be parallel, which means their vectors have to be constant multiples of one another. It's easy to see by inspection that $\langle1,1,0 \rangle$ is not parallel to $\langle1,-2,1\rangle$. This shows that the lines aren't parallel. To prove they don't intersect, set the two lines' equations equal to one another. This gives
$$
\langle3,-1,2\rangle + t\langle1,1,0 \rangle = \langle0,5,2\rangle+s\langle1,-2,1\rangle
$$
Solving this component-wise leads to the system
$$
\begin{align}
3&=s-t\\
-6&=-2s-t\\
0&=s
\end{align},
$$
which is clearly not possible because substituting $s=0$ as given by the last equation makes the first two equations $t=-3$ and $t=6$, which is a contradiction, so the system has no solution, so the lines don't intersect.
For the third part, the distance, there are two simple ways to solve the problem. First, you can use the distance formula and then minimize it. Let $D(s,t)$, where $s$ is the parameter for $l_2$ and $t$ is the parameter for $l_1$, be the square of the distance between arbitrary points on $l_1$ and $l_2$. Then minimizing
$$ D(s,t) = ((3+t)-s)^2 + ((-1+t)-(5-2s))^2 + (2-(2+s))^2$$
gives the minimum distance between the two lines.
Alternatively, and much more easily, you can find a vector $\vec n$ that is perpendicular to both $l_1$ and $l_2$. Then, the distance between $l_1$ and $l_2$ is just the projection of $\overrightarrow{Q_1Q_2}$ onto $\vec n$, where $Q_1$ is a point on $l_1$ and $Q_2$ is a point on $l_2$. You can find $\vec n$ by taking the cross product of $\vec{v_1} = \langle1,1,0 \rangle$ and $\vec{v_2} = \langle1,-2,1\rangle$. That is,
$$
\vec n = \vec v_1 \times \vec v_2.
$$
Using that value, the minimum distance between the lines is
$$
\mbox {proj}_\overrightarrow{Q_1Q_2}\vec n = \frac{\|\overrightarrow{Q_1Q_2}\cdot\vec n\|}{\|\vec n\|}
$$
First the equations for all vectors $x$ on line $g$ and all vectors $y$ on line $h$:
$$
\begin{align}
g: x &= a + \lambda b \quad \\
h: y &= c + \mu d
\end{align}
\quad (*)
$$
The difference vector between two of those vectors is
$$
D = y - x = c - a + \mu d - \lambda b
$$
The length of $D$ squared is:
\begin{align}
q(\lambda, \mu)
&= D \cdot D \\
&= (c - a)^2 + (\mu d-\lambda b)^2 + 2 (c-a)\cdot(\mu d - \lambda b) \\
&= (c - a)^2 + \mu^2 d^2 + \lambda^2 b^2 - 2 \mu \lambda (d\cdot b)
+ 2 \mu ((c-a)\cdot d) - 2 \lambda ((c-a)\cdot b) \\
\end{align}
The gradient of $q$ is:
$$
q_\lambda = 2\lambda b^2 - 2\mu (d\cdot b)-2((c-a)\cdot b) \\
q_\mu = 2\mu d^2 - 2\lambda (d\cdot b)+2((c-a)\cdot d)
$$
It should vanish for local extrema of $q$ which leads to the system
$$
A u = v
$$
with
$$
A =
\left(
\begin{matrix}
b^2 & - d\cdot b \\
- d \cdot b & d^2
\end{matrix}
\right)
\quad
u =
\left(
\begin{matrix}
\lambda \\
\mu
\end{matrix}
\right)
\quad
v =
\left(
\begin{matrix}
(c-a) \cdot b \\
-(c-a) \cdot d
\end{matrix}
\right)
$$
A unique solution exists for
$$
0 \ne \mbox{det A} = b^2 d^2 - (d \cdot b)^2
$$
Note that the dot operator stands for the scalar product.
That solution is
$$
u = A^{-1} v
$$
with
$$
A^{-1} = \frac{1}{\mbox{det } A}
\left(
\begin{matrix}
d^2 & d\cdot b \\
d \cdot b & b^2
\end{matrix}
\right)
$$
Inserting the found values for $\lambda$ and $\mu$ into equations $(*)$ will provide you the two vectors, whose points are closest to each other.
Example:
$$
a = (2,0,0), b=(1,1,1), c=(0,1,-1), d=(-1,0,-1)
$$
leads to
$$
\lambda = 1, \mu = -2.5, x_\min = (3,1,1), y_\min=(2.5,1,1.5), D=(-0.5,0,0.5)
$$
The image renders the $x$ values in green, the $y$ values in purple, and $D$ in red.
Best Answer
Take the common normal direction.
$$\mathbf{n} = \pmatrix{1\\2\\1} \times \pmatrix{-1\\-1\\1} = \pmatrix{3 \\-2 \\1 } $$
Now project any point from the lines onto this direction. Their difference is the distance between the lines
$$ d = \frac{ \mathbf{n} \cdot ( \mathbf{r}_1 - \mathbf{r}_2 )}{\| \mathbf{n} \|} $$
$$ d = \frac{ \pmatrix{3\\-2\\1} \cdot \left( \pmatrix{1\\-1\\0} - \pmatrix{1\\2\\1} \right) }{ \| \pmatrix{3\\-2\\1} \|} = \frac{ \pmatrix{3\\-2\\1} \cdot \pmatrix{0\\-4\\0} } {\sqrt{14}} = \frac{8}{\sqrt{14}} = 2.1380899352993950$$
NOTE: The $\cdot$ is the vector inner product, and $\times$ is the cross product