匈牙利算法 模版

poj 3692

//POJ3692
//匈牙利算法,二分图的最大匹配
//最大完全数：最大完全子图中顶点的个数   最大完全数=原图的补图的最大独立数

#include <iostream>
#include <cstdio>
#include <cstring>

using namespace std;

#define MAXN 210

int nx,ny;
bool vis[MAXN];
int map[MAXN][MAXN];
int link[MAXN];//第y个点与第x相连

bool dfs(int x)//寻找增广路
{
    for(int y=1;y<=ny;y++ )//遍历右边点集
    {
        if(!vis[y] && map[x][y])
        {
            vis[y]=1;
            if(link[y]==-1 || dfs(link[y]))//无配对 或 构成增广序列
            {
                link[y]=x;
                return true;
            }
        }
    }
    return false;
}

void maxmatch()
{
    int ans=0;
    memset(link,-1,sizeof(link));
    for(int x=1;x<=nx;x++)//遍历左边点集
    {
        memset(vis,0,sizeof(vis));
        if(dfs(x))
            ans++;
    }
    printf("%d\n",nx+ny-ans);
}

int main()
{
    int a,b,m;
    int cas=1;
    while(scanf("%d%d%d",&nx,&ny,&m) && (nx+ny+m))
    {
        //memset(map,0,sizeof(map));//原图
        memset(map,1,sizeof(map));//补图
        while(m--)
        {
            scanf("%d%d",&a,&b);
            //map[a][b]=1;
            map[a][b]=0;
        }
        printf("Case %d: ",cas++);
        maxmatch();
    }
    return 0;
}

poj 3041(dfs 模版)

//最坏时间复杂度是o(n^3)
#include<iostream>
#include<cstdio>
#include<cstring>

using namespace std;

#define MAXN 505

int N,M;
bool vis[MAXN];
int map[MAXN][MAXN];
int link[MAXN];
bool can(int t)
{
    for(int i=1;i<=N;i++)//对于集合N循环遍历
    {
        if(!vis[i] && map[t][i])//如果i与t有边&&i没被访问过
        {
            vis[i]=1;
            if(link[i]==-1 || can(link[i]))//如果i与点集合gn中的点没有边，没有参与匹配的边则找到一条增广路,rpos到i存在没有参与匹配的边，而i到link[i]存在参与匹配的边然后往后找直到找到增广路，然后用M‘来代替M增加一条边
            {
                link[i]=t;
                return true;
            }
        }
    }
    return false;
}

void maxmatch()
{
    int ans=0;
    memset(link,-1,sizeof(link));
    for(int i=1;i<=N;i++)//遍历一个点集
    {
        memset(vis,0,sizeof(vis));//每次都要将上一次走过的路径清空
        if(can(i))
            ans++;
    }
    printf("%d\n",ans);
}

int main()
{
    int a,b;
    while(scanf("%d%d",&N,&M)==2)
    {
        memset(map,0,sizeof(map));
        for(int i=0;i<M;i++)
        {
            scanf("%d%d",&a,&b);
            map[a][b]=1;
        }
        maxmatch();
    }
    return 0;
}