HDU2063-过山车[Hopcroft-Carp]二分图匹配

过山车

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

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

6 3 31 11 21 32 12 33 10

 

Sample Output

3

 

Author

PrincessSnow

 

Source

RPG专场练习赛

 

Recommend

lcy

二分图匹配

#include<bits/stdc++.h>
using namespace std;
const int maxn = 1005;
const int INF = 0x3f3f3f3f;
vector<int> G[maxn];
/// mx  是 左边点集所匹配的右边的点 ，同理 my 是 右边点集所匹配的左边的点集，dx是源点到x点集的距离，dy是源点到y点集的距离。
int mx[maxn],my[maxn],dx[maxn],dy[maxn];
bool vis[maxn];
int cntx,cnty;
int dis;
bool SearchPath() {
    //bool flag = false;
    dis = INF;
    queue<int> qu;
    memset(dx,-1,sizeof(dx));
    memset(dy,-1,sizeof(dy));
    for(int i=1;i<=cntx;i++)if(mx[i] == -1) { /// 以左点集(未匹配的点)作为源点集，
        qu.push(i);dx[i] = 0;
    }
    while(!qu.empty())
    {
        int u = qu.front();qu.pop();
        if(dx[u] > dis) break;
        for(int i=0;i<G[u].size();i++) {
            int to = G[u][i];
            if(dy[to] == -1) { /// 如果to还没有访问到
                dy[to] = dx[u] + 1;
                if(my[to] != -1) {  ///如果to已经有了匹配的对象
                    dx[my[to]] = dy[to] + 1; /// 距离加1
                    qu.push(my[to]); /// 将to的对象加入队列.
                }
                //else flag = true; ///to没有匹配的对象，则找到了增广路.
                else dis = dy[to];
            }
        }
        //if(flag) break;
    }
    return dis != INF;
}
bool dfs(int u) { /// 类似与匈牙利的find函数
    for(int i=0;i<G[u].size();i++) {
        int to = G[u][i];
        if(!vis[to] && dy[to] == dx[u] + 1) ///这里不同，只需要距离为1的匹配。
        {
            vis[to] = true;
            if(my[to]!=-1&&dy[to]==dis) continue;
            if(my[to]==-1||dfs(my[to])) { /// 如果to没有匹配，或者to的匹配节点还能找到另一个匹配.
                mx[u] = to;               /// 改变匹配关系。
                my[to] = u;
                return true;
            }
        }
    }
    return false;
}
int MaxMatch() {
    int ans = 0;
    memset(mx,-1,sizeof(mx));
    memset(my,-1,sizeof(my));
    while(SearchPath()) /// search the path
    {
        memset(vis,0,sizeof(vis));
        for(int i=1;i<=cntx;i++) {   /// dfs is match
            if(mx[i] == -1 && dfs(i)) ans++;
        }
    }
    return ans;
}
void init()
{
    for(int i=0;i<maxn;i++) G[i].clear();
}
int main()
{
    int k;
    while(~scanf("%d",&k),k) {
        init();
        scanf("%d%d",&cntx,&cnty);
        for(int i=0,u,v;i<k;i++) {
            scanf("%d%d",&u,&v);
            G[u].push_back(v);
        }
        printf("%d\n",MaxMatch());
    }
    return 0;
}