HDU 4398 whosyourdaddy 精确覆盖,允许重复覆盖

本文介绍了一种解决领主攻击问题的算法实现,通过最小覆盖的方法计算领主最少需要攻击多少次才能确保所有点至少被攻击一次。该算法允许重复覆盖,并采用特殊的删除策略仅移除同列点。

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

题目大意:有n个点,其中一些点是相连的。领主的攻击具有溅射,即攻击一个点,此点相邻的点也会收到攻击。问,领主最少攻击多少次,使得每个点都至少被攻击一次。


与传统精确覆盖相比,此题允许重复覆盖。那么,我们删除的时候,只把同列的点删掉,而不是把同列的点所在的行删掉。

需要注意的是,点要算上重边以及表头,共有55*(4+1)+56个点。我最开始写的300个点,TLE无数次。。。


//#pragma comment(linker, "/STACK:102400000,102400000")
#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
#include<cmath>
#include<cctype>
#include<string>
#include<algorithm>
#include<iostream>
#include<ctime>
#include<map>
#include<set>
using namespace std;
#define MP(x,y) make_pair((x),(y))
#define PB(x) push_back(x)
typedef long long LL;
//typedef unsigned __int64 ULL;
/* ****************** */
const int INF=100011122;
const double INFF=1e100;
const double eps=1e-8;
//const LL mod=1000000007;
const int NN=100010;
const int MM=2000010;
/* ****************** */

const int maxn=60;
const int maxr=60;
const int maxnode=360;
int sz,ansd;
int S[maxn],H[maxr];
int row[maxnode],col[maxnode];
int L[maxnode],D[maxnode],U[maxnode],R[maxnode];
bool mat[60][60];
bool vis[60];

void dl_init(int n)
{
    memset(H,-1,sizeof(H));
    int i;
    for(i=0;i<=n;i++)
    {
        U[i]=D[i]=i;
        L[i]=i-1;
        R[i]=i+1;
        S[i]=0;
    }
    L[0]=n;
    R[n]=0;
    sz=n+1;
}

void dl_add(int r,int c)
{
    row[sz]=r,col[sz]=c;
    S[c]++;

    U[sz]=U[c];
    D[sz]=c;

    D[ U[c] ]=sz;
    U[c]=sz;

    if(H[r]==-1)
    {
        H[r]=L[sz]=R[sz]=sz;
    }
    else
    {
        L[sz]=L[H[r]];
        R[sz]=H[r];
        R[ L[sz] ]=sz;
        L[ H[r] ]=sz;
    }
    sz++;
}

void remove(int id)
{
    int i;
    for(i=D[id];i!=id;i=D[i])
    {
        L[R[i]]=L[i];
        R[L[i]]=R[i];
    }
}
void restore(int id)
{
    int i;
    for(i=U[id];i!=id;i=U[i])
    {
        R[L[i]]=i;
        L[R[i]]=i;
    }
}

int h()
{
    int i,j,k,ans=0;
    memset(vis,false,sizeof(vis));
    for(i=R[0];i!=0;i=R[i])
    {
        if(!vis[i])
        {
            ans++;
            vis[i]=true;
            for(j=D[i];j!=i;j=D[j])
                for(k=R[j];k!=j;k=R[k])
                    vis[ col[k] ]=true;
        }
    }
    return ans;
}

void dl_dfs(int d)
{
    if(d+h()>=ansd)
        return;
    if(R[0]==0)
    {
        ansd=d;
        return;
    }

    int i,j,c=R[0];

    for(i=R[0];i!=0;i=R[i])
        if(S[i]<S[c])
            c=i;

    for(i=D[c];i!=c;i=D[i])
    {
        remove(i);
        for(j=R[i];j!=i;j=R[j])
            remove(j);
        dl_dfs(d+1);
        for(j=L[i];j!=i;j=L[j])
            restore(j);
        restore(i);
    }
}
int main()
{
    int n,m,u,v,i,j;
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        memset(mat,false,sizeof(mat));
        dl_init(n);
        for(i=0;i<m;i++)
        {
            scanf("%d%d",&u,&v);
            mat[u][v]=true;
            mat[v][u]=true;
        }
        for(i=1;i<=n;i++)
        {
            for(j=1;j<=n;j++)
            {
                if(mat[i][j] || i==j)
                {
                    dl_add(i,j);
                }
            }
        }
        ansd=n;
        dl_dfs(0);
        printf("%d\n",ansd);
    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值