HDOJ 2063 过山车【匈牙利算法求二分图最大匹配 DFS增广 BFS增广】

过山车

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

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 3
1 1
1 2
1 3
2 1
2 3
3 1
0
  
  

 

Sample Output

  
  

   
   3
  
  

 

Author

PrincessSnow

 

Source

RPG专场练习赛

 

二分图匹配裸题。匈牙利算法，DFS增广  和 BFS增广

#include <iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
#include<algorithm>
#define maxn 1010
using namespace std;
int g[maxn][maxn];
int k,nx,ny;
int cx[maxn],cy[maxn];
int a[maxn],b[maxn];
int mark[maxn];
int find(int u)//这个较好理解，不理解的话，拿笔在本上自己画画就懂了
{
    for(int v=1;v<=k;v++)
    {
        int t=b[v];
        if(g[u][t]&&!mark[t])//相连接并且在本次查找当中未使用过，mark数组每次查找都会重新赋值为0
        {
            mark[t]=1;
            if(cy[t]==-1||find(cy[t]))//要么当前的t没用过，或者可以将t腾出来
            {
                cx[u]=t;
                cy[t]=u;
                return 1;
            }
        }
    }
    return 0;
}
int maxmatch()
{
    int res=0;
    memset(cx,-1,sizeof(cx));
    memset(cy,-1,sizeof(cy));
    for(int i=1;i<=k;++i)
    {
        if(cx[a[i]]==-1)
        {
            memset(mark,0,sizeof(mark));
            res+=find(a[i]);
        }
    }
    return res;
}
int main()
{
    while(scanf("%d",&k),k)
    {
        scanf("%d%d",&nx,&ny);
        memset(g,0,sizeof(g));
        memset(a,0,sizeof(a));
        memset(b,0,sizeof(b));
        for(int i=1;i<=k;++i)
        {
            scanf("%d%d",&a[i],&b[i]);
            g[a[i]][b[i]]=1;
        }
        int ans=maxmatch();
        printf("%d\n",ans);
    }
    return 0;
}

这个是BFS没上面的容易理解

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
#define maxn 1010
using namespace std;
int g[maxn][maxn];
int pred[maxn];
int a[maxn],b[maxn];
int cx[maxn],cy[maxn];
int k;
int mark[maxn];
int maxmatch()
{
    int res=0,tx,ty,y;
    memset(cx,-1,sizeof(cx));
    memset(cy,-1,sizeof(cy));
    for(int i=1;i<=k;++i)
    {
        tx=a[i];
        if(cx[tx]!=-1)
            continue;
        queue<int>q;
        memset(mark,0,sizeof(mark));
        for(int j=1;j<=k;++j)
        {
            ty=b[j];
            pred[ty]=-2;//初始化标记为-2
        }
        for(int j=1;j<=k;++j)
        {
            ty=b[j];
            if(g[tx][ty]&&!mark[ty])//这里mark 数组是为了一个值只入队列一次
            {
                pred[ty]=-1;//与tx相连标记为-1并入队列
                q.push(ty);
                mark[ty]=1;
            }
        }
        while(!q.empty())
        {
            y=q.front();
            if(cy[y]==-1)
                break;
            q.pop();//这个放到后面是为了之后判断是否找到增广路
            for(int j=1;j<=k;++j)
            {
                ty=b[j];
                if(pred[ty]==-2&&g[cy[y]][ty]&&!mark[ty])//寻找交替增广路，
                {
                    pred[ty]=y;
                    q.push(ty);
                    mark[ty]=1;
                }
            }
        }
        if(q.empty())
            continue;
        while(pred[y]>-1)//将找到的交替增广路反转
        {
            cx[cy[pred[y]]]=y;
            cy[y]=cy[pred[y]];
            y=pred[y];
        }
        cy[y]=tx;
        cx[tx]=y;
        res++;
    }
    return res;
}
int main()
{
    int m,n;
    while(scanf("%d",&k),k)
    {
        scanf("%d%d",&m,&n);
        memset(g,0,sizeof(g));
        memset(a,0,sizeof(a));
        memset(b,0,sizeof(b));
        for(int i=1;i<=k;++i)
        {
            scanf("%d%d",&a[i],&b[i]);
            g[a[i]][b[i]]=1;
        }
        int ans=maxmatch();
        printf("%d\n",ans);
    }
    return 0;
}