Suppose we have a well shuffled deck of 52 cards and we want to divide those cards evenly among four non distinct players. We're asked to calculate the probability that all of them receive exactly one ace.
Without the added constraint of four of the players being non distinct this chance would have been easily calculable by using the multinomial theorem like so:
$$n_a = \frac{4!\cdot48!}{(12!)^4},\, N=\frac{52!}{(13!)^4},\, P=\frac{n_a}{N}=\frac{4!\cdot48!\cdot(13!)^4}{52!\cdot(12!)^4}$$
where $n_a$ are the outcomes we're interested in and $N$ is the complete outcome space but with the added constraint I'm stuck. Initially, I thought that this could be solved by multiplying $n_a$ with $4!$ but this won't fly since this would result in a higher probability than 1. Could someone point me in the right direction?