poj 3041 Asteroids 二分图最大匹配

题目

题目链接：http://poj.org/problem?id=3041

题目来源：《挑战》例题。

简要题意：给定矩阵，及上面一些东西，每次可以搞掉一行或一列的东西，问最少几次可以杀光。

题解

应该算是非常经典的二分图最大匹配的题目了，也能从中一窥图论的魅力。

将每行和每列作为点，每个输入的点行和列作为二分图的两边去连边。

可以发现要求的就是最小点覆盖，即求最小的边集令每条边至少一个端点在边集构成的点集中。

由于二分图的最小点覆盖等于最大匹配，所以可以去求匹配，需要两次转化。

二分图匹配的话可以匈牙利或者最大流去弄。

代码

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <cstring>
#include <stack>
#include <queue>
#include <string>
#include <vector>
#include <set>
#include <map>
#define fi first
#define se second
using namespace std;
typedef long long LL;
typedef pair<int,int> PII;
// head
const int N = 505;
const int M = 1e4+5;

struct Edge {
    int to, nxt;
    Edge(int to, int nxt) : to(to), nxt(nxt) {}
    Edge() {}
};
int head[N*2];
Edge e[M*2];
void addEdge(int from, int to, int cnt) {
    e[cnt] = Edge(to, head[from]);
    head[from] = cnt;
}
void init(int n) {
    for (int i = 0; i <= n; i++) head[i] = -1;
}

int match[N*2];
bool vis[N*2];

bool dfs(int x) {
    vis[x] = true;
    for (int i = head[x]; ~i; i = e[i].nxt) {
        int to = e[i].to, w = match[to];
        if (w == -1 || (!vis[w] && dfs(w))) {
            match[x] = to;
            match[to] = x;
            return true;
        }
    }
    return false;
}

int bipatiteMatch(int n) {
    memset(match, -1, sizeof match);
    int ans = 0;
    for (int i = 0; i < n; i++) {
        if (match[i] == -1) {
            memset(vis, 0, sizeof vis);
            if (dfs(i)) ans++;
        }
    }
    return ans;
}

int main() {
    int n, k, r, c;
    while (scanf("%d%d", &n, &k) == 2) {
        init(2*n);
        int ec = 0;
        for (int i = 0; i < k; i++) {
            scanf("%d%d", &r, &c);
            r--, c--;
            addEdge(r, c+n, ec++);
            addEdge(c+n, r, ec++);
        }
        printf("%d\n", bipatiteMatch(2*n));
    }
    return 0;
}