匈牙利算法之hdu1150

//题目大意：有AB两种机器，分别有很多种工作状态.如果从一种状态换到另外一种状态的话.

//需要重启机器..现在给一个关系...表示任务i可以用A中的一个机器或者是B中的一个机器来完成.

//求最小的重启次数..完成所有任务..

//题目连接:http://acm.hdu.edu.cn/showproblem.php?pid=1150

//思路：由于每个任务都可以用A或B机器的某种工作状态来完成，

//因此，对于每个任务，可以把A，B的工作状态作为顶点，
//在对应的Ａ和Ｂ的工作状态间连线，则对于每条边的两个顶点，
//因为每个任务都要完成，故至少要有一个顶点（某个机器的一种状态）存在，
//这样这个问题就转化成了在这些顶点中选取尽量少的顶点（尽量少的工作状态，尽量少的重启次数），
//使得每一条边都至少和其中一个点相关联，即最小点覆盖问题。
//由于二分图的最小点覆盖数=最大匹配数，故可用匈牙利算法解决。
#include <iostream>
#include <math.h>
#include <string.h>
using namespace std;
#define MAX 110
int map[MAX][MAX];
int m,n,k,z,x,y;
int mk[MAX];
//从X集合中的顶点u出发用深度优先的策略寻找增广路
//(这种增广路只能使当前的匹配数增加1)
int nx,ny; //X和Y集合中顶点的个数
int cx[MAX],cy[MAX];
//cx[i]表示最终求得的最大匹配中与Xi匹配的Y顶点, cy[i]同理
int path(int u)
{
    for(int v=1; v<=m; v++) //考虑所有Yi顶点v
    {
        if(map[u][v]&&!mk[v])
        {
            mk[v]=1;
            if(cy[v]==-1|| path(cy[v]))//如果v没有匹配,
            {//或者如果v已经匹配了，但从y[v]出发可以找到一条增广路
                cx[u] = v; //把v匹配给u
                cy[v] = u; //把u匹配给v
                return 1; //找到可增广路
            }
        }
    }
    return 0 ; //如果不存在从u出发的增广路
}
int MaxMatch() //求二部图最大匹配的匈牙利算法
{
    int res=0;
    memset(cx,0xff,sizeof(cx)); //从0匹配开始增广
    memset(cy,0xff,sizeof(cy));
    for(int i=1; i<=n; i++)
    {
        if(cx[i]==-1) //从每个未盖点出发进行寻找增广路
        {
            memset(mk,0,sizeof(mk));
            res+=path(i); //每找到一条增广路，可使得匹配数加1
        }
    }
    return res;
}
int main()
{
    int i;
    while(cin>>n>>m>>k)
    {
        if(n==0)break;
        nx=n;
        ny=m;
        memset(map,0,sizeof(map));
        for(i=0; i<k; i++)
        {
            cin>>z>>x>>y;
            map[x][y]=1;
        }
        int max=MaxMatch();
        cout<<max<<endl;
    }
    return 0;