Stable Matching - Gale Shapely Algorithm

基本定义

Matching: A matching S = { ( m , w ) ∣ m ∈ M  and  w ∈ W } S = \{(m, w) | m \in M \text{ and } w \in W\} S={(m,w)∣m∈M and w∈W}, where each $m$ and $w$ appears in at most one pair in $S$.

Perfect Matching: A perfect matching S ′ = { ( m , w ) ∣ m ∈ M  and  w ∈ W } S' = \{(m, w) | m \in M \text{ and } w \in W\} S′={(m,w)∣m∈M and w∈W}, where each $m$ and $w$ appears in exactly one pair in $S^{\prime }$.

Preference: Each man $m\in M$ ranks all the women. We say "$m$ prefers $w$ to $w^{\prime }$" if $m$ ranks $w$ higher than $w^{\prime }$. Similarly, each woman $w\in W$ ranks all the men.

Instability: An instability with respect to $S$: $(m,w^{\prime })\notin S$, but $m$ prefers $w^{\prime }$ to his partner and $w^{\prime }$ prefers $m$ to her partner.

Stable Matching: A matching $S$ is stable if
\quad\quad\quad\quad\quad\quad 1. $S$ is perfect
\quad\quad\quad\quad\quad\quad 2. no instability with respect to $S$

Valid Partner: A woman $w$ is a valid partner of a man $m$ if there is a stable matching that contains pair $(m,w)$.

Best Valid Partner: $w$ is the best valid partner of $m$ if
\quad\quad\quad\quad\quad\quad 1. $w$ is a valid of $m$
\quad\quad\quad\quad\quad\quad 2. $\forall w^{\prime }>w$, $m$ and $w^{\prime }$ is not a valid partner
Denote $best(m)$

$w^{\prime }>w$ indicates $m$ ranks $w^{\prime }$ higher than $w$.

伪代码

1. Initially all $m\in M$ and $w\in W$ are free
2. While $\exists m$ who is free and hasn’t proposed to any woman
3. \quad\quad Choose such a man $m$
4. \quad\quad Let $w$ be the highest-ranked woman in $m$'s preference list to whom $m$ has not yet proposed
5. \quad\quad If $w$ is free then
6. \quad\quad\quad\quad $(m,w)$ become engaged
7. \quad\quad Else ($w$ is currently engaged to $m^{\prime }$)
8. \quad\quad\quad\quad If $w$ prefers $m^{\prime }$ to $m$ then
9. \quad\quad\quad\quad\quad\quad $m$ remains free
10. \quad\quad\quad\quad Else ($w$ prefers $m^{\prime }$ to $m$)
11. \quad\quad\quad\quad\quad\quad $(m,w)$ become engaged
12. \quad\quad\quad\quad\quad\quad $m^{\prime }$ becomes free
13. \quad\quad\quad\quad Endif
14. \quad\quad Endif
15. Endwhile
16. Return the set $S$ of engaged pair

GS算法定理

定理1：在 $w$ 获取第一个 proposal 之后，$w$ 将一直有配偶。她得到的配偶会越来越好 (in terms of her preference list)。

定理2：$m$ propose 的配偶越来越差 (in terms of his preference list)。

定理3：GS 算法在至多 $n^2$ 次 iterations 之后终止。
证明：
令 $P(t)$ 为 { ( m , w ) ∣ m  has proposed to  w } \{(m,w)| m \text{ has proposed to } w \} {(m,w)∣m has proposed to w}
$\because$ 每次循环男人向不同的女人求婚
$\therefore |P(t+1)|>|P(t)|$
一共最多有 $n^2$ 个男女配对。

定理4：如果在算法运行过程中 $m$ 是 free 的，那么存在一个 $m$ 还没有propose 过的 $w$ 。
证明：（反证法）
假设在一个时刻 $m$ is free but has already proposed to every women. 由定理1可知每一个女人都有配偶。
$\because$ the engaged pairs form a matching $n\times n$
$\therefore$ there are $n$ engaged men at this point
矛盾

定理5：程序终止时 return 的 set $S$ 是一个 perfect matching。
证明：
The set of engaged pairs always forms a matching.
假设这个算法在终止的时候存在一个 free man $m$. 那么，这个男人 $m$ 肯定向所有女人都求过婚，否则 while loop 不会结束。
与定理4矛盾

