[二分图最大独立集] bzoj3175: [Tjoi2013]攻击装置

本文提供了一道名为[Tjoi2013]攻击装置的编程题解决方案,该题通过构建图论模型解决攻击装置布局问题。文章详细介绍了使用邻接表存储图数据结构的方法,并采用匈牙利算法求解最大匹配问题,最终计算出攻击装置的最佳布局方案。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

bzoj 3175: [Tjoi2013]攻击装置
同 bzoj4808 马 http://blog.youkuaiyun.com/qq_36038511/article/details/79601632

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
struct node
{
    int x,y,next;
}a[410000];
int last[41000],len;
int match[81000],v[81000],c[210][210];
int t;
void build(int x,int y)
{
    len++;
    a[len].x=x;a[len].y=y;a[len].next=last[x];last[x]=len;
}
bool find(int x)
{
    for (int k=last[x];k;k=a[k].next)
    {
        int y=a[k].y;
        if (v[y]!=t)
        {
            v[y]=t;
            if (match[y]==0||find(match[y])==1) 
            {
                match[y]=x;
                return 1;
            }
        }
    }
    return 0;
}
char ch[210];
int main()
{
    int n,sum=0;
    scanf("%d",&n);
    memset(c,-1,sizeof(c));
    for (int i=1;i<=n;i++)
    {
        scanf("%s",ch+1);
        for (int j=1;j<=n;j++) 
        {
            c[i][j]=ch[j]-'0';
            if (c[i][j]==0) sum++;
        }
    }
    for (int i=1;i<=n;i++)
    {
        for (int j=1;j<=n;j++)
        {
            if (c[i][j]==0)
            {
                int x=(i-1)*n+j;
                if (i-1>=1&&j-2>=1&&c[i-1][j-2]==0) build(x,(i-2)*n+j-2+sum);
                if (i+1<=n&&j-2>=1&&c[i+1][j-2]==0) build(x,i*n+j-2+sum);
                if (i-1>=1&&j+2<=n&&c[i-1][j+2]==0) build(x,(i-2)*n+j+2+sum);
                if (i+1<=n&&j+2<=n&&c[i+1][j+2]==0) build(x,i*n+j+2+sum);
                if (i-2>=1&&j-1>=1&&c[i-2][j-1]==0) build(x,(i-3)*n+j-1+sum);
                if (i-2>=1&&j+1<=n&&c[i-2][j+1]==0) build(x,(i-3)*n+j+1+sum);
                if (i+2<=n&&j-1>=1&&c[i+2][j-1]==0) build(x,(i+1)*n+j-1+sum);
                if (i+2<=n&&j+1<=n&&c[i+2][j+1]==0) build(x,(i+1)*n+j+1+sum);

            }
        }
    }
    t=0;
    int ans=0;
    for (int i=1;i<=n;i++)
    {
        for (int j=1;j<=n;j++)
        {
            if (c[i][j]==0)
            {
                t++;
                if (find((i-1)*n+j)) ans++;
            }
        }
    }
    printf("%d\n",sum-ans/2);
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值