Dancing Links_模版

这个模版是我从SuDoKu_DLX上面分离出来的，我也是刚接触DLX，总的来说只要能转换为精确覆盖问题的问题，都可以用DLX来解决，而且速度绝对是数一数二的。。。

const int SIZE = 16;
const int MAX_C = SIZE * SIZE * 4;                    //最大列 
const int MAX_R = SIZE * SIZE * SIZE;                 //最大行
const int MAX_LINK = MAX_C * MAX_R;                   //链表最大范围

int L[MAX_LINK],R[MAX_LINK],U[MAX_LINK],D[MAX_LINK];  //抽象链表
int C[MAX_LINK],O[MAX_LINK],S[MAX_C],H[MAX_R];        //C&O代表列&行，S每一列的节点数，H每一行的第一个节点
int NodeNumber;                                       //用来指向节点 
//////////////////////Dancing Links模版//////////////////////
void init(void);        //Dancing Links的抽象链表初始化 
void insert(int,int);   //在链表的一个位置中添加标记 
void remove(int);       //删除一列，同时删除这一列中的行 
void resume(int);       //恢复一列，同时恢复这一列中的行
//////////////////////Dancing Links模版//////////////////////

void init(void)
{
   for(int i=0;i<=MAX_C;i++)
      {
         L[i]=i-1;
         R[i]=i+1;
         U[i]=i;
         D[i]=i;
         C[i]=i;
         O[i]=0;   
      }
   L[0]=MAX_C;
   R[MAX_C]=0;
   NodeNumber=MAX_C+1;
   memset(S,0,sizeof(S));
   memset(H,0,sizeof(H));
}

void insert(int i,int j)   
{
   if(H[i])
      {
         L[NodeNumber]=L[H[i]];
         R[NodeNumber]=H[i];
         L[R[NodeNumber]]=NodeNumber;
         R[L[NodeNumber]]=NodeNumber;
      }
   else
      {
         L[NodeNumber]=NodeNumber;
         R[NodeNumber]=NodeNumber;
         H[i]=NodeNumber;
      }
   U[NodeNumber]=U[j];
   D[NodeNumber]=j;
   U[D[NodeNumber]]=NodeNumber;
   D[U[NodeNumber]]=NodeNumber;
   C[NodeNumber]=j;
   O[NodeNumber]=i;
   S[j]++;
   NodeNumber++;
}

void remove(int c)
{
   L[R[c]]=L[c];
   R[L[c]]=R[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)
{
   for(int i=U[c];i!=c;i=U[i])
      {
         for(int j=L[i];j!=i;j=L[j])
            {
               U[D[j]]=j;
               D[U[j]]=j;
               S[C[j]]++;
            } 
      }
   L[R[c]]=c;
   R[L[c]]=c;   
}