[最小割]BZOJ 3175—— [Tjoi2013]攻击装置

题目描述

给定一个矩阵，这个矩阵上有一些位置可以放置装置。

装置会日字形互相攻击，求最多可放置多少个装置。

解题思路

因为日字形相互攻击的两个点肯定满足一个坐标和是奇数，另一个是偶数。

将矩阵看成零一奇偶矩阵，那么互相攻击的装置可以形成一个二分图。

答案显然是总数-最小割。

怎么思考？最小割不存在同时选择两个装置且互相攻击（因为不满足最大流）。

#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;
}