POJ 3074 Sudoku （Dancing Links）

#include  
#include  
#include  
#include  
#include  
#include  
#include  
  
using namespace std;  
//   列：（行+列+块）*9种可能+9*9个格子  
//   行： 9*9*9  表示第i行第j列填k  
const int MAXN=(9+9+9)*9+9*9+9*9*9*9*9*4+10;  
#define INF 0xFFFFFF  
int size;  
int head,sz;  
int U[MAXN],D[MAXN],L[MAXN],R[MAXN];  
int H[MAXN],ROW[MAXN],C[MAXN],S[MAXN],O[MAXN];  
  
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])  
        {  
            ++S[C[j]];  
            U[D[j]]=j;  
            D[U[j]]=j;  
          }  
     }  
     L[R[c]]=c;  
     R[L[c]]=c;  
}  
  
bool dfs(int k)  
{  
    if(R[head]==head)  
    {  
        sort(O,O+9*9);  
        int p=0;  
        for(int i=0;i<9;i++)  
        {  
            for(int j=0;j<9;j++)  
            {  
                int num=O[p++];  
                //cout<