[POJ2594]Treasure Exploration（最小路径覆盖变种，floyd算法，匈牙利算法）

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

题意：给一张图，有单向边。现在要往图上某几个点放几个机器人，机器人要把所有的点都走过，问最少需要放多少个机器人，路径可以重复走。

路径可以重复走，则可以用floyd求出所有可达的点，再做最小路径覆盖，根据公式 最小路径覆盖=节点数-最大匹配直接做二分图最大匹配就可以。

/*
━━━━━┒ギリギリ♂ eye！
┓┏┓┏┓┃キリキリ♂ mind！
┛┗┛┗┛┃＼○／
┓┏┓┏┓┃ /
┛┗┛┗┛┃ノ)
┓┏┓┏┓┃
┛┗┛┗┛┃
┓┏┓┏┓┃
┛┗┛┗┛┃
┓┏┓┏┓┃
┛┗┛┗┛┃
┓┏┓┏┓┃
┃┃┃┃┃┃
┻┻┻┻┻┻
*/
#include <algorithm>
#include <iostream>
#include <iomanip>
#include <cstring>
#include <climits>
#include <complex>
#include <fstream>
#include <cassert>
#include <cstdio>
#include <bitset>
#include <vector>
#include <deque>
#include <queue>
#include <stack>
#include <ctime>
#include <set>
#include <map>
#include <cmath>
using namespace std;
#define fr first
#define sc second
#define cl clear
#define BUG puts("here!!!")
#define W(a) while(a--)
#define pb(a) push_back(a)
#define Rint(a) scanf("%d", &a)
#define Rs(a) scanf("%s", a)
#define Cin(a) cin >> a
#define FRead() freopen("in", "r", stdin)
#define FWrite() freopen("out", "w", stdout)
#define Rep(i, len) for(int i = 0; i < (len); i++)
#define For(i, a, len) for(int i = (a); i < (len); i++)
#define Cls(a) memset((a), 0, sizeof(a))
#define Clr(a, x) memset((a), (x), sizeof(a))
#define Full(a) memset((a), 0x7f7f7f, sizeof(a))
#define lrt rt << 1
#define rrt rt << 1 | 1
#define pi 3.14159265359
#define RT return
#define lowbit(x) x & (-x)
#define onecnt(x) __builtin_popcount(x)
typedef long long LL;
typedef long double LD;
typedef unsigned long long ULL;
typedef pair<int, int> pii;
typedef pair<string, int> psi;
typedef pair<LL, LL> pll;
typedef map<string, int> msi;
typedef vector<int> vi;
typedef vector<LL> vl;
typedef vector<vl> vvl;
typedef vector<bool> vb;

const int maxn = 555;
int nu, nv;
int G[maxn][maxn];
int linker[maxn];
bool vis[maxn];

bool dfs(int u) {
For(v, 1, nv+1) {
if(G[u][v] && !vis[v]) {
vis[v] = 1;
if(linker[v] == -1 || dfs(linker[v])) {
linker[v] = u;
return 1;
}
}
}
return 0;
}

int hungary() {
int ret = 0;
Clr(linker, -1);
For(u, 1, nu+1) {
Cls(vis);
if(dfs(u)) ret++;
}
return ret;
}

int n, m;

int main() {
// FRead();
int u, v;
while(~Rint(n) && ~Rint(m) && n + m) {
Cls(G);
nu = nv = n;
Rep(i, m) {
Rint(u); Rint(v);
G[u][v] = 1;
}
For(k, 1, n+1) {
For(i, 1, n+1) {
For(j, 1, n+1) {
if(G[i][k] + G[k][j] == 2) {
G[i][j] = 1;
}
}
}
}
printf("%d\n", nu - hungary());
}
RT 0;
}