Задача об устойчивом паросочетании

Определение:

Устойчивое паросочетание (англ. stable matching) — паросочетание без неустойчивых пар.

Задача:

Найти полное устойчивое паросочетание между элементами двух множеств размера [math]N[/math], имеющими свои предпочтения.

Основная задача

Есть [math]n[/math] мужчин и [math]n[/math] женщин. Они обладают следующими особенностями:

1. Каждый человек оценивает лишь людей противоположного пола (все гетеросексуальны);
2. Каждый мужчина может отсортировать женщин от наименее привлекательной к наиболее привлекательной, причем его предпочтения не меняются (у каждого мужчины своя функция оценки);
3. Каждая женщина может отсортировать мужчин от наименее привлекательного к наиболее привлекательному, причем её предпочтения не меняются (у каждой женщины своя функция оценки).

Очевидным образом по такому определению строится полный двудольный граф (левая доля — мужчины, правая — женщины), назовем его МЖ.

Рассмотрим некоторое паросочетание в МЖ.

Определение:

Пара [math]\langle A, b\rangle[/math] называется неустойчивой (англ. unstable pair), если: В паросочетании есть пары [math]\langle A, a\rangle[/math] и [math]\langle B, b\rangle[/math] ([math]A[/math] женат на [math]a[/math], [math]B[/math] женат на [math]b[/math]) [math]A[/math] считает [math]b[/math] привлекательней, чем [math]a[/math] [math]b[/math] считает [math]A[/math] привлекательней, чем [math]B[/math]

Алгоритм Гейла-Шепли

Решение задачи было описано в [math]1962[/math] году математиками Девидом Гейлом (Университета Брауна) и Ллойдом Шепли (Принстонский университет) в статье «Поступление в колледж и стабильность браков» (College admissions and the stability of marriage) в журнале American Mathematical Monthly . Набор правил, следование которым всегда приводит к образованию стабильных пар, получил название алгоритма Гейла-Шепли или «алгоритма отложенного согласия» (алгоритм предложи-и-откажи).

Интуитивное описание

1. Мужчины делают предложение наиболее предпочитаемой женщине.
2. Каждая женщина из всех поступивших предложений выбирает наилучшее и отвечает на него «может быть» (помолвка), на все остальные отвечает «нет» (отказ).
3. Мужчины, получившие отказ, обращаются к следующей женщине из своего списка предпочтений, мужчины, получившие ответ «может быть», ничего не делают.
4. Если женщине пришло предложение лучше предыдущего, то она прежнему претенденту (которому ранее сказала «может быть») говорит «нет», а новому претенденту говорит «может быть».
5. Шаги [math]1[/math]-[math]4[/math] повторяются, пока у всех мужчин не исчерпается список предложений, в этот момент женщины отвечают «да» на те предложения «может быть», которые у них есть в настоящий момент.

Описание в псевдокоде

  // Изначально все мужчины не женаты и все женщины незамужние.
  k = 0 // количество сформированных пар
  while Существует M <- некоторый свободный мужчина, не делавший предложения всем женщинам
     w <- первая женщина из текущего списка M
     if w свободна
        помечаем M и w помолвленными
     else if w предпочитает M своему "текущему" жениху M'
        помечаем M и w помолвленными, M' — свободным
     else
        w отказывает M

Время работы алгоритма — [math]O(n^2)[/math], так как количество итераций цикла [math]\mathrm {while}[/math] не превосходит [math]O(n^2)[/math].

Доказательство корректности

Утверждение (Наблюдение [math]1[/math]):

Мужчины делают предложения женщинам в порядке убывания симпатии.

Утверждение (Наблюдение [math]2[/math]):

Как только женщина была помолвлена, она не может стать непомолвленной, она может только улучшить свой выбор (сказать «может быть» более предпочтительному кандидату).

Для начала покажем, что алгоритм завершит свою работу.

Лемма (Лемма [math]1[/math]):

Алгоритм завершается после максимум [math]n^2[/math] итераций цикла [math]\mathrm{\mathbf{while}}[/math].

Доказательство:

[math]\triangleright[/math]

На каждой итерации мужчина делает предложение очередной женщине. Но всего может быть не более [math]n^2[/math] предложений.

[math]\triangleleft[/math]

Теперь покажем, что по завершении алгоритма задача будет решена.

Лемма (Лемма [math]2[/math]):

Все мужчины и женщины будут заняты.

Доказательство:

[math]\triangleright[/math]

Предположим, что некоторый мужчина ([math]A[/math]) не женат по завершении алгоритма. Тогда и некоторая женщина ([math]a[/math]) незамужняя. По наблюдению [math]2[/math], [math]a[/math] не получала предложений. Но [math]A[/math] сделал предложения всем женщинам, так как он остался не женат. Получаем противоречие. Таким образом, все мужчины заняты. Аналогичные рассуждения применяем для женщин. Пусть какая-то женщина незамужняя. Значит, есть мужчина, который остался не женат. Но мы доказали, что по завершении алгоритма все мужчины заняты. Снова пришли к противоречию.

[math]\triangleleft[/math]

Лемма (Лемма [math]3[/math]):

После завершения алгоритма не будет неустойчивых пар.

Доказательство:

[math]\triangleright[/math]

