【模板整合】匈牙利算法和Hopcroft-Karp算法

匈牙利的DFS实现

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MAXN 1010
using namespace std;
int n,m;
int map[MAXN][MAXN],num[MAXN];
bool vis[MAXN];
int ans;
bool find(int x)
{
    for (int i=1;i<=n;i++)
    {
        if (!vis[i]&&map[x][i])
        {
            vis[i]=1;
            if (!num[i]||find(num[i]))
            {
                num[i]=x;
                return 1;
            }
        }
    }
    return 0;
}
int main()
{
    scanf("%d%d",&n,&m);
    /*for (int i=1;i<=m;i++)
    {
        int x,y;
        scanf("%d%d",&x,&y);x++;y++;
        map[i][x]=map[i][y]=1;
    }建图请自动忽略*/
    for (int i=1;i<=m;i++)
    {
        memset(vis,0,sizeof(vis));
        if (find(i))    ans++;
        else    break;
    }
    printf("%d\n",ans);
}

Hopcroft-Karp算法

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
#define MAXN 1010
#define MAXINT 0x7fffffff
using namespace std;
int n,m,top,x,y;
int ans;
int disx[MAXN],disy[MAXN],matx[MAXN],maty[MAXN];//x,y,分别为二分图的两个点集,mat为每个点在对侧集合的匹配点,如果当前没有匹配点则为-1 
struct edge
{
    int to;
    edge *next;
}e[MAXN*MAXN],*prev[MAXN];
void insert(int u,int v)
{
    e[++top].to=v;e[top].next=prev[u];prev[u]=&e[top];
}
bool bfs()//寻找最短增广路 
{
    bool ret=0;
    queue<int> q;
    memset(disx,0,sizeof(disx));memset(disy,0,sizeof(disy));
    for (int i=1;i<=n;i++)
        if (matx[i]==-1)    q.push(i);//找到未盖点,入队 
    while (!q.empty())//在二分图另一个点集的非盖点中寻找增广路 
    {
        int x=q.front();q.pop();
        for (edge *i=prev[x];i;i=i->next)
            if (!disy[i->to])
            {
                disy[i->to]=disx[x]+1;
                if (maty[i->to]==-1)    ret=1;//找到增广路 
                else    disx[maty[i->to]]=disy[i->to]+1,q.push(maty[i->to]);    
            }
    }
    return ret;
}
bool dfs(int x)//沿增广路增广 
{
    for (edge *i=prev[x];i;i=i->next)
    {
        if (disy[i->to]==disx[x]+1)
        {
            disy[i->to]=0;
            if (maty[i->to]==-1||dfs(maty[i->to]))
            {
                matx[x]=i->to;maty[i->to]=x;return 1;
            }
        }
    }
    return 0;
}
int main()
{
    scanf("%d%d",&n,&m);
    memset(matx,-1,sizeof(matx));memset(maty,-1,sizeof(maty));
    /*for (int i=1;i<=n;i++)
    {
        scanf("%d%d",&x,&y);x++;y++;
        insert(i,x);insert(i,y);
    }建图请自动忽略*/
    while (bfs())
    {
        for (int i=1;i<=m;i++)
            if (matx[i]==-1&&dfs(i))    ans++;
    }
    cout<<ans<<endl;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值