又一个二分图最大匹配问题。匈牙利算法总是很好写的,就是有时候求的不是匹配数,比如最大独立集等。
做了七八个图论的题,其他都交了一次,只有这个交了两次才过,原因是没注意intersection的编号是从1开始,还按照0算了,于是贡献一次wa!
最近做了好几个这种题,直接贴代码:
/*
Hdoj: 1151 Air Raid
*/
#include <iostream>
#include <cstdio>
#include <cstring>
#define MaxV 125
using namespace std;
bool map[MaxV][MaxV];
bool checked[MaxV];
int xM[MaxV], yM[MaxV];
int interNo, streetNo;
bool search(int i)
{
for(int j = 0; j < interNo; j++) {
if(map[i][j] && !checked[j]) {
checked[j] = true;
if(yM[j] == -1 || search(yM[j])) {
xM[i] = j;
yM[j] = i;
return true;
}
}
}
return false;
}
int solve()
{
int res = 0;
memset(xM, -1, sizeof(xM));
memset(yM, -1, sizeof(yM));
for(int i = 0; i < interNo; i++) {
if(xM[i] == -1) {
memset(checked, false, sizeof(checked));
if(search(i))
res++;
}
}
return res;
}
int main()
{
//freopen("data.in", "rb", stdin);
int t;
scanf("%d", &t);
while(t--) {
memset(map, false, sizeof(map));
scanf("%d%d", &interNo, &streetNo);
for(int i = 0; i < streetNo; i++) {
int a, b;
scanf("%d%d", &a, &b);
map[a-1][b-1] = true;
}
int result = solve();
printf("%d\n", interNo - result);
}
return 0;
}