Division algorithm for multivariate polynomials?
Solution 1:
No, the polynomial division algorithm does not immediately generalize to multivariate rings. Here is a simple proof. Just as for $\rm\:\Bbb Z,\:$ a domain having an algorithm for division with smaller remainder, also enjoys Euclid's algorithm for gcds, which, in extended form, yields Bezout's identity. Therefore gcds have linear representation $\rm\:gcd(a,b) = r a + s b\:$ (i.e. Euclidean domains are Bezout). But this fails in multivariate polynomial rings $\rm\:F[x_1,\ldots, x_n],\ n \ge 2,\:$ since $\rm\:gcd(x_1,x_2) = 1\:$ but there is no Bezout equation $\rm\: 1 = x_1 f + x_2 g\ $ (evaluating at $\rm\:x_1 = 0 = x_2\,$ $\Rightarrow$ $\rm\,1 = 0\:$ in $\rm\,F,\,$ contra field axioms).
But all Euclidean is not lost, since one can generalize the polynomial division algorithm in a way that recovers many of the important properties. For this, look up the Grobner basis algorithm, which may be viewed as a nonlinear generalization of the Euclidean/division algorithm and Gaussian elimination.
Solution 2:
If you can do Euclidean division in a reasonable sense in a domain $R$, then you can prove that $R$ is a principal ideal domain. This is not the case for the $F[x_{1}, \dots, x_{n}]$ for $n > 1$.
Please see the Wikipedia page for a clear and precise definition.