二分图匹配

hdu 3360 National Treasures http://acm.hdu.edu.cn/showproblem.php?pid=3360　　　　2013-2-3

最少的保安数=n*m-最多的物品数

最多的物品数就是最大独立集 物品的奇偶不同才可能有边 是二分图

#include <cstdio>
#include <cstring>
#include <vector>
using namespace std;

const int N=2600;
struct BGMatch{
    int n;
    vector<int> G[N];
    int mat[N];
    bool vis[N];
    void init(int tn){
        n=tn;
        for(int i=0;i<n;i++) G[i].clear();
    }
    void addMatch(int u,int v){
        G[u].push_back(v);
    }
    bool find(int u){
        int sz=G[u].size();
        for(int i=0;i<sz;i++){
            int v=G[u][i];
            if(vis[v]) continue;
            vis[v]=true;
            if(mat[v]==-1 || find(mat[v])){
                mat[v]=u;
                return true;
            }
        }
        return false;
    }
    int maxMatch(){
        int cnt=0;
        memset(mat,-1,sizeof(mat));
        for(int u=0;u<n;u++){
            memset(vis,0,sizeof(vis));
            if(find(u)) cnt++;
        }
        return cnt;
    }
};

BGMatch g;
const int dx[]={-1,-2,-2,-1,1,2,2,1,-1,0,1,0};
const int dy[]={-2,-1,1,2,2,1,-1,-2,0,1,0,-1};
int mz[55][55];
int main(){
    int n,m,C=0;
    while(scanf("%d%d",&n,&m),n||m){
        g.init(n*m);
        for(int i=0;i<n;i++)
            for(int j=0;j<m;j++) scanf("%d",&mz[i][j]);
        for(int x=0;x<n;x++)
            for(int y=0;y<m;y++) if(mz[x][y]!=-1){
                int t=mz[x][y];
                for(int i=0;i<12;i++) if((t>>i)&1){
                    int nx=x+dx[i], ny=y+dy[i];
                    if(nx>=0 && nx<n && ny>=0 && ny<m && mz[nx][ny]!=-1){
                        g.addMatch(x*m+y,nx*m+ny);
                        g.addMatch(nx*m+ny,x*m+y);
                    }
                }
            }
        int ans=g.maxMatch()/2;
        printf("%d. %d\n",++C,ans);
    }
    return 0;
}