HDU 3111 Sudoku（精确覆盖）

数独问题，输入谜题，输出解

既然都把重复覆盖的给写成模板了，就顺便把精确覆盖的模板也写好看点吧。。。赤裸裸的精确覆盖啊~~~水一水~~~然后继续去搞有点难度的题了。。。

 #include <cstdio>
 #include <cstring>
 #include <iostream>
 #include <algorithm>
 #include <cmath>
 #include <string>
 #include <vector>
 #include <queue>
 #include <set>
 using namespace std;

 #define ll long long
 #define eps 1e-8
 #define mod 21092013

 #define inf 0x3f3f3f3f
 #define maxr 777
 #define maxn (maxr*maxr)
 int n,m;
 int L[maxn],R[maxn],U[maxn],D[maxn],cnt;
 int row[maxn],col[maxn];
 int N[maxr],use[maxr],head[maxr];
 void init(){
     memset(head,-,sizeof(head));
     memset(N,,sizeof(N));
     for(int i=;i<=m;++i){
         L[i]=i-,R[i]=i+;
         U[i]=D[i]=i;
         row[i]=,col[i]=i;
     }
     L[]=m,R[m]=;
     cnt=m;
 }
 void remove(int c){// 删除列以及所在列含有1的行
     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],--N[col[j]];
 }
 void resume(int c){// 恢复列以及所在列含有1的行
     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,++N[col[j]];
     L[R[c]]=R[L[c]]=c;
 }
 int low(){
     int mi=maxr,idx=;
     for(int i=R[];i;i=R[i])if(N[i]<mi&&N[i])mi=N[i],idx=i;
     return idx;
 }
 void link(int r,int c){
     ++N[c],++cnt;
     row[cnt]=r,col[cnt]=c;
     U[cnt]=U[c],D[cnt]=c;
     U[D[cnt]]=D[U[cnt]]=cnt;
     if(head[r]==-)
         head[r]=L[cnt]=R[cnt]=cnt;
     else {
         L[cnt]=L[head[r]];
         R[cnt]=head[r];
         L[R[cnt]]=R[L[cnt]]=cnt;
     }
 }
 int xy[maxr][];char num[maxr];
 char ch[][];
 bool dance(int dep){
     if(R[]==)return true;
     int c=low();
     if(c==)return false;
     for(int i=D[c];i!=c;i=D[i]){
         int r=row[i];
         ch[xy[r][]][xy[r][]]=num[r];
         use[dep]=i;
         remove(col[i]);
         for(int j=R[i];j!=i;j=R[j])remove(col[j]);
         if(dance(dep+))return true;
         for(int j=L[i];j!=i;j=L[j])resume(col[j]);
         resume(col[i]);
     }
     return false;
 }

 int main(){
     int t,ca=;
     scanf("%d",&t);
     while(t--){
         if(ca)scanf("%s",ch[]);
         if(ca)puts("---");
         ca=;
         n=,m=;
         init();
         for(int i=;i<=;++i)scanf("%s",ch[i]+);
         for(int i=;i<=;++i){
             for(int j=;j<=;++j){
                 if(ch[i][j]>=''&&ch[i][j]<=''){
                     int val=ch[i][j]-'';
                     ++n;
                     link(n,(i-)*+val);
                     link(n,+(j-)*+val);
                     link(n,+(i-)*+j);
                     link(n,+((i-)/*+(j+)/-)*+val );
                     xy[n][]=i,xy[n][]=j;
                     num[n]=val+'';
                 }
                 else {
                     for(int k=;k<=;++k){
                         int val=k;
                         ++n;
                         link(n,(i-)*+val);
                         link(n,+(j-)*+val);
                         link(n,+(i-)*+j);
                         link(n,+((i-)/*+(j+)/-)*+val );
                         xy[n][]=i,xy[n][]=j;
                         num[n]=val+'';
                     }
                 }
             }
         }
         if(dance())for(int i=;i<=;++i)printf("%s\n",ch[i]+);
         else puts("impossible");
     }
     return ;
 }