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,-1,sizeof(head));
    memset(N,0,sizeof(N));
    for(int i=0;i<=m;++i){
        L[i]=i-1,R[i]=i+1;
        U[i]=D[i]=i;
        row[i]=0,col[i]=i;
    }
    L[0]=m,R[m]=0;
    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=0;
    for(int i=R[0];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]==-1)
        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][2];char num[maxr];
char ch[12][12];
bool dance(int dep){
    if(R[0]==0)return true;
    int c=low();
    if(c==0)return false;
    for(int i=D[c];i!=c;i=D[i]){
        int r=row[i];
        ch[xy[r][0]][xy[r][1]]=num[r];
        use[dep]=i;
        remove(col[i]);
        for(int j=R[i];j!=i;j=R[j])remove(col[j]);
        if(dance(dep+1))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=0;
    scanf("%d",&t);
    while(t--){
        if(ca)scanf("%s",ch[0]);
        if(ca)puts("---");
        ca=1;
        n=0,m=324;
        init();
        for(int i=1;i<=9;++i)scanf("%s",ch[i]+1);
        for(int i=1;i<=9;++i){
            for(int j=1;j<=9;++j){
                if(ch[i][j]>='0'&&ch[i][j]<='9'){
                    int val=ch[i][j]-'0';
                    ++n;
                    link(n,(i-1)*9+val);
                    link(n,81+(j-1)*9+val);
                    link(n,162+(i-1)*9+j);
                    link(n,243+((i-1)/3*3+(j+2)/3-1)*9+val );
                    xy[n][0]=i,xy[n][1]=j;
                    num[n]=val+'0';
                }
                else {
                    for(int k=1;k<=9;++k){
                        int val=k;
                        ++n;
                        link(n,(i-1)*9+val);
                        link(n,81+(j-1)*9+val);
                        link(n,162+(i-1)*9+j);
                        link(n,243+((i-1)/3*3+(j+2)/3-1)*9+val );
                        xy[n][0]=i,xy[n][1]=j;
                        num[n]=val+'0';
                    }
                }
            }
        }
        if(dance(0))for(int i=1;i<=9;++i)printf("%s
",ch[i]+1);
        else puts("impossible");
    }
    return 0;
}