BZOJ 1741: [Usaco2005 nov]Asteroids 穿越小行星群

Description

贝茜想驾驶她的飞船穿过危险的小行星群．小行星群是一个NxN的网格(1≤N≤500)，在网格内有K个小行星(1≤K≤10000)． 幸运地是贝茜有一个很强大的武器，一次可以消除所有在一行或一列中的小行星，这种武器很贵，所以她希望尽量地少用．给出所有的小行星的位置，算出贝茜最少需要多少次射击就能消除所有的小行星．

Input

第1行：两个整数N和K，用一个空格隔开．

第2行至K+1行：每一行有两个空格隔开的整数R，C(1≤R，C≤N)，分别表示小行星所在的行和列．

Output

一个整数表示贝茜需要的最少射击次数，可以消除所有的小行星

题解：

二分图，两列点分别表示行、列。

对于一个点(x,y) 由 x向y连边。

最小覆盖即是答案。

代码：

#include<cstdio>

#include<algorithm>

#include<cstring>

//by zrt

//problem:

using namespace std;

typedef long long ll;

const double eps=1e-9;

int H[505],P[10005],X[10005],tot;

int n,k;

inline void add(int x,int y){

    P[++tot]=y;X[tot]=H[x];H[x]=tot;

}

bool cover[505];

int link[505];

int ans;

bool find(int k){

    for(int i=H[k];i;i=X[i]){

        if(!cover[P[i]]){

            int q=link[P[i]];

            link[P[i]]=k;

            cover[P[i]]=1;

            if(q==-1){

                ans++;return true;

            }

            if(find(q)) return true;

            link[P[i]]=q;

        }

    }

    return false;

}

int main(){

    #ifdef LOCAL

    freopen("in.txt","r",stdin);

    freopen("out.txt","w",stdout);

    #endif

    scanf("%d%d",&n,&k);

    for(int i=0,x,y;i<k;i++){

        scanf("%d%d",&x,&y);

        add(x,y);

    }

    memset(link,-1,sizeof link);

    for(int i=1;i<=n;i++){

        memset(cover,0,sizeof cover);

        find(i);

    }

    printf("%d\n",ans);

    return 0;

}