SuDoKu_DLX_9*9

比DFS快了4倍只用10MS就能得出结果。。

#include <cstdio>
#include <cstring>
#include <ctime>
const int XSIZE = 3; 
const int SIZE = XSIZE * XSIZE;
const int MAX_C = SIZE * SIZE * 4;                    //最大列 
const int MAX_R = SIZE * SIZE * SIZE;                 //最大行
const int MAX_SUDOKU = 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,RecordNumber;                          //用来指向节点
int state[MAX_SUDOKU],ans[MAX_SUDOKU],record[MAX_SUDOKU];
//////////////////////Dancing Links模版//////////////////////
void init(void);        //Dancing Links的抽象链表初始化 
void insert(int,int);   //在链表的一个位置中添加标记 
void remove(int);       //删除一列，同时删除这一列中的行 
void resume(int);       //恢复一列，同时恢复这一列中的行
//////////////////////Dancing Links模版//////////////////////
bool input(void)
{
   char buffer[SIZE+1][SIZE+1];//留一个空间 
   if(scanf("%s",buffer[0])==EOF)   return false;
   for(int i=1;i<SIZE;i++)
      scanf("%s",buffer[i]);
   memset(state,0,sizeof(state));
   for(int i=0;i<SIZE;i++)
      for(int j=0;j<SIZE;j++)
         {
            if(buffer[i][j]!='0')
               state[i*SIZE+j]=buffer[i][j]-'0';
         }
   return true;
}

void add(int i,int j,int k)
{
   int row=i*MAX_SUDOKU+j*SIZE+k;
   insert(row,i*SIZE+j+1);//注意 
   insert(row,i*SIZE+k+MAX_SUDOKU);
   insert(row,j*SIZE+MAX_SUDOKU*2+k);
   insert(row,(i/XSIZE*XSIZE +j/XSIZE)*SIZE+MAX_SUDOKU*3+k);
}

void build(void)
{
   int pos,row;
   for(int i=0;i<SIZE;i++)
      for(int j=0;j<SIZE;j++)
         {
            pos=i*SIZE+j;
            if(state[pos]!=0)
               {
                add(i,j,state[pos]);
               }
            else if(state[pos]==0)
               {
                  for(int k=1;k<=SIZE;k++)
                     {
                        add(i,j,k);
                     }   
               }
         }   
}

bool dfs(int k)
{
   if(!R[0])
      {
         RecordNumber=k;
         return true;   
      }
   int count=~(1<<31),c;
   for(int i=R[0];i;i=R[i])
      {
         if(S[i]<count)
            {
               count=S[i];
               c=i;
               if(count==1)
                  break;   
            }  
      }
   remove(c);
   for(int i=D[c];i!=c;i=D[i])
      {
         for(int j=R[i];j!=i;j=R[j])
            {
               remove(C[j]);   
            }
         record[k]=O[i];
         if(dfs(k+1))
            return true;
         for(int j=L[i];j!=i;j=L[j])
            {
               resume(C[j]);   
            }  
      }
   resume(c);
   return false;
}

void output(void)
{
   for(int i=0;i<RecordNumber;i++)
      {
         ans[(record[i]-1)/SIZE]=(record[i]-1)%SIZE+1;
      }
   for(int i=0;i<SIZE;i++)
      {
         for(int j=0;j<SIZE;j++)
            printf("%-2d",ans[i*SIZE+j]);
         printf("\n");  
      }
}

int main()
{
   while(input())
      {
         time_t start=clock();
         init();
         build();
         dfs(0);
         output();
         printf("Time:%ldMs\n",clock()-start);
      }
   return 0;   
}

//////////////////////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;   
}
/*
070680050
000000001
409507000
560000340
000000000
034000017
000703804
300000000
090015030
*/