If $P(x)=c$ for a constant $c,$ then $\dfrac{c(x+y+z)}{xyz}=3c$ . We have $2c=3c.$ Therefore $c=0.$
Now consider the case of non-constant polynomials.
First we have $$xP(x)+yP(y)+zP(z)=xyz(P(x-y)+P(y-z)+P(z-x))$$ for all nonzero real numbers $x,y,z$ satisfying $2xyz=x+y+z$ . Both sides of the equality are polynomials (of $x,y,z$ ). They have the same values on the 2-dimensional surface $2xyz=x+y+z$ , except for some 1-dimensional curves in it. By continuity, the equality holds for all points on the surface, including those with $z=0.$ Let $z=0,$ we have $y=-x$ and $x(P(x)-P(-x))=0.$ Therefore $P$ is an even function.
(Here is a sketch of an elementary proof. Let $z=\dfrac{x+y}{2xy-1}.$ We have $$xP(x)+yP(y)+\dfrac{x+y}{2xy-1}P(\dfrac{x+y}{2xy-1})=xy\dfrac{x+y}{2xy-1}(P(x-y)+P(y-\dfrac{x+y}{2xy-1})+P(\dfrac{x+y}{2xy-1}-x)).$$ This is an equality of rational expressions. By multiplying $(2xy-1)^N$ on both sides for a sufficiently large $N$ , they become polynomials, say $A(x,y)=B(x,y)$ for all real $x, y$ with $x\ne 0, y\ne 0, x+y\ne 0$ and $2xy-1\ne 0.$ For a fixed $x,$ we have two polynomials (of $y$ ) having same values for infinitely many $y$ . They must be identical. Let $y=0,$ we have $x^{N+1}(P(x)-P(-x))=0.$ )
Notice that if $P(x)$ is a solution, then is $cP(x)$ for any constant $c.$ For simplicity, we assume the leading coefficient of $P$ is $1$ : $$P(x)=x^n+a_{n-2}x^{n-2}+\cdots +a_2x^2+a_0,$$ where $n$ is a positive even number.
Let $y=\dfrac{1}{x}$ , $z=x+\dfrac{1}{x}.$ we have $$xP(x)+\dfrac{1}{x}P\left (\dfrac{1}{x}\right )+\left ( x+\dfrac{1}{x}\right ) P\left ( x+\dfrac{1}{x}\right ) =\left (x+\dfrac{1}{x}\right )\left ( P\left (x-\dfrac{1}{x}\right )+P(-x)+P\left (\dfrac{1}{x}\right )\right ).$$
Simplify using $P(x)=P(-x),$ $$\left (x+\dfrac{1}{x}\right ) \left (P\left (x+\dfrac{1}{x}\right )-P\left (x-\dfrac{1}{x}\right )\right )=\dfrac{1}{x}P(x)+xP\left (\dfrac{1}{x}\right ).$$
Expand and combine like terms, both sides are of the form $$c_{n-1}x^{n-1}+c_{n-3}x^{n-3}+\cdots+c_1x+c_{-1}x^{-1}+\cdots+c_{-n+1}x^{-n+1}.$$
They have the same values for infinitely many $x.$ They must be identical. We just compare their leading terms. On the left hand side it is $2nx^{n-1}$ . There are two cases for the right hand sides: If $n>2$ , it is $x^{n-1}$ ; If $n=2$ , it is $(1+a_0)x.$ It does not work for $n>2.$ When $n=2,$ we have $4=1+a_0.$ therefore $a_0=3.$
The solution: $P(x)=c(x^2+3)$ for any constant $c.$
-JZ