[step:Prove linear independence of $\mathcal{B}_S$ via the polynomial ring isomorphism]We establish linear independence of $\mathcal{B}_S$ in arbitrary characteristic by constructing an explicit isomorphism between $\operatorname{Sym}^r(V)$ and the space of homogeneous polynomials of degree $r$.
Let $\{e_1, \ldots, e_n\}$ be the chosen basis for $V$, and let $k[x_1, \ldots, x_n]_r$ denote the $k$-vector space of homogeneous polynomials of degree $r$ in $n$ variables over $k$. Define the $k$-linear map
\begin{align*}
\Phi: V^{\otimes r} &\to k[x_1, \ldots, x_n]_r \\
e_{i_1} \otimes \cdots \otimes e_{i_r} &\mapsto x_{i_1} \cdots x_{i_r}
\end{align*}
extended by linearity. Since polynomial multiplication is commutative, $\Phi(v_1 \otimes \cdots \otimes v_r)$ is unchanged under any permutation of the factors $v_1, \ldots, v_r$. Therefore $\Phi$ vanishes on all elements of the form $v_1 \otimes \cdots \otimes v_r - v_{\sigma(1)} \otimes \cdots \otimes v_{\sigma(r)}$ for $\sigma \in S_r$, which generate the kernel of the quotient map $\pi: V^{\otimes r} \to \operatorname{Sym}^r(V)$. By the universal property of quotients, $\Phi$ descends to a well-defined $k$-linear map
\begin{align*}
\bar{\Phi}: \operatorname{Sym}^r(V) &\to k[x_1, \ldots, x_n]_r \\
e_{i_1} \cdots e_{i_r} &\mapsto x_{i_1} \cdots x_{i_r}.
\end{align*}
**Surjectivity.** Every monomial $x_{i_1} \cdots x_{i_r} \in k[x_1, \ldots, x_n]_r$ is the image of $e_{i_1} \cdots e_{i_r} \in \operatorname{Sym}^r(V)$. Since the monomials span $k[x_1, \ldots, x_n]_r$, the map $\bar{\Phi}$ is surjective.
**Dimension count on the target.** The monomials $x_1^{\mu_1} \cdots x_n^{\mu_n}$ with $\mu_k \ge 0$ and $\mu_1 + \cdots + \mu_n = r$ form a basis for $k[x_1, \ldots, x_n]_r$ over $k$. This is a standard fact: distinct monomials are linearly independent because they are distinguished by their exponent vectors (the coefficient of any given monomial in a polynomial is well-defined). The number of such monomials is the number of solutions to $\mu_1 + \cdots + \mu_n = r$ in non-negative integers, which by stars-and-bars equals $\binom{n+r-1}{r}$. Therefore $\dim_k k[x_1, \ldots, x_n]_r = \binom{n+r-1}{r}$.
**Injectivity from dimension comparison.** Since $\mathcal{B}_S$ spans $\operatorname{Sym}^r(V)$ (established in Step 1) and $|\mathcal{B}_S| = \binom{n+r-1}{r}$ (the stars-and-bars count), we have $\dim \operatorname{Sym}^r(V) \le \binom{n+r-1}{r} = \dim k[x_1, \ldots, x_n]_r$. Since $\bar{\Phi}$ is a surjection between finite-dimensional spaces with $\dim \operatorname{Sym}^r(V) \le \dim k[x_1, \ldots, x_n]_r$, and surjectivity forces $\dim \operatorname{Sym}^r(V) \ge \dim k[x_1, \ldots, x_n]_r$, we conclude $\dim \operatorname{Sym}^r(V) = \binom{n+r-1}{r}$ and $\bar{\Phi}$ is an isomorphism.
In particular, since $\bar{\Phi}$ maps the spanning set $\mathcal{B}_S$ onto the basis of monomials bijectively, $\mathcal{B}_S$ is linearly independent.[/step]