Hopcroft-Karp算法模板

http://blog.csdn.net/u011742541/article/details/14104573

//Hopcroft-Karp算法 有点类似dinic 都是先对图BFS分层再沿层数DFS找增广路    
//***************************************************************************    
bool searchPath()        //BFS 对二分图分层    
{  
    dist = inf;  
    queue<int>que;  
    memset( dx,-1,sizeof(dx) );  
    memset( dy,-1,sizeof(dy) );  
    for( int i = 1; i <= nx; i ++ ){  //找到x集合所有未被匹配的点压入队列中  
        if( cx[i] == -1 ){  
            que.push(i);  
            dx[i] = 0;  
        }  
    }  
    while( !que.empty() ){  
        int u = que.front(); que.pop();  
        if( dx[u] > dist )    //只分层到第一个找个的可匹配点层数  
            break;  
        for( int v = 1; v <= ny; v ++ ) {  
            if( map[u][v] && dy[v] == -1 ){  
                dy[v] = dx[u] + 1;  
                if( cy[v] == -1 )      //找个一个可以匹配点 标记分层的层数   
                    dist = dy[v];  
                else{  
                    dx[cy[v]] = dy[v] + 1;  
                    que.push( cy[v] );  
                }  
            }  
        }  
    }  
    return dist != inf;          
}  
bool findPath( int u )    //沿着层数DFS    
{  
    for( int v = 1; v <= ny; v ++ ){  
        if( map[u][v] && !vis[v] && dy[v] == dx[u] + 1 ){  
            vis[v] = true;  
            if( cy[v] != -1 && dy[v] == dist  )    //如果v已经有匹配了且v的层数为dist（ 最大层数为dist 所以v原来匹配的不可能再匹配 ）  
                continue;  
            if( cy[v] == -1 || findPath( cy[v] ) ){   //如果v未匹配就跟u匹配v 否则看v原来匹配的是否还能跟其他的匹配 能就跟u匹配 不能就不匹配   
                cy[v] = u;  
                cx[u] = v;  
                return true;  
            }  
        }  
    }  
    return false;  
}  
int HK_MaxMatch()  
{  
    int ans = 0;  
    memset( cx,-1,sizeof(cx) );  
    memset( cy,-1,sizeof(cy) );  
    while( searchPath() ){   //分层 + 判断是否还有未匹配点  
        memset( vis,0,sizeof(vis) );  
        for( int i = 1; i <= nx; i ++ ){    
            if( cx[i] == -1 )  
                ans += findPath(i);  
        }  
    }  
    return ans;  
}  
//***************************************************************************