定理6：GS 算法 return a set of pairs $S$. 这个集合 $S$ 是 stable matching。
证明：（反证法）
假设集合 $S$ 中存在 instability，即 $\exists (m,w),(m^{\prime },w^{\prime })\in S$ $m$ 更喜欢 $w^{\prime }$ 且 $w^{\prime }$ 更喜欢 $m$。
由算法代码可知，$m$ 最后 propose 的 $w$。
若 $m$ 在之前没有向 $w^{\prime }$ propose 过，那么 $m$ 更喜欢 $w$。（矛盾）
若 $m$ 在之前有向 $w^{\prime }$ propose 过，那么 $m$ 被 $w^{\prime }$ 拒绝过，因为 $w^{\prime }$ 更喜欢 $m^{\prime \prime }$。
$\because m^{\prime }$ 是 $w^{\prime }$ 最终的配偶
$\therefore m^{\prime \prime }=m^{\prime }$ 或 $m^{\prime }>m^{\prime \prime }>m$ （矛盾）
因此得出结论，$S$ 是 stable matching。

定理7：令 S ∗ = { ( m , b e s t ( m ) ) : m ∈ M } S^* = \{ (m, best(m)) : m\in M\} S∗={(m,best(m)):m∈M} ，every execution of the GS Algorithm results in the set $S^{\ast }$。【所有execution（choose free man的顺序不同）得到同一个配对】
证明：（反证法）
假设一个G-S算法的 execution $ϵ$ 得到一个 matching $S$。
其中$\exists m\in S$, $m$ 的配偶不是 $best(m)$ [not best valid pair]。
这说明，$m$ 在之前被一个 valid partner $w$ 拒绝过。
假设 $m$ 是第一个被 valid partner 拒绝的男人。
因为 $m$ propose in decreasing order of preference, $m$ 被 $w=best(m)$ 拒绝。
$w$ 拒绝 $m$ 有两种可能：
（1）$w$ 更喜欢现在的配偶
（2）$w$ 跟 $m$ 配对之后有更喜欢的配偶出现
均形成配对$(w,m^{\prime })$ ，其中$m^{\prime }>m$
$\because$ $(m,w)$是 valid pair
$\therefore$ 存在一个 stable matching $(m,w)\in S^{\prime }$。假设在$S^{\prime }$中，$(m^{\prime },w^{\prime })$。
$\because m$ 是第一个被valid partner $w$ 拒绝的男人
$\therefore m^{\prime }$ 在这之前没有被valid partner 拒绝过
$m^{\prime }$ propose in decreasing order 且 $w^{\prime }$ 是 $m^{\prime }$ 的一个valid partner。
因此对于 $m^{\prime }$，$w>w^{\prime }$。
但是前面我们已经证明：对于 $w$, $m^{\prime }>m$。
$\therefore$ $(m^{\prime },w^{\prime })$ 是一个 unstable pair
这与前面假设 $S^{\prime }$ 矛盾

定理8：在一个 stable matching $S^{\ast }$ 中，所有女人的配偶都是 worst valid partner。
证明：（反证法）
假设 $(m,w)\in S^{\ast }$ 其中 $m$ 不是 $w$ 的 worst valid partner。
那么存在一个 stable matching $(m^{\prime },w)\in S^{\prime }$，其中 $m^{\prime }<m$。
假设在 $S^{\prime }$ 中 $(m,w^{\prime })$.
$\because w$ 是 $m$ 的best valid partner, $w^{\prime }$ 是 $m$ 的 valid partner
$\therefore$ 对于 $m$, $w>w^{\prime }$
发现 $(m,w)\in S^{\prime }$ 是一个unstable pair
与前面假设 $S^{\prime }$ stable 矛盾

代码实现思路

算法时间复杂度 $O(n^2)$，以下是代码思路：

parameter:
MenPref[m, i] denote $i^{th}$ woman on man $m$'s preference list
WomenPref[w, i] denote $i^{th}$ man on woman $w$'s preference list

以下四个操作 $O(1)$:

1. 找到一个 free man
2. 对于一个男人 $m$，找到一个他没有求过婚的女人
3. 知道一个女人 $w$ 是否已经有配对，且这个配偶是谁
4. 对于一个女人 $w$ 和两个男人 $m，$$m^{\prime }$，知道女人更喜欢谁

需要用到的数据结构：

1. $\forall m$, an array Next[m] is the position of the next woman he will propose to on his list. (初始：Next[m]=1)

   *每次 propose MenPrefs[m, Next[m]]
   Next[m]++

2. $\forall w$，an array Current[w] 是女人 $w$ 现在的配偶。

3. 生成一个 $n\times n$ array Ranking[w, m]

   row 是女人（按序号），column 是男人（按序号），entry是此男 人在女人那里的 ranking。
   *每次只需要compare Rank[w, m] 与 Rank[w, m’]即可