题目描述
给定一个矩阵,这个矩阵上有一些位置可以放置装置。
装置会日字形互相攻击,求最多可放置多少个装置。
解题思路
因为日字形相互攻击的两个点肯定满足一个坐标和是奇数,另一个是偶数。
将矩阵看成零一奇偶矩阵,那么互相攻击的装置可以形成一个二分图。
答案显然是总数-最小割。
怎么思考?最小割不存在同时选择两个装置且互相攻击(因为不满足最大流)。
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=40005,maxm=800050,INF=2147483647;
const int flg[8][2]={{-1,-2},{-2,-1},{1,-2},{2,-1},{-1,2},{-2,1},{1,2},{2,1}};
int tot=-1,lnk[maxn],son[maxm],cap[maxm],flw[maxm],nxt[maxm];
int n,B[maxn],que[maxn],dis[maxn],hed,til,s,t,ans;
bool vis[maxn];
char a[205][205];
inline char _read(){
char ch=getchar();
while(ch!='0'&&ch!='1') ch=getchar();
return ch;
}
void add(int x,int y,int c){
nxt[++tot]=lnk[x];lnk[x]=tot;son[tot]=y;cap[tot]=c;flw[tot]=0;
nxt[++tot]=lnk[y];lnk[y]=tot;son[tot]=x;cap[tot]=0;flw[tot]=0;
}
int min(int x,int y){if (x<y) return x;return y;}
bool BFS(){
memset(vis,0,sizeof(vis));
hed=0,til=1;que[1]=s;vis[s]=1;
while(hed!=til){
int x=que[++hed];
for (int j=lnk[x];j!=-1;j=nxt[j])
if (!vis[son[j]]&&cap[j]>flw[j]){
que[++til]=son[j];vis[son[j]]=1;
dis[son[j]]=dis[x]+1;
}
}
return vis[t];
}
int DFS(int x,int mi){
if (x==t||mi==0) return mi;
int num=0;
for (int &j=B[x];j!=-1;j=nxt[j])
if (dis[x]+1==dis[son[j]]){
int now=DFS(son[j],min(mi,cap[j]-flw[j]));
if (now){
flw[j]+=now;flw[j^1]-=now;
num+=now;mi-=now;
if (!mi) break;
}
}
return num;
}
int Dinic(){
int num=0;
while(BFS()){
memcpy(B,lnk,sizeof(B));
num+=DFS(s,INF);
}
return num;
}
int id(int x,int y){return (x-1)*n+y;}
bool check(int x,int y){
if (x<1||y<1||x>n||y>n) return 0;
if (a[x][y]=='1') return 0;
return 1;
}
int main(){
freopen("exam.in","r",stdin);
freopen("exam.out","w",stdout);
memset(lnk,255,sizeof(lnk));
scanf("%d",&n);s=n*n+2;t=n*n+1;
for (int i=1;i<=n;i++)
for (int j=1;j<=n;j++) a[i][j]=_read();
for (int i=1;i<=n;i++)
for (int j=1;j<=n;j++)
if (a[i][j]=='0'){
if ((i+j)%2){
add(s,id(i,j),1);
for (int k=0;k<8;k++) if (check(i+flg[k][0],j+flg[k][1])) add(id(i,j),id(i+flg[k][0],j+flg[k][1]),1);
}else add(id(i,j),t,1);
ans++;
}
printf("%d\n",ans-Dinic());
return 0;
}