Чарльз Мостеллер: Пятьдесят занимательных вероятностных задач с решениями

Итак, оценкой числа π в этом случае является 4∙5.2/(62/10) ≈ 3.35 вместо 3.14. При другом опыте, состоящем также из 10 подбрасываний, было получено 67 пересечений, что дает оценку 3.10.

Е. Молина предложил эту задачу, чтобы дать формулировку знаменитой проблемы Ферма на вероятностном языке.

Пусть z обозначает число белых шаров в первой урне, x — число белых шаров и y — число черных шаров во второй урне. Тогда задача состоит в том, чтобы найти целые числа n , x , y и z такие, что

или

z n = x n + y n .

Хотя для многих значений n известно, что это уравнение не имеет корней, но не установлено, так ли это при всех n ≥ 3. Доказано, однако, что целочисленных решений нет при n < 2000.

Из таблиц или из непосредственного расчета нетрудно выписать распределения числа простых делителей для небольших значений N .

В таблице 1 приведены результаты для N = 100 и N = 1000 вместе со средними x̅ и дисперсиями s ².

Из этой таблицы, например, видно, что среди первых 100 натуральных чисел ровно 26 простых, у 34 чисел два простых делителя и только у двух шесть простых делителей.

Распределение числа делителей при N = 100 напоминает выборку из закона Пуассона. Для пуассоновских распределений среднее равно дисперсии. Из таблицы видно, что для N = 100 среднее несколько больше дисперсии. Если рассмотреть величину x − 1 вместо x , то новое среднее будет равно 1.40, а дисперсия, равная 1.48, не изменится. Полезно сравнить полученные результаты с табличными вероятностями для закона Пуассона. (Сумма элементов последней строки первой половины табл. 2 не равна 100 из-за округления значений.)

Видно, что при N = 100 совпадение лучше, нежели при N = 1000. Для N = 1000 более точная аппроксимация при небольших значениях x − 1 может быть получена за счет выбора меньшего математического ожидания пуассоновского распределения.

Таблица 2 подтверждает предположение о пуассоновости распределения числа простых делителей, однако картина слишком сложна, чтобы можно было угадать вид параметра этого закона для больших N .

Мы знаем, что вероятность отсутствия простых делителей, т. е. того, что само число просто, равна приближенно 1/ln( N ). Для закона Пуассона вероятность появления 0 равна e −m , где m — математическое ожидание этого распределения (см. задачу 29). Отсюда выводим:

и

− m = −ln(ln( N )),

или

m = ln(ln( N )).

Любопытно сравнить эту формулу с полученными ранее результатами.

Имеем

ln(ln(100)) = 1.53,

что надо сравнить со средним 1.4 при N = 100. Для N = 1000 среднее равнялось 1.88, а