HDU 1426 dancing links解决数独问题

题目大意：

这是一个最简单的数独填充题目，题目保证只能产生一种数独，所以这里的初始9宫格较为稠密，可以直接dfs也没有问题

但最近练习dancing links，这类数据结构解决数独无疑效率会高很多

dancing links的数独限制条件是：

1.每行有9个元素，共9行 对应dlx81列

2.每列有9个元素，共9行 对应dlx81列

3.每个九宫格有9个元素，共9行 对应dlx81列

4.81个格子，每个格子最多有一个数

 #include <iostream>
 #include <cstring>
 #include <cstdio>
 #include <algorithm>
 #include <queue>
 #include <climits>
 #include <cmath>

 using namespace std;
 #define N 1000
 #define MAXNODE 1000000
 const int INF = INT_MAX;
 const double eps = 1e-;

 int a[][];
 char str[];

 int belong[][] =  {
                         {,,,,,,,,},
                         {,,,,,,,,},
                         {,,,,,,,,},
                         {,,,,,,,,},
                         {,,,,,,,,},
                         {,,,,,,,,},
                         {,,,,,,,,},
                         {,,,,,,,,},
                         {,,,,,,,,},
                     };
 struct DLX{
     int n ,m , size;
     int col[MAXNODE] , row[MAXNODE];
     int U[MAXNODE] , D[MAXNODE] , L[MAXNODE] , R[MAXNODE];
     int cnt_col[N] , first[N];
     int ans[];

     void init(int _n , int _m)
     {
         n = _n , m = _m;
         size= m ;
         for(int i= ; i<=m ; i++){
             L[i] = i- , R[i] = i+;
             U[i] = D[i] = i;
         }
         L[] = m , R[m] = ;
         for(int i= ; i<=m ; i++) cnt_col[i] = ;
         for(int i= ; i<=n ; i++) first[i] = -;
     }

     void link(int r , int c)
     {
         ++size;
         U[D[c]] = size , D[size] = D[c];
         U[size] = c , D[c] = size;

         if(first[r]<) L[size]=R[size]=first[r] = size;
         else{
             L[R[first[r]]] = size , R[size] = R[first[r]];
             L[size] = first[r] , R[first[r]] = size;
         }
         row[size] = r , col[size] = c , cnt_col[c]++;
     }

     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];
                 cnt_col[col[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]] = D[U[j]] = j;
                 cnt_col[col[j]]++;
             }
         }
         L[R[c]] = R[L[c]] = c;
     }

     bool Dance(int d)
     {
         if(!R[]){
             for(int i= ; i<d ; i++){
                     int r = (ans[i]-)/;
                     int c = ((ans[i]-)%)/;
                     a[r][c] = ((ans[i]-)%)+;
             }
             return true;
         }
         int st=R[];
         for(int i=R[] ; i!= ; i=R[i])
             if(cnt_col[i]<cnt_col[st])
                 st = i;
         Remove(st);
         for(int i=D[st] ; i!=st ; i=D[i]){
             ans[d] = row[i];
             for(int j=R[i] ; j!=i ; j=R[j]) Remove(col[j]);
             if(Dance(d+)) return true;
             for(int j=L[i] ; j!=i ; j=L[j]) Resume(col[j]);
         }
         Resume(st);
         return false;
     }

 }dlx;

 void printM()
 {
     for(int i= ; i< ; i++)
         for(int j= ; j< ; j++){
            if(j<) printf("%d " , a[i][j]);
            else printf("%d\n" , a[i][j]);
         }
 }

 int main()
 {
   //  freopen("a.in" , "r" , stdin);
     int index=;
     bool flag = false;
     while(~scanf("%s" , str))
     {
         if(flag) puts("");
         flag = true;
         if(str[] == '?') a[index/][index%] = ;
         else a[index/][index%] = str[]-'';
         index++;
         while(index<){
             scanf("%s" , str);
             if(str[] == '?') a[index/][index%] = ;
             else a[index/][index%] = str[]-'';
             index++;
         }
        // printM();
         dlx.init( , );
         for(int i= ; i< ; i++)
             for(int j= ; j< ; j++)
             {
                 if(a[i][j]){
                     dlx.link((i*+j)*+a[i][j] , i*+a[i][j]);
                     dlx.link((i*+j)*+a[i][j] , +j*+a[i][j]);
                     dlx.link((i*+j)*+a[i][j] , +(belong[i][j]-)*+a[i][j]);
                     dlx.link((i*+j)*+a[i][j] , +i*+j+);
                 }
                 else{
                     for(int k= ; k<= ; k++){
                         dlx.link((i*+j)*+k , i*+k);
                         dlx.link((i*+j)*+k , +j*+k);
                         dlx.link((i*+j)*+k , +(belong[i][j]-)*+k);
                         dlx.link((i*+j)*+k , +i*+j+);
                     }
                 }
             }
         dlx.Dance();

         printM();
         index=;
     }

     return ;
 }