BZOJ3175 攻击装置

本文介绍了一道名为BZOJ3175的算法题——攻击装置的解题思路及代码实现。通过奇偶性构造二分图并运用最大匹配算法求解最大独立集问题。

BZOJ3175 攻击装置

题目传送门

题解

emm...似乎说这个是最大独立集,可能之前并不知道有这个东西。由于攻击装置是攻击日字形的,所以我们可以根据奇偶性构造二分图,跑出二分图的最大匹配之后,用总的点数减去匹配个数就是答案了(因为这所有的匹配中,每一对只能够选择一个)。

code

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
bool Finish_read;
template<class T>inline void read(T &x){Finish_read=0;x=0;int f=1;char ch=getchar();while(!isdigit(ch)){if(ch=='-')f=-1;if(ch==EOF)return;ch=getchar();}while(isdigit(ch))x=x*10+ch-'0',ch=getchar();x*=f;Finish_read=1;}
template<class T>inline void print(T x){if(x/10!=0)print(x/10);putchar(x%10+'0');}
template<class T>inline void writeln(T x){if(x<0)putchar('-');x=abs(x);print(x);putchar('\n');}
template<class T>inline void write(T x){if(x<0)putchar('-');x=abs(x);print(x);}
/*================Header Template==============*/
const int maxn=5e5+500;;
int atkx[8]={-1,-2,1,2,-1,-2,1,2};
int atky[8]={-2,-1,-2,-1,2,1,2,1};
int n,m,tot,c,tim;
struct edge {
    int to,nxt;
}E[maxn<<2];
int head[maxn];
char s[maxn];
int num[300][300],belong[maxn],vis[maxn],st[maxn];
/*==================Define Area================*/
void addedge(int u,int v) {
    E[++tot].to=v;E[tot].nxt=head[u];head[u]=tot;
}

bool Judge(int x,int y) {
    return (x>=1&&x<=n&&y>=1&&y<=n&&num[x][y]!=-1);
}

bool Solve(int u) {
    for(int i=head[u];~i;i=E[i].nxt) {
        int to=E[i].to;
        if(vis[to]==tim) continue;
        vis[to]=tim;
        if(belong[to]==-1||Solve(belong[to])) return belong[to]=u;
    }
    return 0;
}

int main() {
    read(n);
    memset(head,-1,sizeof head);
    int cnt=0;
    for(int i=1;i<=n;i++) {
        scanf("%s",s);
        for(int j=0;j<n;j++) {
            if(s[j]=='1') num[i][j+1]=-1;
            else num[i][j+1]=++cnt,belong[cnt]=-1;
        }
    }
    for(int i=1;i<=n;i++) {
        for(int j=1;j<=n;j++) {
            if(num[i][j]!=-1&&!((i+j)&1)) {
                for(int k=0;k<8;k++) {
                    int nx=i+atkx[k],ny=j+atky[k];
                    if(Judge(nx,ny)) {
                        addedge(num[i][j],num[nx][ny]);
                    }
                }
                st[++c]=num[i][j];
            }
        }
    }
    int ans=0;
    tim=0;
    for(int i=1;i<=c;i++) {
        ++tim;
        ans+=Solve(st[i]);
    }
    printf("%d\n",cnt-ans);
    return 0;
}

转载于:https://www.cnblogs.com/Apocrypha/p/9432755.html

评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符  | 博主筛选后可见
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值