Suppose we have $k$ buckets, with an infinite number of balls in each bucket. The balls within one bucket are indistinguishable, those between buckets are. We assign to bucket $i$ a probability $p_i$ (so that $\sum_{i=1}^k p_i=1$). If we draw a ball from one of these buckets with probability $p_i$ without replacing it, and repeat this $N$ times, we are sampling from a multinomial distribution with $N$ trials and $k$ categories. The fact that there is an infinite number of balls in each bucket makes replacement superfluous.
Now, what happens if bucket $i$ only contains $c_i$ balls? The $c_i$ are defined arbitrarily, but such that $\sum_{i=1}^k c_i\ge N$. When drawing from the buckets with probability $p_i$, the number of balls in the bucket has no influence on the outcome as long as there is a ball in the bucket. That is, if a bucket $i$ is depleted, you cannot draw from that bucket anymore. You therefore have to draw from another bucket, with probability $p_j/(1-\sum_{i\in depleted} p_i)$. What is this distribution? Can its pmf be given? Can it be related to the multivariate hypergeometric distribution? Are there approximations for it?