二分图匹配(最小顶点覆盖,多重二分图匹配)

二分图匹配:

最基础的一题:http://acm.hdu.edu.cn/showproblem.php?pid=2063

匈牙利算法的核心是使用used[],每次完成一个选取者的搜索前需要重置used[]

利用used[i]=1和found(boy[i]),来让已经选择了的点绕开这次需要选择的点

 input:

4 3 3

1 1

1 2

2 3

3 1

模拟算法:

1选1    boy[1]=1;

2选3    boy[3]=2;

3想选1    used[1]=1;发现boy[1]=1,尝试found(boy[1]),

found(1),想选1,但是used[1]=1;所以只能绕开它,选择其它的,发现可以选择2(used[2]=0)且boy[2]=0;

boy[2]=1;return 1,然后可以boy[3]=1;

#include<stdio.h>
#include<string.h>
int line[510][510];///矩阵保存关系
int boy[510];
int used[510];
int n,m;
                                      ///*2和*3模拟能否绕开选取
int found(int x)
{
    int i,j;
    for(i=1;i<=m;i++)                 ///遍历所有被选者的编号
    {
        if(line[x][i]==1&&used[i]==0) ///如果 x对i有好感且在这一个阶段选取阶段没有被选取(||假设选取)*1
        {                             
            used[i]=1;                  ///标记被选取(看其它的能否绕开这个选择)*2
            if(boy[i]==0||found(boy[i]))///如果被选者没有归属或他的归属着可以调换
            {                          
                boy[i]=x;                ///将归属定为 x
                return 1;
            }
        }
    }
    return 0;
}
 
int main()
{
    int i,j,k,x,y,sum;
    while(scanf("%d %d %d",&k,&n,&m),k!=0)
    {
       memset(line,0,sizeof(line));
       memset(boy,0,sizeof(boy));
       memset(used,0,sizeof(used));
       for(i=0;i<k;i++)
       {
          scanf("%d %d",&x,&y);
          line[x][y]=1;
       }
       sum=0;          
       for(i=1;i<=n;i++)///遍历选取者
       {
           memset(used,0,sizeof(used));///每次都要清 0
           if(found(i)) sum++;         
       }
       printf("%d\n",sum);
    }
    return 0;
}

邻接表模板:

#include<bits/stdc++.h>
#define MAXN 100005
using namespace std;
vector<int> vc[MAXN];
int boy[510];
int used[510];
int n,m;

int found(int x)
{
    for(int i=0;i<vc[x].size();i++){
        int v=vc[x][i];
        if(used[v]==0){
            used[v]=1;
            if(boy[v]==0||found(boy[v])){
                boy[v]=x;
                return 1;
            }
        }
    }
    return 0;
}

int main()
{
    int k,x,y,sum;
    while(scanf("%d %d %d",&k,&n,&m),k!=0)
    {
        for(int i=0;i<MAXN;i++)
            vc[i].clear();
        memset(boy,0,sizeof(boy));
        memset(used,0,sizeof(used));
        for(int i=0;i<k;i++)
        {
            scanf("%d%d",&x,&y);
            vc[x].push_back(y);
        }
        sum=0;
        for(int i=1;i<=n;i++)///遍历选取者
        {
            memset(used,0,sizeof(used));///每次都要清 0
            if(found(i))
            sum++;
        }
        printf("%d\n",sum);
    }
    return 0;
}

最小顶点覆盖数:

Muddy Fields

用宽度为1长度不限的木板将水洼‘*’盖住而不盖住草‘.'

用匈牙利算法求最大匹配数,需要自己建图

最小覆盖点数=最大匹配数

横着盖:(图中数字表示用的是第几块木板),rx为至少的横板数

1.2.

.333

444.

. . 5.

竖着盖:ry为至少的竖板数

1.2.

.324

532.

. . 2.

rx与ry匹配

AC:

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define MAXN 3005
using namespace std;
char s[MAXN][MAXN];
int r[MAXN][MAXN],c[MAXN][MAXN];
int mp[MAXN][MAXN];
int boy[MAXN],used[MAXN];
int rx,ry;

int found(int x)
{
    for(int i=1;i<=ry;i++){//ry第二种匹配物的个数
        if(mp[x][i]==1&&used[i]==0){
            used[i]=1;
            if(boy[i]==0||found(boy[i])){
                boy[i]=x;
                return 1;
            }
        }
    }
    return 0;
}

int main()
{
    int n,m;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
        scanf("%s",s[i]+1);
    rx=0,ry=0;
    for(int i=1;i<=n;i++){//建图
        for(int j=1;j<=m;j++){
            if(s[i][j]=='*'){
                ++rx;
                while(j<=m&&s[i][j]=='*')
                    r[i][j]=rx,j++;
                j=j-1;
            }
        }
    }
    for(int j=1;j<=m;j++){
        for(int i=1;i<=n;i++){
            if(s[i][j]=='*'){
                ++ry;
                while(i<=n&&s[i][j]=='*')
                    c[i][j]=ry,i++;
                i=i-1;
            }
        }
    }
    for(int i=1;i<=n;i++){//建边
        for(int j=1;j<=m;j++){
            if(s[i][j]=='*'){
                mp[r[i][j]][c[i][j]]=1;
            }
        }
    }
    int ans=0;
    for(int i=1;i<=rx;i++)//第一种匹配物个数
    {
        memset(used,0,sizeof(used));
        if(found(i))
            ans++;
    }
    printf("%d\n",ans);
    return 0;
}

POJ - 3041

