HD-2063 匈牙利算法解决二分图问题

最新推荐文章于 2020-12-23 17:43:47 发布

JobsandCzj

最新推荐文章于 2020-12-23 17:43:47 发布

阅读量662

点赞数 1

分类专栏： ACM-算法与刷题文章标签： ACM算法匈牙利算法二分图 C++

ACM-算法与刷题专栏收录该内容

45 篇文章 0 订阅

订阅专栏

过山车

Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 6950 Accepted Submission(s): 3019

Problem Description

RPG girls今天和大家一起去游乐场玩，终于可以坐上梦寐以求的过山车了。可是，过山车的每一排只有两个座位，而且还有条不成文的规矩，就是每个女生必须找个个男生做partner和她同坐。但是，每个女孩都有各自的想法，举个例子把，Rabbit只愿意和XHD或PQK做partner，Grass只愿意和linle或LL做partner，PrincessSnow愿意和水域浪子或伪酷儿做partner。考虑到经费问题，boss刘决定只让找到partner的人去坐过山车，其他的人，嘿嘿，就站在下面看着吧。聪明的Acmer，你可以帮忙算算最多有多少对组合可以坐上过山车吗？

Input

输入数据的第一行是三个整数K , M , N，分别表示可能的组合数目，女生的人数，男生的人数。0<K<=1000
1<=N 和M<=500.接下来的K行，每行有两个数，分别表示女生Ai愿意和男生Bj做partner。最后一个0结束输入。

Output

对于每组数据，输出一个整数，表示可以坐上过山车的最多组合数。

Sample Input

3 3

1

2

3

1

3

1

0

Sample Output

3

    题意： 
  

    1.有n个女生和m个男生去，做过山车 
  

    2.女生一定要搭一个男生 
  

    3.女生有要求对于男生，有K个情况 
  

    4.请问最多可以有几对可以玩过 
   山车 
  

    思路： 
  

    1.这明显是一个要建图的题目，我一开始想用回溯法，列出所有项，然后找到一个最多的情况，然而其时间复杂度O(m*n*n)，显然当数据量过大时，不太适用 
  

    2.想一想，这题讲的就是女生和男生的关系，而女女与男男之间是没多大关系，因而我觉得可以用 匈牙利算法来求解，时间复杂度O（m*n），显然适合。 
  

    匈牙利算法（ 
   用增广路径求二分图最大匹配的算法 
   ） 
  

 
通过两者的关系，从头开始找，如果本身没有匹配过，直接匹配 
  

 
如果之前匹配过，然后站在后者的角度，去找可以匹配的，然后进行互换 
  

 
然后从头到尾一遍搜索，就可以找到你想要的答案 
  

 
   例（案例）： 
  

3 3

1

2

3

1

3

3 1

我们先来讲讲匈牙利算法：

【书本上的算法往往讲得非常复杂，我和我的朋友计划用一些简单通俗的例子来描述算法的流程】

匈牙利算法是由匈牙利数学家Edmonds于1965年提出，因而得名。匈牙利算法是基于Hall定理中充分性证明的思想，它是部图匹配最常见的算法，该算法的核心就是寻找增广路径，它是一种用增广路径求二分图最大匹配的算法。

-------等等，看得头大？那么请看下面的版本：

通过数代人的努力，你终于赶上了剩男剩女的大潮，假设你是一位光荣的新世纪媒人，在你的手上有N个剩男，M个剩女，每个人都可能对多名异性有好感（-_-||暂时不考虑特殊的性取向），如果一对男女互有好感，那么你就可以把这一对撮合在一起，现在让我们无视掉所有的单相思（好忧伤的感觉快哭了），你拥有的大概就是下面这样一张关系图，每一条连线都表示互有好感。

本着救人一命，胜造七级浮屠的原则，你想要尽可能地撮合更多的情侣，匈牙利算法的工作模式会教你这样做：

===============================================================================

一：先试着给1号男生找妹子，发现第一个和他相连的1号女生还名花无主，got it，连上一条蓝线

