bzoj3175 [Tjoi2013]攻击装置

最新推荐文章于 2019-01-10 20:33:00 发布

原创最新推荐文章于 2019-01-10 20:33:00 发布 · 180 阅读

0 ·

CC 4.0 BY-SA版权

网络流专栏收录该内容

134 篇文章

订阅专栏

本文介绍了一种解决01矩阵中放置攻击装置问题的算法。该算法利用二分图匹配原理，在保证装置间不互相攻击的前提下，计算出最多能放置多少个装置。通过构建特定的图模型并运用图算法来求解。

http://www.elijahqi.win/archives/3129
Description

给定一个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)
求在装置互不攻击的情况下，最多可以放置多少个装置。

Input

第一行一个整数N，表示矩阵大小为N*N。接下来N行每一行一个长度N的01串，表示矩阵。
Output

一个整数，表示在装置互不攻击的情况下最多可以放置多少个装置。
Sample Input

3
010
000
100
Sample Output

4
HINT

100%数据 N<=200

Source

二分图黑白染色求原图最大独立集等于总点数-（最小点覆盖边集=最大匹配）

#include<queue>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=220;
const int inf=0x3f3f3f3f;
char s[N][N];
int dx[]={-2,-1,1,2,2,1,-1,-2},dy[]={1,2,2,1,-1,-2,-2,-1};
int cnt,n,id[N][N],num=1,sum;
int h[N*N],level[N*N],T,cur[N*N];
struct node{
    int y,next,z;
}data[N*N*16];
inline void insert1(int x,int y,int z){
    data[++num].y=y;data[num].next=h[x];h[x]=num;data[num].z=z;
    data[++num].y=x;data[num].next=h[y];h[y]=num;data[num].z=0;
}
inline bool bfs(){
    memset(level,0,sizeof(level));queue<int>q;q.push(0);level[0]=1;
    while(!q.empty()){
        int x=q.front();q.pop();
        for (int i=h[x];i;i=data[i].next){
            int y=data[i].y,z=data[i].z;
            if (level[y]||!z) continue;level[y]=level[x]+1;if (y==T) return 1;q.push(y);
        }
    }return 0;
}
inline int dfs(int x,int s){
    if (x==T) return s;int ss=s;
    for (int &i=cur[x];i;i=data[i].next){
        int y=data[i].y,z=data[i].z;
        if(level[x]+1==level[y]&&z){
            int xx=dfs(y,min(z,s));if (!xx) level[y]=0;
            s-=xx;data[i].z-=xx;data[i^1].z+=xx;if (!s) return ss;
        }
    }return ss-s;
}
int main(){
    freopen("bzoj3175.in","r",stdin);
    scanf("%d",&n);
    for (int i=1;i<=n;++i) scanf("%s",s[i]+1);
    for (int i=1;i<=n;++i) 
        for (int j=1;j<=n;++j) if(s[i][j]=='0') id[i][j]=++cnt;
    T=cnt+1;
    for (int i=1;i<=n;++i){
        for (int j=1;j<=n;++j){
            if(s[i][j]=='0') ++sum;else continue;
            if (i+j&1) {insert1(id[i][j],T,1);continue;}insert1(0,id[i][j],1);
            for (int k=0;k<8;++k){
                int x=i+dx[k],y=j+dy[k];
                if(x<1||x>n||y<1||y>n||s[x][y]=='1') continue;
                insert1(id[i][j],id[x][y],1);
            }
        }
    }while(bfs()) memcpy(cur,h,sizeof(h)),sum-=dfs(0,inf);
    printf("%d\n",sum);
    return 0;
}