（算法）稳定婚姻匹配

题目：

婚介所登记了N位男孩和N位女孩，每个男孩都对N个女孩的喜欢程度做了排序，每个女孩都对N个男孩的喜欢程度做了排序，你作为月老，能否给出稳定的牵手方案？

稳定的定义：如果男孩i和女孩a牵手，但男孩i对女孩b更喜欢，而女孩b的男朋友j拼不过男孩i，则没有力量阻碍男孩i和女孩b的私奔，这即是不稳定的。

思路：

   1962 年，美国数学家 David Gale 和 Lloyd Shapley 发明了一种寻找稳定婚姻的策略。不管男女各有多少人，不管他们各自的偏好如何，应用这种策略后总能得到一个稳定的婚姻搭配。换句话说，他们证明了稳定的婚姻搭配总是存在的。有趣的是，这种策略反映了现实生活中的很多真实情况。

    算法中采用了男生主动追求女孩的形式。

算法步骤描述：

    第一轮，每个男人都选择自己名单上排在首位的女人，并向她表白。这种时候会出现两种情况：

　（1）该女士还没有被男生追求过，则该女士接受该男生的请求。

　（2）若该女生已经接受过其他男生的追求，那么该女生会将该男士与她的现任男友进行比较，若更喜欢她的男友，那么拒绝这个人的追求，否则，抛弃其男友

    第一轮结束后，有些男人已经有女朋友了，有些男人仍然是单身。

    在第二轮追女行动中，每个单身男都从所有还没拒绝过他的女孩中选出自己最中意的那一个，并向她表白，不管她现在是否是单身。这种时候还是会遇到上面所说的两种情况，还是同样的解决方案。直到所有人都不再是单身。

怎么证明这个算法肯定能够得到稳定的婚姻？

（1）随着轮数的增加，总有一个时候所有人都能配上对。因为男生根据自己心目中的排名依次对女士进行表白，假如有一个人没有配上对，那么这个人必定是向所有的女孩进行表白了。但是女孩只要被表白过一次，就不可能是单身，也就是说此时所有的女生都不是单身的，这与有一个人没有配上对是相悖的。所以假设不成立。该算法一定会使得所有人都能够配对成功。

（2）随着轮数的增加，男士追求的对象越来越糟，而女士的男友则可能变得越来越好。假设男A和女1各有各自的对象，但是比起现在的对象，男A更喜欢女1，所以，在此之前男A肯定已经跟女1表白过的，并且女1拒绝了男A，也就是女1有了比男A更好的男友，不会出现私奔的情况……。

代码：

#include<iostream>
using namespace std;

const int N=4;

void GaleShapley(const int (&man)[N][N],const int (&woman)[N][N],int (&match)[N]){
    int wm[N][N];    // wm[i][j]: rank from girl i to boy j
    int choose[N];    // choose[i]: current boyfriend of girl i
    int manIndex[N]; //    manIndex[i]: how many girls that have rejected boy i
    int i,j;
    int w,m;
    for(i=0;i<N;i++){
        match[i]=-1;
        choose[i]=-1;
        manIndex[i]=0;
        for(j=0;j<N;j++)
            wm[i][woman[i][j]]=j;
    }

    bool bSingle=false;
    while(!bSingle){
        bSingle=true;
        for(i=0;i<N;i++){
            if(match[i]!=-1) // boy i already have a girlfriend
                continue;
            bSingle=false;
            j=manIndex[i]++; // the jth girl that boy i like most
            w=man[i][j];    
            m=choose[w];    // current girl w's boyfriend
            if(m==-1 || wm[w][i]<wm[w][m]){ // if girl w prefer boy i
                match[i]=w;
                choose[w]=i;
                if(m!=-1)
                    match[m]=-1;
            }
        }
    }
}

void Print(const int (&match)[N],int N){
    for(int i=0;i<N;i++)
        cout<<"Boy "<<i<<" matches "<<"Girl "<<match[i]<<endl;
    cout<<endl;
}

int main(){
    int man[N][N]={
        {2,3,1,0},
        {2,1,3,0},
        {0,2,3,1},
        {1,3,2,0},
    };
    int woman[N][N]={
        {0,3,2,1},
        {0,1,2,3},
        {0,2,3,1},
        {1,0,3,2},
    };

    int match[N];
    GaleShapley(man,woman,match);
    Print(match,N);

    return 0;
}