Esta cuestión ya ha sido vinculado a una pregunta acerca de la generalización para el caso de los $M$ personas con el mismo cumpleaños (la paradoja de Cumpleaños con M compartido cumpleaños), me gustaría añadir, la respuesta general. Comenzando con la expresión dada por Michael Lugo (una derivación se puede encontrar aquí), es más fácil obtener la solución general. En primer lugar, la expresión se puede simplificar mediante el uso de la caída de factoriales:
\begin{ecuación}
\displaystyle P(W\ge1) = 1 - \sum_{i=0}^{\lfloor n/2 \rfloor} \frac{(n)_{2i}(m)_{n-i}}{i! 2^i m^n}
\end{ecuación}
donde $(n)_k := n\cdot(n-1)\cdot\ldots\cdot(n-k+1)$. El $i$ésimo término en la suma representa la probabilidad de tener $i$ distintos pares de personas.
Para generalizar este resultado observamos que el número 2 es sólo el tamaño de un par y puede ser sustituido por cualquier grupo mayor tamaño. Sin embargo, debemos ser cuidadosos, ya que los $2^i$ en el denominador es realmente $(2!)^i$ (véase la derivación para más detalles).
Por lo tanto, si queremos saber la probabilidad de que al menos $M$ personas de un grupo de $n$ comparten un cumpleaños, se suma sobre todos los tamaños de grupo menor que $M$, en sustitución de $2i$ en el numerador por $ki$ (donde $k$ es la suma del índice) y $2^i$ en el denominador por $(k!)^i$:
\begin{ecuación}
\displaystyle P(W\ge1) = 1 - \sum_{k=2}^{M-1}\sum_{i=0}^{\lfloor n/k \rfloor} \frac{(n)_{ki}(m)_{n-i}}{i! (k!)^i m^n}
\end{ecuación}