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);
}
}