In general, synthetic division in multiple variables is not well defined since the notion of remainder is not unique. In Gröbner basis theory, one defines for that purpose a "monomial order" that then allows to reduce multivariate polynomials like univariate polynomials and to determine minimal generating sets for ideals.
If one knows the remainder, which is zero here, and the final degree structure, as in this case, one can use one of Schönhages tricks, division by multiplication. To divide $f$ by $g=1+x^ay^bz^c$, multiply successively f by $1-x^ay^bz^c$, $1+(x^ay^bz^c)^2$, $1+(x^ay^bz^c)^4$, $1+(x^ay^bz^c)^8$ etc. until the degree of $(x^ay^bz^c)^{2^k}$ is outside of the degree range of the quotient. Then the lower degree terms of this product form the quotient.
Best Answer
Now, suppose $G(x)=x^3-19x+30$ and note that $G(2)=0$.
So $(x-2)$ is a factor of $G(x)$ by the factor theorem.
By synthetic division, we get $$\color{red}{G(x)=(x-2)(x^2+2x-15)}$$
Again, suppose $H(x)=x^2+2x-15$ and note that $G(3)=0$.
So $(x-3)$ is a factor of $G(x)$ by the factor theorem.
By synthetic division, we get $$\color{orange}{H(x)=(x-3)(x+5)}$$
So we can write that $$\color{blue}{x^4-x^3-19x^2+49x-30=(x-1)(x-2)(x-3)(x+5)}$$