===============================================================================

二：接着给2号男生找妹子，发现第一个和他相连的2号女生名花无主，got it

===============================================================================

三：接下来是3号男生，很遗憾1号女生已经有主了，怎么办呢？

我们试着给之前1号女生匹配的男生（也就是1号男生）另外分配一个妹子。

(黄色表示这条边被临时拆掉)

与1号男生相连的第二个女生是2号女生，但是2号女生也有主了，怎么办呢？我们再试着给2号女生的原配()重新找个妹子(注意这个步骤和上面是一样的，这是一个递归的过程)

此时发现2号男生还能找到3号女生，那么之前的问题迎刃而解了，回溯回去

2号男生可以找3号妹子~~~ 1号男生可以找2号妹子了~~~ 3号男生可以找1号妹子

所以第三步最后的结果就是：

===============================================================================

四：接下来是4号男生，很遗憾，按照第三步的节奏我们没法给4号男生腾出来一个妹子，我们实在是无能为力了……香吉士同学走好。

===============================================================================

这就是匈牙利算法的流程，其中找妹子是个递归的过程，最最关键的字就是“ 腾”字

其原则大概是：有机会上，没机会创造机会也要上

【code】

[cpp]view plain copy 
     
 bool find(int x){  
     int i,j;  
     for (j=1;j<=m;j++){    //扫描每个妹子  
         if (line[x][j]==true && used[j]==false)        
         //如果有暧昧并且还没有标记过(这里标记的意思是这次查找曾试图改变过该妹子的归属问题，但是没有成功，所以就不用瞎费工夫了）  
         {  
             used[j]=1;  
             if (girl[j]==0 || find(girl[j])) {   
                 //名花无主或者能腾出个位置来，这里使用递归  
                 girl[j]=x;  
                 return true;  
             }  
         }  
     }  
     return false;  
 }  

在主程序我们这样做：每一步相当于我们上面描述的一二三四中的一步

[cpp]view plain copy 
     
 for (i=1;i<=n;i++)  
 {  
     memset(used,0,sizeof(used));    //这个在每一步中清空  
     if find(i) all+=1;  
 }  

我们再回到这题。

    在此，我就用图来表示 
  

    初始化： 
  

杭电 <wbr>acm <wbr>2063 <wbr>( <wbr>过山车 <wbr>)

0对，无匹配

第一遍

1对， a--A

第二遍

2对， a--B ，b--A

第三遍

3对， a--B ，b--C ，c--A

整个流程就是这样子的

现在应该思路清晰了，但是如何编程呢

编程（匈牙利算法）：

bool 寻找从k出发的对应项出的可增广路

{

while (从邻接表中列举k能关联到顶点j)

{

if (j不在增广路上)

{

把j加入增广路;

if (j是未盖点或者从j的对应项出发有可增广路)

{

修改j的对应项为k;

则从k的对应项出有可增广路,返回true;

}

则从k的对应项出没有可增广路,返回false;

}

void 匈牙利hungary()

{

for i->1 to n

{

if (则从i的对应项出有可增广路)

匹配数++;

}

输出匹配数;

}

#include <iostream>
using namespace std;
bool map[505][505],d[505];
int n,m;
int s[505];
bool find(int a)
{
    int i;
    for(i=1;i<=m;i++)
    {
        if(!d[i]&&map[a][i])
        {
            d[i]=1;
            if(!s[i]||find(s[i]))
            {
                s[i]=a;
                return true;
            }
        }
    }
    return false;
}
int main()
{
    int t,a,b,i,ans;
    while(cin>>t,t)
    {
        ans=0;
        memset(map,0,sizeof(map));
        memset(s,0,sizeof(s));
        cin>>n>>m;
        while(t--)
        {
            cin>>a>>b;
            map[a][b]=true;
        }
        for(i=1;i<=n;i++)
        {
            memset(d,0,sizeof(d));
            if(find(i))
                ans++;
        }
        cout<<ans<<endl;
    }
    return 0;
}