Эта задача была на старом нашем форуме несколько лет назад. Привожу свое решение вместе с кусочком дискуссии.

Задача такая:
Из колоды в 52 карты наудачу берется 6 карт. Какова вероятность того, что среди них будут представители всех четырех мастей.



Я уже писал, что мое решение неверно. В нем действительно (как заметил Кэро)некоторые комбинации при вычислении m учитываются несколько раз. Но задача на удивление интересна, т.к. по постановке она классическая (а потому кажется, что решение простое), но решение ее требует достаточно большого умственного напряжения (по крайней мере для меня). Я хочу предложить свое решение для этой задачи, хотя мой ответ совпадает с ответом Кэро (я проверил окончательные цифры - он молодец!), да и решение у Кэро короче. Но его решение написано так кратко, что я его не понимаю. Поэтому и привожу свое решение , хотя оно, по-видимому, длиннее.
Сначала пришлось вывести формулу для вероятности суммы четырех событий (хотя это и несложно):
P(A+B+C+D)=P(A)+P(B )+P(С )+P(D) – P(AB) – P(AC) – P(AD) – P(CD) – P(BC) – P(BD) + P(ABC) + P(ABD) + P(ACD) + P(BCD) – P(ABCD) .
Обозначим масти номерами (так удобнее) : 1 – пики, 2 – крести, 3 – бубни, 4 – черви.
Обозначим события:
А – среди выданных 6 карт есть представители всех мастей.
А1 – среди выданных карт есть пиковые, … , А4 – среди выданных есть червовые.
Обозначим через В, В1, …, В4 события, противоположные событиям А, А1, …,А4.
Тогда Р(А) = 1 – Р(В). Найдем Р(В). Ясно, что В=В1+В2+В3+В4. Поэтому сейчас для вычисления Р(В) надо записать формулу выше
P(В1+B2+В3+В4)=P(В1)+P(В2)+P(В3)+P(В4) – P(В1B2) – P(В1В3) – P(В1В4) – P(В3В4) – P(B2В3) – P(B2В4) + P(В1В2В3) + P(В1В2В4) + P(В1В3В4) + P(B2В3В4) – P(В1В2В3В4) .
Ясно, что последнее слагаемое = 0. Из симметрии мастей все слагаемые в каждой из трех групп слагаемых в этой формуле одинаковы, а потому
P(В)=4*P(В1) – 6*P(В1B2) + 4*P(В1В2В3) . Теперь уже проще. Пусть n = С(52,6) – как и было раньше – общее число различных шестерок карт. Тогда ясно, что
Р(В1)=С(39,6)/n, Р(В1В2)=С(26,6)/n, а Р(В1В2В3)=С(13,6)/n . Подставляя это в формулу выше, получим Р(В)=0.573517917. А потому
Р(А) = 0.426482082 .