Dancing links

Dancing links是一种能高效实现Knuth的X算法的技术，它可以使很多搜索问题得到极大的优化。

假设x是一个双向链表中的一个节点，L[x]表示X的前驱，R[x]表示x的后继，

则R[L[x]] = R[x], L[R[x]] = L[x]这一操作可以把x从链表中移除，这是众所周知的，

当然，一个细致的程序员还会用 L[x] = R[x] = x或 L[x] = R[x] = NULL这样的操作来清除x，

以免发生内存泄露，所以只有很少的程序员意识到，若不清除x则 R[L[x]] = x, L[R[x]] = x

这个操作可以把刚才移除的x恢复到之前的链表中，而这样的恢复功能在有些时候是很有用的，

例如，在一个交互性的程序中，用户很可能想撤销他所做的一个或一系列操作，恢复到先前的状态。

#define head 100
int U[N],D[N],L[N],R[N];
int C[N],H[N],ans[N];//C[N]表示N的列标，H[N]表示N的行标，ans[]用来储存结果
bool dance(int k)
{
     int c = R[head];
     if(c==head) {
           Output();//Output the solution
           return true;
     }
     remove(c);     
     for(int i=D[c]; i!=c; i=D[i])
     {
           ans[k] = H[i];
           for(int j=R[i]; j!=i; j=R[j]) remove(C[j]);
           if(dance(k+1)) return true;
           for(int j=L[i]; j!=i; j=L[j]) resume(C[j]);
     }     
     resume(c);
     return false;
}

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

void resume(int c)
{
     L[R[c]] = c;
     R[L[c]] = c;
     for(int i=U[c]; i!=c; i=U[i]) {
           for(int j=R[i]; j!=i; j=R[j]) {
                 U[D[j]] = j;
                 D[U[j]] = j;
           }
     }
}