HDU 3498 舞蹈链精确覆盖

题意:n个敌人,和m个敌人之间的关系,u与v是朋友,当攻击敌人I时,他的朋友也可以被打死,问最少攻击几次可以将所有的敌人消灭

思路:比较裸的一道舞蹈链题目,今天也看了看算法的思想,比较好理解,但是自己写代码的话还是不可能的(毕竟太弱~~~),然后在网上找了几个模版,调了一个自己喜欢的照着谢了这个模版题,决定今天就刷舞蹈链了,舞蹈链精确覆盖推荐这篇博客,大神写的真是详细,懂了思想找个模版这题就可以过了

#include <vector>  
#include <stdio.h>  
#include <string.h>  
#include <stdlib.h>  
#include <iostream>  
#include <algorithm>  
using namespace std;  
typedef long long ll;  
typedef unsigned long long ull;  
const int inf=0x3f3f3f3f;  
const ll INF=0x3f3f3f3f3f3f3f3fll;  
const int maxn=60;  
int L[maxn*maxn],R[maxn*maxn],U[maxn*maxn],D[maxn*maxn];//节点的上下左右四个方向的链表  
int C[maxn*maxn],H[maxn],cnt[maxn],vis[maxn];//C列H行cnt列链表中元素个数  
int n,m,id,fans;  
vector<int>G[maxn];  
void init(){  
    for(int i=0;i<=n;i++){  
        G[i].clear();  
        G[i].push_back(i);  
        cnt[i]=0;U[i]=D[i]=i;  
        L[i+1]=i;R[i]=i+1;  
    }  
    R[n]=0;id=n+1;  
    memset(H,-1,sizeof(H));  
}  
void Link(int r,int c){  
    cnt[c]++;C[id]=c;  
    U[id]=U[c];D[U[c]]=id;  
    D[id]=c;U[c]=id;  
    if(H[r]==-1) H[r]=L[id]=R[id]=id;  
    else{  
        L[id]=L[H[r]];R[L[H[r]]]=id;  
        R[id]=H[r];L[H[r]]=id;  
    }  
    id++;  
}  
void Remove(int Size){  
    for(int j=D[Size];j!=Size;j=D[j])  
        L[R[j]]=L[j],R[L[j]]=R[j];  
}  
void Resume(int Size){  
    for(int j=D[Size];j!=Size;j=D[j])  
        L[R[j]]=R[L[j]]=j;  
}  
int h()  
{  
    int sum=0;  
    memset(vis,0,sizeof(vis));  
    for(int i=R[0];i;i=R[i]){  
        if(vis[i]) continue;  
        sum++;  
        for(int j=D[i];j!=i;j=D[j]){  
            for(int k=R[j];k!=j;k=R[k])  
                vis[C[k]]=1;  
        }  
    }  
    return sum;  
}  
void Dance(int k){  
    int mm=maxn,pos;  
    if(k+h()>=fans) return;  
    if(!R[0]){  
        if(k<fans) fans=k;  
        return;  
    }  
    for(int i=R[0];i;i=R[i])  
        if(mm>cnt[i]) mm=cnt[i],pos=i;  
    for(int i=D[pos];i!=pos;i=D[i])  
    {  
        Remove(i);  
        for(int j=R[i];j!=i;j=R[j]) Remove(j);  
        Dance(k+1);  
        for(int j=R[i];j!=i;j=R[j]) Resume(j);  
        Resume(i);  
    }  
}  
int main(){  
    int u,v;  
    while(scanf("%d%d",&n,&m)!=-1){  
        init();  
        for(int i=0;i<m;i++){  
            scanf("%d%d",&u,&v);  
            G[u].push_back(v);  
            G[v].push_back(u);  
        }  
        for(int i=1;i<=n;i++){  
            for(unsigned int j=0;j<G[i].size();j++){  
                int t=G[i][j];  
                Link(i,t);  
            }  
        }  
        fans=n;  
        Dance(0);  
        printf("%d\n",fans);  
    }  
    return 0;  
}  


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值