[bzoj3175][Tjoi2013]攻击装置——最大独立集

题目大意：

给定一个01矩阵，其中你可以在0的位置放置攻击装置。每一个攻击装置（x,y）都可以按照“日”字攻击其周围的 8个位置（x-1,y-2）,(x-2,y-1),(x+1,y-2),(x+2,y-1),(x-1,y+2),(x-2,y+1), (x+1,y+2),(x+2,y+1)
求在装置互不攻击的情况下，最多可以放置多少个装置。

思路：

要求不能互相攻击， 那么我们把可以互相攻击到的点之间都连一条边，这样便形成了一个图，然后即对这个图求最大独立集。发现马如果走x步可以回到原地，所有的x都是偶数，即这个图中所有的环都是偶环，便可以二分图染色构造二分图。
在二分图中，最大独立集=点的总数-最小点覆盖=点的总数-最匹配数。

/*===============================
 * Author : ylsoi
 * Problem : bzoj3175
 * Algorithm : Graph Matching
 * Time : 2018.5.26
 * ============================*/
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<climits>
#include<queue>
using namespace std;
void File(){
    freopen("bzoj3175.in","r",stdin);
    freopen("bzoj3175.out","w",stdout);
}
#define REP(i,a,b) for(register int i=a;i<=b;++i)
#define DREP(i,a,b) for(register int i=a;i>=b;--i)
#define MREP(i,x) for(register int i=beg[x];i;i=E[i].last)
#define mem(a) memset(a,0,sizeof(a))
#define ll long long
#define inf INT_MAX
#define Num(x,y) ((x-1)*n+y)
const int maxn=200+10;
int n,cnt,beg[maxn*maxn],ans,tot,vis[maxn*maxn];
bool mp[maxn][maxn],c[maxn*maxn];
char s[maxn][maxn];
int dx[9]={0,1,1,2,2,-1,-1,-2,-2};
int dy[9]={0,2,-2,1,-1,2,-2,1,-1};
struct edge{
    int to;
    int last;
}E[maxn*maxn*16];
void add(int u,int v){
    ++cnt;
    E[cnt].to=v;
    E[cnt].last=beg[u];
    beg[u]=cnt;
}
bool pd(int x,int y){return x>=1 && x<=n && y>=1 && y<=n && (!mp[x][y]);}
void color(int u){
    vis[u]=1;
    MREP(i,u){
        int v=E[i].to;
        if(vis[v])continue;
        c[v]=c[u]^1;
        color(v);
    }
}
int b[maxn*maxn];
bool dfs(int u,int mark){
    MREP(i,u){
        int v=E[i].to;
        if(vis[v]==mark)continue;
        vis[v]=mark;
        if(!b[v] || dfs(b[v],mark)){
            b[v]=u;
            return true;
        }
    }
    return false;
}
int main(){
    File();
    scanf("%d",&n);
    REP(i,1,n)scanf("%s",s[i]+1);
    REP(i,1,n)REP(j,1,n)
        mp[i][j]= s[i][j]=='0' ? 0 : 1;
    REP(i,1,n)REP(j,1,n){
        if(mp[i][j])continue;
        REP(k,1,8){
            int xx=i+dx[k],yy=j+dy[k],num=Num(xx,yy);
            if(!pd(xx,yy))continue;
            add(Num(i,j),num);
        }
    }
    REP(i,1,n)REP(j,1,n){
        if(mp[i][j])continue;
        if(!vis[Num(i,j)])
            color(Num(i,j));
    }
    mem(vis);
    REP(i,1,n)REP(j,1,n){
        if(mp[i][j])continue;
        ++tot;
        if(!c[Num(i,j)])
            ans+=dfs(Num(i,j),Num(i,j));
    }
    printf("%d\n",tot-ans);
    return 0;
}