hdu 1068 Girls and Boys

最新推荐文章于 2013-05-17 18:02:57 发布

KubeSpace

最新推荐文章于 2013-05-17 18:02:57 发布

阅读量612

点赞数

分类专栏：二分图匹配文章标签： hdu

本文链接：https://blog.csdn.net/u010087314/article/details/8922249

版权

二分图匹配专栏收录该内容

12 篇文章 0 订阅

订阅专栏

博客探讨了HDU 1068问题，重点在于解决最大二分匹配问题。通过Hopcroft-Karp算法进行求解，该算法在寻找图的最大独立集时应用，其大小等于图中边数减去最大匹配的数量。

摘要由CSDN通过智能技术生成

最大二分匹配。

题目所求为最大独立集=n-最大匹配

Hopcroft-krap()算法：

#include <iostream>
#include <cstring>
#include <cstdlib>
#include <queue>

using namespace std;
#define INF 10000000

queue<int> q;
int g[505][505];
int vis[505];
int distx[505],disty[505];
int linkx[505],linky[505];
int n;
int dis;

bool bfs()
{
    while( !q.empty() ) q.pop();
    memset(distx,-1,sizeof(distx));
    memset(disty,-1,sizeof(disty));
    for(int i=0;i<n;i++)
        if(linkx[i]==-1)
        {
            q.push(i);
            distx[i]=0;
        }
    dis=INF;
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        if(distx[u]>dis) break;
        for(int i=0;i<n;i++)
            if( disty[i]==-1 && g[u][i] )
            {
                disty[i]=distx[u]+1;
                if(linky[i]==-1) dis=disty[i];
                else
                {
                    distx[linky[i]]=disty[i]+1;
                    q.push(linky[i]);
                }
            }
    }
    return dis<INF;
}

int dfs(int u)
{
    for(int i=0;i<n;i++)
        if( !vis[i] && g[u][i] && disty[i]==distx[u]+1 )
        {
            vis[i]=1;
            if( disty[i]==dis && linky[i]!=-1 )
                continue;
            if( linky[i]==-1 || dfs(linky[i]) )
            {
                linky[i]=u;
                linkx[u]=i;
                return 1;
            }
        }
    return 0;
}

int main()
{
    while(cin>>n)
    {
        memset(g,0,sizeof(g));
        for(int i=0;i<n;i++)
        {
            int a,b;
            char c;
            cin>>a>>c;
            cin>>c;
            cin>>b>>c;
            for(int i=0;i<b;i++)
            {
                int c;
                cin>>c;
  //              if(!g[c][a])                   //为什么加了这个判断就错了，只需要求一次的啊。。。

		      g[a][c]=1;
            }
        }
        int ans=0;
        memset(linkx,-1,sizeof(linkx));
        memset(linky,-1,sizeof(linky));
        while( bfs() )
        {
            memset(vis,0,sizeof(vis));
            for(int i=0;i<n;i++)
                if(linkx[i]==-1)
                    ans+=dfs(i);
        }
        cout<<n-ans/2<<endl;
    }
    return 0;
}