Prove this inequality with $xyz\le 1$ 
If $x,y,z\ge 0$ and $\color{red}{xyz\le 1}$, show that
$$\color{blue}{\dfrac{x^2-x+1}{x^2+y^2+1}+\dfrac{y^2-y+1}{y^2+z^2+1}
+\dfrac{z^2-z+1}{z^2+x^2+1}\ge 1}.$$


Multiplying and dividing each of the expressions numerator by $(1 + 
   x)$, $(1 + y)$ and $(1 + z)$ respectively:
$$
\frac{z^3+1}{(z+1) \left(x^2+z^2+1\right)}+\frac{x^3+1}{(x+1) \left(x^2+y^2+1\right)}+\frac{y^3+1}{(y+1) \left(y^2+z^2+1\right)}
$$
Further simplifying:
$$
1+ \frac{x^2-x+1}{x^2+y^2+1}-\frac{x^2+z}{x^2+z^2+1}+\frac{y^2-y+1}{y^2+z^2+1}
$$
Giving similar treatment to second and fourth terms:
$$
1+\frac{x^3+1}{(x+1) \left(x^2+y^2+1\right)}+\frac{y^3+1}{(y+1) \left(y^2+z^2+1\right)}-\frac{x^2+z}{x^2+z^2+1}
$$
Note that the $1$ in the above expression suggests that this expression could be $\geq1$ as long as the sum of second, third and fourth term is $\geq0$. Simplifying further:
$$
1+\frac{y^2-y+1}{y^2+z^2+1}-\frac{x^4+\left(y^2+z\right) x^2+x+y^2 z+z-1}{\left(x^2+y^2+1\right) \left(x^2+z^2+1\right)}
$$
$$
\Rightarrow 1+\frac{y^3+1}{(y+1) \left(y^2+z^2+1\right)}+\frac{-x^4-x^2 y^2-x^2 z-x-y^2 z-z+1}{\left(x^2+y^2+1\right) \left(x^2+z^2+1\right)}
$$
Note that since the second term is always positive, this means first and second terms will always be greater than 1. Now let's concentrate on third term to see if that is negative and if so, whether it is big enough to make the whole term negative.
We are given that $xyz\le1$, let us consider the upper bound of this inequality and assume $xyz=1$, this implies, $z=\frac{1}{xz}$, which upon substitute in above expression, we get:
$$
\frac{x^2+y^2+1+x y \left(-x^2+\left(2 x^3 y-x^2 y\right)+x y^3+x y-y^2-1\right)}{\left(x^2+y^2+1\right) \left(x^2 \left(x^2+1\right) y^2+1\right)}
$$
How do we find out if the above term is positive or not? Let's try proof by contradiction. The denominator is always positive, so let's concentrate on the numerator and assume that the term is negative.
The approach I took was to consider only the nominator and try to solve the following equations simulataneously:
$$
2 x^4 y^2-x^3 y^2-x^3 y+x^2 y^4+x^2 y^2+x^2-x y^3-x y+y^2+1 < 0
$$
$$x\gt0$$
$$y\gt0$$
And since there does not exist any solution so that means that numerator will be non-negative.
Hence we can conclude that the term will be $\ge1$
