HUST 1017 Dancing Links 模板

#include<cstdio>
#define N 2011
#define M 102011
int U[M],D[M],L[M],R[M],C[M],X[M];//x代表M所属的行，c代表M所属的列，U，D，L，R为一个元素的上下左右指针
int H[N],S[N],Q[N];//H是水平循环链表的头指针，S代表每一列的元素个数，Q存储一个可能结果

int size,n,m;
void remove(int c)
{
    R[L[c]]=R[c],L[R[c]]=L[c];
    for(int i=D[c]; i!=c; i=D[i])
        for(int j=R[i]; j!=i; j=R[j])
            U[D[j]]=U[j],D[U[j]]=D[j],--S[C[j]];
}
void resume(int c)
{
    R[L[c]]=L[R[c]]=c;
    for(int i=U[c]; i!=c; i=U[i])
        for(int j=L[i]; j!=i; j=L[j])
            ++S[C[U[D[j]]=D[U[j]]=j]];
}
bool Dance(int k)
{
    int i,j,tmp,c;
    if(!R[0])
    {
        printf("%d",k);
        for(i=0; i<k; ++i)printf(" %d",X[Q[i]]);
        puts("");
        return 1;
    }
    for(tmp=N,i=R[0]; i; i=R[i])
        if(S[i]<tmp)tmp=S[c=i];
    remove(c);
    for(i=D[c]; i!=c; i=D[i])
    {
        Q[k]=i;
        for(j=R[i]; j!=i; j=R[j])remove(C[j]);
        if(Dance(k+1))return 1;
        for(j=L[i]; j!=i; j=L[j])resume(C[j]);//在这里我遵循的原则是先删除的后还原，后删除的先还原。
    }
    resume(c);
    return 0;
}
void Link(int r,int c)
{
    ++S[C[size]=c];
    D[size]=D[c];
    U[D[c]]=size;
    U[size]=c;
    D[c]=size;
    if(H[r]<0)H[r]=L[size]=R[size]=size;
    else
    {
        R[size]=R[H[r]];
        L[R[H[r]]]=size;
        L[size]=H[r];
        R[H[r]]=size;
    }
    X[size++]=r;
}
int main()
{
    int i,j,num;
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        for(i=0; i<=m; ++i)
        {
            S[i]=0;
            D[i]=U[i]=i;
            L[i+1]=i;
            R[i]=i+1;
        }
        R[m]=0;
		L[0]=m;
        size=m+1;
        for(i=1; i<=n; ++i)
        {
            H[i]=-1;
            scanf("%d",&num);
            while(num--)scanf("%d",&j),Link(i,j);
        }
        if(!Dance(0))puts("NO");
    }
    return 0;
}

 

 在网上搜了好多资料，论文，终于参考别人代码看懂了dancing links。。。。。。。不过对于remove和resume顺序相反还是有点不理解。。。