Ira Gessel has very kindly explained how this can be solved by counting nonintersecting lattice paths. Again, working through the example of $(3,2,1,1)$ will explain the general approach.
First, convert the problem to counting distinct part partitions contained in $(3+3,2+2,1+1,1+0) = (6,4,2,1)$. Now a collection of four nonintersecting paths from $(0,0), (1,0), (2,0), (3,0)$ to $(1,6), (2,4), (3,2), (4,1)$ gives a subpartition by looking at the height of each column of boxes. E.g.,

corresponds to $(5,3,2,0)$ in $(6,4,2,1)$ and the partition $(5-3,3-2,2-1,0-0) = (2,1,1)$ in $(3,2,1,1)$. Using the Gessel-Viennot Lemma, the number of such quartets of nonintersecting paths is given by
$$ \begin{vmatrix} \binom{7}{1} & \binom{6}{2} & \binom{5}{3} & \binom{5}{4} \\
\binom{6}{0} & \binom{5}{1} & \binom{4}{2} & \binom{4}{3} \\
\binom{7}{-1} & \binom{4}{0} & \binom{3}{1} & \binom{3}{2} \\
\binom{8}{-2} & \binom{5}{-1} & \binom{2}{0} & \binom{2}{1} \end{vmatrix}
= \begin{vmatrix} 7 & 15 & 10 & 5 \\
1 & 5 & 6 & 4 \\
0 & 1 & 3 & 3 \\
0 & 0 & 1 & 2 \end{vmatrix}
= 19.$$
Note that moving to distinct part partitions is necessary for this approach since, for instance, the empty partition $(0,0,0,0)$ in $(3,2,1,1)$ would correspond to paths that intersect other dots along the bottom. In the equivalent distinct part subpartition problem, the empty partition corresponds to $(3,2,1,0)$.
Let me close with a historical note. In the context of Young diagrams, this determinant result is usually attributed to Kreweras 1965. In the context of partitions, MacMahon's solution is in his 1915 collection and may date to even earlier. (The MacMahon and Kreweras solutions are very similar, and Gessel and Viennot connect to Kreweras as an application of their results.)