HDU1914（稳定婚姻）

http://acm.hdu.edu.cn/showproblem.php?pid=1914

思路：Gale-Shapley算法。算法过程是男士不停地求婚，女士不停地拒绝。在每一轮中，每个尚未订婚的男士在他还没有求过婚的女士中选一个自己最喜欢的求婚（不管她有没有订婚）。然后每个女士在向她求婚的人之中选择她最喜欢的一个订婚，并且拒绝其他人。注意，这些向她求婚的人中包含她的未婚夫，因此她可以选择另一个自己更喜欢的人订婚，而抛弃自己的现任未婚夫。这个算法的结果是使男士都能娶到自己有可能娶到的最好的妻子，所以是对男士的理想配对。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<string>
#include<vector>
#include<map>
#include<set>
#include<queue>
using namespace std;
struct male
{
    int f,rev[],tag;
} m[];
struct female
{
    int tag,temp,val,wait[];
} f[];
int _,n,t,k,mf[][],fm[][];
char ch[];
bool ok()
{
    int i;
    for (i=;i<=;i++)
        if (m[i].f==&&m[i].tag>) return true;
    return false;
}
int main()
{
    scanf("%d",&_);
    while (_--)
    {
        scanf("%d",&n);
        int i,j;
        for (i=;i<=;i++)
        {
            f[i].tag=;
            m[i].tag=;
        }
        for (i=;i<=n;i++)
        {
            scanf("%s",ch);
            t=ch[]-'a'+;
            m[t].f=;
            m[t].tag=t;
            memset(m[t].rev,,sizeof(m[t].rev));
        }
        for (i=;i<=n;i++)
        {
            scanf("%s",ch);
            t=ch[]-'A'+;
            f[t].tag=t;
            f[t].temp=;
            f[t].val=;
            memset(f[t].wait,,sizeof(f[i].wait));
        }
        for (i=;i<=n;i++)
        {
            scanf("%s",ch);
            t=ch[]-'a'+;
            for (j=;j<=n+;j++)
                mf[t][j-]=ch[j]-'A'+;
        }
        for (i=;i<=n;i++)
        {
            //cout<<1;
            scanf("%s",ch);
            t=ch[]-'A'+;
            for (j=;j<=n+;j++)
                fm[t][j-]=ch[j]-'a'+;
        }
        //cout<<"hhhhhhh"<<endl;
        while (ok())
        {
            for (i=;i<=;i++)
            {
                if (m[i].f==&&m[i].tag>)
                {
                    for (j=;j<=n;j++)
                    {
                        t=mf[i][j];
                        if (m[i].rev[t]==)
                        {
                            m[i].rev[t]=;
                            m[i].f=;
                            k=++f[t].wait[];
                            f[t].wait[k]=i;
                            break;
                        }
                    }
                }
            }
            for (i=;i<=;i++)
            {
                if (f[i].tag>)
                {
                    for (j=;j<=f[i].wait[];j++)
                    {
                        t=f[i].wait[j];
                        for (k=;k<=n;k++)
                            if (fm[i][k]==t) break;
                        if (f[i].val>k)
                        {
                            m[f[i].temp].f=;
                            f[i].temp=t;
                            f[i].val=k;
                        }
                        else m[t].f=;
                    }
                    f[i].wait[]=;
                }
            }
        }
        int out[];
        memset(out,,sizeof(out));
        for (i=;i<=;i++)
            if (f[i].tag>)
            {
                j=f[i].temp;
                out[j]=i;
            }
        for (i=;i<=;i++)
            if (out[i]) printf("%c %c\n",i-+'a',out[i]-+'A');
        if (_) printf("\n");
    }
    return ;
}