舞蹈链 —— 重复覆盖 模板

D a n c i n g L i n k s ( 重 复 覆 盖 模 板 ) : Dancing Links(重复覆盖模板): DancingLinks
本模板取自 H D U HDU HDU 3498 3498 3498

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;

struct DancingLink{
	const static int N = 500010;
	const static int M = 1010;
    int n, s, ansd;	// 列数 节点总数 
    int S[M], A[M], H[M];	// S[]该列节点总数  A[]答案  H[]行首指针 
    int L[N], R[N], U[N], D[N];	// L[],R[],U[],D[] 上下左右 
    int X[N], C[N], vis[M];	// X[] C[] 行列编号 
    void init(int n){	// 初始化 
        this->n = n;
        for(int i=0; i<=n; i++)
            U[i]=i, D[i]=i, L[i]=i-1, R[i]=i+1;
        R[n]=0, L[0]=n; s=n+1;
        memset(S, 0, sizeof(S));
        memset(H, -1, sizeof(H));
    }
    void DelCol(int c){	// 删除列 
        for(int i=D[c]; i!=c; i=D[i])
            L[R[i]]=L[i], R[L[i]]=R[i];
    }
    void ResCol(int c){	// 恢复列 
        for(int i=U[c]; i!=c; i=U[i])
            L[R[i]]=R[L[i]]=i;
    }
    void AddNode(int r,int c){	// 添加节点 
        ++S[c], C[++s]=c, X[s]=r;
        D[s]=D[c], U[D[c]]=s, U[s]=c, D[c]=s;
        if(H[r]<0) H[r]=L[s]=R[s]=s;	// 行首节点
        else R[s]=R[H[r]], L[R[H[r]]]=s, L[s]=H[r], R[H[r]]=s;
    }
    int f(){
        int ret=0;
        memset(vis, 0, sizeof(vis));
        for(int i=R[0]; i; i=R[i])
            if(!vis[i]){
                ret++; vis[i]=1;
                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 ret;
    }
    void dfs(int d){	// 深度,深搜遍历 
        if(d+f() >= ansd) return;
        if(!R[0]){
            if(d<ansd) ansd=d;
            return;
        }
        int c=R[0];
        for(int i=R[0]; i; i=R[i]) if(S[i]<S[c]) c=i;
        for(int i=D[c]; i!=c; i=D[i]){
            DelCol(i); A[d]=X[i];
            for(int j=R[i]; j!=i; j=R[j]) DelCol(j);
            dfs(d+1);
            for(int j=L[i]; j!=i; j=L[j]) ResCol(j);
            ResCol(i);
        }
    }
} dlx; 

int main(){
    int n, m;
    while(~scanf("%d%d", &n, &m)) {
        dlx.init(n);
        for(int i=1; i<=m; i++) {
            int a, b;
            scanf("%d%d", &a, &b);	// 第 i行有 num个 1
            dlx.AddNode(a, b);
            dlx.AddNode(b, a);
        }
        for(int i=1; i<=n; i++)
            dlx.AddNode(i, i);
        dlx.ansd = 0x3f3f3f;
        dlx.dfs(0);
        printf("%d\n", dlx.ansd);
    }
}
  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值