Пусть у нас есть
брачных пар и
лодок. Будем считать, что лодки различимы (пронумерованы), а места в лодках нет. Тогда количество всех рассаживаний людей по лодкам (исходов) равно
.
Подсчитаем количество благоприятных исходов, используя принцип включений-исключений. Пусть свойство
, где
и
, состоит в том, что в лодку номер
попадает брачная пара номер
. При этом рассаживание будет благоприятным тогда и только тогда, когда оно не удовлетворяет ни одному из свойств
. Пусть
- это множество всех свойств, и для любого подмножества
число рассаживаний удовлетворяющих свойствам из
равно
. Тогда по принципу включений-исключений число рассаживаний, не удовлетворяющий никакому свойству из
равно:
Попробуем представить эту сумму в удобоваримом виде.
Прежде всего заметим, что если в
больше двух свойств вида
для фиксированного
, то
(в лодку не могут сесть больше 2-х пар). Аналогично, если в
больше одного
для фиксированного
, то
(каждая пара может сесть только в одну лодку). Для остальных
определим сигнатуру как двухкомпонентный вектор
, где
и
, то есть
- это число лодок куда согласно
должно сесть как минимум
брачных пар. Понятно, что если
имеет сигнатуру
, то
.
Вычислим
для подмножества
заданной сигнатуры
. Понятно, что элементы
однозначно определяют
лодок и
брачных пар, которые в них сидят. Для оставшихся
пар остается
пустых лодок и
лодок, заполненных наполовину (т.е. с одной уже сидящей парой). Поэтому
Чтобы дать окончательный ответ, осталось подсчитать, а сколько у нас есть подмножеств
заданной сигнатуры
. Так как
, то чтобы определить
, можно сначала выбрать брачные пары, которым
предписывает сидеть в лодках, потом выбрать
лодок, и наконец распределить выбранные пары по выбранным лодкам. Отсюда мы немедленно получаем, что количество различных
заданной сигнатуры
равно
Собирая все вместе, получаем количество благоприятных исходов:
Численные значения для
:
Код:
? f(n)=sum(m1=0,n,sum(m2=0,n-m1,(-1)^m1*(2*n)!*n!*(2*n-m1-2*m2)!/(n-m1-m2)!/m1!/m2!/2^(2*n-2*m1-m2)))
? vector(10,n,f(n))
%1 = [0, 6, 900, 748440, 1559930400, 6928346502000, 58160619655538400, 845986566719614320000, 19957466912796971445888000, 724891264860942581350908960000]
Соответственно, для
получаем искомую вероятность: