poj 3487 zoj 3120 The Stable Marriage Problem 延迟认可算法（Gale-Shapley算法）

/*
传说中的Gale-Shapley算法。。。
第一次听说这个算法
百度一下才知道此算法已经有其实际应用了。。。
此题从早上八点开始看
先是用了n多队列写
测试倒是通过了
就是提交时老SF
然后就看到了网上说的Gale-Shapley算法
从学习算法然后写代码调试
最后到提交AC时已经是下午4点半
本来以为今天肯定过不了
没先到提交时一次通过
喜悦之情顿生。。。 
不说废话了
下面讲一下思路：
该算法的总体思路是
让男生根据他们对女生的喜欢程度对女生进行追求
女生则根据前来追求的对象决定是接受还是拒绝
接受还是拒绝取决于该女生对该男生的喜欢程度 
具体实现是用两个矩阵
1. libman[][]记录男生喜欢的女生列表
从喜欢到不喜欢降序排列
2.ladyscore[][]记录女生对各个男生所打的分数
分数越高则表明女生对该男生的心仪程度越高 
男生有两种状态
自由男和约会男
我们的目标是让所有的自由男变成约会男
所以用一个栈freeman(队列亦可)装自由男
当栈空时表明约会结束
用数组man[]来记录对应编号的男生下一个约会目标
用数组lady[]来记录对应编号的女生现任男友的编号
为了便于处理
虚构一个不存在的男友，编号最大，分数最低
意即此人为某女生现任男友时只要有人追求
该女生就会毫不留情地抛弃他 
*/
#define LOCAL
#include<iostream>
#include<cstdio>
#include<cstring>
#include<stack>
#include<algorithm> //sort()包含在这个头文件中 ps:在大多数情况下sort()比qsort()要快 
#define N 26
using namespace std;
int main()
{
#ifdef LOCAL
       freopen("input.txt","r",stdin);
       freopen("output.txt","w",stdout);
#endif

    short ncase,couple,libman[N][N],man[N],ladyscore[N][N],lady[N],male,female,i,j,first=1,a[N];
    stack<int> freeman;
    string list;
    char temp;  
    cin>>ncase;
    while(ncase--)
    {
            if(first) first=0; 
            else cout<<endl;                 //打印实例间的空格 
            while(!freeman.empty())  freeman.pop(); //清空栈 
            cin>>couple;
            //得到数据 
            for(i=0;i<couple;i++)
            {
                   cin>>temp;
                   a[i]=temp-'a';   //用a[]记录男生列表 
                   freeman.push(a[i]);   //所有男生(一开始都是自由男)入栈                  
            }  
            sort(a,a+couple);//对所有输入的男生进行排序一边与后面按字典序输出 
            for(i=0;i<couple;i++) cin>>temp;
            //接收女生列表(由于女生和男生是一一对应的所以这里只是将其接收而并不记录) 
            for(i=0;i<couple;i++)//记录男生的心仪女生 
            {
                   cin>>list;
                   for(j=0;j<couple;j++)
                   {
                           libman[i][j]=list[j+2]-'A';                     
                   }                     
            }            
            for(i=0;i<couple;i++)//女生对他们的追求者一一打分 
            {
                   cin>>list;
                   for(j=0;j<couple;j++)
                          ladyscore[i][list[j+2]-'a']=couple-j;              
            }            
            for(i=0;i<couple;i++)
                   ladyscore[i][couple]=0;//对"不存在"的男生打最低分 
            memset(man,0,sizeof(man));
            for(i=0;i<N;i++)  lady[i]=couple;//这里不能用memset()因为couple是变量而不是常数 
            while(!freeman.empty())//进行匹配 
            {
                     male=freeman.top();//用male记录当前追求者 
                     female=libman[male][man[male]];//用female记录male所追求的女生 
                     if(ladyscore[female][male]>ladyscore[female][lady[female]]) //若当前追求者比其现任男友优秀 
                     {
                              freeman.pop();//该男生约会成功跳出自由男行列 
                              if(lady[female]!=couple)//如果该女生有前男友 
                              {freeman.push(lady[female]);man[lady[female]]++;}
                              //将其前男友打入自由男行列，该男生刻意追求他的次喜欢女生 
                              lady[female]=male;//追求者成为该女生现任男友 
                     }                      
                     else if(ladyscore[female][male]<ladyscore[female][lady[female]]) //若当前追求者没有其现任男友优秀 
                              man[male]++;//该男生可以追求其次喜欢女生 
            }
            for(i=0;i<couple;i++) //输出约会结果 
                     printf("%c %c\n",a[i]+'a',libman[a[i]][man[a[i]]]+'A');
                               
    }
    return 0;
}