Предположим [math]\langle A, b\rangle[/math] (где [math]A[/math], [math]B[/math] — мужчины; [math]a[/math], [math]b[/math] — женщины; [math]A[/math] женат на [math]a[/math], [math]B[/math] женат на [math]b[/math]) — нестабильная пара в паросочетании, найденном алгоритмом Гейла-Шепли. Возможны два случая: [math]A[/math] не делал предложения [math]b[/math]. Значит, [math]A[/math] находит [math]a[/math] более привлекательной, чем [math]b[/math]. Но чтобы рассматриваемая пара была неустойчивой, необходимо, чтобы [math]b[/math] для [math]A[/math] была более привлекательна, чем [math]a[/math]. Следовательно, [math]\langle A, b\rangle[/math] — устойчивая пара. [math]A[/math] делал предложение [math]b[/math]. Тогда был такой момент, когда [math]b[/math] отказала [math]A[/math], значит, [math]b[/math] находит [math]B[/math] более привлекательным, чем [math]A[/math]. Снова получается, что [math]\langle A, b\rangle[/math] — устойчивая пара.

[math]\triangleleft[/math]

Анализ полученного алгоритмом паросочетания

Алгоритм Гейла-Шепли гарантирует, что будет найдено некоторое решение задачи. Но решений может быть более одного. Зададимся вопросом, какими свойствами обладает решение, найденное алгоритмом.

Лемма (man-optimality):

Из всех возможных решений алгоритмом Гейла-Шепли будет найдено решение, наилучшее для мужчин (каждый мужчина получает в жены женщину, наилучшую из всех возможных при условии корректности решения).

Доказательство:

[math]\triangleright[/math]

Докажем от противного, что для каждого мужчины не существует устойчивого паросочетания, в котором его супругой была бы более желанная для него женщина. Предположим, для мужчины [math]A[/math] это свойство не выполняется. Так как он оказался женат не на лучшей из кандидатур, то существует женщина [math]a[/math], которая предпочла ему другого, более привлекательного мужчину [math]B[/math], при этом женщина [math]a[/math] для мужчины [math]B[/math] стоит на первом месте в его текущем списке. Предположим, существует устойчивое паросочетание, содержащее [math]\langle A, a\rangle[/math]. По определению, в устойчивом паросочетании нет неустойчивых пар. Пара [math]\langle B, a\rangle[/math] станет неустойчивой, если [math]B[/math] будет предпочитать [math]a[/math] своей супруге. Значит, [math]B[/math] женат на ком-то, кто лучше, чем [math]a[/math]. Но такое невозможно, так как [math]a[/math] стоит для него на первом месте. Таким образом, если женщина [math]a[/math] вычёркивается из списка предпочтений мужчины [math]A[/math], то любое паросочетание, содержащее [math]\langle A, a\rangle[/math], неустойчиво.

[math]\triangleleft[/math]

Лемма (woman-pessimality):

Из всех возможных решений алгоритмом Гейла-Шепли будет найдено решение, наихудшее для женщин.

Доказательство:

[math]\triangleright[/math]

Пусть [math]A[/math], [math]B[/math] — мужчины; [math]c[/math], [math]d[/math] — женщины; [math]A[/math] женат на [math]d[/math], [math]B[/math] женат на [math]c[/math]. Предположим, [math]\langle A, c\rangle[/math] — стабильная пара в паросочетании [math]S'[/math], найденном алгоритмом Гейла-Шепли, но [math]A[/math] не самый худший выбор для [math]c[/math]. Тогда существует стабильная пара в паросочетании [math]S[/math], в которой [math]c[/math] замужем за [math]B[/math], который менее привлекателен, чем [math]A[/math]. Тогда пусть мужем [math]d[/math] будет [math]A[/math] в паросочетании [math]S[/math]. Получается [math]A[/math] считает [math]c[/math] более привлекательной, чем [math]d[/math]. Соответственно [math]\langle A, c\rangle[/math] — нестабильная пара в паросочетании [math]S[/math]. То есть для [math]c[/math] есть мужчина, который более привлекателен, чем её муж.

[math]\triangleleft[/math]

Обобщения задачи

Интересно, что данная задача не всегда имеет решение, если допустить однополые пары (устойчивого паросочетания может не быть) .

Случай же, когда у нас есть [math]N[/math] мужчин и [math]M[/math] женщин ([math]N \neq M[/math]) легко сводится к описанной выше задаче. Рассмотрим [math]M \gt N[/math] ([math]M \lt N[/math] аналогично). Добавим [math]M - N[/math] фиктивных мужчин, которые являются наименее привлекательными с точки зрения каждой из женщин. Тогда если в найденном алгоритмом Гейла-Шепли паросочетании некоторая женщина будет замужем за таким фиктивным мужчиной, то это будет означать, что она на самом деле осталась без пары.

Также интересна задача о выборе учебного заведения: вместо множества мужчин введем множество университетов, а вместо множества женщин — множество кандидатов, подающих заявления на поступление. Причем в каждом университете есть квота на количество студентов, которое университет может принять. Задача очевидно сводится к основной добавлением [math](K-1)[/math] "филиалов" для каждого университета ([math]K[/math] — квота). И добавлением фиктивного университета (поступление в который означает, что кандидату придется попробовать поступить через год).

Применения в реальной жизни

Задача о нахождении устойчивого паросочетания и её решение имеют множество применений в реальной жизни, лишь некоторые из них:

• Распределение студентов по коллеждам в США
• Распределение интернов по больницам
• Распределение донорских органов по нуждающимся в них людям

Решение данной задачи было отмечено при вручении Нобелевской премии по экономике в [math]2012[/math] году за «теорию стабильного распределения и практическое применение рыночных моделей». Её получили один из создателей алгоритма, Ллойд Шепли, а также Элвин Рот, во многом развивший исследования Ллойда Шепли и Дэвида Гейла. Сам Гейл не был удостоен премии, вероятно, лишь в силу того, что умер в [math]2008[/math] году.

См. также

• Паросочетания

Источники информации

• Stable matching, Prinston lecture's presentation
• Stable marriage problem
• Задача о марьяже
• Устойчивость супружеских пар и другие комбинаторные задачи (Статья Дональда Кнута)
• Устойчивость супружеских пар и другие комбинаторные задачи - Казанский государственный университет