hdoj 1151 Air Raid

又一个二分图最大匹配问题。匈牙利算法总是很好写的,就是有时候求的不是匹配数,比如最大独立集等。

做了七八个图论的题,其他都交了一次,只有这个交了两次才过,原因是没注意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;
}


你可能感兴趣的:(hdoj 1151 Air Raid)