题目大意:
给定一个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;
}