HDU 3360National Treasures(行列匹配+最大匹配)

题意：在一个n*m的图中每个点有宝贝要雇佣警卫来保护，-1是已经有的了，其他的整数，代表这个位置周围至少有的警卫布局，参考题目的图片。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<set>
#include<map>
#include<string>
#include<cstring>
#include<stack>
#include<queue>
#include<vector>
#include<cstdlib>
#define lson (rt<<1),L,M
#define rson (rt<<1|1),M+1,R
#define M ((L+R)>>1)
#define cl(a,b) memset(a,b,sizeof(a));
#define LL long long
#define P pair<int,int>
#define X first
#define Y second
#define pb push_back
#define fread (zcc) freopen(zcc,"r",stdin);
#define fwrite(zcc) freopen(zcc,"w",stdout);
using namespace std;
const int maxn=2505;
const int inf=999999;

vector<int> G[maxn];
int matching[maxn];
bool vis[maxn];
int num;
bool dfs(int u){
    int N=G[u].size();
    for(int i=0;i<N;i++){
        int v=G[u][i];
        if(vis[v])continue;
        vis[v]=true;
        if(matching[v]==-1||dfs(matching[v])){
            matching[v]=u;
            return true;
        }
    }
    return false;
}
int hungar(){
    int ans=0;
    cl(matching,-1);
    for(int i=0;i<num;i++){
        cl(vis,false);
        if(dfs(i))ans++;
    }
    return ans;
}
int  dir[][2]={-1,-2,-2,-1,-2,1,-1,2,1,2,2,1,2,-1,1,-2,-1,0,0,1,1,0,0,-1};
int a[maxn][maxn];
int main(){
    int n,m;
    int cas=1;
    while(~scanf("%d%d",&n,&m)&&(m||n)){
        for(int i=0;i<=n*m;i++)G[i].clear();
        num=0;
        for(int i=0;i<n;i++){
            for(int j=0;j<m;j++){
                scanf("%d",&a[i][j]);
            }
        }
        for(int i=0;i<n;i++){
            for(int j=0;j<m;j++){
                int x=a[i][j];
                if(x==-1)continue;
                for(int s=0;s<12;s++)if(x&(1<<s)){
                    int xx=dir[s][0]+i;
                    int yy=dir[s][1]+j;
                    if(xx>=0&&xx<n&&yy>=0&&yy<m&&a[xx][yy]!=-1){
                        G[i*m+j].pb(xx*m+yy);
                        G[xx*m+yy].pb(i*m+j);
                    }
                }
            }
        }
        num=n*m;
        printf("%d. %d\n",cas++,hungar()/2);

    }
    return 0;
}