In Stanley's Enumerative Combinatorics, there is the following identity
$$\sum_{n \geq k}S(n,k) x^n = \frac{x^k}{(1-x)(1-2x) \dots (1-kx)}$$
where $S(n,k)$ denotes the number of partitions of an $n$-set into $k$ blocks. (A partition of an $n$-set $N$ into $k$ blocks is a collection of $k$ non-empty and disjoint of subsets of $N$ whose union is $N$).
What I have tried, following Stanley's perceived advice is expanding the denominator terms
$$\frac{1}{1-ax} = \sum_{n \geq 0} a^nx^n$$
so that
$$\sum_{n \geq k}S(n,k) x^n = x^k \left(\sum_{n \geq 0} x^n\right)\left(\sum_{n \geq 0} 2^nx^n\right)\dots\left(\sum_{n \geq 0} k^nx^n\right)$$
and then comparing coefficients. What I obtained is that for $n = k + q$, $q \geq 0$, then
$$\Large \sum_{\substack{(q_1,q_2, \dots, q_k) \\q_1,q_2,\dots,q_k \geq 0 \\ \sum_{i=1}^n = q}} 1^k1^{q_1}2^{q_2}\dots k^{q_k} = S(n,k)$$
but I fail to grasp how this sum represents $S(n,k)$. I understand the summation goes over $k$-compositions of $n$ but I can't seem to relate it to $S(n,k)$, neither combinatorially nor by the closed form expression
$$S(n,k) = \frac{1}{k!}\sum_{i=0}^k(-1)^{k-i}\binom{k}{i}i^n$$
Maybe I'm entirely on the wrong track. Any ideas?