poj3041——Asteroids

题目大意：在N*N的格子中，小行星占有K个格子，发射一道光束可以消灭一行或一列的所有小行星，问最少发射多少道光束可以将小行星全部消灭

输入：N  K(1 <= N <= 500) (1 <= K <= 10,000)

           第i个小行星的行数R  列数C（共K行）

输出：最少光束个数

分析：每行每列作为结点，每个小行星作为边。

           例如对于输入：（1，1）、（1，3）、（2，2）、（3，2） 四颗小行星构造图G：

          我们要消灭所有小行星，就是消灭图中所有边。由于要射击最少的行列数，所以问题转化为在图中选择最少的点，使得它们邻接到所有边，也就是最小点覆盖数问题，根据定理：最小点覆盖数=二分图最大匹配边数（选择最多的边，这些边的端点都不能重复选择，比如上图最多选择边1-1、2-2，再选的话端点就重复匹配了）。所以我们用匈牙利算法解决这道题，它是一种用增广路求最大匹配的算法。

           匈牙利算法详解参照https://www.cnblogs.com/shenben/p/5573788.html

代码：

#include <iostream>  
#include <vector>
  
using namespace std;  
  
int n, k;
int v1, v2; //二分图顶点集，都等于n  
vector<vector<bool> > map;
vector<int> link; //记录v2中的点y在v1中所匹配的点x的编号  
  
bool dfs(int x, vector<bool>& vis){
    for (int y = 1; y <= v2; y++){
        if (map[x][y] && !vis[y]){
            vis[y] = true;
            if (link[y] == 0 || dfs(link[y], vis)){
                link[y] = x;  
                return true;  
            }  
        }  
    }  
    return false;  
}  
  
int search(){ //匈牙利算法hungary algorithm 
    int result = 0; //最大匹配数
    for (int x = 1; x <= v1; x++){
        vector<bool> vis(n + 1);
        if (dfs(x, vis)) //从v1中的节点x开始寻找增广路径p  
            result++;  
    }  
    return result;
}  
  
int main(){
    cin >> n >> k;  
    v1 = v2 = n;  
    int x, y;  
    map = vector<vector<bool> >(n + 1,vector<bool>(n + 1)); 
    link = vector<int>(n + 1);
    for (int i = 1; i <= k; i++){
        cin >> x >> y;  
        map[x][y] = true;
    }  
    cout << search() << endl;  
    return 0;  
}