求最少需要几枪可以击碎所以小行星->求该图的最小顶点覆盖数

3 4               out:2
1 1
1 3
2 2
3 2

由König定理定理可知,二分图的最小顶点覆盖数等于二分图的最大匹配数

ac:

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<vector>
#define MAXN 100005
using namespace std;
vector<int> vc[MAXN];
int boy[510];
int used[510];
int n,m;

int found(int x)
{
    for(int i=0;i<vc[x].size();i++){
        int v=vc[x][i];
        if(used[v]==0){
            used[v]=1;
            if(boy[v]==0||found(boy[v])){
                boy[v]=x;
                return 1;
            }
        }
    }
    return 0;
}

int main()
{
    int k,x,y,sum;
    while(scanf("%d %d",&n,&k)!=EOF)
    {
        for(int i=0;i<MAXN;i++)
            vc[i].clear();
        memset(boy,0,sizeof(boy));
        memset(used,0,sizeof(used));
        for(int i=0;i<k;i++)
        {
            scanf("%d%d",&x,&y);
            vc[x].push_back(y);
        }
        sum=0;
        for(int i=1;i<=n;i++)///遍历选取者
        {
            memset(used,0,sizeof(used));///每次都要清 0
            if(found(i))
            sum++;
        }
        printf("%d\n",sum);
    }
    return 0;
}

多重二分图匹配:

Jamie's Contact Groups

题意:

n个人,m个通讯录

给定每个人可以去若干个通讯录

求所以通讯录中人数的最大值的最小情况

解析:

最大求最小

二分+二分图多重匹配

ac:

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <vector>
#include <string>
#define MAXN 1005
using namespace std;
int n,m,mid,x;
char str[20];
int vis[MAXN];
vector<int> vc[MAXN],isin[MAXN];
 
bool Find(int u)
{
    int len=vc[u].size();
    for(int i=0;i<len;i++)
    {
        int v=vc[u][i];
        if(!vis[v])
        {
            vis[v]=1;
            int blen=isin[v].size();
            if(blen<mid)//可以装
            {
                isin[v].push_back(u);
                return 1;
            }
            for(int j=0;j<blen;j++)
            {
                if(Find(isin[v][j]))
                {
                    isin[v][j]=u;
                    return 1;
                }
            }
        }
    }
    return 0;
}
 
bool maxmatch()
{
    for(int i=0;i<m;i++)
        isin[i].clear();
    for(int i=0;i<n;i++)
    {
        memset(vis,0,sizeof(vis));
        if(!Find(i))
            return 0;//匹配不了,出结果
    }
    return 1;
}
 
int main()
{
    while(~scanf("%d%d",&n,&m)&&(n+m))
    {
        for(int i=0;i<=n;i++)
            vc[i].clear();
        for(int i=0;i<n;i++)
        {
            scanf("%s",str);
            while(scanf("%d",&x)==1)
            {
                vc[i].push_back(x);
                if(getchar()=='\n') break;
            }
        }
        int l=1,r=n,ans=n;
        while(l<=r)//二分找最小
        {
            mid=(l+r)>>1;
            if(maxmatch())
            {
                r=mid-1;
                ans=mid;
            }
            else l=mid+1;
        }
        printf("%d\n",ans);
    }
    return 0;
}

输入数据包含多组测试用例,每组数据的第一行输入n,表示男生的数量,接下来有n行,每行n个数表示第i个女生对第j个男生的好感度(n<=300)

模板代码:

#include<iostream>//slack:匹配
#include<algorithm>//match:松弛
#include<cstring>//gap:差距
#include<cstdio>
#define inf 0x3f3f3f3f
using namespace std;
 
int n;
int love[301][301];
int vis_girl[301],vis_boy[301];
int ex_girl[301],ex_boy[301];
int match[301],slack[301];
 
int dfs(int x)
{
    vis_girl[x]=1;
    for(int i=0;i<n;i++)
    {
        if(vis_boy[i]) continue;
        int gap=ex_girl[x]+ex_boy[i]-love[x][i];
        if(!gap){
            vis_boy[i]=1;
            if(match[i]==-1||dfs(match[i])){
                match[i]=x;
                return 1;
            }
        }
        else slack[i]=min(slack[i],gap);
    }
    return 0;
}
 
void km()
{
    memset(match,-1,sizeof(match));
    memset(ex_boy,0,sizeof(ex_boy));
    for(int i=0;i<n;i++)
    {
        ex_girl[i]=love[i][0];
        for(int j=1;j<n;j++)
            ex_girl[i]=max(ex_girl[i],love[i][j]);
    }
 
    for(int i=0;i<n;i++)
    {
        fill(slack,slack+n,inf);
        while(1){
            memset(vis_boy,0,sizeof(vis_boy));
            memset(vis_girl,0,sizeof(vis_girl));
            if(dfs(i)) break;
            int d=inf;
            for(int j=0;j<n;j++)
                if(!vis_boy[j]) d=min(d,slack[j]);
            for(int j=0;j<n;j++){
                if(vis_girl[j]) ex_girl[j]-=d;
                if(vis_boy[j]) ex_boy[j]+=d;
                else slack[j]-=d;
            }
        }
    }
    int res=0;
    for(int i=0;i<n;i++)
        res+=love[match[i]][i];
    printf("%d\n",res);
}
 
int main()
{
    while(scanf("%d",&n)!=EOF)
    {
        for(int i=0;i<n;i++)
            for(int j=0;j<n;j++)
                scanf("%d",&love[i][j]);
        km();
    }
    return 0;
}

 

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值